public static CellArray Inverse(CellArray A) { if (!CellArray.IsMatrix(A)) { throw new Exception("A must be a matrix"); } int A_Rows = A.Count; int A_Cols = CellArray.ColumnCount(A); if (A_Rows != A_Cols) { throw new Exception(string.Format("Matrix must be square {0}:{1}", A_Rows, A_Cols)); } LUDecomposition engine = new LUDecomposition(A); CellArray B = CellArray.Identity(A_Rows); return(engine.solve(B)); }
public static CellArray Identity(int Dimension) { return(CellArray.Identity(Dimension, CellAffinity.DOUBLE)); }