Elliptic integrals for Rust

Ellip is a pure-Rust implementation of elliptic integrals. Ellip also provides less common functions like Bulirsch's cel and el. Some applications of the elliptic integrals include computing the lengths of plane curves, magnetism, astrophysics, and string theory.

Quick Start

Start by installing Ellip.

>> cargo add ellip

Let's compute the circumference of an ellipse.

use ellip::*;

fn ellipse_length(a: f64, b: f64) -> Result<f64, StrErr> {
    Ok(8.0 * elliprg(0.0, a * a, b * b)?)
}

let ans = ellipse_length(5.0, 3.0).unwrap();
ellip::util::assert_close(ans, 25.526998863398124, 1e-15);

Learn more at doc.rs.

Features

• Legendre's complete integrals
  • ellipk: Complete elliptic integral of the first kind (K).
  • ellipe: Complete elliptic integral of the second kind (E).
  • ellippi: Complete elliptic integral of the third kind (Π).
  • ellipd: Complete elliptic integral of Legendre's type (D).
• Legendre's incomplete integrals
  • ellipf: Incomplete elliptic integral of the first kind (F).
  • ellipeinc: Incomplete elliptic integral of the second kind (E).
  • ellippiinc: Incomplete elliptic integral of the third kind (Π).
  • ellipdinc: Incomplete elliptic integral of Legendre's type (D).
• Bulirsch's integrals
  • cel: General complete elliptic integral in Bulirsch's form.
  • cel1: Complete elliptic integral of the first kind in Bulirsch's form.
  • cel2: Complete elliptic integral of the second kind in Bulirsch's form.
  • el1: Incomplete elliptic integral of the first kind in Bulirsch's form.
  • el2: Incomplete elliptic integral of the second kind in Bulirsch's form.
  • el3: Incomplete elliptic integral of the third kind in Bulirsch's form.
• Carlson's symmetric integrals
  • elliprf: Symmetric elliptic integral of the first kind (RF).
  • elliprg: Symmetric elliptic integral of the second kind (RG).
  • elliprj: Symmetric elliptic integral of the third kind (RJ).
  • elliprc: Degenerate elliptic integral of RF (RC).
  • elliprd: Degenerate elliptic integral of the third kind (RD).
• Miscellaneous functions
  • jacobi_zeta: Jacobi Zeta function (Z).
  • heuman_lambda: Heuman Lambda function (Λ0).

Testing

In the unit tests, the functions are tested against the Boost Math and Wolfram test data. Since Ellip accepts the argument m (parameter) instead of k (modulus) to allow larger domain support, the full accuracy report uses exclusively the Wolfram data. The full accuracy report, test data, and test generation scripts can be found here. The performance benchmark is presented to provide comparison between functions in Ellip. Comparing performance with other libraries is non-trivial, since they accept different domains.

Benchmark on AMD Ryzen 5 4600H with Radeon Graphics running x86_64-unknown-linux-gnu rustc 1.89.0 using ellip v0.5.0 at f64 precision (ε≈2.22e-16).

Legendre's Elliptic Integrals

Bulirsch's Elliptic Integrals

Carlson's Symmetric Integrals

Miscellaneous Functions

Learn more at docs.rs.