Crate gooding_lambert 

gooding-lambert

Gooding’s method for solving Lambert’s orbital boundary-value problem.

Given two position vectors and a time of flight, computes the velocity vectors at each endpoint of the connecting Keplerian arc. Handles all conic sections (elliptic, parabolic, hyperbolic), all transfer angles, and multi-revolution solutions.

Quick Start

use gooding_lambert::{lambert, Direction};

let mu = 398600.4418_f64; // Earth GM (km³/s²)
let r1 = [6678.0, 0.0, 0.0]; // km
let r2 = [0.0, 42164.0, 0.0]; // km (GEO)
let tof = 5.0 * 3600.0; // seconds

let sol = lambert(mu, r1, r2, tof, 0, Direction::Prograde).unwrap();
let speed = (sol.v1[0].powi(2) + sol.v1[1].powi(2) + sol.v1[2].powi(2)).sqrt();
assert!(speed > 5.0 && speed < 15.0); // km/s at LEO departure

References

Gooding, R. H. (1990). “A procedure for the solution of Lambert’s orbital boundary-value problem.” Celestial Mechanics and Dynamical Astronomy, 48(2), 145–165.

Structs

LambertSolution

Velocity solution returned by lambert.

Enums

Direction

Transfer direction for Lambert’s problem.

LambertError

Errors returned by lambert.

Functions

lambert

Solve Lambert’s problem using Gooding’s (1990) method.