|
From: | Rafael Gonzalez |
Subject: | Re: eig function project |
Date: | Fri, 19 Apr 2013 12:16:59 -0300 |
On Thu, Apr 18, 2013 at 11:46 AM, Rafael Gonzalez <address@hidden> wrote:
I did some research about the way Matlab does the eig function and I got this table:It seems to me that matching this logic with the proper LAPACK subroutine call in the init function of the EIG and FloatEIG classes a Matlab-complaint eig function can be achieved. Any comments on this?A reference on why balancing is sometimes harmful:And yikes! I just discovered that dgeev does scaling and permutation ('B' arg to dgebal)
unlike dggev which only does permutations ('P' arg to dggbal). I'm going to ask
the lapack folks about this - I gave them a case years ago which failed due to
balancing.--
Ed Meyer
[Prev in Thread] | Current Thread | [Next in Thread] |