Symmetric covariance matrix causes exception in multi_normal_lpdf for being non-symmetric

[jsocolar]

I don’t know for certain that this will fix the problem (I suspect it will), but you could check if tcrossprod is more stable. See here: 5.5 Dot Products and Specialized Products | Stan Functions Reference