| Re: [eigen] inverse unit test |
[ Thread Index |
Date Index
| More lists.tuxfamily.org/eigen Archives
]
- To: eigen@xxxxxxxxxxxxxxxxxxx
- Subject: Re: [eigen] inverse unit test
- From: Manuel Yguel <manuel.yguel@xxxxxxxxx>
- Date: Thu, 10 Jun 2010 09:52:32 +0200
- Dkim-signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=gmail.com; s=gamma; h=domainkey-signature:received:mime-version:received:in-reply-to :references:from:date:message-id:subject:to:content-type :content-transfer-encoding; bh=wbYvZt+Zvk+2xlTKQhct/ujl1GLgrCqWxias/NE166Y=; b=gPfytRWmXHUwc2rSqcT/mvR/ehCnA6XunsBIaStxvzyZdVVu3unQBE2pn/qnXWfhU+ cC3gGQliGl0GrNTYdbcc54ZTDCjqcv1PDKURlSGcX98b2TI16Li7asfhLzrVgt4Bj8GA fGa68FGg1tVNu9LrqhBLhSvs6hmKGcMDCi9DU=
- Domainkey-signature: a=rsa-sha1; c=nofws; d=gmail.com; s=gamma; h=mime-version:in-reply-to:references:from:date:message-id:subject:to :content-type:content-transfer-encoding; b=jYMIAYJXHQMK/+/4nVnDCJjLA1BPnjv5htSTGu+tMclTiunkBj9rSghNLEa+9KuUST eCC8ozrjIlRwBOYjhekHy8albgyY2EE46WfIdd6wiYfuYh1zDWhMh1wImc96/lGxn+U/ 0Ooo32TLbQzWHbvYTL9rEyRxF5lwXWmA1ZmL4=
my last e-mail was not clear,
what I wanted to explain is that, the function in MathFunction.h is
not sufficient (in fact it has not the good signature).
However somewhere else, (I guess in Fuzzy.h)
Manu
On Thu, Jun 10, 2010 at 9:42 AM, Manuel Yguel <manuel.yguel@xxxxxxxxx> wrote:
> I hit the same problem several times and the implementation of
> isMuchSmallerThan for float is:
>
> inline bool ei_isMuchSmallerThan(float a, float b, float prec =
> NumTraits<float>::dummy_precision())
> {
> return ei_abs(a) <= ei_abs(b) * prec;
> }
>
> if b is zero, then whatever prec is, it will failed if a is not zero,
> so I do agree with Hauke that his implementation is needed.
> So in fact in my own code I finally used the function he proposed.
>
>
> - cheers,
>
> Manu
>