Batched Gaussian Processes

Have you look into using MatrixNormal or KroneckerNormal?