[Octave-bug-tracker] [bug #54167] eigs called with user defined function gives incorrect results if B is not the identity matrix

[Marco Caliari]

Follow-up Comment #2, bug #54167 (project octave):

I started to work on this. Eigs...Func were not able to manage the generalized
case. In the end, they should manage it, both with mode=1 and with mode=3. For
the matrix case, the two modes are managed by different functions
(Eigs...Matrix and Eigs...MatrixShift). For the moment, I implemented just
mode=1 for EigsRealNonSymmetricFunc. I am not sure that a different function
is needed for mode=3. But of course, some if-else statements are necessary. I
added a simple test to eigs.m which fails with current octave and passes with
the attached patch. The original submission require mode=3, so it is not ready
yet. Comments are welcome.

(file #44442)
    _______________________________________________________

Additional Item Attachment:

File name: eigs.diff                      Size:11 KB


    _______________________________________________________

Reply to this item at:

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

_______________________________________________
  Message sent via Savannah
  https://savannah.gnu.org/