[Octave-bug-tracker] [bug #59488] eigs error in dsaupd for zero matrix

[John W. Eaton]

Follow-up Comment #7, bug #59488 (project octave):

In addition to matrices of zeros, do we also need to have special cases for
matrices that contain Inf or NaN?  For example, I see the following:


octave> x = sprand (20, 20, 0.1); x(10) = NaN; eigs (x, 1)
DLASCL: parameter number 4 is invalid
error: Fortran procedure terminated by call to XERBLA
error: called from
    eigs at line 285 column 20

octave> x = sprand (20, 20, 0.1); x(10) = Inf; eigs (x, 1)
error: __eigs__: eigs: error in dneupd: DNAUPD did not find any eigenvalues to
sufficient accuracy.
This should not happen.  Please, see
https://www.gnu.org/software/octave/bugs.html, and file a bug report
error: called from
    eigs at line 285 column 20



    _______________________________________________________

Reply to this item at:

  <https://savannah.gnu.org/bugs/?59488>

_______________________________________________
  Message sent via Savannah
  https://savannah.gnu.org/