@interstellar-tools/equations 0.61.1 → 0.63.0

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

@@ -5,7 +5,8 @@ import type { OrbitInterface, Radians, TimeStepInterface } from '@interstellar-t
  *
  * 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.
+ * `angle` to move the body over time. To find a body at a given date instead,
+ * use {@link trueAnomalyAt}.
  *
  * **Step 1: Mean anomaly after the time step ($M$)**
  *
@@ -15,9 +16,9 @@ import type { OrbitInterface, Radians, TimeStepInterface } from '@interstellar-t
  * - $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$).
+ *   $P$ is ignored, so a positive step always moves the body forward.
+ * - $\Delta t$ is `timeStep.value` in days. It can be any length, and negative
+ *   to go back in time.
  *
  * The result is wrapped to $[0, 2\pi)$.
  *

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,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"}
+{"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;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GAuDG;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

@@ -7,7 +7,8 @@ import { solveKepler } from '../kepler/solve-kepler.js';
  *
  * 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.
+ * `angle` to move the body over time. To find a body at a given date instead,
+ * use {@link trueAnomalyAt}.
  *
  * **Step 1: Mean anomaly after the time step ($M$)**
  *
@@ -17,9 +18,9 @@ import { solveKepler } from '../kepler/solve-kepler.js';
  * - $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$).
+ *   $P$ is ignored, so a positive step always moves the body forward.
+ * - $\Delta t$ is `timeStep.value` in days. It can be any length, and negative
+ *   to go back in time.
  *
  * The result is wrapped to $[0, 2\pi)$.
  *

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,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"}
+{"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;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GAuDG;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

@@ -3,26 +3,29 @@ import type { OrbitInterface, Radians, TemporalInterface } from '@interstellar-t
  * 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}.
+ * anomaly from the result, pass it to {@link solveKepler}. To find a body at a
+ * given date, use {@link meanAnomalyAt} or {@link trueAnomalyAt}.
  *
  * **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:
+ * elliptical orbit. It grows linearly with time:
  *
- * $$ M = M_0 + n \cdot \Delta T $$
+ * $$ M = M_0 + n \cdot \Delta T, \qquad n = \frac{2\pi}{|P|} $$
  *
  * 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.
+ * - $M_0$ is the **initial mean anomaly**, converted from the true anomaly `angle`.
+ * - $n$ is the **mean motion**, from the period $P$ in days. The sign of $P$ is
+ *   ignored, so a positive step always moves the body forward along its orbit.
  * - $\Delta T$ is the time step, converted to days from its unit.
  *
+ * The step can be any length, from a fraction of a second to many orbits, and
+ * negative to go back in time. The result is wrapped to $[0, 2\pi)$.
+ *
  * @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.
+ * @param {TemporalInterface} timeStep - How far to advance, in any temporal unit; negative goes back in time.
+ * @returns {Radians} The mean anomaly after the time step, in $[0, 2\pi)$.
  *
  * @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$.
  *
@@ -36,9 +39,9 @@ import type { OrbitInterface, Radians, TemporalInterface } from '@interstellar-t
  *   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
