I have soon finished the first version of Jacobi.jl: Jacobi theta functions and related functions. I still have to implement the Weierstrass sigma-function and the Weierstrass zeta-function, and the first version will be done.

The package also includes a couple of elliptic integrals, because they are needed for other functions, and in the Elliptic package they are not available for complex values of the variable (I don’t know why because there’s a simple implementation).

Thanks. Are you sure that the Julia mod function is the same as “my” modulo function? Mine can take negative values, it mimics the Fortran function. Anyway I will test.