Re: [eigen] ordering of eigenvalues of EigenSolver |
[ Thread Index |
Date Index
| More lists.tuxfamily.org/eigen Archives
]
- To: eigen@xxxxxxxxxxxxxxxxxxx
- Subject: Re: [eigen] ordering of eigenvalues of EigenSolver
- From: Susanne Suter <susanne.suter@xxxxxxxxx>
- Date: Tue, 22 Feb 2011 08:57:16 +0100
- Dkim-signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=gmail.com; s=gamma; h=domainkey-signature:mime-version:in-reply-to:references:date :message-id:subject:from:to:content-type:content-transfer-encoding; bh=1CPkfcJSzvJB77GqnURC//DtYzHJPGpCE52hfrZM+kw=; b=dhLHdBK5Mr5DM1VySCZHojmj9Q6BmKfNG/gyK4xvQK+wTQqfbMLuP11dpCwWDmQ1rl kjNRMdnMWSZbmkq0MIcDa+zGtV7eVzaI5wUhdbotRHrGLPp7WvsLWZMRE/DAxccHXjfY cR1V29FuxUEb57jE8B2bl7tIgPskVhsXk6G88=
- Domainkey-signature: a=rsa-sha1; c=nofws; d=gmail.com; s=gamma; h=mime-version:in-reply-to:references:date:message-id:subject:from:to :content-type:content-transfer-encoding; b=WMf2yM3s691temgNe0TUv6E19e5LbsCyrMO0gic2afjJrTyHpqQSm5plNkS5jYIPT8 J1tmFdNP1EWpKQr0AVLNnicGhJ0DcrYO0MK4IMg4XH64Xu5jNnVfPwzHzEahzP6OLGY3 YF+q8nKSPEywRqkkyRinrW+OFETIchmyEOhDE=
hi
thank you for taking the time to check that issue.
I noticed now that I use instead of EIG/SSYEV the function EIGS in
MatLab and the function SSYEVX in Lapack, which both return the x most
significant eigenvalues. so I've learned by now that there are two (or
more) types of eigenvalue decomposition methods provided in Matlab and
Lapack while eigen supports one.
thx again,
susanne
On Tue, Feb 22, 2011 at 8:35 AM, Gael Guennebaud
<gael.guennebaud@xxxxxxxxx> wrote:
> On Tue, Feb 22, 2011 at 1:18 AM, Benoit Jacob <jacob.benoit.1@xxxxxxxxx> wrote:
>>>
>>> Furthermore, other eigenvalue decomposition libraries seem to order
>>> eigenvalues, too. E.g., MatLab (based on ARPACK)
>>
>> Ah, I didn't know about Matlab.
>
> actually MatLab does as us and Lapack do, for instance:
>
> $ eig( diag([1 -1 2 -3 4]))'
>
> ans =
>
> -3 -1 1 2 4
>
> Same for Octave or Scilab.
>
> I don't see any reason to change for a convention that nobody follow
> while our current behavior seems to be quite standard.
>
> On the other hand we could allow for sorting capability such that user
> can easily sort rows or column as they like.
>
> gael
>
>
>