Так что я знаю, как решить систему, используя не-CSR способ, но я не могу обернуться вокруг этого метода. Я уже рассчитал матрицу итераций G_J и вектор итераций c_J, и я не знаю, что с ними делать. Я уже искал inte rnet и не нашел ничего полезного. У меня также реализован метод массива CSR_matrix *. Не могли бы вы объяснить мне, что мне нужно сделать? Вот мои аргументы функции:
function [x] = Jacobi_sparse(G_values, G_colind, G_rowptr, c, tol)