public Matrix GetP() // Function returns permutation matrix "P" due to permutation vector "pi" { if (L == null) MakeLU(); Matrix matrix = ZeroMatrix(rows, cols); for (int i = 0; i < rows; i++) matrix[pi[i], i] = 1; return matrix; }