@interstellar-tools/equations 0.60.1 → 0.61.0

package/dist/categories/angle/compute-angle.d.ts

@@ -1,62 +1,51 @@
-import type { CelestialBodyType, Radians, TimeStepInterface } from '@interstellar-tools/types';
+import type { OrbitInterface, Radians, TimeStepInterface } from '@interstellar-tools/types';
 /**
- * Computes the orbital angle (true anomaly, $ν$) of a celestial body for a given date and time step.
+ * Advances a body along its orbit by one time step and returns its new
+ * **true anomaly** ($\nu$).
  *
- * **Mathematical Explanation:**
+ * The orbit's `angle` is the starting true anomaly. There is no date input:
+ * call this once per frame or step, and feed the result back in as the next
+ * `angle` to move the body over time.
  *
- * This function determines the **true anomaly** ($ν$) based on the time elapsed since the
- * **J2000 epoch** (January 1, 2000) and the provided time step. The calculation follows these steps:
+ * **Step 1: Mean anomaly after the time step ($M$)**
  *
- * **Step 1: Compute Time Elapsed Since J2000**
- *
- * The number of days since **J2000** is computed as:
- * $$ \Delta T = \frac{t - t_{J2000}}{\text{MILLISECONDS\_PER\_DAY}} $$
+ * Done by {@link advanceMeanAnomaly}:
+ * $$ M = M_0 + n \cdot \Delta t, \qquad n = \frac{2\pi}{|P|} $$
  * where:
- * - $t$ is the current date in milliseconds.
- * - $t_{J2000}$ is **J2000** (2000-01-01T00:00:00Z).
- * - **MILLISECONDS_PER_DAY** is the number of milliseconds in one day.
- *
- * **Step 2: Compute the Mean Anomaly ($M$)**
+ * - $M_0$ is the mean anomaly at the start, converted from the true anomaly
+ *   `angle` (equal to it when $e = 0$).
+ * - $n$ is the mean motion, from the period $P$ converted to days. The sign of
+ *   $P$ is ignored, so the angle always increases.
+ * - $\Delta t$ is `timeStep.value` in days, clamped to at least $10^{-5}$ days
+ *   and at most $10\,|P|$ ($|P|$ when $e \geq 0.9$).
  *
- * The **mean anomaly** is calculated as:
- * $$ M = M_0 + n \cdot (\Delta T + \text{timeStep}) $$
- * where:
- * - $M_0$ is the initial mean anomaly at **J2000**.
- * - $n$ is the mean motion (orbital angular velocity).
- * - $\text{timeStep}$ is the time step in days.
+ * The result is wrapped to $[0, 2\pi)$.
  *
- * **Step 3: Solve Kepler’s Equation for Eccentric Anomaly ($E$)**
+ * **Step 2: Solve Kepler's equation for the eccentric anomaly ($E$)**
  *
- * Kepler’s equation:
  * $$ M = E - e \sin(E) $$
- * is solved numerically to obtain the **eccentric anomaly** ($E$).
- *
- * **Step 4: Convert Eccentric Anomaly ($E$) to True Anomaly ($ν$)**
+ * is solved numerically by {@link solveKepler}.
  *
- * Using the relation:
- * $$ ν = 2 \tan^{-1} \left( \sqrt{\frac{1+e}{1-e}} \tan\left(\frac{E}{2}\right) \right) $$
+ * **Step 3: Convert the eccentric anomaly to the true anomaly ($\nu$)**
  *
- * **Step 5: Adjust the True Anomaly for Certain Bodies**
- * - For **comets**, the final **true anomaly** is adjusted by subtracting the **eccentricity** ($e$).
- * - For **stars and other celestial objects**, the computed **true anomaly** is used directly.
+ * Done by {@link eccentricToTrueAnomaly}:
+ * $$ \nu = 2 \tan^{-1} \left( \sqrt{\frac{1+e}{1-e}} \tan\left(\frac{E}{2}\right) \right) $$
  *
- * @param {CelestialBodyType} body - The celestial body for which to compute the true anomaly.
- * @param {TimeStepInterface} timeStep - The time step in days since the last frame.
- * @returns {Radians} The computed true anomaly in radians.
+ * @param {Pick<OrbitInterface, 'e' | 'angle' | 'period'>} orbit - The body's orbit, for example `body.orbit` from a dataset.
+ * @param {TimeStepInterface} timeStep - How far to advance, in days.
+ * @returns {Radians} The true anomaly after the time step, in radians.
+ * @throws {@link https://developer.mozilla.org/en-US/docs/Web/JavaScript/Reference/Global_Objects/RangeError | RangeError} If the orbit's eccentricity is outside $0 \leq e < 1$.
  *
  * @example
  * ```ts
+ * import { PLANETS } from '@interstellar-tools/constants';
  * import { computeAngle } from '@interstellar-tools/equations';
  *
- * const earth: CelestialBodyType = {
- *   name: 'Earth',
- *   e: 0.0167, // Orbital eccentricity
- *   angle: 0,  // Initial angle
- * };
+ * const earth = PLANETS.find((planet) => planet.id === 'earth')!;
+ * const timeStep: TimeStepInterface = { value: 1, unit: 'd' };
  *
- * const timeStep: TimeStepInterface = { days: 1 }; // Time step of 1 day
- * const angle: Radians = computeAngle(earth, timeStep);
- * console.log(angle); // Output: Computed true anomaly in radians
+ * // Earth's true anomaly one day after its dataset position
+ * const angle: Radians = computeAngle(earth.orbit!, timeStep);
  * ```
  *
  * @see [Kepler's Equation (Wikipedia)](https://en.wikipedia.org/wiki/Kepler%27s_equation)
@@ -65,5 +54,5 @@ import type { CelestialBodyType, Radians, TimeStepInterface } from '@interstella
  * @see [Orbital Mechanics (NASA)](https://solarsystem.nasa.gov/basics/chapter2-2/)
  * @group Angles
  */
-export declare const computeAngle: (body: CelestialBodyType, timeStep: TimeStepInterface) => Radians;
+export declare const computeAngle: (orbit: Pick<OrbitInterface, "e" | "angle" | "period">, timeStep: TimeStepInterface) => Radians;
 //# sourceMappingURL=compute-angle.d.ts.map

package/dist/categories/angle/compute-angle.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"compute-angle.d.ts","sourceRoot":"","sources":["../../../src/categories/angle/compute-angle.ts"],"names":[],"mappings":"AAAA,OAAO,KAAK,EACV,iBAAiB,EACjB,OAAO,EACP,iBAAiB,EAClB,MAAM,2BAA2B,CAAC;AAMnC;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GAiEG;AACH,eAAO,MAAM,YAAY,GACvB,MAAM,iBAAiB,EACvB,UAAU,iBAAiB,KAC1B,OASF,CAAC"}
+{"version":3,"file":"compute-angle.d.ts","sourceRoot":"","sources":["../../../src/categories/angle/compute-angle.ts"],"names":[],"mappings":"AAAA,OAAO,KAAK,EACV,cAAc,EACd,OAAO,EACP,iBAAiB,EAClB,MAAM,2BAA2B,CAAC;AAMnC;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GAsDG;AACH,eAAO,MAAM,YAAY,GACvB,OAAO,IAAI,CAAC,cAAc,EAAE,GAAG,GAAG,OAAO,GAAG,QAAQ,CAAC,EACrD,UAAU,iBAAiB,KAC1B,OASF,CAAC"}

