Hi Ricardo, Thank you for the quick answer!
using → shape = xys.shape
Leads to the same error
xyh = EulerMaruyama(“xyh”, dt, osc_sde, (τh, ah), shape= xys.shape, testval=xys)
52 zh = pm.Normal(“zh”, mu=mh * xyh[:, 0] + (1 - mh) * xyh[:, 1], sigma=0.1, observed=zs)
53
8 frames
/usr/local/lib/python3.10/dist-packages/pytensor/tensor/var.py in getitem(self, args)
496 # Check if the number of dimensions isn’t too large.
497 if self.ndim < index_dim_count:
→ 498 raise IndexError(“too many indices for array”)
499
500 # Convert an Ellipsis if provided into an appropriate number of
IndexError: too many indices for array
Am I missing something?