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

[Rik]

Update of bug #49523 (project octave):

                  Status:                   Fixed => In Progress            

    _______________________________________________________

Follow-up Comment #22:

Have to re-open the bug temporarily.  The BIST tests are creating intermittent
failures on the automated buildbot testing system.  Dmitri first reported
this.  See log files at http://buildbot.octave.org:8010/waterfall.

Some sample errors are:


***** testif HAVE_ARPACK
 A = toeplitz ([-2, 1, zeros(1, 8)]);
 A = kron (A, eye (10)) + kron (eye (10), A);
 opts.v0 = (1:100)';
 opts.maxit = 3;
 warning ("off", "Octave:eigs:UnconvergedEigenvalues", "local");
 d = eigs (A, 4, "lm", opts);
 assert (d(2:4), [NaN; NaN; NaN]);
!!!!! test failed
ASSERT errors for:  assert (d (2:4),[NaN; NaN; NaN])

  Location  |  Observed  |  Expected  |  Reason
    (1)         -7.365        NaN        'NaN' mismatch
***** testif HAVE_ARPACK
 A = toeplitz ([-2, 1, zeros(1, 8)]);
 A = kron (A, eye (10)) + kron (eye (10), A);
 Afun = @(x) A * x;
 opts.v0 = (1:100)';
 opts.maxit = 3;
 opts.issym = true;
 warning ("off", "Octave:eigs:UnconvergedEigenvalues", "local");
 d = eigs (Afun, 100, 4, "sm", opts);
 assert (d(2:4), [NaN; NaN; NaN]);
!!!!! test failed
ASSERT errors for:  assert (d (2:4),[NaN; NaN; NaN])

  Location  |  Observed  |  Expected  |  Reason
    (1)        -0.13578       NaN        'NaN' mismatch
***** testif HAVE_ARPACK
 A = magic (100);
 opts.v0 = (1:100)';
 opts.maxit = 1;
 warning ("off", "Octave:eigs:UnconvergedEigenvalues", "local");
 d = eigs (A, 4, "lm", opts);
 assert (d(4), NaN);
!!!!! test failed
ASSERT errors for:  assert (d (4),NaN)

  Location  |  Observed  |  Expected  |  Reason
     ()       3.8548e-11      NaN        'NaN' mismatch
***** testif HAVE_ARPACK
 A = magic (100) / 100 + eye (100);
 opts.v0 = (1:100)';
 opts.maxit = 1;
 warning ("off", "Octave:eigs:UnconvergedEigenvalues", "local");
 d = eigs (A, 4, "sm", opts);
 assert (d(3:4), [NaN+1i*NaN; NaN+1i*NaN]);
!!!!! test failed
eigs: error -14 in dneupd
***** testif HAVE_ARPACK
 A = 1i * magic (100) / 100 + eye (100);
 opts.v0 = (1:100)';
 opts.maxit = 1;
 warning ("off", "Octave:eigs:UnconvergedEigenvalues", "local");
 d = eigs (A, 4, "sm", opts);
 assert (d(2:4), [NaN+1i*NaN; NaN+1i*NaN; NaN+1i*NaN]);
!!!!! test failed
eigs: error -14 in zneupd


This is stressing the ARPACK library.  And maybe it indicates the ARPACK
library used by the buildbots is bad.


    _______________________________________________________

Reply to this item at:

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

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