At Mon, 25 Aug 2008 17:02:49 +0200, Thies Heidecke wrote: > line 489 reads: > const double tb = -log(y) * (b[0] * y*(b[1] + y*b[2])); > > but it should be > const double tb = -log(y) * (b[0] + y*(b[1] + y*b[2])); > > after recompiling the function now gives reasonable values. >
Thanks for the excellent bug report. I have committed your correction to the source repository for the next release. http://git.savannah.gnu.org/gitweb/?p=gsl.git;a=commit;h=163bd7269d3c74af6b0a81ac82d5396a2afd3b3e -- Brian Gough GNU Scientific Library - http://www.gnu.org/software/gsl/ _______________________________________________ Bug-gsl mailing list [email protected] http://lists.gnu.org/mailman/listinfo/bug-gsl
