[Octave-bug-tracker] [bug #53390] lu factorization of large sparse matrices often fails

[Marco Caliari]

Follow-up Comment #8, bug #53390 (project octave):

I think we have to disable for the moment the calls [L,U]=lu(S) and
[L,U,P]=lu(S) for sparse matrices. Even if umfpack does not fail, the result
is unpredictable. I do not think it is a big issue: I cannot imagine how to
use an LU factorization of a sparse matrix without pivoting and reduction of
fill-in ([L,U,P,Q]=lu(S) syntax). Maybe for didactic purposes. Then, there are
some options:

1) we implement it from scratch
2) GSOC
3) we try to adapt ilu (as far as I know, ilu is lu with an if clause)

    _______________________________________________________

Reply to this item at:

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

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