+ * advanceMeanAnomaly(marsOrbit, { value: 1, unit: 'd' }); // one day on
+ * advanceMeanAnomaly(marsOrbit, { value: 10, unit: 'yr' }); // ten years on
+ * advanceMeanAnomaly(marsOrbit, { value: -1, unit: 'd' }); // one day back
  * ```
  *
  * @see {@link solveKepler} to convert the mean anomaly to the eccentric anomaly.

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

@@ -1 +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"}
+{"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;AAOnC;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GAiDG;AAEH,eAAO,MAAM,kBAAkB,GAC7B,OAAO,IAAI,CAAC,cAAc,EAAE,GAAG,GAAG,OAAO,GAAG,QAAQ,CAAC,EACrD,UAAU,iBAAiB,KAC1B,OAiBF,CAAC"}

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

@@ -1,33 +1,35 @@
 import { TWO_PI } from '@interstellar-tools/constants';
 import { convertTemporalUnit } from '@interstellar-tools/temporal';
-import { wrapAngle } from '../angle/wrap-angle.js';
+import { norm2pi } from '../helpers/misc.js';
 import { trueToMeanAnomaly } from './true-to-mean-anomaly.js';
-import { areEqual } from './utils/are-equal.js';
 /* node:coverage disable */
 /**
  * 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}.
+ * anomaly from the result, pass it to {@link solveKepler}. To find a body at a
+ * given date, use {@link meanAnomalyAt} or {@link trueAnomalyAt}.
  *
  * **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:
+ * elliptical orbit. It grows linearly with time:
  *
- * $$ M = M_0 + n \cdot \Delta T $$
+ * $$ M = M_0 + n \cdot \Delta T, \qquad n = \frac{2\pi}{|P|} $$
  *
  * 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.
+ * - $M_0$ is the **initial mean anomaly**, converted from the true anomaly `angle`.
+ * - $n$ is the **mean motion**, from the period $P$ in days. The sign of $P$ is
+ *   ignored, so a positive step always moves the body forward along its orbit.
  * - $\Delta T$ is the time step, converted to days from its unit.
  *
+ * The step can be any length, from a fraction of a second to many orbits, and
+ * negative to go back in time. The result is wrapped to $[0, 2\pi)$.
+ *
  * @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.
+ * @param {TemporalInterface} timeStep - How far to advance, in any temporal unit; negative goes back in time.
+ * @returns {Radians} The mean anomaly after the time step, in $[0, 2\pi)$.
  *
  * @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$.
  *
@@ -41,9 +43,9 @@ import { areEqual } from './utils/are-equal.js';
  *   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
+ * advanceMeanAnomaly(marsOrbit, { value: 1, unit: 'd' }); // one day on
+ * advanceMeanAnomaly(marsOrbit, { value: 10, unit: 'yr' }); // ten years on
+ * advanceMeanAnomaly(marsOrbit, { value: -1, unit: 'd' }); // one day back
  * ```
  *
  * @see {@link solveKepler} to convert the mean anomaly to the eccentric anomaly.
@@ -61,21 +63,8 @@ export const advanceMeanAnomaly = (orbit, timeStep) => {
     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
+    // Convert true anomaly (V) to mean anomaly (M0)
     const M0 = e === 0 ? angle : trueToMeanAnomaly(angle, e);
-    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(timeStepInDays, maxTimeStepFactor * Math.abs(periodInDays)));
-    // Compute new mean anomaly with hysteresis filter and sign preservation
-    const prevM = M0;
-    let newM = wrapAngle(M0 + meanMotion * clampedDeltaT);
-    if (areEqual(newM, prevM)) {
-        newM = prevM;
-    }
-    return newM;
+    return norm2pi(M0 + meanMotion * timeStepInDays);
 };
 //# sourceMappingURL=advance-mean-anomaly.js.map

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

@@ -1 +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"}
+{"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,OAAO,EAAE,MAAM,oBAAoB,CAAC;AAE7C,OAAO,EAAE,iBAAiB,EAAE,MAAM,2BAA2B,CAAC;AAE9D,2BAA2B;AAC3B;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GAiDG;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,gDAAgD;IAChD,MAAM,EAAE,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,iBAAiB,CAAC,KAAK,EAAE,CAAC,CAAC,CAAC;IAEzD,OAAO,OAAO,CAAC,EAAE,GAAG,UAAU,GAAG,cAAc,CAAY,CAAC;AAC9D,CAAC,CAAC"}

package/dist/categories/anomalies/anomaly-at.d.ts

