[Octave-bug-tracker] [bug #50546] wrong result from eigs function

[Marco Caliari]

Follow-up Comment #6, bug #50546 (project octave):

The first method involves B\A, which is not symmetric. It is reasonable that
the problem is more difficult. The proper way to rewrite the problem is


R = chol(B);
issymmetric (R' \ A / R, 1e-14)
[w5,d5] = eigs(R' \ A / R, 5, 'lm');
v5 = R \ A;


If you try this example with 'la' instead of 'lm', eigs will tell you that the
matrix is unsymmetric. Probably the test for symmetry in eigs is more
stringent than my issymmetric(...,1e-14).
Now I move to #39573.


    _______________________________________________________

Reply to this item at:

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

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