This is the mail archive of the
libc-alpha@sources.redhat.com
mailing list for the glibc project.
Re: Many new math failures
- From: Stephen L Moshier <steve at moshier dot net>
- To: Andreas Jaeger <aj at suse dot de>
- Cc: <libc-alpha at sources dot redhat dot com>
- Date: Thu, 27 Jun 2002 08:29:39 -0400 (EDT)
- Subject: Re: Many new math failures
> Test "Real part of: cacosh (0.75 + 1.25 i) ==
>1.13239363160530819522266333696834467 +
>1.11752014915610270578240049553777969 i":
>@@ -310,8 +308,8 @@
> float: 3
> idouble: 1
> ifloat: 3
>-ildouble: 37
-ldouble: 37
>+ildouble: 370
>+ldouble: 370
This is due to arithmetic cancellation error inside the function routine.
It might be better to use the identity cacosh(x) = I * cacos(x) here.