/// <summary> /// X = A^T * A; /// </summary> public MatrixXD TMult() { double[] outMatrix = new double[Cols * Cols]; EigenDenseUtilities.TMult(GetValues(), Rows, Cols, outMatrix); return(new MatrixXD(outMatrix, Cols, Cols)); }