[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: |
Marco Caliari |
Subject: |
[Octave-bug-tracker] [bug #49523] eigs (R2013b and later) uses NaNs, not zeros, for unconverged Ritz values |
Date: |
Tue, 11 Apr 2017 03:52:57 -0400 (EDT) |
User-agent: |
Mozilla/5.0 (X11; Ubuntu; Linux x86_64; rv:52.0) Gecko/20100101 Firefox/52.0 |
Follow-up Comment #9, bug #49523 (project octave):
In the nonsymmetric case, the nonconverged Ritz values are set to NaN+NaNi in
Matlab. You can try
A = toeplitz ([-2, 1, zeros(1, 8)]) + magic (10) * 1e-10;
A = kron (A, eye (10)) + kron (eye (10), A);
opts.v0 = (1:100)';
opts.maxit = 1;
eigs (A, 8, 'sm', opts)
The converged values are complex with imaginary part zero. If you omit opts,
you see that all the converged eigenvalues are double. If we follow Matlab and
comment #5, in case of nonconverged Ritz values we should set them and the
corresponding eigenvectors to NaN+NaNi. An alternative would be: if the
converged Ritz values are double, set the nonconverged and the corresponding
eigenvectors to (double) NaN. This would save some memory.
_______________________________________________________
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, Marco Caliari, 2017/04/07
- [Octave-bug-tracker] [bug #49523] eigs (R2013b and later) uses NaNs, not zeros, for unconverged Ritz values, Marco Caliari, 2017/04/07
- [Octave-bug-tracker] [bug #49523] eigs (R2013b and later) uses NaNs, not zeros, for unconverged Ritz values, Marco Caliari, 2017/04/10
- [Octave-bug-tracker] [bug #49523] eigs (R2013b and later) uses NaNs, not zeros, for unconverged Ritz values, Rik, 2017/04/10
- [Octave-bug-tracker] [bug #49523] eigs (R2013b and later) uses NaNs, not zeros, for unconverged Ritz values, anonymous, 2017/04/10
- [Octave-bug-tracker] [bug #49523] eigs (R2013b and later) uses NaNs, not zeros, for unconverged Ritz values, Rik, 2017/04/10
- [Octave-bug-tracker] [bug #49523] eigs (R2013b and later) uses NaNs, not zeros, for unconverged Ritz values, anonymous, 2017/04/10
- [Octave-bug-tracker] [bug #49523] eigs (R2013b and later) uses NaNs, not zeros, for unconverged Ritz values,
Marco Caliari <=
- [Octave-bug-tracker] [bug #49523] eigs (R2013b and later) uses NaNs, not zeros, for unconverged Ritz values, Rik, 2017/04/11
- [Octave-bug-tracker] [bug #49523] eigs (R2013b and later) uses NaNs, not zeros, for unconverged Ritz values, Marco Caliari, 2017/04/12
- [Octave-bug-tracker] [bug #49523] eigs (R2013b and later) uses NaNs, not zeros, for unconverged Ritz values, Rik, 2017/04/12
- [Octave-bug-tracker] [bug #49523] eigs (R2013b and later) uses NaNs, not zeros, for unconverged Ritz values, Rik, 2017/04/16