| Re: [eigen] Plans for EigenSolver |
[ Thread Index |
Date Index
| More lists.tuxfamily.org/eigen Archives
]
- To: eigen@xxxxxxxxxxxxxxxxxxx
- Subject: Re: [eigen] Plans for EigenSolver
- From: Gael Guennebaud <gael.guennebaud@xxxxxxxxx>
- Date: Tue, 18 Aug 2009 16:40:50 +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 :content-transfer-encoding; bh=JJjyS84X7KGCBEOEms54kqMjp3igaVR4nF2ZeyzvkHo=; b=vflI9oRxGF0cZPbZlAMCWUE6xfC4iS6WYZ7IPt+WiqQoq9h6A4P2xRTxuVdBnv07KX RVyGTmRKDPEriD4VkmtWKPXGxZwe2+wimjvz8ctxu7mF7srkmRIJVhTKwjIQvCmjG+ef 7VTXGJ3A7iMVprGlwvPN5rZLLSC9zJ/TZtaNY=
- 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=QOQsm/KUpIK3PcQ0S2UbFoXsT9mKqYR+O/iPk6RjKEdX2UGwyewtbpK9bGSo+Utdds 8UhrpSt3RBCSuRSiQIraGKrvPYBx+sz0pX0+wUp5HgbJTeUZnzHr6SSfG2h8lncREnn1 H2AxtNPiQQ2pPQwWuskYPTIF2kw3ba+w6oghk=
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.
gael.
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.
>
> Benoit
>
>
>