|Re: [eigen] Progress and plans|
[ Thread Index |
| More lists.tuxfamily.org/eigen Archives
- To: eigen@xxxxxxxxxxxxxxxxxxx
- Subject: Re: [eigen] Progress and plans
- From: Benoit Jacob <jacob.benoit.1@xxxxxxxxx>
- Date: Mon, 22 Feb 2010 16:03:18 -0500
- 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=TdlGA0ZADDuxcFxTvSA0iEDN+2dfH3pF/mfI7ii6lfw=; b=AspPXedi3Mv+79iPmBl60g6JshAfFgt8N2FtF2szFRsXIcFy21NFel4WtD9zuQkcSY fcAB19ILcrTSDuGY8i3wMJSIcrPszes9nFlKudCuFhCEm+vsYTkGXAm0NzvHCFvRFHB9 x+8uN7glJVoUcdbFIPGsfxhCDsmKT77Gq/7Dw=
- 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=ZknVOol77xOMQ4zIPiFE0UAjWXstqfShItzwsEb4VDmVQKqw6VvbM/5o2jyZjFMrdK s3y36LCRG+8ZMZJkekcl7lC6e0wR89eB7SH5+/NjOrRY+dQpAxRf3KRdJVno9mptlZnG fFNBhn7HuoKPaI6uklSTx/TnHwWIGJxrLZcjo=
2010/2/22 Jitse Niesen <jitse@xxxxxxxxxxxxxxxxx>:
> On Thu, 18 Feb 2010, Benoit Jacob wrote:
>>> I think the module is usable at the moment. However, there are some
>>> if the matrices have clustered eigenvalues (this is the most difficult
>>> case). One issue is that ComplexSchur often reaches the maximum number of
>>> iterations. I don't know whether that's an inherent weakness in the
>>> or only in the implementation - it may be as simple as increasing the
>>> maximum number of iterations.
>> It'd be interesting to know if increasing the number of iterations helps..
> It helps but not a lot. The algorithm starts fine but gets stuck before it
> reaches the termination condition. It seems that either the condition is too
> strict, or there is an instability in some part of the computation. I'll
> look further into it.
ok. there are many examples of iterative algorithms where the last bit
of precision is hard to get. for example in jacobisvd, for some
matrices, it is, which is why i currently terminate when the relative
error is < 2*epsilon.
>>> There seems to be another issue in matrix_function_3, which I haven't
>>> investigated yet.
> I solved this. The mistake was not in the computation of the matrix
> function, but in the computation of the supposedly exact value I compared
> the result with. :-/
>> That's certainly the most elegant possible syntax!
> OK, I can agree with that. The only worry is that array.exp() means
> something else: apply the exponential function entry-wise. So there is a
> source of confusion. On the other hand, array * array and matrix * matrix
> already have different meanings.
Exactly. That is exactly the point of having Array vs. Matrix classes.
So, that a function has a different meaning on either side, is
> I will change the code to the matrix.exp() syntax soon if nobody argues
> against it.
>>> * related, what needs to be done to get the module out of unstable/ ?
>> You're doing pretty much exactly what you have to :)
>> 1) say that you want it out of unstable/ ---> check
>> 2) make sure the API is ready for prime time ---> in progress
>> 3) tests, documentation ----> check
>> 4) commit to maintaining that stuff in the future ? ;)
> Ad 4) Of course. Otherwise it's not worth writing in the first place.