Как решить линейную систему с использованием трехзначного выходного вектора представления данных в качестве входных данных для Cholesky?

Question

Я пытаюсь изучить некоторые методы для решения линейных систем с Python. Я реализовал некоторые из этих методов. Теперь я бы хотел проверить их с помощью больших и разреженных матриц. Чтобы сделать это, я начал узнавать о методе представления трех столбцов, потому что заметил, что я должен уменьшить мою разреженную матрицу, прежде чем вводить ее в свой метод. Три представления столбцов кажутся простыми, но я не могу понять, как использовать его вывод как вход моего метода Холески (например). Как использовать его вывод (массив из трех столбцов с, значения и ссылки) в качестве ввода моего метода? Нужно ли переписывать мой метод Холески?

Вот мой Cholesky метод: https://raw.githubusercontent.com/angellicacardozo/linear-algebra-methods/master/P03CHOLE.py

Спасибо

Answer 1

Может быть, это может помочь вам:

For i = 1 To n 
    For j = 1 To i 
     Sum = a(i, j) 
     For k = 1 To j-1 
      Sum = Sum - a(i, k) * a(j, k) 
     If i > j Then 
      a(i, j) = Sum/a(j, j) 
     Else If Sum > 0 Then 
      a(i, i) = Sqrt(Sum) 
     Else 
      ERROR