[Octave-bug-tracker] [bug #47236] output of eigs() with 'sa' option differs from Matlab

[Doug Stewart]

Follow-up Comment #5, bug #47236 (project octave):

I have narrowed it down some more.


clear
load yaledata;

[v, m] = eigs (A, 3, "sa");
m
L= diag(sum(A))-A;
[v, d] = eigs (L, 3, "sa");
d
gplot (A, v(:,[2 1]))


m shows that it is sorted from small to big
but d is wrong.

why the eigs of L are sorted different than the eigs of A is the problem.




(file #36443)
    _______________________________________________________

Additional Item Attachment:

File name: yaledata                       Size:61 KB


    _______________________________________________________

Reply to this item at:

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

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