protected static double[,] chol(double[,] a) { var chol = new CholeskyDecomposition(a); return chol.LeftTriangularFactor; } #endregion