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");
}
}