@@ -0,0 +1,61 @@
+import type { OrbitInterface, Radians } from '@interstellar-tools/types';
+/**
+ * The **mean anomaly** ($M$) of an orbit at a given date.
+ *
+ * Every `orbit.angle` in the datasets is the true anomaly at **J2000**
+ * (2000-01-01 12:00 TT), so this advances the orbit from J2000 to `date` with
+ * {@link advanceMeanAnomaly}. Dates before J2000 work too.
+ *
+ * $$ M = M_0 + \frac{2\pi}{|P|} \, (t - t_{J2000}) $$
+ *
+ * The time is measured in UTC, as JavaScript dates are, against
+ * {@link J2000_UTC}. Leap seconds since 2000 are ignored, which moves a body by
+ * a few seconds of its orbit at most.
+ *
+ * @param {Pick<OrbitInterface, 'e' | 'angle' | 'period'>} orbit - An orbit whose `angle` is the true anomaly at J2000, for example `body.orbit` from a dataset.
+ * @param {Date} date - The date to find the body at.
+ * @returns {Radians} The mean anomaly at `date`, in $[0, 2\pi)$.
+ * @throws {TypeError} If `date` is not a valid `Date`.
+ * @throws {RangeError} If the orbit's eccentricity is outside $0 \leq e < 1$.
+ *
+ * @example
+ * ```ts
+ * import { PLANETS } from '@interstellar-tools/constants';
+ * import { meanAnomalyAt } from '@interstellar-tools/equations';
+ *
+ * const earth = PLANETS.find((planet) => planet.id === 'earth')!;
+ *
+ * meanAnomalyAt(earth.orbit!, new Date('2026-10-17T12:00:00Z')); // ≈ 4.934 rad (282.7°)
+ * ```
+ *
+ * @see {@link trueAnomalyAt} for where the body actually is.
+ * @group Anomalies
+ */
+export declare const meanAnomalyAt: (orbit: Pick<OrbitInterface, "e" | "angle" | "period">, date: Date) => Radians;
+/**
+ * The **true anomaly** ($\nu$) of an orbit at a given date: where the body is.
+ *
+ * Takes the mean anomaly at `date` from {@link meanAnomalyAt}, solves Kepler's
+ * equation with {@link solveKepler}, and converts the result with
+ * {@link eccentricToTrueAnomaly}.
+ *
+ * @param {Pick<OrbitInterface, 'e' | 'angle' | 'period'>} orbit - An orbit whose `angle` is the true anomaly at J2000, for example `body.orbit` from a dataset.
+ * @param {Date} date - The date to find the body at.
+ * @returns {Radians} The true anomaly at `date`, in $[0, 2\pi)$.
+ * @throws {TypeError} If `date` is not a valid `Date`.
+ * @throws {RangeError} If the orbit's eccentricity is outside $0 \leq e < 1$.
+ *
+ * @example
+ * ```ts
+ * import { PLANETS } from '@interstellar-tools/constants';
+ * import { trueAnomalyAt } from '@interstellar-tools/equations';
+ *
+ * const earth = PLANETS.find((planet) => planet.id === 'earth')!;
+ *
+ * trueAnomalyAt(earth.orbit!, new Date('2026-10-17T12:00:00Z')); // ≈ 4.901 rad (280.8°)
+ * ```
+ *
+ * @group Anomalies
+ */
+export declare const trueAnomalyAt: (orbit: Pick<OrbitInterface, "e" | "angle" | "period">, date: Date) => Radians;
+//# sourceMappingURL=anomaly-at.d.ts.map

package/dist/categories/anomalies/anomaly-at.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"anomaly-at.d.ts","sourceRoot":"","sources":["../../../src/categories/anomalies/anomaly-at.ts"],"names":[],"mappings":"AACA,OAAO,KAAK,EAAE,cAAc,EAAE,OAAO,EAAE,MAAM,2BAA2B,CAAC;AAkBzE;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GA+BG;AAEH,eAAO,MAAM,aAAa,GACxB,OAAO,IAAI,CAAC,cAAc,EAAE,GAAG,GAAG,OAAO,GAAG,QAAQ,CAAC,EACrD,MAAM,IAAI,KACT,OACoE,CAAC;AAGxE;;;;;;;;;;;;;;;;;;;;;;;;GAwBG;AAEH,eAAO,MAAM,aAAa,GACxB,OAAO,IAAI,CAAC,cAAc,EAAE,GAAG,GAAG,OAAO,GAAG,QAAQ,CAAC,EACrD,MAAM,IAAI,KACT,OAMF,CAAC"}

package/dist/categories/anomalies/anomaly-at.js

