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

[Marco Caliari]

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

Hi,

your example is not correct, since B is not positive definite. On the other
hand, I think you are right and line 1275 (and 2238 and 3259) of eigs-base.cc
should be changed to


vector_product (b, workd+iptr(0)-1, dtmp);


as for the case ido == 2.

Marco

    _______________________________________________________

Reply to this item at:

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

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