Of the 5 (the last one that I checked of course), CEPHES has just raw source, like LAPACK, and its in C. That is much better as we will eventually need to vectorize this code, it is a kernel.
On Sat, Oct 22, 2016 at 10:24 AM, Mark Adams <[email protected]> wrote: > > > On Sat, Oct 22, 2016 at 8:07 AM, Matthew Knepley <[email protected]> > wrote: > >> Boost is an unmitigated nightmare. What exactly do you need? >> > > 7.13.3 Legendre Form of Complete Elliptic IntegralsFunction: *double* > *gsl_sf_ellint_Kcomp* *(double k, gsl_mode_t mode)*Function: *int* > *gsl_sf_ellint_Kcomp_e* *(double k, gsl_mode_t mode, gsl_sf_result > * result)* > > These routines compute the complete elliptic integral *K(k)* to the > accuracy specified by the mode variable mode. Note that Abramowitz & > Stegun define this function in terms of the parameter *m = k^2*. > Function: *double* *gsl_sf_ellint_Ecomp* *(double k, gsl_mode_t mode)* > Function: *int* *gsl_sf_ellint_Ecomp_e* *(double k, gsl_mode_t mode, > gsl_sf_result * result)* > > >> DLMF is a good reference: >> >> http://dlmf.nist.gov/19.39 >> > > Yes, that is where I got a pointer to GSL. THere are 5 other open source > options: > >
