[Octave-bug-tracker] [bug #39573] Wrong result from eigs in generalized shift-invert mode

[Mike Miller]

Update of bug #39573 (project octave):

                  Status:                    None => Fixed                  
             Open/Closed:                    Open => Closed                 

    _______________________________________________________

Follow-up Comment #4:

The original example here gives a different result with the current version of
Octave


>> A = [1 0 1; 0 1 0; 1 0 0];
>> B = [1 0 0 ; 0 1 0 ; 0 0 0];
>> [V, D] = eigs (A, B, 1, 0);
>> residual_norm = norm (A*V - B*V*D)
residual_norm = 0


So at least something has changed with respect to this example since 3.6.4. If
eigs is still not correct for some cases, please feel free to comment on this
bug report or open a new bug report with details about what still is not
working correctly.

    _______________________________________________________

Reply to this item at:

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

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