[Top][All Lists]

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

[Octave-bug-tracker] [bug #44004] eig and eigs routine

From: Lachlan Andrew
Subject: [Octave-bug-tracker] [bug #44004] eig and eigs routine
Date: Sat, 24 Jan 2015 05:41:13 +0000
User-agent: Mozilla/5.0 (X11; Linux x86_64; rv:35.0) Gecko/20100101 Firefox/35.0

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

Thanks, Jordi, for pointing out that the initialization is random.  However,
that is not causing the problem here.

It seems that the call to the   dnaupd   function in ARPACK works reliably
(and that is where the randomness in the algorithm is).

The problem is in the call to  dneupd.  Most of the time, it correctly returns
the first  k  values of workl[iptr(5)] as expected.  However sometimes it
returns the *last* two (which are not valid eigenvalues) instead of the first

Sometimes it replaces three or four, depending on the value of  k.

Note that the main difference between
  D = eigs(...)
  [V,D] = eigs(...)
is in the first parameter passed to  dneupd.

bug #41305 also suggests that there is a (different) bug in  dneupd, so I
wouldn't rule out that possibility.


Reply to this item at:


  Message sent via/by Savannah

reply via email to

[Prev in Thread] Current Thread [Next in Thread]