[Octave-bug-tracker] [bug #36211] bug: crashes and endless loops

[John W. Eaton]

Follow-up Comment #3, bug #36211 (project octave):

I'm attaching a very simple Fortran program that calls dgelsd in a way similar
to what Octave does from C++ using the values from your test.  It shows the
same problems, either an infinite loop or an error in calling dlascl.  The
call to dgelsd seems correct to me in this case, and works for other (non-Inf)
values of A\b.  So unless someone can show that we are doing the wrong thing
in Matrix::lssolve (dMatrix.cc around line 2240), I think the error is in
lapack, not Octave, and so should be reported to the lapack maintainers.

On my Debian system, I have the following packages installed:


  libatlas3gf-base                          3.8.3-27
  liblapack3gf                              3.3.1-1



(file #25672)
    _______________________________________________________

Additional Item Attachment:

File name: dgelsd-bug.f                   Size:0 KB


    _______________________________________________________

Reply to this item at:

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

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