| Re: [eigen] Status of the eigen solver |
[ Thread Index |
Date Index
| More lists.tuxfamily.org/eigen Archives
]
- To: eigen@xxxxxxxxxxxxxxxxxxx
- Subject: Re: [eigen] Status of the eigen solver
- From: "Timothy Hunter" <tjhunter@xxxxxxxxxxxx>
- Date: Wed, 24 Sep 2008 17:33:31 -0700
- Dkim-signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=gmail.com; s=gamma; h=domainkey-signature:received:received:message-id:date:from:sender :to:subject:in-reply-to:mime-version:content-type :content-transfer-encoding:content-disposition:references :x-google-sender-auth; bh=yf0WeN8R2kZVl9DykUqSxn4ZWAFadIc+hJ49l7bNfIE=; b=Jw+r1IlnESq+b1Np1G3/5SuFfxxe35Nh7FpjznklKkAZ4MzVlTjbNqZgHMHftOjd/T KCwNF7hFT1yvVT+5tsNSeRuDaY8FlUBj/ZYd2AcDpfMkQnVwCo/Mly0KGhdxjJBH6Tmt sX07yUoVGsEL9xfr/Ox5Ufe0aVzj4Nr5KyITM=
- Domainkey-signature: a=rsa-sha1; c=nofws; d=gmail.com; s=gamma; h=message-id:date:from:sender:to:subject:in-reply-to:mime-version :content-type:content-transfer-encoding:content-disposition :references:x-google-sender-auth; b=JZBIci07nxcfRKmiX3qGZrN2EOg3LZ2Yz6iTRsGxvqtw9OjdAgCHMU6bKdT3xSTgwg x0MZAcFBySWEaTSMl8oaYZDu/XwipnDCkgHkixFeWREMjAoWO16qHe2q7Jl6cwYFqShM V4H7/yTeh0euBzKYtFSA0rlk/Gb4Tp6NEGFLE=
For your information, the eigen values returned seem correct (as
checked with random matrices against matlab's computations). Also, the
computed eigen vectors have no imaginary part and are different.
On Wed, Sep 24, 2008 at 1:47 PM, Gael Guennebaud
<gael.guennebaud@xxxxxxxxx> wrote:
>
> ah, thanks for the fix. Actually it is supposed to work for non symmetric
> matrices. At least it worked a while ago. I'll check that when I'll have
> time....
>
> gael.
>
> On Wed, Sep 24, 2008 at 10:41 PM, Benoît Jacob <jacob@xxxxxxxxxxxxxxx>
> wrote:
>>
>> On Wednesday 24 September 2008 22:14:12 Timothy Hunter wrote:
>> > Hello,
>> > I would like to compute eigen values for non-symmetric matrices. I tried
>> > to
>> > use to the EigenSolver but it does not compile because it cannot find
>> > the
>> > function .block(int, int) (see attached log).
>>
>> OK, this function was recently renamed to segment(int,int).
>> Fixed in r864468.
>>
>> > Also, the
>> > tests for non self-adjoint matrices are disabled in test/eigensolver.cpp
>> > .
>> > Does it mean this class is not ready yet?
>>
>> I tried re-enabling this test, and fixed it so it compiled, and the result
>> is
>> that with a real, non-selfadjoint matrix, the unit-test failed at line
>> 117.
>> So please consider it as not yet ready for now.
>>
>> Gael, if it's not yet ready maybe #ifdef it out for now.
>>
>> Cheers,
>> Benoit
>>
>>
>
>