octave-bug-tracker
[Top][All Lists]
Advanced

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

[Octave-bug-tracker] [bug #54179] sparse/eigs.m test failure


From: Dmitri A. Sergatskov
Subject: [Octave-bug-tracker] [bug #54179] sparse/eigs.m test failure
Date: Tue, 3 Jul 2018 07:21:36 -0400 (EDT)
User-agent: Mozilla/5.0 (X11; Fedora; Linux x86_64; rv:61.0) Gecko/20100101 Firefox/61.0

Follow-up Comment #10, bug #54179 (project octave):

I have two tests failing:


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

  Location  |  Observed  |  Expected  |  Reason
    (1)      -2.4903e-71+3.6015i   NaN+NaNi     'NaN' mismatch
    (2)      -2.4903e-71-3.6015i   NaN+NaNi     'NaN' mismatch
<...>
***** testif HAVE_ARPACK
 A = toeplitz ([0, 1, zeros(1, 8)], [0, -1, zeros(1, 8)]);
 A = kron (A, eye (10)) + kron (eye (10), A);
 Afun = @(x) A * x;
 opts.v0 = (1:100)';
 opts.maxit = 5;
 warning ("off", "Octave:eigs:UnconvergedEigenvalues", "local");
 d = eigs (Afun, 100, 4, "lm", opts);
 assert (d(3:4), [NaN+1i*NaN; NaN+1i*NaN]);
!!!!! test failed
ASSERT errors for:  assert (d (3:4),[NaN + 1i * NaN; NaN + 1i * NaN])

  Location  |  Observed  |  Expected  |  Reason
    (1)      -2.4903e-71+3.6015i   NaN+NaNi     'NaN' mismatch
    (2)      -2.4903e-71-3.6015i   NaN+NaNi     'NaN' mismatch
<...>
n =  180
nmax =  182



and yes, changing opts.maxit to 4 in both cases allow them to pass.

Dmitri.
-- 


    _______________________________________________________

Reply to this item at:

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

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




reply via email to

[Prev in Thread] Current Thread [Next in Thread]