package/dist/categories/angle/compute-angle.js

@@ -1,64 +1,53 @@
+import { advanceMeanAnomaly } from '../anomalies/advance-mean-anomaly.js';
 import { eccentricToTrueAnomaly } from '../anomalies/eccentric-to-true-anomaly.js';
-import { meanToEccentricAnomaly } from '../anomalies/mean-to-eccentric-anomaly.js';
 import { solveKepler } from '../kepler/solve-kepler.js';
 /**
- * Computes the orbital angle (true anomaly, $ν$) of a celestial body for a given date and time step.
+ * Advances a body along its orbit by one time step and returns its new
+ * **true anomaly** ($\nu$).
  *
- * **Mathematical Explanation:**
+ * The orbit's `angle` is the starting true anomaly. There is no date input:
+ * call this once per frame or step, and feed the result back in as the next
+ * `angle` to move the body over time.
  *
- * This function determines the **true anomaly** ($ν$) based on the time elapsed since the
- * **J2000 epoch** (January 1, 2000) and the provided time step. The calculation follows these steps:
+ * **Step 1: Mean anomaly after the time step ($M$)**
  *
- * **Step 1: Compute Time Elapsed Since J2000**
- *
- * The number of days since **J2000** is computed as:
- * $$ \Delta T = \frac{t - t_{J2000}}{\text{MILLISECONDS\_PER\_DAY}} $$
+ * Done by {@link advanceMeanAnomaly}:
+ * $$ M = M_0 + n \cdot \Delta t, \qquad n = \frac{2\pi}{|P|} $$
  * where:
- * - $t$ is the current date in milliseconds.
- * - $t_{J2000}$ is **J2000** (2000-01-01T00:00:00Z).
- * - **MILLISECONDS_PER_DAY** is the number of milliseconds in one day.
- *
- * **Step 2: Compute the Mean Anomaly ($M$)**
+ * - $M_0$ is the mean anomaly at the start, converted from the true anomaly
+ *   `angle` (equal to it when $e = 0$).
+ * - $n$ is the mean motion, from the period $P$ converted to days. The sign of
+ *   $P$ is ignored, so the angle always increases.
+ * - $\Delta t$ is `timeStep.value` in days, clamped to at least $10^{-5}$ days
+ *   and at most $10\,|P|$ ($|P|$ when $e \geq 0.9$).
  *
- * The **mean anomaly** is calculated as:
- * $$ M = M_0 + n \cdot (\Delta T + \text{timeStep}) $$
- * where:
- * - $M_0$ is the initial mean anomaly at **J2000**.
- * - $n$ is the mean motion (orbital angular velocity).
- * - $\text{timeStep}$ is the time step in days.
+ * The result is wrapped to $[0, 2\pi)$.
  *
- * **Step 3: Solve Kepler’s Equation for Eccentric Anomaly ($E$)**
+ * **Step 2: Solve Kepler's equation for the eccentric anomaly ($E$)**
  *
- * Kepler’s equation:
  * $$ M = E - e \sin(E) $$
- * is solved numerically to obtain the **eccentric anomaly** ($E$).
- *
- * **Step 4: Convert Eccentric Anomaly ($E$) to True Anomaly ($ν$)**
+ * is solved numerically by {@link solveKepler}.
  *
- * Using the relation:
- * $$ ν = 2 \tan^{-1} \left( \sqrt{\frac{1+e}{1-e}} \tan\left(\frac{E}{2}\right) \right) $$
+ * **Step 3: Convert the eccentric anomaly to the true anomaly ($\nu$)**
  *
- * **Step 5: Adjust the True Anomaly for Certain Bodies**
- * - For **comets**, the final **true anomaly** is adjusted by subtracting the **eccentricity** ($e$).
- * - For **stars and other celestial objects**, the computed **true anomaly** is used directly.
+ * Done by {@link eccentricToTrueAnomaly}:
+ * $$ \nu = 2 \tan^{-1} \left( \sqrt{\frac{1+e}{1-e}} \tan\left(\frac{E}{2}\right) \right) $$
  *
- * @param {CelestialBodyType} body - The celestial body for which to compute the true anomaly.
- * @param {TimeStepInterface} timeStep - The time step in days since the last frame.
- * @returns {Radians} The computed true anomaly in radians.
+ * @param {Pick<OrbitInterface, 'e' | 'angle' | 'period'>} orbit - The body's orbit, for example `body.orbit` from a dataset.
+ * @param {TimeStepInterface} timeStep - How far to advance, in days.
+ * @returns {Radians} The true anomaly after the time step, in radians.
+ * @throws {@link https://developer.mozilla.org/en-US/docs/Web/JavaScript/Reference/Global_Objects/RangeError | RangeError} If the orbit's eccentricity is outside $0 \leq e < 1$.
  *
  * @example
  * ```ts
+ * import { PLANETS } from '@interstellar-tools/constants';
  * import { computeAngle } from '@interstellar-tools/equations';
  *
- * const earth: CelestialBodyType = {
- *   name: 'Earth',
- *   e: 0.0167, // Orbital eccentricity
- *   angle: 0,  // Initial angle
- * };
+ * const earth = PLANETS.find((planet) => planet.id === 'earth')!;
+ * const timeStep: TimeStepInterface = { value: 1, unit: 'd' };
  *
- * const timeStep: TimeStepInterface = { days: 1 }; // Time step of 1 day
- * const angle: Radians = computeAngle(earth, timeStep);
- * console.log(angle); // Output: Computed true anomaly in radians
+ * // Earth's true anomaly one day after its dataset position
+ * const angle: Radians = computeAngle(earth.orbit!, timeStep);
  * ```
  *
  * @see [Kepler's Equation (Wikipedia)](https://en.wikipedia.org/wiki/Kepler%27s_equation)
@@ -67,13 +56,13 @@ import { solveKepler } from '../kepler/solve-kepler.js';
  * @see [Orbital Mechanics (NASA)](https://solarsystem.nasa.gov/basics/chapter2-2/)
  * @group Angles
  */
