Weierstrass P function and its derivative implemented in C by Lattès iteration (a variant of the Coquereaux-Grossmann-Lautrup method [1]).
- GNU Scientific Library (GSL), http://www.gnu.org/software/gsl/ (for complex number types, arithmetic, and complex trig functions)
- R. Coquereaux, A. Grossmann, and B. Lautrup. Iterative method for calculation of the Weierstrass elliptic function. IMA J. Numer. Anal. 10 (1990), no. 1, 119–128. http://dx.doi.org/10.1093/imanum/10.1.119