Any suggestions on continous error "ValueError: shape mismatch: value array of shape (62,) could not be broadcast to indexing result of shape (0,)"

10

I think the slow sampling is related to
U = pm.MvNormal("U", mu = th_t, cov = ct, shape = (tt,specise))
This 62 dimensional normal distribution accounts the most time spend, which is understandable.