[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
[Octave-bug-tracker] [bug #49523] eigs (R2013b and later) uses NaNs, not
From: |
Dmitri A. Sergatskov |
Subject: |
[Octave-bug-tracker] [bug #49523] eigs (R2013b and later) uses NaNs, not zeros, for unconverged Ritz values |
Date: |
Thu, 30 Nov 2017 00:33:44 -0500 (EST) |
User-agent: |
Mozilla/5.0 (X11; Fedora; Linux x86_64; rv:57.0) Gecko/20100101 Firefox/57.0 |
Follow-up Comment #23, bug #49523 (project octave):
I also get different test failures depending on the blas version used:
LD_PRELOAD=/usr/lib64/libopenblas.so.0 ./run-octave -f -q --no-gui
octave:1> test eigs
***** 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
LD_PRELOAD=/usr/lib64/atlas/libsatlas.so ./run-octave -f -q --no-gui
octave:1> test eigs
***** 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
ASSERT errors for: assert (d (3:4),[NaN + 1i * NaN; NaN + 1i * NaN])
Location | Observed | Expected | Reason
(1) 1+1.783e-14i NaN+NaNi 'NaN' mismatch
LD_PRELOAD=/usr/lib64/libblas.so.3 ./run-octave -f -q --no-gui
octave:1> test eigs
***** testif HAVE_ARPACK
A = toeplitz ([0, 1, zeros(1, 8)], [0, -1, zeros(1, 8)]);
A(1, 1) = 1;
A = kron (A, eye (10)) + kron (eye (10), A);
opts.v0 = (1:100)';
opts.maxit = 1;
warning ("off", "Octave:eigs:UnconvergedEigenvalues", "local");
d = eigs (A, 4, "sm", opts);
assert (d(3:4), [NaN; NaN]);
assert (imag (d(3:4)), [0; 0]);
!!!!! test failed
ASSERT errors for: assert (d (3:4),[NaN; NaN])
Location | Observed | Expected | Reason
(1) 0.067993 NaN 'NaN' mismatch
Those are 100% reproducible.
I use
arpack-devel-3.5.0-4.fc27.x86_64
arpack-3.5.0-4.fc27.x86_64
(same as fedora buildbots)
Dmitri.
_______________________________________________________
Reply to this item at:
<http://savannah.gnu.org/bugs/?49523>
_______________________________________________
Message sent via/by Savannah
http://savannah.gnu.org/
- [Octave-bug-tracker] [bug #49523] eigs (R2013b and later) uses NaNs, not zeros, for unconverged Ritz values, Rik, 2017/11/14
- [Octave-bug-tracker] [bug #49523] eigs (R2013b and later) uses NaNs, not zeros, for unconverged Ritz values, Rik, 2017/11/23
- [Octave-bug-tracker] [bug #49523] eigs (R2013b and later) uses NaNs, not zeros, for unconverged Ritz values, Marco Caliari, 2017/11/23
- [Octave-bug-tracker] [bug #49523] eigs (R2013b and later) uses NaNs, not zeros, for unconverged Ritz values, Rik, 2017/11/23
- [Octave-bug-tracker] [bug #49523] eigs (R2013b and later) uses NaNs, not zeros, for unconverged Ritz values, Marco Caliari, 2017/11/24
- [Octave-bug-tracker] [bug #49523] eigs (R2013b and later) uses NaNs, not zeros, for unconverged Ritz values, Marco Caliari, 2017/11/27
- [Octave-bug-tracker] [bug #49523] eigs (R2013b and later) uses NaNs, not zeros, for unconverged Ritz values, Marco Caliari, 2017/11/27
- [Octave-bug-tracker] [bug #49523] eigs (R2013b and later) uses NaNs, not zeros, for unconverged Ritz values, anonymous, 2017/11/27
- [Octave-bug-tracker] [bug #49523] eigs (R2013b and later) uses NaNs, not zeros, for unconverged Ritz values, Rik, 2017/11/29
- [Octave-bug-tracker] [bug #49523] eigs (R2013b and later) uses NaNs, not zeros, for unconverged Ritz values, Rik, 2017/11/30
- [Octave-bug-tracker] [bug #49523] eigs (R2013b and later) uses NaNs, not zeros, for unconverged Ritz values,
Dmitri A. Sergatskov <=
- [Octave-bug-tracker] [bug #49523] eigs (R2013b and later) uses NaNs, not zeros, for unconverged Ritz values, Rik, 2017/11/30
- [Octave-bug-tracker] [bug #49523] eigs (R2013b and later) uses NaNs, not zeros, for unconverged Ritz values, Dmitri A. Sergatskov, 2017/11/30
- [Octave-bug-tracker] [bug #49523] eigs (R2013b and later) uses NaNs, not zeros, for unconverged Ritz values, Marco Caliari, 2017/11/30
- [Octave-bug-tracker] [bug #49523] eigs (R2013b and later) uses NaNs, not zeros, for unconverged Ritz values, Dmitri A. Sergatskov, 2017/11/30
- [Octave-bug-tracker] [bug #49523] eigs (R2013b and later) uses NaNs, not zeros, for unconverged Ritz values, Marco Caliari, 2017/11/30
- [Octave-bug-tracker] [bug #49523] eigs (R2013b and later) uses NaNs, not zeros, for unconverged Ritz values, Dmitri A. Sergatskov, 2017/11/30
- [Octave-bug-tracker] [bug #49523] eigs (R2013b and later) uses NaNs, not zeros, for unconverged Ritz values, Marco Caliari, 2017/11/30
- [Octave-bug-tracker] [bug #49523] eigs (R2013b and later) uses NaNs, not zeros, for unconverged Ritz values, Dmitri A. Sergatskov, 2017/11/30
- [Octave-bug-tracker] [bug #49523] eigs (R2013b and later) uses NaNs, not zeros, for unconverged Ritz values, Rik, 2017/11/30
- [Octave-bug-tracker] [bug #49523] eigs (R2013b and later) uses NaNs, not zeros, for unconverged Ritz values, Rik, 2017/11/30