[Octave-bug-tracker] [bug #53499] Factorization of symmetric positive definite sparse matrix fails whereas Cholesky factorization works

[anonymous]

Follow-up Comment #5, bug #53499 (project octave):

clear all;
load linsolve_test_i.dat
load linsolve_test_j.dat
load linsolve_test_v.dat
load linsolve_test_b.dat
A = sparse(i, j, v);
tic
disp("chol(A)");
[L, P, Q] = chol(A, "lower", "vector");
disp("L^-T L-^1 b");
x(Q, :) = L.' \ (L \ b(Q, :));
toc
max(max(abs(A * x - b)))
disp("A\b");
tic;x=A\b;toc
max(max(abs(A * x - b)))


(file #43719, file #43720, file #43721, file #43722)
    _______________________________________________________

Additional Item Attachment:

File name: linsolve_test_i.7z             Size:514 KB
File name: linsolve_test_v.7z             Size:1785 KB
File name: linsolve_test_j.7z             Size:64 KB
File name: linsolve_test_b.7z             Size:7 KB


    _______________________________________________________

Reply to this item at:

  <http://savannah.gnu.org/bugs/?53499>

_______________________________________________
  Message sent via/by Savannah
  http://savannah.gnu.org/