[petsc-dev] elliptic functions

[Mark Adams]

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 <mfadams at lbl.gov> wrote:

>
>
> On Sat, Oct 22, 2016 at 8:07 AM, Matthew Knepley <knepley at gmail.com>
> 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:
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20161022/abfc94da/attachment.html>