-export const computeAngle = (body, timeStep) => {
+export const computeAngle = (orbit, timeStep) => {
     // Compute mean anomaly (M) with time step
-    const M = meanToEccentricAnomaly(body, timeStep);
+    const M = advanceMeanAnomaly(orbit, timeStep);
     // Solve Kepler's equation for eccentric anomaly (E)
-    const E = solveKepler(M, body.e);
+    const E = solveKepler(M, orbit.e);
     // Convert eccentric anomaly (E) to true anomaly (ν)
-    const V = eccentricToTrueAnomaly(E, body.e);
+    const V = eccentricToTrueAnomaly(E, orbit.e);
     return V;
 };
 //# sourceMappingURL=compute-angle.js.map

package/dist/categories/angle/compute-angle.js.map

@@ -1 +1 @@
-{"version":3,"file":"compute-angle.js","sourceRoot":"","sources":["../../../src/categories/angle/compute-angle.ts"],"names":[],"mappings":"AAMA,OAAO,EAAE,sBAAsB,EAAE,MAAM,2CAA2C,CAAC;AACnF,OAAO,EAAE,sBAAsB,EAAE,MAAM,2CAA2C,CAAC;AACnF,OAAO,EAAE,WAAW,EAAE,MAAM,2BAA2B,CAAC;AAExD;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GAiEG;AACH,MAAM,CAAC,MAAM,YAAY,GAAG,CAC1B,IAAuB,EACvB,QAA2B,EAClB,EAAE;IACX,0CAA0C;IAC1C,MAAM,CAAC,GAAY,sBAAsB,CAAC,IAAI,EAAE,QAAQ,CAAC,CAAC;IAC1D,oDAAoD;IACpD,MAAM,CAAC,GAAY,WAAW,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,CAAC,CAAC;IAC1C,oDAAoD;IACpD,MAAM,CAAC,GAAY,sBAAsB,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,CAAC,CAAC;IAErD,OAAO,CAAC,CAAC;AACX,CAAC,CAAC"}
+{"version":3,"file":"compute-angle.js","sourceRoot":"","sources":["../../../src/categories/angle/compute-angle.ts"],"names":[],"mappings":"AAMA,OAAO,EAAE,kBAAkB,EAAE,MAAM,sCAAsC,CAAC;AAC1E,OAAO,EAAE,sBAAsB,EAAE,MAAM,2CAA2C,CAAC;AACnF,OAAO,EAAE,WAAW,EAAE,MAAM,2BAA2B,CAAC;AAExD;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GAsDG;AACH,MAAM,CAAC,MAAM,YAAY,GAAG,CAC1B,KAAqD,EACrD,QAA2B,EAClB,EAAE;IACX,0CAA0C;IAC1C,MAAM,CAAC,GAAY,kBAAkB,CAAC,KAAK,EAAE,QAAQ,CAAC,CAAC;IACvD,oDAAoD;IACpD,MAAM,CAAC,GAAY,WAAW,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,CAAC;IAC3C,oDAAoD;IACpD,MAAM,CAAC,GAAY,sBAAsB,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,CAAC;IAEtD,OAAO,CAAC,CAAC;AACX,CAAC,CAAC"}

package/dist/categories/anomalies/advance-mean-anomaly.d.ts

@@ -0,0 +1,50 @@
+import type { OrbitInterface, Radians, TemporalInterface } from '@interstellar-tools/types';
+/**
+ * Advances an orbit by a time step and returns its new **mean anomaly** ($M$).
+ *
+ * The orbit's `angle` (true anomaly) is the starting point. To get the eccentric
+ * anomaly from the result, pass it to {@link solveKepler}.
+ *
+ * **Mathematical Explanation:**
+ *
+ * The **mean anomaly** is a measure of the position of an orbiting body as if it moved
+ * with uniform angular motion along a circular orbit with the same period as the actual
+ * elliptical orbit. It is computed as:
+ *
+ * $$ M = M_0 + n \cdot \Delta T $$
+ *
+ * where:
+ * - $M_0$ is the **initial mean anomaly** (converted from true anomaly if necessary).
+ * - $n$ is the **mean motion**, given by:
+ *   $$ n = \frac{2\pi}{P} $$
+ *   where $P$ is the orbital period in days.
+ * - $\Delta T$ is the time step, converted to days from its unit.
+ *
+ * @param {Pick<OrbitInterface, 'e' | 'angle' | 'period'>} orbit - The body's orbit, for example `body.orbit` from a dataset.
+ * @param {TemporalInterface} timeStep - The time step over which to compute the change, in any temporal unit.
+ * @returns {Radians} The computed mean anomaly in radians.
+ *
+ * @throws {@link https://developer.mozilla.org/en-US/docs/Web/JavaScript/Reference/Global_Objects/RangeError | RangeError} If the orbit's eccentricity is outside the range $0 \leq e < 1$.
+ *
+ * @example
+ * ```ts
+ * import { advanceMeanAnomaly } from '@interstellar-tools/equations';
+ *
+ * const marsOrbit = {
+ *   e: 0.0934, // Eccentricity of Mars' orbit
+ *   angle: 1.047 as Radians, // Initial true anomaly (in radians)
+ *   period: { value: 687, unit: 'd' } // Orbital period in days
+ * };
+ *
+ * const timeStep: TemporalInterface = { value: 1, unit: 'd' }; // 1-day step
+ * const meanAnomaly = advanceMeanAnomaly(marsOrbit, timeStep);
+ * console.log(meanAnomaly); // Output: Computed mean anomaly in radians
+ * ```
+ *
+ * @see {@link solveKepler} to convert the mean anomaly to the eccentric anomaly.
+ * @see [Mean Anomaly (Wikipedia)](https://en.wikipedia.org/wiki/Mean_anomaly)
+ * @see [Orbital Mechanics (NASA)](https://solarsystem.nasa.gov/basics/chapter2-2/)
+ * @group Anomalies
+ */
+export declare const advanceMeanAnomaly: (orbit: Pick<OrbitInterface, "e" | "angle" | "period">, timeStep: TemporalInterface) => Radians;
+//# sourceMappingURL=advance-mean-anomaly.d.ts.map

package/dist/categories/anomalies/advance-mean-anomaly.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"advance-mean-anomaly.d.ts","sourceRoot":"","sources":["../../../src/categories/anomalies/advance-mean-anomaly.ts"],"names":[],"mappings":"AAEA,OAAO,KAAK,EACV,cAAc,EACd,OAAO,EACP,iBAAiB,EAClB,MAAM,2BAA2B,CAAC;AAQnC;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GA8CG;AAEH,eAAO,MAAM,kBAAkB,GAC7B,OAAO,IAAI,CAAC,cAAc,EAAE,GAAG,GAAG,OAAO,GAAG,QAAQ,CAAC,EACrD,UAAU,iBAAiB,KAC1B,OAqCF,CAAC"}

package/dist/categories/anomalies/{mean-to-eccentric-anomaly.js → advance-mean-anomaly.js}

