Bug corrected in matrix multiplication.
Added simple tests for multivariate normal sampling.