static void Main(string[] args) { Console.WriteLine("Number of rows/columns:"); int size = int.Parse(Console.ReadLine()); double[,] matrixA = new double[size, size]; matrixA = MultidimensionalArrayOperations.Read(size); Console.WriteLine(); Console.Write("Matrix A:"); MultidimensionalArrayOperations.Write(matrixA, size); new Cholesky(matrixA, size); }
public Cholesky(double[,] matrixA, int size) { if (IsSymmetric(matrixA, size)) { double[,] matrixL = decomposition(matrixA, size); double[,] transposedMatrixL = MultidimensionalArrayOperations.Transpose(matrixL, size); Console.WriteLine(); Console.Write("Matrix L:"); MultidimensionalArrayOperations.Write(matrixL, size); Console.WriteLine(); Console.Write("Transposed matrix L:"); MultidimensionalArrayOperations.Write(transposedMatrixL, size); } else { Console.WriteLine(); Console.WriteLine("The matrix is not symmetric"); } }