@@ -0,0 +1,81 @@
+import { J2000_UTC, MILLISECONDS_PER_DAY } from '@interstellar-tools/constants';
+import { norm2pi } from '../helpers/misc.js';
+import { solveKepler } from '../kepler/solve-kepler.js';
+import { advanceMeanAnomaly } from './advance-mean-anomaly.js';
+import { eccentricToTrueAnomaly } from './eccentric-to-true-anomaly.js';
+/** Days from J2000 to `date`, after checking it's a valid date. */
+const daysSinceJ2000 = (date) => {
+    if (!(date instanceof Date) || !Number.isFinite(date.getTime())) {
+        throw new TypeError(`date must be a valid Date. Received: ${String(date)}`);
+    }
+    return (date.getTime() - J2000_UTC.getTime()) / MILLISECONDS_PER_DAY;
+};
+/* node:coverage disable */
+/**
+ * The **mean anomaly** ($M$) of an orbit at a given date.
+ *
+ * Every `orbit.angle` in the datasets is the true anomaly at **J2000**
+ * (2000-01-01 12:00 TT), so this advances the orbit from J2000 to `date` with
+ * {@link advanceMeanAnomaly}. Dates before J2000 work too.
+ *
+ * $$ M = M_0 + \frac{2\pi}{|P|} \, (t - t_{J2000}) $$
+ *
+ * The time is measured in UTC, as JavaScript dates are, against
+ * {@link J2000_UTC}. Leap seconds since 2000 are ignored, which moves a body by
+ * a few seconds of its orbit at most.
+ *
+ * @param {Pick<OrbitInterface, 'e' | 'angle' | 'period'>} orbit - An orbit whose `angle` is the true anomaly at J2000, for example `body.orbit` from a dataset.
+ * @param {Date} date - The date to find the body at.
+ * @returns {Radians} The mean anomaly at `date`, in $[0, 2\pi)$.
+ * @throws {TypeError} If `date` is not a valid `Date`.
+ * @throws {RangeError} If the orbit's eccentricity is outside $0 \leq e < 1$.
+ *
+ * @example
+ * ```ts
+ * import { PLANETS } from '@interstellar-tools/constants';
+ * import { meanAnomalyAt } from '@interstellar-tools/equations';
+ *
+ * const earth = PLANETS.find((planet) => planet.id === 'earth')!;
+ *
+ * meanAnomalyAt(earth.orbit!, new Date('2026-10-17T12:00:00Z')); // ≈ 4.934 rad (282.7°)
+ * ```
+ *
+ * @see {@link trueAnomalyAt} for where the body actually is.
+ * @group Anomalies
+ */
+/* node:coverage enable */
+export const meanAnomalyAt = (orbit, date) => advanceMeanAnomaly(orbit, { value: daysSinceJ2000(date), unit: 'd' });
+/* node:coverage disable */
+/**
+ * The **true anomaly** ($\nu$) of an orbit at a given date: where the body is.
+ *
+ * Takes the mean anomaly at `date` from {@link meanAnomalyAt}, solves Kepler's
+ * equation with {@link solveKepler}, and converts the result with
+ * {@link eccentricToTrueAnomaly}.
+ *
+ * @param {Pick<OrbitInterface, 'e' | 'angle' | 'period'>} orbit - An orbit whose `angle` is the true anomaly at J2000, for example `body.orbit` from a dataset.
+ * @param {Date} date - The date to find the body at.
+ * @returns {Radians} The true anomaly at `date`, in $[0, 2\pi)$.
+ * @throws {TypeError} If `date` is not a valid `Date`.
+ * @throws {RangeError} If the orbit's eccentricity is outside $0 \leq e < 1$.
+ *
+ * @example
+ * ```ts
+ * import { PLANETS } from '@interstellar-tools/constants';
+ * import { trueAnomalyAt } from '@interstellar-tools/equations';
+ *
+ * const earth = PLANETS.find((planet) => planet.id === 'earth')!;
+ *
+ * trueAnomalyAt(earth.orbit!, new Date('2026-10-17T12:00:00Z')); // ≈ 4.901 rad (280.8°)
+ * ```
+ *
+ * @group Anomalies
+ */
+/* node:coverage enable */
+export const trueAnomalyAt = (orbit, date) => {
+    const M = meanAnomalyAt(orbit, date);
+    const E = solveKepler(M, orbit.e);
+    const nu = eccentricToTrueAnomaly(E, orbit.e);
+    return norm2pi(nu);
+};
+//# sourceMappingURL=anomaly-at.js.map

