| Subject: | Error in t/slatec |
My config: perl-5.8.3 (shared libperl, ithreads, Solaris 8, gcc-3.2.3)
Pdl-2.4.1
I get "not ok 37" in t/slatec; I added some debugging output, see
here:
ok 33
DEBUG: x=[ -0.3 0.7 1.7 2.7 3.7 4.7 5.7
6.7 7.7 8.7 9.7]
DEBUG: f=[ 0.09 0.49 2.89 7.29 13.69 22.09 32.49
44.89 59.29 75.69 94.09]
DEBUG: d=[0 0.685714 3.10588 5.21482 7.26486 9.29362 11.3123 13.3254
15.3351 17.3425 19.4], err=0
DEBUG: ans=[243 170.33333]
DEBUG: int=[ 243.03072 170.29642], err=[0 0]
ok 34
ok 35
DEBUG: hi=[ 8.6999998 6.6999998], lo=[-0.30000001 0.69999999]
ok 36
DEBUG: int=[3.67884 3.39096], err=[0 0], ans=[ 219.50999 100.13999],
abs=[ 215.83115 96.749028], all=1
not ok 37
To be honest, I do not understand the calculations, but hopefully
you know what is wrong here.
Let me know if I should try any additional tests
For comparison, on Linux (RedHat 7.2, gcc-3.2.3, same perl) I get this
correct result:
DEBUG: int=[219.565 100.087], err=[0 0], ans=[ 219.50999 100.13999],
abs=[0.054787995 0.053306263], all=1
ok 37
Best regards,
Marek