|Re: [eigen] Plans for EigenSolver|
[ Thread Index |
| More lists.tuxfamily.org/eigen Archives
- To: eigen@xxxxxxxxxxxxxxxxxxx
- Subject: Re: [eigen] Plans for EigenSolver
- From: Gael Guennebaud <gael.guennebaud@xxxxxxxxx>
- Date: Fri, 4 Sep 2009 17:43:46 +0200
- Dkim-signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=gmail.com; s=gamma; h=domainkey-signature:mime-version:received:in-reply-to:references :date:message-id:subject:from:to:content-type; bh=iQdp9miTK8BmOcOwyqQzZZNpsP3kAqibNISFB6ugvH8=; b=iW14H4K4JBrMLx3vcipF3g1FX7nv3/Hg8hsNUM0Pic5X/EbOcoaLUHRI0bVMW82KpD dbnV8VCmD+PJ7H9xRvDLRb0wzLDA1QhxgMRElPxbMwHj7jDyXaZiwtZSknxYxKmEL1Ai A/DuFlAINnMdmvearAY+zPf9XB6ReowZ1kA0s=
- 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; b=fliThP5QarifOC8MWlVoxu6arQwUup0AB1EARFuL1XfrosLc7M9VvGs5GBUZy+9Fbu o12YF0FXVb+xn7fvcXExkYzYuiiQX5I8/Ak6477v+FQubvglgzry1TBE0Bk1qly28hjz S73KPbm2s3Jzxl07lInbKuAMnNEHvsLx0JaL8=
to update the list:
- I've split the QR module into a pure QR one and an Eigenvalues
module for all eigen values related stuff
- I've recently merged this patch into the main devel branch: it
provides a ComplexSchur decomposition and a ComplexEigenSolver for
- Actually the real and complex Schur decompositions are too different
to factorize them, and the way you use them is also very different. So
I plan to have a RealSchur decomposition but no generic Schur
- On the other hand once we'll have a clean RealSchur class, I plane
to have a unique EigenSolver and remove the ComplexEigenSolver class.
if anybody have some idea to build tricky matrices for our eigen solvers...
On Sun, Aug 23, 2009 at 6:29 PM, Gael
> here is a general eigen solver working for complexes:
> this a very rough implementation ported from JAMAPACK, but at least it
> passes all the unit tests and unlike the current EigenSolver the
> algorithm is easy to follow.
> Any help to make it better is more than welcome, some ideas:
> - make it suitable for real (it is currently specialized for complexes)
> - more unit tests to check the accuracy, etc.
> these 2 todo are mandatory to include it in Eigen, and then:
> - offer the possibility not to compute the eigenvectors
> - make the internal Shur decomposition available
> - more cleaning
> - improve speed,
> - improve accuracy,
> On Tue, Aug 18, 2009 at 4:40 PM, Gael
> Guennebaud<gael.guennebaud@xxxxxxxxx> wrote:
>> actually someone told me yesterday she managed to hack an EigenSolver
>> for complexes on top of our HessenbergDecomposition and mimicking the
>> implementation of the Shur decomposition of JAMPACK. I waiting for her
>> reply about its integration in Eigen... She acknowledged her
>> implementation was certainly not perfect, but I'm sure this is already
>> a good step in the right direction and that would be a very nice
>> contribution. I'll let you now as soon as I have more info.
>> On Tue, Aug 18, 2009 at 3:22 PM, Benoit Jacob<jacob.benoit.1@xxxxxxxxx> wrote:
>>> 2009/8/18 Jitse Niesen <jitse@xxxxxxxxxxxxxxxxx>:
>>>> 2. I rewrite EigenSolver along the lines of SelfAdjointEigenSolver, possibly
>>>> working with Benoit (and others) to achieve this. This may be more useful
>>>> (there is more interest in eigenvalue problem solvers than matrix function),
>>>> but I'm not sure it's efficient because Benoit works much faster than me.
>>> .... BUT if you want to rewrite EigenSolver on your own, by all means
>>> do it ! ;) it's just that i'm not able to help at the moment.
>>> I think that Gael just made Hessenberg use the new Householder
>>> transformations, which is a good first step, you'll have to do on your
>>> own the rest.