[Octave-bug-tracker] [bug #49523] eigs (R2013b and later) uses NaNs, not zeros, for unconverged Ritz values

[Marco Caliari]

Follow-up Comment #25, bug #49523 (project octave):

@Rik, Dmitri: are the results really intermittent or 100% reproducible?
In the first case, no idea, since in the tests I initialized the initial
vector which should be the only otherwise random choice.
I tried standard blas, atlas and openblas on a Ubuntu systems, I was able to
reproduce some failures. Convergence clearly depends on the blas/lapack
libraries (and the very stringent tolerance, which is eps by default). I
relaxed the tests and now they succeed with all my libraries. Please try the
patch.

(file #42528)
    _______________________________________________________

Additional Item Attachment:

File name: eigs2.diff                     Size:2 KB


    _______________________________________________________

Reply to this item at:

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

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