package/dist/categories/anomalies/anomaly-at.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"anomaly-at.js","sourceRoot":"","sources":["../../../src/categories/anomalies/anomaly-at.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,SAAS,EAAE,oBAAoB,EAAE,MAAM,+BAA+B,CAAC;AAGhF,OAAO,EAAE,OAAO,EAAE,MAAM,oBAAoB,CAAC;AAC7C,OAAO,EAAE,WAAW,EAAE,MAAM,2BAA2B,CAAC;AAExD,OAAO,EAAE,kBAAkB,EAAE,MAAM,2BAA2B,CAAC;AAC/D,OAAO,EAAE,sBAAsB,EAAE,MAAM,gCAAgC,CAAC;AAExE,mEAAmE;AACnE,MAAM,cAAc,GAAG,CAAC,IAAU,EAAU,EAAE;IAC5C,IAAI,CAAC,CAAC,IAAI,YAAY,IAAI,CAAC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,IAAI,CAAC,OAAO,EAAE,CAAC,EAAE,CAAC;QAChE,MAAM,IAAI,SAAS,CAAC,wCAAwC,MAAM,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC;IAC9E,CAAC;IAED,OAAO,CAAC,IAAI,CAAC,OAAO,EAAE,GAAG,SAAS,CAAC,OAAO,EAAE,CAAC,GAAG,oBAAoB,CAAC;AACvE,CAAC,CAAC;AAEF,2BAA2B;AAC3B;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GA+BG;AACH,0BAA0B;AAC1B,MAAM,CAAC,MAAM,aAAa,GAAG,CAC3B,KAAqD,EACrD,IAAU,EACD,EAAE,CACX,kBAAkB,CAAC,KAAK,EAAE,EAAE,KAAK,EAAE,cAAc,CAAC,IAAI,CAAC,EAAE,IAAI,EAAE,GAAG,EAAE,CAAC,CAAC;AAExE,2BAA2B;AAC3B;;;;;;;;;;;;;;;;;;;;;;;;GAwBG;AACH,0BAA0B;AAC1B,MAAM,CAAC,MAAM,aAAa,GAAG,CAC3B,KAAqD,EACrD,IAAU,EACD,EAAE;IACX,MAAM,CAAC,GAAG,aAAa,CAAC,KAAK,EAAE,IAAI,CAAC,CAAC;IACrC,MAAM,CAAC,GAAG,WAAW,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,CAAC;IAClC,MAAM,EAAE,GAAG,sBAAsB,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,CAAC;IAE9C,OAAO,OAAO,CAAC,EAAE,CAAY,CAAC;AAChC,CAAC,CAAC"}

package/dist/categories/orbits/cw-hill-derivatives.d.ts

@@ -56,7 +56,7 @@ import { CwState, CwStateDerivative } from '@interstellar-tools/types';
  *
  * const mu = 398600.4418;  // km³/s²
  * const a  = 7000;         // km
- * const n  = meanMotion(mu, a); // rad/s
+ * const n  = meanMotion(a, mu); // rad/s
  *
  * const s: CwState = [100, 0, 0, 0, 0.01, 0];
  * const ds = cwHillDerivatives(s, n);

package/dist/categories/orbits/cw-hill-derivatives.js

@@ -55,7 +55,7 @@
  *
  * const mu = 398600.4418;  // km³/s²
  * const a  = 7000;         // km
- * const n  = meanMotion(mu, a); // rad/s
+ * const n  = meanMotion(a, mu); // rad/s
  *
  * const s: CwState = [100, 0, 0, 0, 0.01, 0];
  * const ds = cwHillDerivatives(s, n);

package/dist/categories/orbits/hyperbolic-periapsis-speed.d.ts

@@ -16,19 +16,19 @@
  * - If $ \mu $ is in km³/s², then $r_p$ must be in km and speeds are in km/s.
  *
  * @param vInfinity Hyperbolic excess speed $ v_{\infty} $ (must be finite, >= 0).
- * @param mu Standard gravitational parameter $ \mu $ (must be finite, > 0).
  * @param rp Periapsis radius $ r_p $ from the central body's center (must be finite, > 0).
+ * @param mu Standard gravitational parameter $ \mu $ (must be finite, > 0).
  * @returns Periapsis speed $ v_p $ in the same speed units as `vInfinity`.
  * @throws {TypeError} If any input is not finite.
- * @throws {RangeError} If `vInfinity < 0`, `mu <= 0`, or `rp <= 0`.
+ * @throws {RangeError} If `vInfinity < 0`, `rp <= 0`, or `mu <= 0`.
  * @group Orbits
  *
  * @example
  * ```ts
  * // Earth example (km-based):
  * // mu = 398600.4418 km³/s², rp ~ 6678 km (LEO-ish), v∞ = 3.2 km/s
- * const vp = hyperbolicPeriapsisSpeed(3.2, 398600.4418, 6678);
+ * const vp = hyperbolicPeriapsisSpeed(3.2, 6678, 398600.4418);
  * ```
  */