@@ -5,7 +5,10 @@ import { trueToMeanAnomaly } from './true-to-mean-anomaly.js';
 import { areEqual } from './utils/are-equal.js';
 /* node:coverage disable */
 /**
- * Computes the **mean anomaly** ($M$) of a celestial body for a given time step.
+ * Advances an orbit by a time step and returns its new **mean anomaly** ($M$).
+ *
+ * The orbit's `angle` (true anomaly) is the starting point. To get the eccentric
+ * anomaly from the result, pass it to {@link solveKepler}.
  *
  * **Mathematical Explanation:**
  *
@@ -20,52 +23,53 @@ import { areEqual } from './utils/are-equal.js';
  * - $n$ is the **mean motion**, given by:
  *   $$ n = \frac{2\pi}{P} $$
  *   where $P$ is the orbital period in days.
- * - $\Delta T$ is the time step.
+ * - $\Delta T$ is the time step, converted to days from its unit.
  *
- * @param {CelestialBodyType} body - The celestial body for which to compute the mean anomaly.
- * @param {TemporalInterface} timeStep - The time step over which to compute the change.
+ * @param {Pick<OrbitInterface, 'e' | 'angle' | 'period'>} orbit - The body's orbit, for example `body.orbit` from a dataset.
+ * @param {TemporalInterface} timeStep - The time step over which to compute the change, in any temporal unit.
  * @returns {Radians} The computed mean anomaly in radians.
  *
- * @throws {@link https://developer.mozilla.org/en-US/docs/Web/JavaScript/Reference/Global_Objects/RangeError | RangeError} If the body's eccentricity is outside the range $0 \leq e < 1$.
+ * @throws {@link https://developer.mozilla.org/en-US/docs/Web/JavaScript/Reference/Global_Objects/RangeError | RangeError} If the orbit's eccentricity is outside the range $0 \leq e < 1$.
  *
  * @example
  * ```ts
- * import { meanToEccentricAnomaly } from '@interstellar-tools/equations';
+ * import { advanceMeanAnomaly } from '@interstellar-tools/equations';
  *
- * const mars: CelestialBodyType = {
- *   name: 'Mars',
+ * const marsOrbit = {
  *   e: 0.0934, // Eccentricity of Mars' orbit
- *   angle: 1.047, // Initial true anomaly (in radians)
- *   period: { value: 687, unit: 'day' }, // Orbital period in days
+ *   angle: 1.047 as Radians, // Initial true anomaly (in radians)
+ *   period: { value: 687, unit: 'd' } // Orbital period in days
  * };
  *
- * const timeStep: TemporalInterface = { value: 1, unit: 'day' }; // 1-day step
- * const meanAnomaly = meanToEccentricAnomaly(mars, timeStep);
+ * const timeStep: TemporalInterface = { value: 1, unit: 'd' }; // 1-day step
+ * const meanAnomaly = advanceMeanAnomaly(marsOrbit, timeStep);
  * console.log(meanAnomaly); // Output: Computed mean anomaly in radians
  * ```
  *
+ * @see {@link solveKepler} to convert the mean anomaly to the eccentric anomaly.
  * @see [Mean Anomaly (Wikipedia)](https://en.wikipedia.org/wiki/Mean_anomaly)
  * @see [Orbital Mechanics (NASA)](https://solarsystem.nasa.gov/basics/chapter2-2/)
  * @group Anomalies
  */
 /* node:coverage enable */
