I just got the GSL book today and immediately
tried the sample program on page 9.
I did not get the exact same answer shown in the book:
JO(5) = -1.775967713143382920e-01
I got J)(5) = -1.775967713143382642e-01
^^^
I'm running 64-bit AMD OpenBSD 4.1 on AM2 Athlon 64.
gcc (GCC) 3.3.5 (propolice)
Copyright (C) 2003 Free Software Foundation, Inc.
This is free software; see the source for copying conditions. There is NO
warranty; not even for MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.
FP computation is done using MMX registers except for long double which
is done using the 8087 registers. I am now working on a version of GSL
compiled using long double to force use of the 80-bit FP data format
and see if I get the result listed in the book.
Comments?
Thanks,
Dave Feustel
--
_______________________________________________
Help-gsl mailing list
[email protected]
http://lists.gnu.org/mailman/listinfo/help-gsl