-export declare const hyperbolicPeriapsisSpeed: (vInfinity: number, mu: number, rp: number) => number;
+export declare const hyperbolicPeriapsisSpeed: (vInfinity: number, rp: number, mu: number) => number;
 //# sourceMappingURL=hyperbolic-periapsis-speed.d.ts.map

package/dist/categories/orbits/hyperbolic-periapsis-speed.js

@@ -17,40 +17,40 @@
  * - If $ \mu $ is in km³/s², then $r_p$ must be in km and speeds are in km/s.
  *
  * @param vInfinity Hyperbolic excess speed $ v_{\infty} $ (must be finite, >= 0).
- * @param mu Standard gravitational parameter $ \mu $ (must be finite, > 0).
  * @param rp Periapsis radius $ r_p $ from the central body's center (must be finite, > 0).
+ * @param mu Standard gravitational parameter $ \mu $ (must be finite, > 0).
  * @returns Periapsis speed $ v_p $ in the same speed units as `vInfinity`.
  * @throws {TypeError} If any input is not finite.
- * @throws {RangeError} If `vInfinity < 0`, `mu <= 0`, or `rp <= 0`.
+ * @throws {RangeError} If `vInfinity < 0`, `rp <= 0`, or `mu <= 0`.
  * @group Orbits
  *
  * @example
  * ```ts
  * // Earth example (km-based):
  * // mu = 398600.4418 km³/s², rp ~ 6678 km (LEO-ish), v∞ = 3.2 km/s
- * const vp = hyperbolicPeriapsisSpeed(3.2, 398600.4418, 6678);
+ * const vp = hyperbolicPeriapsisSpeed(3.2, 6678, 398600.4418);
  * ```
  */
 /* node:coverage enable */