-export const meanToEccentricAnomaly = (body, timeStep) => {
-    const { e, angle, period } = body;
+export const advanceMeanAnomaly = (orbit, timeStep) => {
+    const { e, angle, period } = orbit;
     if (e < 0 || e >= 1) {
         throw new RangeError(`Invalid eccentricity: ${e}. Eccentricity must be in the range [0, 1].`);
     }
-    // Convert period to days and compute mean motion
+    // Convert period and time step to days and compute mean motion
     const periodInDays = convertTemporalUnit(period, 'd').value;
+    const timeStepInDays = convertTemporalUnit(timeStep, 'd').value;
     const meanMotion = TWO_PI / Math.abs(periodInDays);
     // Convert true anomaly (V) to mean anomaly (M0) if necessary
     const M0 = e === 0 ? angle : trueToMeanAnomaly(angle, e);
-    if (timeStep.value === 0) {
+    if (timeStepInDays === 0) {
         return M0;
     }
     // Clamping small time steps to prevent unnecessary updates
     const minTimeStep = 1e-5; // Minimum allowed time step
     const maxTimeStepFactor = e >= 0.9 ? 1 : 10; // Stricter limit for highly eccentric orbits
-    const clampedDeltaT = Math.max(minTimeStep, Math.min(timeStep.value, maxTimeStepFactor * Math.abs(periodInDays)));
+    const clampedDeltaT = Math.max(minTimeStep, Math.min(timeStepInDays, maxTimeStepFactor * Math.abs(periodInDays)));
     // Compute new mean anomaly with hysteresis filter and sign preservation
     const prevM = M0;
     let newM = wrapAngle(M0 + meanMotion * clampedDeltaT);
@@ -74,4 +78,4 @@ export const meanToEccentricAnomaly = (body, timeStep) => {
     }
     return newM;
 };
-//# sourceMappingURL=mean-to-eccentric-anomaly.js.map
+//# sourceMappingURL=advance-mean-anomaly.js.map

package/dist/categories/anomalies/advance-mean-anomaly.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"advance-mean-anomaly.js","sourceRoot":"","sources":["../../../src/categories/anomalies/advance-mean-anomaly.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,MAAM,EAAE,MAAM,+BAA+B,CAAC;AACvD,OAAO,EAAE,mBAAmB,EAAE,MAAM,8BAA8B,CAAC;AAOnE,OAAO,EAAE,SAAS,EAAE,MAAM,wBAAwB,CAAC;AAEnD,OAAO,EAAE,iBAAiB,EAAE,MAAM,2BAA2B,CAAC;AAC9D,OAAO,EAAE,QAAQ,EAAE,MAAM,sBAAsB,CAAC;AAEhD,2BAA2B;AAC3B;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GA8CG;AACH,0BAA0B;AAC1B,MAAM,CAAC,MAAM,kBAAkB,GAAG,CAChC,KAAqD,EACrD,QAA2B,EAClB,EAAE;IACX,MAAM,EAAE,CAAC,EAAE,KAAK,EAAE,MAAM,EAAE,GAAG,KAAK,CAAC;IAEnC,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,IAAI,CAAC,EAAE,CAAC;QACpB,MAAM,IAAI,UAAU,CAClB,yBAAyB,CAAC,6CAA6C,CACxE,CAAC;IACJ,CAAC;IAED,+DAA+D;IAC/D,MAAM,YAAY,GAAG,mBAAmB,CAAC,MAAM,EAAE,GAAG,CAAC,CAAC,KAAK,CAAC;IAC5D,MAAM,cAAc,GAAG,mBAAmB,CAAC,QAAQ,EAAE,GAAG,CAAC,CAAC,KAAK,CAAC;IAChE,MAAM,UAAU,GAAG,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,YAAY,CAAC,CAAC;IACnD,6DAA6D;IAC7D,MAAM,EAAE,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,iBAAiB,CAAC,KAAK,EAAE,CAAC,CAAC,CAAC;IAEzD,IAAI,cAAc,KAAK,CAAC,EAAE,CAAC;QACzB,OAAO,EAAa,CAAC;IACvB,CAAC;IAED,2DAA2D;IAC3D,MAAM,WAAW,GAAG,IAAI,CAAC,CAAC,4BAA4B;IACtD,MAAM,iBAAiB,GAAG,CAAC,IAAI,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,6CAA6C;IAC1F,MAAM,aAAa,GAAG,IAAI,CAAC,GAAG,CAC5B,WAAW,EACX,IAAI,CAAC,GAAG,CAAC,cAAc,EAAE,iBAAiB,GAAG,IAAI,CAAC,GAAG,CAAC,YAAY,CAAC,CAAC,CACrE,CAAC;IACF,wEAAwE;IACxE,MAAM,KAAK,GAAG,EAAE,CAAC;IAEjB,IAAI,IAAI,GAAG,SAAS,CAAC,EAAE,GAAG,UAAU,GAAG,aAAa,CAAC,CAAC;IAEtD,IAAI,QAAQ,CAAC,IAAI,EAAE,KAAK,CAAC,EAAE,CAAC;QAC1B,IAAI,GAAG,KAAK,CAAC;IACf,CAAC;IAED,OAAO,IAAe,CAAC;AACzB,CAAC,CAAC"}

package/dist/categories/kepler/solve-kepler.d.ts

@@ -6,6 +6,9 @@ import type { Radians } from '@interstellar-tools/types';
  * - **High-eccentricity solver** for extreme orbits ($e > 0.9$).
  * - **Bisection fallback** if the selected solver doesn't return a root.
  *
+ * This is the conversion from **mean anomaly** to **eccentric anomaly**. Use
+ * {@link eccentricToTrueAnomaly} on the result to get the true anomaly.
+ *
  * ---
  *
  * **Mathematical Explanation:**

package/dist/categories/kepler/solve-kepler.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"solve-kepler.d.ts","sourceRoot":"","sources":["../../../src/categories/kepler/solve-kepler.ts"],"names":[],"mappings":"AACA,OAAO,KAAK,EAAE,OAAO,EAAE,MAAM,2BAA2B,CAAC;AAQzD;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GAoFG;AACH,eAAO,MAAM,WAAW,GACtB,GAAG,OAAO,EACV,GAAG,MAAM,EACT,gBAAY,EACZ,kBAAgB,KACf,OAuBF,CAAC"}
+{"version":3,"file":"solve-kepler.d.ts","sourceRoot":"","sources":["../../../src/categories/kepler/solve-kepler.ts"],"names":[],"mappings":"AACA,OAAO,KAAK,EAAE,OAAO,EAAE,MAAM,2BAA2B,CAAC;AAQzD;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GAuFG;AACH,eAAO,MAAM,WAAW,GACtB,GAAG,OAAO,EACV,GAAG,MAAM,EACT,gBAAY,EACZ,kBAAgB,KACf,OAuBF,CAAC"}

package/dist/categories/kepler/solve-kepler.js

@@ -10,6 +10,9 @@ import { solveKeplerNewtonRaphson } from './solve-kepler-newton-raphson.js';
  * - **High-eccentricity solver** for extreme orbits ($e > 0.9$).
  * - **Bisection fallback** if the selected solver doesn't return a root.
  *
+ * This is the conversion from **mean anomaly** to **eccentric anomaly**. Use
+ * {@link eccentricToTrueAnomaly} on the result to get the true anomaly.
+ *
  * ---
  *
  * **Mathematical Explanation:**

package/dist/categories/kepler/solve-kepler.js.map

@@ -1 +1 @@
-{"version":3,"file":"solve-kepler.js","sourceRoot":"","sources":["../../../src/categories/kepler/solve-kepler.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,MAAM,EAAE,MAAM,+BAA+B,CAAC;AAGvD,OAAO,EAAE,OAAO,EAAE,MAAM,oBAAoB,CAAC;AAE7C,OAAO,EAAE,oBAAoB,EAAE,MAAM,6BAA6B,CAAC;AACnE,OAAO,EAAE,2BAA2B,EAAE,MAAM,qCAAqC,CAAC;AAClF,OAAO,EAAE,wBAAwB,EAAE,MAAM,kCAAkC,CAAC;AAE5E;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GAoFG;AACH,MAAM,CAAC,MAAM,WAAW,GAAG,CACzB,CAAU,EACV,CAAS,EACT,OAAO,GAAG,EAAE,EACZ,SAAS,GAAG,IAAI,EACP,EAAE;IACX,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,IAAI,CAAC,EAAE,CAAC;QACpB,MAAM,IAAI,UAAU,CAAC,yBAAyB,CAAC,2BAA2B,CAAC,CAAC;IAC9E,CAAC;IAED,8DAA8D;IAC9D,MAAM,EAAE,GAAG,OAAO,CAAC,CAAC,CAAY,CAAC;IAEjC,wDAAwD;IACxD,IAAI,CAAC,GACH,CAAC,GAAG,GAAG;QACL,CAAC,CAAC,2BAA2B,CAAC,EAAE,EAAE,CAAC,EAAE,OAAO,EAAE,SAAS,CAAC;QACxD,CAAC,CAAC,wBAAwB,CAAC,EAAE,EAAE,CAAC,EAAE,OAAO,EAAE,SAAS,CAAC,CAAC;IAE1D,6EAA6E;IAC7E,sEAAsE;IACtE,MAAM,QAAQ,GAAG,OAAO,CAAC,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC;IAEnD,IAAI,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,QAAQ,EAAE,MAAM,GAAG,QAAQ,CAAC,IAAI,SAAS,CAAC,EAAE,CAAC;QAC1D,CAAC,GAAG,oBAAoB,CAAC,EAAE,EAAE,CAAC,EAAE,OAAO,EAAE,SAAS,CAAC,CAAC;IACtD,CAAC;IAED,OAAO,OAAO,CAAC,CAAC,CAAY,CAAC;AAC/B,CAAC,CAAC"}
+{"version":3,"file":"solve-kepler.js","sourceRoot":"","sources":["../../../src/categories/kepler/solve-kepler.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,MAAM,EAAE,MAAM,+BAA+B,CAAC;AAGvD,OAAO,EAAE,OAAO,EAAE,MAAM,oBAAoB,CAAC;AAE7C,OAAO,EAAE,oBAAoB,EAAE,MAAM,6BAA6B,CAAC;AACnE,OAAO,EAAE,2BAA2B,EAAE,MAAM,qCAAqC,CAAC;AAClF,OAAO,EAAE,wBAAwB,EAAE,MAAM,kCAAkC,CAAC;AAE5E;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GAuFG;AACH,MAAM,CAAC,MAAM,WAAW,GAAG,CACzB,CAAU,EACV,CAAS,EACT,OAAO,GAAG,EAAE,EACZ,SAAS,GAAG,IAAI,EACP,EAAE;IACX,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,IAAI,CAAC,EAAE,CAAC;QACpB,MAAM,IAAI,UAAU,CAAC,yBAAyB,CAAC,2BAA2B,CAAC,CAAC;IAC9E,CAAC;IAED,8DAA8D;IAC9D,MAAM,EAAE,GAAG,OAAO,CAAC,CAAC,CAAY,CAAC;IAEjC,wDAAwD;IACxD,IAAI,CAAC,GACH,CAAC,GAAG,GAAG;QACL,CAAC,CAAC,2BAA2B,CAAC,EAAE,EAAE,CAAC,EAAE,OAAO,EAAE,SAAS,CAAC;QACxD,CAAC,CAAC,wBAAwB,CAAC,EAAE,EAAE,CAAC,EAAE,OAAO,EAAE,SAAS,CAAC,CAAC;IAE1D,6EAA6E;IAC7E,sEAAsE;IACtE,MAAM,QAAQ,GAAG,OAAO,CAAC,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC;IAEnD,IAAI,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,QAAQ,EAAE,MAAM,GAAG,QAAQ,CAAC,IAAI,SAAS,CAAC,EAAE,CAAC;QAC1D,CAAC,GAAG,oBAAoB,CAAC,EAAE,EAAE,CAAC,EAAE,OAAO,EAAE,SAAS,CAAC,CAAC;IACtD,CAAC;IAED,OAAO,OAAO,CAAC,CAAC,CAAY,CAAC;AAC/B,CAAC,CAAC"}

package/dist/index.d.ts

@@ -19,8 +19,8 @@ export * from './categories/helpers/radians.js';
 export * from './categories/helpers/rot-1.js';
 export * from './categories/helpers/rot-3.js';
 export * from './categories/helpers/transpose-3.js';
+export * from './categories/anomalies/advance-mean-anomaly.js';
 export * from './categories/anomalies/eccentric-to-true-anomaly.js';
-export * from './categories/anomalies/mean-to-eccentric-anomaly.js';
 export * from './categories/anomalies/true-to-mean-anomaly.js';
 export * from './categories/kepler/solve-kepler-bisection.js';
 export * from './categories/kepler/solve-kepler-high-eccentricity.js';

package/dist/index.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"index.d.ts","sourceRoot":"","sources":["../src/index.ts"],"names":[],"mappings":"AAAA;;;;GAIG;AAGH,cAAc,qCAAqC,CAAC;AAEpD,cAAc,kCAAkC,CAAC;AAIjD,cAAc,8DAA8D,CAAC;AAE7E,cAAc,kEAAkE,CAAC;AAEjF,cAAc,wDAAwD,CAAC;AAEvE,cAAc,8DAA8D,CAAC;AAE7E,cAAc,qDAAqD,CAAC;AAEpE,cAAc,4DAA4D,CAAC;AAI3E,cAAc,wCAAwC,CAAC;AAEvD,cAAc,+BAA+B,CAAC;AAE9C,cAAc,kCAAkC,CAAC;AAEjD,cAAc,8BAA8B,CAAC;AAE7C,cAAc,iCAAiC,CAAC;AAEhD,cAAc,+BAA+B,CAAC;AAE9C,cAAc,+BAA+B,CAAC;AAE9C,cAAc,qCAAqC,CAAC;AAIpD,cAAc,qDAAqD,CAAC;AAEpE,cAAc,qDAAqD,CAAC;AAEpE,cAAc,gDAAgD,CAAC;AAI/D,cAAc,+CAA+C,CAAC;AAE9D,cAAc,uDAAuD,CAAC;AAEtE,cAAc,oDAAoD,CAAC;AAEnE,cAAc,qCAAqC,CAAC;AAIpD,cAAc,8CAA8C,CAAC;AAE7D,cAAc,uCAAuC,CAAC;AAEtD,cAAc,6CAA6C,CAAC;AAE5D,cAAc,iDAAiD,CAAC;AAIhE,cAAc,sDAAsD,CAAC;AAErE,cAAc,iDAAiD,CAAC;AAEhE,cAAc,uCAAuC,CAAC;AAEtD,cAAc,4CAA4C,CAAC;AAE3D,cAAc,qCAAqC,CAAC;AAEpD,cAAc,4DAA4D,CAAC;AAE3E,cAAc,mDAAmD,CAAC;AAElE,cAAc,oDAAoD,CAAC;AAEnE,cAAc,iDAAiD,CAAC;AAEhE,cAAc,sCAAsC,CAAC;AAErD,cAAc,oCAAoC,CAAC;AAEnD,cAAc,4CAA4C,CAAC;AAE3D,cAAc,gEAAgE,CAAC;AAE/E,cAAc,kDAAkD,CAAC;AAEjE,cAAc,mDAAmD,CAAC;AAElE,cAAc,mDAAmD,CAAC;AAElE,cAAc,uCAAuC,CAAC;AAItD,cAAc,kDAAkD,CAAC;AAEjE,cAAc,yDAAyD,CAAC;AAExE,cAAc,6CAA6C,CAAC;AAE5D,cAAc,+CAA+C,CAAC;AAE9D,cAAc,iDAAiD,CAAC;AAEhE,cAAc,oDAAoD,CAAC;AAEnE,cAAc,mDAAmD,CAAC"}
+{"version":3,"file":"index.d.ts","sourceRoot":"","sources":["../src/index.ts"],"names":[],"mappings":"AAAA;;;;GAIG;AAGH,cAAc,qCAAqC,CAAC;AAEpD,cAAc,kCAAkC,CAAC;AAIjD,cAAc,8DAA8D,CAAC;AAE7E,cAAc,kEAAkE,CAAC;AAEjF,cAAc,wDAAwD,CAAC;AAEvE,cAAc,8DAA8D,CAAC;AAE7E,cAAc,qDAAqD,CAAC;AAEpE,cAAc,4DAA4D,CAAC;AAI3E,cAAc,wCAAwC,CAAC;AAEvD,cAAc,+BAA+B,CAAC;AAE9C,cAAc,kCAAkC,CAAC;AAEjD,cAAc,8BAA8B,CAAC;AAE7C,cAAc,iCAAiC,CAAC;AAEhD,cAAc,+BAA+B,CAAC;AAE9C,cAAc,+BAA+B,CAAC;AAE9C,cAAc,qCAAqC,CAAC;AAIpD,cAAc,gDAAgD,CAAC;AAE/D,cAAc,qDAAqD,CAAC;AAEpE,cAAc,gDAAgD,CAAC;AAI/D,cAAc,+CAA+C,CAAC;AAE9D,cAAc,uDAAuD,CAAC;AAEtE,cAAc,oDAAoD,CAAC;AAEnE,cAAc,qCAAqC,CAAC;AAIpD,cAAc,8CAA8C,CAAC;AAE7D,cAAc,uCAAuC,CAAC;AAEtD,cAAc,6CAA6C,CAAC;AAE5D,cAAc,iDAAiD,CAAC;AAIhE,cAAc,sDAAsD,CAAC;AAErE,cAAc,iDAAiD,CAAC;AAEhE,cAAc,uCAAuC,CAAC;AAEtD,cAAc,4CAA4C,CAAC;AAE3D,cAAc,qCAAqC,CAAC;AAEpD,cAAc,4DAA4D,CAAC;AAE3E,cAAc,mDAAmD,CAAC;AAElE,cAAc,oDAAoD,CAAC;AAEnE,cAAc,iDAAiD,CAAC;AAEhE,cAAc,sCAAsC,CAAC;AAErD,cAAc,oCAAoC,CAAC;AAEnD,cAAc,4CAA4C,CAAC;AAE3D,cAAc,gEAAgE,CAAC;AAE/E,cAAc,kDAAkD,CAAC;AAEjE,cAAc,mDAAmD,CAAC;AAElE,cAAc,mDAAmD,CAAC;AAElE,cAAc,uCAAuC,CAAC;AAItD,cAAc,kDAAkD,CAAC;AAEjE,cAAc,yDAAyD,CAAC;AAExE,cAAc,6CAA6C,CAAC;AAE5D,cAAc,+CAA+C,CAAC;AAE9D,cAAc,iDAAiD,CAAC;AAEhE,cAAc,oDAAoD,CAAC;AAEnE,cAAc,mDAAmD,CAAC"}

package/dist/index.js

@@ -23,8 +23,8 @@ export * from './categories/helpers/rot-1.js';
 export * from './categories/helpers/rot-3.js';
 export * from './categories/helpers/transpose-3.js';
 // categories/anomalies
+export * from './categories/anomalies/advance-mean-anomaly.js';
 export * from './categories/anomalies/eccentric-to-true-anomaly.js';
-export * from './categories/anomalies/mean-to-eccentric-anomaly.js';
 export * from './categories/anomalies/true-to-mean-anomaly.js';
 // categories/kepler (solvers)
 export * from './categories/kepler/solve-kepler-bisection.js';

package/dist/index.js.map

@@ -1 +1 @@
-{"version":3,"file":"index.js","sourceRoot":"","sources":["../src/index.ts"],"names":[],"mappings":"AAAA;;;;GAIG;AAEH,mBAAmB;AACnB,cAAc,qCAAqC,CAAC;AAEpD,cAAc,kCAAkC,CAAC;AAEjD,yBAAyB;AAEzB,cAAc,8DAA8D,CAAC;AAE7E,cAAc,kEAAkE,CAAC;AAEjF,cAAc,wDAAwD,CAAC;AAEvE,cAAc,8DAA8D,CAAC;AAE7E,cAAc,qDAAqD,CAAC;AAEpE,cAAc,4DAA4D,CAAC;AAE3E,qBAAqB;AAErB,cAAc,wCAAwC,CAAC;AAEvD,cAAc,+BAA+B,CAAC;AAE9C,cAAc,kCAAkC,CAAC;AAEjD,cAAc,8BAA8B,CAAC;AAE7C,cAAc,iCAAiC,CAAC;AAEhD,cAAc,+BAA+B,CAAC;AAE9C,cAAc,+BAA+B,CAAC;AAE9C,cAAc,qCAAqC,CAAC;AAEpD,uBAAuB;AAEvB,cAAc,qDAAqD,CAAC;AAEpE,cAAc,qDAAqD,CAAC;AAEpE,cAAc,gDAAgD,CAAC;AAE/D,8BAA8B;AAE9B,cAAc,+CAA+C,CAAC;AAE9D,cAAc,uDAAuD,CAAC;AAEtE,cAAc,oDAAoD,CAAC;AAEnE,cAAc,qCAAqC,CAAC;AAEpD,qBAAqB;AAErB,cAAc,8CAA8C,CAAC;AAE7D,cAAc,uCAAuC,CAAC;AAEtD,cAAc,6CAA6C,CAAC;AAE5D,cAAc,iDAAiD,CAAC;AAEhE,oBAAoB;AAEpB,cAAc,sDAAsD,CAAC;AAErE,cAAc,iDAAiD,CAAC;AAEhE,cAAc,uCAAuC,CAAC;AAEtD,cAAc,4CAA4C,CAAC;AAE3D,cAAc,qCAAqC,CAAC;AAEpD,cAAc,4DAA4D,CAAC;AAE3E,cAAc,mDAAmD,CAAC;AAElE,cAAc,oDAAoD,CAAC;AAEnE,cAAc,iDAAiD,CAAC;AAEhE,cAAc,sCAAsC,CAAC;AAErD,cAAc,oCAAoC,CAAC;AAEnD,cAAc,4CAA4C,CAAC;AAE3D,cAAc,gEAAgE,CAAC;AAE/E,cAAc,kDAAkD,CAAC;AAEjE,cAAc,mDAAmD,CAAC;AAElE,cAAc,mDAAmD,CAAC;AAElE,cAAc,uCAAuC,CAAC;AAEtD,wBAAwB;AAExB,cAAc,kDAAkD,CAAC;AAEjE,cAAc,yDAAyD,CAAC;AAExE,cAAc,6CAA6C,CAAC;AAE5D,cAAc,+CAA+C,CAAC;AAE9D,cAAc,iDAAiD,CAAC;AAEhE,cAAc,oDAAoD,CAAC;AAEnE,cAAc,mDAAmD,CAAC"}
+{"version":3,"file":"index.js","sourceRoot":"","sources":["../src/index.ts"],"names":[],"mappings":"AAAA;;;;GAIG;AAEH,mBAAmB;AACnB,cAAc,qCAAqC,CAAC;AAEpD,cAAc,kCAAkC,CAAC;AAEjD,yBAAyB;AAEzB,cAAc,8DAA8D,CAAC;AAE7E,cAAc,kEAAkE,CAAC;AAEjF,cAAc,wDAAwD,CAAC;AAEvE,cAAc,8DAA8D,CAAC;AAE7E,cAAc,qDAAqD,CAAC;AAEpE,cAAc,4DAA4D,CAAC;AAE3E,qBAAqB;AAErB,cAAc,wCAAwC,CAAC;AAEvD,cAAc,+BAA+B,CAAC;AAE9C,cAAc,kCAAkC,CAAC;AAEjD,cAAc,8BAA8B,CAAC;AAE7C,cAAc,iCAAiC,CAAC;AAEhD,cAAc,+BAA+B,CAAC;AAE9C,cAAc,+BAA+B,CAAC;AAE9C,cAAc,qCAAqC,CAAC;AAEpD,uBAAuB;AAEvB,cAAc,gDAAgD,CAAC;AAE/D,cAAc,qDAAqD,CAAC;AAEpE,cAAc,gDAAgD,CAAC;AAE/D,8BAA8B;AAE9B,cAAc,+CAA+C,CAAC;AAE9D,cAAc,uDAAuD,CAAC;AAEtE,cAAc,oDAAoD,CAAC;AAEnE,cAAc,qCAAqC,CAAC;AAEpD,qBAAqB;AAErB,cAAc,8CAA8C,CAAC;AAE7D,cAAc,uCAAuC,CAAC;AAEtD,cAAc,6CAA6C,CAAC;AAE5D,cAAc,iDAAiD,CAAC;AAEhE,oBAAoB;AAEpB,cAAc,sDAAsD,CAAC;AAErE,cAAc,iDAAiD,CAAC;AAEhE,cAAc,uCAAuC,CAAC;AAEtD,cAAc,4CAA4C,CAAC;AAE3D,cAAc,qCAAqC,CAAC;AAEpD,cAAc,4DAA4D,CAAC;AAE3E,cAAc,mDAAmD,CAAC;AAElE,cAAc,oDAAoD,CAAC;AAEnE,cAAc,iDAAiD,CAAC;AAEhE,cAAc,sCAAsC,CAAC;AAErD,cAAc,oCAAoC,CAAC;AAEnD,cAAc,4CAA4C,CAAC;AAE3D,cAAc,gEAAgE,CAAC;AAE/E,cAAc,kDAAkD,CAAC;AAEjE,cAAc,mDAAmD,CAAC;AAElE,cAAc,mDAAmD,CAAC;AAElE,cAAc,uCAAuC,CAAC;AAEtD,wBAAwB;AAExB,cAAc,kDAAkD,CAAC;AAEjE,cAAc,yDAAyD,CAAC;AAExE,cAAc,6CAA6C,CAAC;AAE5D,cAAc,+CAA+C,CAAC;AAE9D,cAAc,iDAAiD,CAAC;AAEhE,cAAc,oDAAoD,CAAC;AAEnE,cAAc,mDAAmD,CAAC"}

package/package.json

@@ -1,6 +1,6 @@
 {
   "name": "@interstellar-tools/equations",
-  "version": "0.60.1",
+  "version": "0.61.0",
   "description": "Curated, well-tested orbital mechanics and dynamics formulas (e.g., Newtonian gravity, vis-viva, Kepler's 3rd law, escape/circular speeds, specific energy, specific angular momentum, peri/apoapsis radii).",
   "keywords": [
     "interstellar",
@@ -77,7 +77,7 @@
     "geo",
     "earth-sun"
   ],
-  "homepage": "https://phun-ky.net/projects/interstellar-tools",
+  "homepage": "https://phun-ky.net/interstellar-tools/",
   "bugs": {
     "url": "https://github.com/phun-ky/interstellar-tools/issues"
   },
@@ -110,12 +110,13 @@
     "test:ci": "glob-bin -c \"node --import tsx --test --no-warnings --experimental-test-coverage --test-coverage-exclude='**/__tests__/helpers/**' --test-reporter=cobertura --test-reporter-destination=coverage/cobertura-coverage.xml --test-reporter=spec --test-reporter-destination=stdout\" \"**/__tests__/**/*.spec.ts\""
   },
   "dependencies": {
-    "@interstellar-tools/constants": "^0.60.0",
+    "@interstellar-tools/constants": "^0.61.0",
     "@interstellar-tools/temporal": "^0.60.0",
-    "@interstellar-tools/types": "^0.60.0"
+    "@interstellar-tools/types": "^0.61.0"
   },
   "devDependencies": {
     "@types/node": "^24.19.2",
+    "fast-check": "4.10.2",
     "prettier": "^3.9.10",
     "typescript": "^6.0.3"
   },

package/dist/categories/anomalies/mean-to-eccentric-anomaly.d.ts

@@ -1,47 +0,0 @@
-import type { CelestialBodyType, Radians, TemporalInterface } from '@interstellar-tools/types';
-/**
- * Computes the **mean anomaly** ($M$) of a celestial body for a given time step.
- *
- * **Mathematical Explanation:**
- *
- * The **mean anomaly** is a measure of the position of an orbiting body as if it moved
- * with uniform angular motion along a circular orbit with the same period as the actual
- * elliptical orbit. It is computed as:
- *
- * $$ M = M_0 + n \cdot \Delta T $$
- *
- * where:
- * - $M_0$ is the **initial mean anomaly** (converted from true anomaly if necessary).
- * - $n$ is the **mean motion**, given by:
- *   $$ n = \frac{2\pi}{P} $$
- *   where $P$ is the orbital period in days.
- * - $\Delta T$ is the time step.
- *
- * @param {CelestialBodyType} body - The celestial body for which to compute the mean anomaly.
- * @param {TemporalInterface} timeStep - The time step over which to compute the change.
- * @returns {Radians} The computed mean anomaly in radians.
- *
- * @throws {@link https://developer.mozilla.org/en-US/docs/Web/JavaScript/Reference/Global_Objects/RangeError | RangeError} If the body's eccentricity is outside the range $0 \leq e < 1$.
- *
- * @example
- * ```ts
- * import { meanToEccentricAnomaly } from '@interstellar-tools/equations';
- *
- * const mars: CelestialBodyType = {
- *   name: 'Mars',
- *   e: 0.0934, // Eccentricity of Mars' orbit
- *   angle: 1.047, // Initial true anomaly (in radians)
- *   period: { value: 687, unit: 'day' }, // Orbital period in days
- * };
- *
- * const timeStep: TemporalInterface = { value: 1, unit: 'day' }; // 1-day step
- * const meanAnomaly = meanToEccentricAnomaly(mars, timeStep);
- * console.log(meanAnomaly); // Output: Computed mean anomaly in radians
- * ```
- *
- * @see [Mean Anomaly (Wikipedia)](https://en.wikipedia.org/wiki/Mean_anomaly)
- * @see [Orbital Mechanics (NASA)](https://solarsystem.nasa.gov/basics/chapter2-2/)
- * @group Anomalies
- */
-export declare const meanToEccentricAnomaly: (body: CelestialBodyType, timeStep: TemporalInterface) => Radians;
-//# sourceMappingURL=mean-to-eccentric-anomaly.d.ts.map

package/dist/categories/anomalies/mean-to-eccentric-anomaly.d.ts.map

@@ -1 +0,0 @@
-{"version":3,"file":"mean-to-eccentric-anomaly.d.ts","sourceRoot":"","sources":["../../../src/categories/anomalies/mean-to-eccentric-anomaly.ts"],"names":[],"mappings":"AAEA,OAAO,KAAK,EACV,iBAAiB,EACjB,OAAO,EACP,iBAAiB,EAClB,MAAM,2BAA2B,CAAC;AAQnC;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GA2CG;AAEH,eAAO,MAAM,sBAAsB,GACjC,MAAM,iBAAiB,EACvB,UAAU,iBAAiB,KAC1B,OAoCF,CAAC"}

package/dist/categories/anomalies/mean-to-eccentric-anomaly.js.map

@@ -1 +0,0 @@
-{"version":3,"file":"mean-to-eccentric-anomaly.js","sourceRoot":"","sources":["../../../src/categories/anomalies/mean-to-eccentric-anomaly.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,MAAM,EAAE,MAAM,+BAA+B,CAAC;AACvD,OAAO,EAAE,mBAAmB,EAAE,MAAM,8BAA8B,CAAC;AAOnE,OAAO,EAAE,SAAS,EAAE,MAAM,wBAAwB,CAAC;AAEnD,OAAO,EAAE,iBAAiB,EAAE,MAAM,2BAA2B,CAAC;AAC9D,OAAO,EAAE,QAAQ,EAAE,MAAM,sBAAsB,CAAC;AAEhD,2BAA2B;AAC3B;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GA2CG;AACH,0BAA0B;AAC1B,MAAM,CAAC,MAAM,sBAAsB,GAAG,CACpC,IAAuB,EACvB,QAA2B,EAClB,EAAE;IACX,MAAM,EAAE,CAAC,EAAE,KAAK,EAAE,MAAM,EAAE,GAAG,IAAI,CAAC;IAElC,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,IAAI,CAAC,EAAE,CAAC;QACpB,MAAM,IAAI,UAAU,CAClB,yBAAyB,CAAC,6CAA6C,CACxE,CAAC;IACJ,CAAC;IAED,iDAAiD;IACjD,MAAM,YAAY,GAAG,mBAAmB,CAAC,MAAM,EAAE,GAAG,CAAC,CAAC,KAAK,CAAC;IAC5D,MAAM,UAAU,GAAG,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,YAAY,CAAC,CAAC;IACnD,6DAA6D;IAC7D,MAAM,EAAE,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,iBAAiB,CAAC,KAAK,EAAE,CAAC,CAAC,CAAC;IAEzD,IAAI,QAAQ,CAAC,KAAK,KAAK,CAAC,EAAE,CAAC;QACzB,OAAO,EAAa,CAAC;IACvB,CAAC;IAED,2DAA2D;IAC3D,MAAM,WAAW,GAAG,IAAI,CAAC,CAAC,4BAA4B;IACtD,MAAM,iBAAiB,GAAG,CAAC,IAAI,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,6CAA6C;IAC1F,MAAM,aAAa,GAAG,IAAI,CAAC,GAAG,CAC5B,WAAW,EACX,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,KAAK,EAAE,iBAAiB,GAAG,IAAI,CAAC,GAAG,CAAC,YAAY,CAAC,CAAC,CACrE,CAAC;IACF,wEAAwE;IACxE,MAAM,KAAK,GAAG,EAAE,CAAC;IAEjB,IAAI,IAAI,GAAG,SAAS,CAAC,EAAE,GAAG,UAAU,GAAG,aAAa,CAAC,CAAC;IAEtD,IAAI,QAAQ,CAAC,IAAI,EAAE,KAAK,CAAC,EAAE,CAAC;QAC1B,IAAI,GAAG,KAAK,CAAC;IACf,CAAC;IAED,OAAO,IAAe,CAAC;AACzB,CAAC,CAAC"}