Many new math failures

Stephen L Moshier steve@moshier.net
Thu Jun 27 09:22:00 GMT 2002


> 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.






More information about the Libc-alpha mailing list