public static MatrixR Minor(MatrixR m, int row, int col) { MatrixR mm = new MatrixR(m.GetRows() - 1, m.GetCols() - 1); int ii = 0, jj = 0; for (int i = 0; i < m.GetRows(); i++) { if (i == row) { continue; } jj = 0; for (int j = 0; j < m.GetCols(); j++) { if (j == col) { continue; } mm[ii, jj] = m[i, j]; jj++; } ii++; } return(mm); }
public static VectorR Transform(MatrixR m, VectorR v) { VectorR result = new VectorR(v.GetSize()); if (!m.IsSquared()) { throw new ArgumentOutOfRangeException( "Dimension", m.GetRows(), "The matrix must be squared!"); } if (m.GetCols() != v.GetSize()) { throw new ArgumentOutOfRangeException( "Size", v.GetSize(), "The size of the vector must be equal" + "to the number of rows of the matrix!"); } for (int i = 0; i < m.GetRows(); i++) { result[i] = 0.0; for (int j = 0; j < m.GetCols(); j++) { result[i] += m[i, j] * v[j]; } } return(result); }
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()); }
public static bool CompareDimension(MatrixR m1, MatrixR m2) { if (m1.GetRows() == m2.GetRows() && m1.GetCols() == m2.GetCols()) { return(true); } else { return(false); } }
public MatrixR(MatrixR m) { Rows = m.GetRows(); Cols = m.GetCols(); matrix = m.matrix; }