I have an application for which extremely long precision seems to be required, and along the way, the GSL integration package would also be helpful. (I have started this process using the quad-double C++ library due to Hida, Li, and Bailey.) However, the coefficients for the Gauss quadrature weights (in integration/qk*.c) are all sized to normal doubles. Can anyone here point me at a place where I could get my hands on the original paper, attributed to L.W. Fullerton, where these are published to 80 decimal digits, or else an algorithm to calculate these to that precision? _______________________________________________ Help-gsl mailing list [email protected] http://lists.gnu.org/mailman/listinfo/help-gsl