-export const hyperbolicPeriapsisSpeed = (vInfinity, mu, rp) => {
+export const hyperbolicPeriapsisSpeed = (vInfinity, rp, mu) => {
     if (!Number.isFinite(vInfinity)) {
         throw new TypeError(`vInfinity must be a finite number. Received: ${vInfinity}`);
     }
-    if (!Number.isFinite(mu)) {
-        throw new TypeError(`mu must be a finite number. Received: ${mu}`);
-    }
     if (!Number.isFinite(rp)) {
         throw new TypeError(`rp must be a finite number. Received: ${rp}`);
     }
+    if (!Number.isFinite(mu)) {
+        throw new TypeError(`mu must be a finite number. Received: ${mu}`);
+    }
     if (vInfinity < 0) {
         throw new RangeError(`vInfinity must be >= 0. Received: ${vInfinity}`);
     }
-    if (mu <= 0) {
-        throw new RangeError(`mu must be > 0. Received: ${mu}`);
-    }
     if (rp <= 0) {
         throw new RangeError(`rp must be > 0. Received: ${rp}`);
     }
+    if (mu <= 0) {
+        throw new RangeError(`mu must be > 0. Received: ${mu}`);
+    }
     const term = vInfinity * vInfinity + (2 * mu) / rp;
     // With valid ranges above, term should be positive, but guard anyway.
     if (term < 0) {

package/dist/categories/orbits/mean-motion.d.ts

@@ -28,16 +28,16 @@
  * orbits per day, the unit used in a satellite's two-line elements (TLEs).
  *
  * ```ts
- * meanMotion(3.986004418e14, 6.371e6 + 420e3); // ≈ 0.001128 rad/s
+ * meanMotion(6.371e6 + 420e3, 3.986004418e14); // ≈ 0.001128 rad/s
  * ```
  *
  * :::
  *
- * @param mu - Standard gravitational parameter $ \mu $ (finite, > 0).
  * @param a - Semi-major axis $ a $ (finite, > 0).
+ * @param mu - Standard gravitational parameter $ \mu $ (finite, > 0).
  * @returns Mean motion $ n $ in **rad/s**.
- * @throws {TypeError} If `mu` or `a` is not a finite number.
- * @throws {RangeError} If `mu <= 0` or `a <= 0`.
+ * @throws {TypeError} If `a` or `mu` is not a finite number.
+ * @throws {RangeError} If `a <= 0` or `mu <= 0`.
  * @group Orbits
  *
  * @example
@@ -45,8 +45,8 @@
  * // Earth example (km-based):
  * const mu = 398600.4418; // km³/s²
  * const a = 7000;         // km
- * const n = meanMotion(mu, a); // rad/s
+ * const n = meanMotion(a, mu); // rad/s
  * ```
  */
-export declare const meanMotion: (mu: number, a: number) => number;
+export declare const meanMotion: (a: number, mu: number) => number;
 //# sourceMappingURL=mean-motion.d.ts.map

package/dist/categories/orbits/mean-motion.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"mean-motion.d.ts","sourceRoot":"","sources":["../../../src/categories/orbits/mean-motion.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GAiDG;AACH,eAAO,MAAM,UAAU,GAAI,IAAI,MAAM,EAAE,GAAG,MAAM,KAAG,MAYlD,CAAC"}
+{"version":3,"file":"mean-motion.d.ts","sourceRoot":"","sources":["../../../src/categories/orbits/mean-motion.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GAiDG;AACH,eAAO,MAAM,UAAU,GAAI,GAAG,MAAM,EAAE,IAAI,MAAM,KAAG,MAYlD,CAAC"}

package/dist/categories/orbits/mean-motion.js

@@ -28,16 +28,16 @@
  * orbits per day, the unit used in a satellite's two-line elements (TLEs).
  *
  * ```ts
- * meanMotion(3.986004418e14, 6.371e6 + 420e3); // ≈ 0.001128 rad/s
+ * meanMotion(6.371e6 + 420e3, 3.986004418e14); // ≈ 0.001128 rad/s
  * ```
  *
  * :::
  *
- * @param mu - Standard gravitational parameter $ \mu $ (finite, > 0).
  * @param a - Semi-major axis $ a $ (finite, > 0).
+ * @param mu - Standard gravitational parameter $ \mu $ (finite, > 0).
  * @returns Mean motion $ n $ in **rad/s**.
- * @throws {TypeError} If `mu` or `a` is not a finite number.
- * @throws {RangeError} If `mu <= 0` or `a <= 0`.
+ * @throws {TypeError} If `a` or `mu` is not a finite number.
+ * @throws {RangeError} If `a <= 0` or `mu <= 0`.
  * @group Orbits
  *
  * @example
@@ -45,18 +45,18 @@
  * // Earth example (km-based):
  * const mu = 398600.4418; // km³/s²
  * const a = 7000;         // km
- * const n = meanMotion(mu, a); // rad/s
+ * const n = meanMotion(a, mu); // rad/s
  * ```
  */
-export const meanMotion = (mu, a) => {
-    if (!Number.isFinite(mu))
-        throw new TypeError(`mu must be finite. Received: ${mu}`);
+export const meanMotion = (a, mu) => {
     if (!Number.isFinite(a))
         throw new TypeError(`a must be finite. Received: ${a}`);
-    if (mu <= 0)
-        throw new RangeError(`mu must be > 0. Received: ${mu}`);
+    if (!Number.isFinite(mu))
+        throw new TypeError(`mu must be finite. Received: ${mu}`);
     if (a <= 0)
         throw new RangeError(`a must be > 0. Received: ${a}`);
+    if (mu <= 0)
+        throw new RangeError(`mu must be > 0. Received: ${mu}`);
     return Math.sqrt(mu / (a * a * a));
 };
 //# sourceMappingURL=mean-motion.js.map

package/dist/categories/orbits/mean-motion.js.map

@@ -1 +1 @@
-{"version":3,"file":"mean-motion.js","sourceRoot":"","sources":["../../../src/categories/orbits/mean-motion.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GAiDG;AACH,MAAM,CAAC,MAAM,UAAU,GAAG,CAAC,EAAU,EAAE,CAAS,EAAU,EAAE;IAC1D,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,EAAE,CAAC;QACtB,MAAM,IAAI,SAAS,CAAC,gCAAgC,EAAE,EAAE,CAAC,CAAC;IAE5D,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC;QACrB,MAAM,IAAI,SAAS,CAAC,+BAA+B,CAAC,EAAE,CAAC,CAAC;IAE1D,IAAI,EAAE,IAAI,CAAC;QAAE,MAAM,IAAI,UAAU,CAAC,6BAA6B,EAAE,EAAE,CAAC,CAAC;IAErE,IAAI,CAAC,IAAI,CAAC;QAAE,MAAM,IAAI,UAAU,CAAC,4BAA4B,CAAC,EAAE,CAAC,CAAC;IAElE,OAAO,IAAI,CAAC,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;AACrC,CAAC,CAAC"}
+{"version":3,"file":"mean-motion.js","sourceRoot":"","sources":["../../../src/categories/orbits/mean-motion.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GAiDG;AACH,MAAM,CAAC,MAAM,UAAU,GAAG,CAAC,CAAS,EAAE,EAAU,EAAU,EAAE;IAC1D,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC;QACrB,MAAM,IAAI,SAAS,CAAC,+BAA+B,CAAC,EAAE,CAAC,CAAC;IAE1D,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,EAAE,CAAC;QACtB,MAAM,IAAI,SAAS,CAAC,gCAAgC,EAAE,EAAE,CAAC,CAAC;IAE5D,IAAI,CAAC,IAAI,CAAC;QAAE,MAAM,IAAI,UAAU,CAAC,4BAA4B,CAAC,EAAE,CAAC,CAAC;IAElE,IAAI,EAAE,IAAI,CAAC;QAAE,MAAM,IAAI,UAAU,CAAC,6BAA6B,EAAE,EAAE,CAAC,CAAC;IAErE,OAAO,IAAI,CAAC,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;AACrC,CAAC,CAAC"}

package/dist/index.d.ts

@@ -20,6 +20,7 @@ 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/anomaly-at.js';
 export * from './categories/anomalies/eccentric-to-true-anomaly.js';
 export * from './categories/anomalies/true-to-mean-anomaly.js';
 export * from './categories/kepler/solve-kepler-bisection.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,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"}
+{"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,sCAAsC,CAAC;AAErD,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

@@ -24,6 +24,7 @@ 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/anomaly-at.js';
 export * from './categories/anomalies/eccentric-to-true-anomaly.js';
 export * from './categories/anomalies/true-to-mean-anomaly.js';
 // categories/kepler (solvers)

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,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"}
+{"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,sCAAsC,CAAC;AAErD,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.61.1",
+  "version": "0.63.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",

package/dist/categories/anomalies/utils/are-equal.d.ts

@@ -1,10 +0,0 @@
-/**
- * Compares two floating-point numbers for equality within a tolerance.
- *
- * @param {number} a - First number to compare.
- * @param {number} b - Second number to compare.
- * @param {number} [epsilon=EPSILON] - Tolerance for floating-point precision.
- * @returns {boolean} `true` if values are approximately equal.
- */
-export declare const areEqual: (a: number, b: number, epsilon?: number) => boolean;
-//# sourceMappingURL=are-equal.d.ts.map

package/dist/categories/anomalies/utils/are-equal.d.ts.map

@@ -1 +0,0 @@
-{"version":3,"file":"are-equal.d.ts","sourceRoot":"","sources":["../../../../src/categories/anomalies/utils/are-equal.ts"],"names":[],"mappings":"AAEA;;;;;;;GAOG;AACH,eAAO,MAAM,QAAQ,GACnB,GAAG,MAAM,EACT,GAAG,MAAM,EACT,UAAS,MAAgB,KACxB,OAEF,CAAC"}

package/dist/categories/anomalies/utils/are-equal.js

@@ -1,13 +0,0 @@
-const EPSILON = 1e-15; // Increased tolerance for high eccentricity
-/**
- * Compares two floating-point numbers for equality within a tolerance.
- *
- * @param {number} a - First number to compare.
- * @param {number} b - Second number to compare.
- * @param {number} [epsilon=EPSILON] - Tolerance for floating-point precision.
- * @returns {boolean} `true` if values are approximately equal.
- */
-export const areEqual = (a, b, epsilon = EPSILON) => {
-    return Math.abs(a - b) < epsilon * Math.max(1, Math.abs(a), Math.abs(b));
-};
-//# sourceMappingURL=are-equal.js.map

package/dist/categories/anomalies/utils/are-equal.js.map

@@ -1 +0,0 @@
-{"version":3,"file":"are-equal.js","sourceRoot":"","sources":["../../../../src/categories/anomalies/utils/are-equal.ts"],"names":[],"mappings":"AAAA,MAAM,OAAO,GAAG,KAAK,CAAC,CAAC,4CAA4C;AAEnE;;;;;;;GAOG;AACH,MAAM,CAAC,MAAM,QAAQ,GAAG,CACtB,CAAS,EACT,CAAS,EACT,UAAkB,OAAO,EAChB,EAAE;IACX,OAAO,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,OAAO,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;AAC3E,CAAC,CAAC"}