Mail Archives: djgpp-workers/1998/09/15/08:23:09
Hello Workers!
I've found a problem in the latest libm.a while testing perl 5.005_02.
Perl fails on two of its tests when I link it with djdev202a (980726+),
but it works with the plain old djdev201. The two failed tests came
from t/lib/complex.t (which contains 1000+ tests):
not ok 975
# 'acoth $z0' expected: '0.14694666622553-0.2318238045004i'
got: '0.14694666622553-0.270209750135292i'
for z = ( 2, 3)
not ok 978
# 'acoth $z0' expected: '0.14694666622553+0.2318238045004i'
got: '0.14694666622553+0.270209750135292i'
for z = ( 2,-3)
The imaginary parts seem to be very unaccurate. Any idea?
Laszlo
ps: Perl uses the following formula to compute acoth:
acoth(z) = 1/2 log((1+z) / (z-1))
- Raw text -