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

[Marco Caliari]

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

I can reproduce the problem with the mex interface of umfpack in matlab.
Therefore, either there is a bug in umfpack or umfpack is not meant to perform
an LU factorization without column pivoting (in fact, [L, U, P, Q] = lu (S)
does not fail). Togheter with this sentence "With three or fewer outputs, lu
uses its own sparse matrix routines" (available in old matlab documentation) I
think that [L, U] = lu (S) in matlab does not use umfpack. I will contact
umfpack's author.

    _______________________________________________________

Reply to this item at:

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

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