public MatrixR LUInverse(MatrixR m) { int n = m.GetRows(); MatrixR u = m.Identity(); LUDecompose(m); VectorR uv = new VectorR(n); for (int i = 0; i < n; i++) { uv = u.GetRowVector(i); LUSubstitute(m, uv); u.ReplaceRow(uv, i); } MatrixR inv = u.GetTranspose(); return(inv); }
public static MatrixR Adjoint(MatrixR m) { if (!m.IsSquared()) { throw new ArgumentOutOfRangeException( "Dimension", m.GetRows(), "The matrix must be squared!"); } MatrixR ma = new MatrixR(m.GetRows(), m.GetCols()); for (int i = 0; i < m.GetRows(); i++) { for (int j = 0; j < m.GetCols(); j++) { ma[i, j] = Math.Pow(-1, i + j) * (Determinant(Minor(m, i, j))); } } return(ma.GetTranspose()); }