Full-precision implementations of elliptic integrals, Jacobi's elliptic and theta functions, Carlson symmetric forms and Weierstrass's ℘ — for MATLAB / Octave and Python (NumPy · PyTorch · JAX). No external toolboxes, tensor support throughout, complex-argument variants included.
Trace Keplerian orbits, solve Kepler's equation, and watch the Carlson duplication algorithm converge.
Real JPL orbital elements on GPU. Throughput benchmarked across NumPy, PyTorch CUDA, and JAX.
Exact period via K(k²), live sn(t|m) animation, phase portrait with separatrix and period divergence.
If you've used any routines in this package — MATLAB, Octave or Python — please cite:
@misc{elliptic,
author = {Moiseev I.},
title = {Elliptic functions for Matlab, Octave and Python},
year = {2008},
publisher = {GitHub},
doi = {10.5281/zenodo.48264},
url = {http://dx.doi.org/10.5281/zenodo.48264}
}
Contributors: @moiseevigor (Igor Moiseev, maintainer) · @drbitboy (Brian Carcich) · @wspr (Will Robertson)