Used above code from Multivariate animal model. The topic advanced recently. There are faster and better ways to calculate the Kronecker product.
@jenvent I don’t mean that the product is faster on vectors than on matrices. I mean sampling on the cholesky-form is much faster and therefore we need the kronecker product based upon cholesky decomposed matrices.