[Octave-bug-tracker] [bug #57196] eigs should detect symmetric/hermitian inputs and return real values only

[David Bateman]

Follow-up Comment #1, bug #57196 (project octave):

Rik,

The basic problem is that ARPACK doesn't supply a ZSAUPD function (or is that
ZHAUPD for hermitian) that is the complex version of the DSAUPD function for
real symmetric problems. So ARPACK doesn't profit from the complex symmetry in
the solution of this type of problem.

The real solution would be to implement the ZSAUPD function,  but that seems
too ambitious. Otherwise in __eigs__ it should be possible to detect hermitian
matrices and at least move the conversion of the eigenvalues to real closer to
the source of the problem

David


    _______________________________________________________

Reply to this item at:

  <https://savannah.gnu.org/bugs/?57196>

_______________________________________________
  Message posté via Savannah
  https://savannah.gnu.org/