>>>>> "Robert" == Robert Dodier <robert_dodier@yahoo.com> writes:
Robert> As Brian Gough points out, there is a version (GPL) of Quadpack
Robert> in the GSL which contains bug fixes. That version is preferable,
Robert> although it is in C which may be harder to translate to Lisp
Robert> than Fortran; perhaps GSL Quadpack can be called through FFI.
A hand-written version in Lisp would allow some neat things to be
done, but I'm converting from Fortran automatically. Maxima doesn't
support an FFI.
Robert> when Quadpack was brought into SLATEC, some clean up was
Robert> done, one bug was fixed (lower limit of integration in
Robert> dqawfe.f and qawfe.f), one was introduced (defabs vs resabs
Robert> in dqagpe.f and qagpe.f), and some decimals were chopped
Robert> off in the Gauss-Kronrod parameters.
Thanks for looking at this! I'll be sure to look into this once
they're in maxima. (I'd like to check in the originals as they are so
we know how they've changed instead of wondering why our version is
different from the SLATEC or netlib version.)
Robert> See also the Quadpack bug list cited by Brian Gough.
Yes, I saw that. I'll be sure to incorporate those as well.
Ray