@interstellar-tools/equations 1.1.0 → 1.3.0

package/dist/categories/anomalies/asymptote-true-anomaly.d.ts

@@ -0,0 +1,50 @@
+import type { Radians } from '@interstellar-tools/types';
+/**
+ * Compute the **true anomaly of a hyperbola's asymptotes**: the limit a body on a hyperbola approaches but never reaches.
+ *
+ * $$
+ * \nu_\infty = \arccos\left(-\frac{1}{e}\right)
+ * $$
+ *
+ * **Units**
+ * - Input: `e` dimensionless.
+ * - Output: $ \nu_\infty $ in **radians**, in (π/2, π).
+ *
+ * ::: info Notes
+ *
+ * - The body comes in from $ -\nu_\infty $ and leaves towards $ +\nu_\infty $.
+ * - The trajectory turns through $ 2\nu_\infty - \pi = 2\arcsin(1/e) $, the turning angle of a flyby; see `gravityAssistTurningAngle`.
+ * - Just above $ e = 1 $ the asymptotes open to almost π (a near-parabola); for large $ e $ they close towards π/2 (a near-straight line).
+ *
+ * :::
+ *
+ * ::: practice Apollo 11, arriving at the Moon
+ *
+ * **134.0°**
+ *
+ * _Asymptote of the approach hyperbola, e = 1.44_
+ *
+ * Without the lunar orbit insertion burn at closest approach, the Moon's gravity would have bent Apollo 11's path through 88° and sent it back out of the Moon's sphere of influence.
+ *
+ * ```ts
+ * asymptoteTrueAnomaly(1.44); // ≈ 2.338 rad, 134.0°
+ * ```
+ *
+ * :::
+ *
+ * @param e - Eccentricity $ e $ (finite, > 1).
+ * @returns Asymptote true anomaly $ \nu_\infty $ in **radians**.
+ * @throws {TypeError} If `e` is not a finite number.
+ * @throws {RangeError} If `e <= 1`.
+ * @group Anomalies
+ *
+ * @example
+ * ```ts
+ * asymptoteTrueAnomaly(2); // 2π/3 (120°)
+ * ```
+ *
+ * @see {@link trueToHyperbolicAnomaly}
+ * @see https://en.wikipedia.org/wiki/Hyperbolic_trajectory Hyperbolic trajectory
+ */
+export declare const asymptoteTrueAnomaly: (e: number) => Radians;
+//# sourceMappingURL=asymptote-true-anomaly.d.ts.map

package/dist/categories/anomalies/asymptote-true-anomaly.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"asymptote-true-anomaly.d.ts","sourceRoot":"","sources":["../../../src/categories/anomalies/asymptote-true-anomaly.ts"],"names":[],"mappings":"AAAA,OAAO,KAAK,EAAE,OAAO,EAAE,MAAM,2BAA2B,CAAC;AAEzD;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GA8CG;AACH,eAAO,MAAM,oBAAoB,GAAI,GAAG,MAAM,KAAG,OAOhD,CAAC"}

package/dist/categories/anomalies/asymptote-true-anomaly.js

@@ -0,0 +1,55 @@
+/**
+ * Compute the **true anomaly of a hyperbola's asymptotes**: the limit a body on a hyperbola approaches but never reaches.
+ *
+ * $$
+ * \nu_\infty = \arccos\left(-\frac{1}{e}\right)
+ * $$
+ *
+ * **Units**
+ * - Input: `e` dimensionless.
+ * - Output: $ \nu_\infty $ in **radians**, in (π/2, π).
+ *
+ * ::: info Notes
+ *
+ * - The body comes in from $ -\nu_\infty $ and leaves towards $ +\nu_\infty $.
+ * - The trajectory turns through $ 2\nu_\infty - \pi = 2\arcsin(1/e) $, the turning angle of a flyby; see `gravityAssistTurningAngle`.
+ * - Just above $ e = 1 $ the asymptotes open to almost π (a near-parabola); for large $ e $ they close towards π/2 (a near-straight line).
+ *
+ * :::
+ *
+ * ::: practice Apollo 11, arriving at the Moon
+ *
+ * **134.0°**
+ *
+ * _Asymptote of the approach hyperbola, e = 1.44_
+ *
+ * Without the lunar orbit insertion burn at closest approach, the Moon's gravity would have bent Apollo 11's path through 88° and sent it back out of the Moon's sphere of influence.
+ *
+ * ```ts
+ * asymptoteTrueAnomaly(1.44); // ≈ 2.338 rad, 134.0°
+ * ```
+ *
+ * :::
+ *
+ * @param e - Eccentricity $ e $ (finite, > 1).
+ * @returns Asymptote true anomaly $ \nu_\infty $ in **radians**.
+ * @throws {TypeError} If `e` is not a finite number.
+ * @throws {RangeError} If `e <= 1`.
+ * @group Anomalies
+ *
+ * @example
+ * ```ts
+ * asymptoteTrueAnomaly(2); // 2π/3 (120°)
+ * ```
+ *
+ * @see {@link trueToHyperbolicAnomaly}
+ * @see https://en.wikipedia.org/wiki/Hyperbolic_trajectory Hyperbolic trajectory
+ */
+export const asymptoteTrueAnomaly = (e) => {
+    if (!Number.isFinite(e))
+        throw new TypeError(`e must be a finite number. Received: ${e}`);
+    if (e <= 1)
+        throw new RangeError(`e must be > 1. Received: ${e}`);
+    return Math.acos(-1 / e);
+};
+//# sourceMappingURL=asymptote-true-anomaly.js.map

package/dist/categories/anomalies/asymptote-true-anomaly.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"asymptote-true-anomaly.js","sourceRoot":"","sources":["../../../src/categories/anomalies/asymptote-true-anomaly.ts"],"names":[],"mappings":"AAEA;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GA8CG;AACH,MAAM,CAAC,MAAM,oBAAoB,GAAG,CAAC,CAAS,EAAW,EAAE;IACzD,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC;QACrB,MAAM,IAAI,SAAS,CAAC,wCAAwC,CAAC,EAAE,CAAC,CAAC;IAEnE,IAAI,CAAC,IAAI,CAAC;QAAE,MAAM,IAAI,UAAU,CAAC,4BAA4B,CAAC,EAAE,CAAC,CAAC;IAElE,OAAO,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,CAAY,CAAC;AACtC,CAAC,CAAC"}

package/dist/categories/anomalies/hyperbolic-to-true-anomaly.d.ts

@@ -0,0 +1,38 @@
+import type { Radians } from '@interstellar-tools/types';
+/**
+ * Convert a **hyperbolic anomaly** $ F $ to the **true anomaly** on a hyperbola.
+ *
+ * $$
+ * \tan\frac{\nu}{2} = \sqrt{\frac{e + 1}{e - 1}}\,\tanh\frac{F}{2}
+ * $$
+ *
+ * **Units**
+ * - Inputs: `F` in **radians**, `e` dimensionless.
+ * - Output: $ \nu $ in **radians**, in $ (-\nu_\infty, \nu_\infty) $.
+ *
+ * ::: info Notes
+ *
+ * - The inverse of {@link trueToHyperbolicAnomaly}.
+ * - As $ F \to \pm\infty $, $ \nu $ approaches the asymptotes $ \pm\nu_\infty $ but never reaches them.
+ * - Returns a signed angle, negative before periapsis. Use `norm2pi` for [0, 2π).
+ *
+ * :::
+ *
+ * @param F - Hyperbolic anomaly $ F $ in **radians** (finite).
+ * @param e - Eccentricity $ e $ (finite, > 1).
+ * @returns True anomaly $ \nu $ in **radians**, in (-π, π).
+ * @throws {TypeError} If any input is not a finite number.
+ * @throws {RangeError} If `e <= 1`.
+ * @group Anomalies
+ *
+ * @example
+ * ```ts
+ * const nu = hyperbolicToTrueAnomaly(Math.acosh(2), 2);
+ * // nu ≈ π/2
+ * ```
+ *
+ * @see {@link trueToHyperbolicAnomaly}
+ * @see {@link solveKeplerHyperbolic}
+ */
+export declare const hyperbolicToTrueAnomaly: (F: number, e: number) => Radians;
+//# sourceMappingURL=hyperbolic-to-true-anomaly.d.ts.map

package/dist/categories/anomalies/hyperbolic-to-true-anomaly.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"hyperbolic-to-true-anomaly.d.ts","sourceRoot":"","sources":["../../../src/categories/anomalies/hyperbolic-to-true-anomaly.ts"],"names":[],"mappings":"AAAA,OAAO,KAAK,EAAE,OAAO,EAAE,MAAM,2BAA2B,CAAC;AAEzD;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GAkCG;AACH,eAAO,MAAM,uBAAuB,GAAI,GAAG,MAAM,EAAE,GAAG,MAAM,KAAG,OAa9D,CAAC"}

package/dist/categories/anomalies/hyperbolic-to-true-anomaly.js

@@ -0,0 +1,49 @@
+/**
+ * Convert a **hyperbolic anomaly** $ F $ to the **true anomaly** on a hyperbola.
+ *
+ * $$
+ * \tan\frac{\nu}{2} = \sqrt{\frac{e + 1}{e - 1}}\,\tanh\frac{F}{2}
+ * $$
+ *
+ * **Units**
+ * - Inputs: `F` in **radians**, `e` dimensionless.
+ * - Output: $ \nu $ in **radians**, in $ (-\nu_\infty, \nu_\infty) $.
+ *
+ * ::: info Notes
+ *
+ * - The inverse of {@link trueToHyperbolicAnomaly}.
+ * - As $ F \to \pm\infty $, $ \nu $ approaches the asymptotes $ \pm\nu_\infty $ but never reaches them.
+ * - Returns a signed angle, negative before periapsis. Use `norm2pi` for [0, 2π).
+ *
+ * :::
+ *
+ * @param F - Hyperbolic anomaly $ F $ in **radians** (finite).
+ * @param e - Eccentricity $ e $ (finite, > 1).
+ * @returns True anomaly $ \nu $ in **radians**, in (-π, π).
+ * @throws {TypeError} If any input is not a finite number.
+ * @throws {RangeError} If `e <= 1`.
+ * @group Anomalies
+ *
+ * @example
+ * ```ts
+ * const nu = hyperbolicToTrueAnomaly(Math.acosh(2), 2);
+ * // nu ≈ π/2
+ * ```
+ *
+ * @see {@link trueToHyperbolicAnomaly}
+ * @see {@link solveKeplerHyperbolic}
+ */
+export const hyperbolicToTrueAnomaly = (F, e) => {
+    for (const [name, x] of [
+        ['F', F],
+        ['e', e]
+    ]) {
+        if (!Number.isFinite(x))
+            throw new TypeError(`${name} must be a finite number. Received: ${x}`);
+    }
+    if (e <= 1)
+        throw new RangeError(`e must be > 1. Received: ${e}`);
+    return (2 *
+        Math.atan(Math.sqrt((e + 1) / (e - 1)) * Math.tanh(F / 2)));
+};
+//# sourceMappingURL=hyperbolic-to-true-anomaly.js.map

package/dist/categories/anomalies/hyperbolic-to-true-anomaly.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"hyperbolic-to-true-anomaly.js","sourceRoot":"","sources":["../../../src/categories/anomalies/hyperbolic-to-true-anomaly.ts"],"names":[],"mappings":"AAEA;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GAkCG;AACH,MAAM,CAAC,MAAM,uBAAuB,GAAG,CAAC,CAAS,EAAE,CAAS,EAAW,EAAE;IACvE,KAAK,MAAM,CAAC,IAAI,EAAE,CAAC,CAAC,IAAI;QACtB,CAAC,GAAG,EAAE,CAAC,CAAC;QACR,CAAC,GAAG,EAAE,CAAC,CAAC;KACA,EAAE,CAAC;QACX,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC;YACrB,MAAM,IAAI,SAAS,CAAC,GAAG,IAAI,uCAAuC,CAAC,EAAE,CAAC,CAAC;IAC3E,CAAC;IAED,IAAI,CAAC,IAAI,CAAC;QAAE,MAAM,IAAI,UAAU,CAAC,4BAA4B,CAAC,EAAE,CAAC,CAAC;IAElE,OAAO,CAAC,CAAC;QACP,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAY,CAAC;AAC3E,CAAC,CAAC"}

package/dist/categories/anomalies/true-to-hyperbolic-anomaly.d.ts

@@ -0,0 +1,40 @@
+import type { Radians } from '@interstellar-tools/types';
+/**
+ * Convert a **true anomaly** on a hyperbola to the **hyperbolic anomaly** $ F $.
+ *
+ * $$
+ * \tanh\frac{F}{2} = \sqrt{\frac{e - 1}{e + 1}}\,\tan\frac{\nu}{2}
+ * $$
+ *
+ * **Units**
+ * - Inputs: `nu` in **radians**, `e` dimensionless.
+ * - Output: $ F $ in **radians** (dimensionless, unbounded).
+ *
+ * ::: info Notes
+ *
+ * - The hyperbolic counterpart of the eccentric anomaly. Feed it to $ M_h = e \sinh F - F $ to get a time from periapsis.
+ * - `nu` is wrapped to (-π, π] first, so the incoming branch can be given either as a negative angle or just below 2π, as {@link stateToElements} returns it.
+ * - Only $ |\nu| < \nu_\infty $ exists on a hyperbola; see {@link asymptoteTrueAnomaly}.
+ *
+ * :::
+ *
+ * @param nu - True anomaly $ \nu $ in **radians** (finite, inside the asymptotes).
+ * @param e - Eccentricity $ e $ (finite, > 1).
+ * @returns Hyperbolic anomaly $ F $ in **radians**: negative before periapsis, positive after.
+ * @throws {TypeError} If any input is not a finite number.
+ * @throws {RangeError} If `e <= 1` or `nu` is at or beyond the asymptotes.
+ * @group Anomalies
+ *
+ * @example
+ * ```ts
+ * import type { Radians } from '@interstellar-tools/types';
+ *
+ * const F = trueToHyperbolicAnomaly((Math.PI / 2) as Radians, 2);
+ * // F ≈ 1.317 (cosh F = e)
+ * ```
+ *
+ * @see {@link hyperbolicToTrueAnomaly}
+ * @see https://en.wikipedia.org/wiki/Hyperbolic_trajectory#Hyperbolic_anomaly Hyperbolic anomaly
+ */
+export declare const trueToHyperbolicAnomaly: (nu: Radians, e: number) => number;
+//# sourceMappingURL=true-to-hyperbolic-anomaly.d.ts.map

package/dist/categories/anomalies/true-to-hyperbolic-anomaly.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"true-to-hyperbolic-anomaly.d.ts","sourceRoot":"","sources":["../../../src/categories/anomalies/true-to-hyperbolic-anomaly.ts"],"names":[],"mappings":"AAAA,OAAO,KAAK,EAAE,OAAO,EAAE,MAAM,2BAA2B,CAAC;AAIzD;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GAoCG;AACH,eAAO,MAAM,uBAAuB,GAAI,IAAI,OAAO,EAAE,GAAG,MAAM,KAAG,MAmBhE,CAAC"}

package/dist/categories/anomalies/true-to-hyperbolic-anomaly.js

@@ -0,0 +1,54 @@
+import { asymptoteTrueAnomaly } from './asymptote-true-anomaly.js';
+/**
+ * Convert a **true anomaly** on a hyperbola to the **hyperbolic anomaly** $ F $.
+ *
+ * $$
+ * \tanh\frac{F}{2} = \sqrt{\frac{e - 1}{e + 1}}\,\tan\frac{\nu}{2}
+ * $$
+ *
+ * **Units**
+ * - Inputs: `nu` in **radians**, `e` dimensionless.
+ * - Output: $ F $ in **radians** (dimensionless, unbounded).
+ *
+ * ::: info Notes
+ *
+ * - The hyperbolic counterpart of the eccentric anomaly. Feed it to $ M_h = e \sinh F - F $ to get a time from periapsis.
+ * - `nu` is wrapped to (-π, π] first, so the incoming branch can be given either as a negative angle or just below 2π, as {@link stateToElements} returns it.
+ * - Only $ |\nu| < \nu_\infty $ exists on a hyperbola; see {@link asymptoteTrueAnomaly}.
+ *
+ * :::
+ *
+ * @param nu - True anomaly $ \nu $ in **radians** (finite, inside the asymptotes).
+ * @param e - Eccentricity $ e $ (finite, > 1).
+ * @returns Hyperbolic anomaly $ F $ in **radians**: negative before periapsis, positive after.
+ * @throws {TypeError} If any input is not a finite number.
+ * @throws {RangeError} If `e <= 1` or `nu` is at or beyond the asymptotes.
+ * @group Anomalies
+ *
+ * @example
+ * ```ts
+ * import type { Radians } from '@interstellar-tools/types';
+ *
+ * const F = trueToHyperbolicAnomaly((Math.PI / 2) as Radians, 2);
+ * // F ≈ 1.317 (cosh F = e)
+ * ```
+ *
+ * @see {@link hyperbolicToTrueAnomaly}
+ * @see https://en.wikipedia.org/wiki/Hyperbolic_trajectory#Hyperbolic_anomaly Hyperbolic anomaly
+ */
+export const trueToHyperbolicAnomaly = (nu, e) => {
+    for (const [name, x] of [
+        ['nu', nu],
+        ['e', e]
+    ]) {
+        if (!Number.isFinite(x))
+            throw new TypeError(`${name} must be a finite number. Received: ${x}`);
+    }
+    if (e <= 1)
+        throw new RangeError(`e must be > 1. Received: ${e}`);
+    const wrapped = Math.atan2(Math.sin(nu), Math.cos(nu));
+    if (Math.abs(wrapped) >= asymptoteTrueAnomaly(e))
+        throw new RangeError(`nu is at or beyond the asymptotes (|nu| >= acos(-1/e)). Received: nu=${nu}, e=${e}`);
+    return 2 * Math.atanh(Math.sqrt((e - 1) / (e + 1)) * Math.tan(wrapped / 2));
+};
+//# sourceMappingURL=true-to-hyperbolic-anomaly.js.map

package/dist/categories/anomalies/true-to-hyperbolic-anomaly.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"true-to-hyperbolic-anomaly.js","sourceRoot":"","sources":["../../../src/categories/anomalies/true-to-hyperbolic-anomaly.ts"],"names":[],"mappings":"AAEA,OAAO,EAAE,oBAAoB,EAAE,MAAM,6BAA6B,CAAC;AAEnE;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GAoCG;AACH,MAAM,CAAC,MAAM,uBAAuB,GAAG,CAAC,EAAW,EAAE,CAAS,EAAU,EAAE;IACxE,KAAK,MAAM,CAAC,IAAI,EAAE,CAAC,CAAC,IAAI;QACtB,CAAC,IAAI,EAAE,EAAE,CAAC;QACV,CAAC,GAAG,EAAE,CAAC,CAAC;KACA,EAAE,CAAC;QACX,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC;YACrB,MAAM,IAAI,SAAS,CAAC,GAAG,IAAI,uCAAuC,CAAC,EAAE,CAAC,CAAC;IAC3E,CAAC;IAED,IAAI,CAAC,IAAI,CAAC;QAAE,MAAM,IAAI,UAAU,CAAC,4BAA4B,CAAC,EAAE,CAAC,CAAC;IAElE,MAAM,OAAO,GAAG,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC;IAEvD,IAAI,IAAI,CAAC,GAAG,CAAC,OAAO,CAAC,IAAI,oBAAoB,CAAC,CAAC,CAAC;QAC9C,MAAM,IAAI,UAAU,CAClB,wEAAwE,EAAE,OAAO,CAAC,EAAE,CACrF,CAAC;IAEJ,OAAO,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,OAAO,GAAG,CAAC,CAAC,CAAC,CAAC;AAC9E,CAAC,CAAC"}

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

@@ -0,0 +1,49 @@
+/**
+ * Solve **Kepler's equation for a hyperbola**: find the hyperbolic anomaly $ F $ from the hyperbolic mean anomaly $ M_h $.
+ *
+ * $$
+ * M_h = e \sinh F - F
+ * $$
+ *
+ * solved with Newton's method:
+ *
+ * $$
+ * F_{n+1} = F_n - \frac{e \sinh F_n - F_n - M_h}{e \cosh F_n - 1}, \qquad F_0 = \min\left(\operatorname{asinh}\frac{|M_h|}{e - 1},\ \sqrt[3]{\frac{6|M_h|}{e}}\right)
+ * $$
+ *
+ * **Units**
+ * - Inputs: `Mh` in **radians**, `e` dimensionless.
+ * - Output: $ F $ in **radians** (dimensionless, unbounded).
+ *
+ * ::: info Notes
+ *
+ * - The equation is odd in $ F $, so the solver works on $ |M_h| $ and restores the sign. Negative values are before periapsis.
+ * - $ F_0 $ is an upper bound on the root, since $ e\sinh F - F $ is at least $ (e-1)\sinh F $ and at least $ (e-1)F + eF^3/6 $. The function is increasing and convex for $ F > 0 $, so Newton's method converges monotonically from there, even close to $ e = 1 $ where the derivative $ e\cosh F - 1 $ is small near periapsis.
+ * - Unlike an ellipse, a hyperbola isn't periodic: $ M_h $ isn't reduced to one turn.
+ * - {@link solveKeplerHighEccentricity} also accepts $ e > 1 $, but this solver throws when it doesn't converge instead of returning its last estimate.
+ *
+ * :::
+ *
+ * @param Mh - Hyperbolic mean anomaly $ M_h $ in **radians** (finite).
+ * @param e - Eccentricity $ e $ (finite, > 1).
+ * @param tolerance - Stop when a step is smaller than `tolerance * max(1, |F|)`. Defaults to `1e-14`.
+ * @param maxIter - Maximum number of Newton steps. Defaults to `50`.
+ * @returns The hyperbolic anomaly $ F $ in **radians**, with the sign of `Mh`.
+ * @throws {TypeError} If any input is not a finite number.
+ * @throws {RangeError} If `e <= 1`, `tolerance <= 0` or `maxIter < 1`.
+ * @throws {Error} If the iteration doesn't converge within `maxIter` steps.
+ * @group Kepler Solvers
+ *
+ * @example
+ * ```ts
+ * import type { Radians } from '@interstellar-tools/types';
+ *
+ * const F = solveKeplerHyperbolic(12.66 as Radians, 1.44);
+ * // F ≈ 3.087, and 1.44 * sinh(F) - F ≈ 12.66
+ * ```
+ *
+ * @see {@link hyperbolicTimeOfFlight}
+ * @see https://en.wikipedia.org/wiki/Hyperbolic_trajectory#Hyperbolic_anomaly Hyperbolic anomaly
+ */
+export declare const solveKeplerHyperbolic: (Mh: number, e: number, tolerance?: number, maxIter?: number) => number;
+//# sourceMappingURL=solve-kepler-hyperbolic.d.ts.map

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

@@ -0,0 +1 @@
+{"version":3,"file":"solve-kepler-hyperbolic.d.ts","sourceRoot":"","sources":["../../../src/categories/kepler/solve-kepler-hyperbolic.ts"],"names":[],"mappings":"AAiBA;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GA8CG;AACH,eAAO,MAAM,qBAAqB,GAChC,IAAI,MAAM,EACV,GAAG,MAAM,EACT,kBAAiB,EACjB,gBAAY,KACX,MAsCF,CAAC"}

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

@@ -0,0 +1,90 @@
+/** sinh x - x, from its Taylor series for small x to avoid cancellation. */
+const sinhMinusX = (x) => {
+    if (Math.abs(x) > 0.5)
+        return Math.sinh(x) - x;
+    const x2 = x * x;
+    let term = (x * x2) / 6;
+    let sum = term;
+    for (let k = 5; Math.abs(term) > 1e-17 * Math.abs(sum); k += 2) {
+        term *= x2 / ((k - 1) * k);
+        sum += term;
+    }
+    return sum;
+};
+/**
+ * Solve **Kepler's equation for a hyperbola**: find the hyperbolic anomaly $ F $ from the hyperbolic mean anomaly $ M_h $.
+ *
+ * $$
+ * M_h = e \sinh F - F
+ * $$
+ *
+ * solved with Newton's method:
+ *
+ * $$
+ * F_{n+1} = F_n - \frac{e \sinh F_n - F_n - M_h}{e \cosh F_n - 1}, \qquad F_0 = \min\left(\operatorname{asinh}\frac{|M_h|}{e - 1},\ \sqrt[3]{\frac{6|M_h|}{e}}\right)
+ * $$
+ *
+ * **Units**
+ * - Inputs: `Mh` in **radians**, `e` dimensionless.
+ * - Output: $ F $ in **radians** (dimensionless, unbounded).
+ *
+ * ::: info Notes
+ *
+ * - The equation is odd in $ F $, so the solver works on $ |M_h| $ and restores the sign. Negative values are before periapsis.
+ * - $ F_0 $ is an upper bound on the root, since $ e\sinh F - F $ is at least $ (e-1)\sinh F $ and at least $ (e-1)F + eF^3/6 $. The function is increasing and convex for $ F > 0 $, so Newton's method converges monotonically from there, even close to $ e = 1 $ where the derivative $ e\cosh F - 1 $ is small near periapsis.
+ * - Unlike an ellipse, a hyperbola isn't periodic: $ M_h $ isn't reduced to one turn.
+ * - {@link solveKeplerHighEccentricity} also accepts $ e > 1 $, but this solver throws when it doesn't converge instead of returning its last estimate.
+ *
+ * :::
+ *
+ * @param Mh - Hyperbolic mean anomaly $ M_h $ in **radians** (finite).
+ * @param e - Eccentricity $ e $ (finite, > 1).
+ * @param tolerance - Stop when a step is smaller than `tolerance * max(1, |F|)`. Defaults to `1e-14`.
+ * @param maxIter - Maximum number of Newton steps. Defaults to `50`.
+ * @returns The hyperbolic anomaly $ F $ in **radians**, with the sign of `Mh`.
+ * @throws {TypeError} If any input is not a finite number.
+ * @throws {RangeError} If `e <= 1`, `tolerance <= 0` or `maxIter < 1`.
+ * @throws {Error} If the iteration doesn't converge within `maxIter` steps.
+ * @group Kepler Solvers
+ *
+ * @example
+ * ```ts
+ * import type { Radians } from '@interstellar-tools/types';
+ *
+ * const F = solveKeplerHyperbolic(12.66 as Radians, 1.44);
+ * // F ≈ 3.087, and 1.44 * sinh(F) - F ≈ 12.66
+ * ```
+ *
+ * @see {@link hyperbolicTimeOfFlight}
+ * @see https://en.wikipedia.org/wiki/Hyperbolic_trajectory#Hyperbolic_anomaly Hyperbolic anomaly
+ */
+export const solveKeplerHyperbolic = (Mh, e, tolerance = 1e-14, maxIter = 50) => {
+    for (const [name, x] of [
+        ['Mh', Mh],
+        ['e', e],
+        ['tolerance', tolerance],
+        ['maxIter', maxIter]
+    ]) {
+        if (!Number.isFinite(x))
+            throw new TypeError(`${name} must be a finite number. Received: ${x}`);
+    }
+    if (e <= 1)
+        throw new RangeError(`e must be > 1. Received: ${e}`);
+    if (tolerance <= 0)
+        throw new RangeError(`tolerance must be > 0. Received: ${tolerance}`);
+    if (maxIter < 1)
+        throw new RangeError(`maxIter must be >= 1. Received: ${maxIter}`);
+    const M = Math.abs(Mh);
+    // Start above the root, where Newton's method can't overshoot
+    let F = Math.min(Math.asinh(M / (e - 1)), Math.cbrt((6 * M) / e));
+    for (let iter = 0; iter < maxIter; iter++) {
+        // e sinh F - F - M, written so nothing cancels when e ≈ 1 and F is small
+        const f = (e - 1) * F + e * sinhMinusX(F) - M;
+        const step = f / (e * Math.cosh(F) - 1);
+        F -= step;
+        if (Math.abs(step) <= tolerance * Math.max(1, Math.abs(F)))
+            return Math.sign(Mh) * F;
+    }
+    throw new Error(`solveKeplerHyperbolic did not converge in ${maxIter} steps. Received: Mh=${Mh}, e=${e}`);
+};
+//# sourceMappingURL=solve-kepler-hyperbolic.js.map

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

@@ -0,0 +1 @@
+{"version":3,"file":"solve-kepler-hyperbolic.js","sourceRoot":"","sources":["../../../src/categories/kepler/solve-kepler-hyperbolic.ts"],"names":[],"mappings":"AAAA,4EAA4E;AAC5E,MAAM,UAAU,GAAG,CAAC,CAAS,EAAU,EAAE;IACvC,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,GAAG;QAAE,OAAO,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;IAE/C,MAAM,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC;IAEjB,IAAI,IAAI,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC;IACxB,IAAI,GAAG,GAAG,IAAI,CAAC;IAEf,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC;QAC/D,IAAI,IAAI,EAAE,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;QAC3B,GAAG,IAAI,IAAI,CAAC;IACd,CAAC;IAED,OAAO,GAAG,CAAC;AACb,CAAC,CAAC;AAEF;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GA8CG;AACH,MAAM,CAAC,MAAM,qBAAqB,GAAG,CACnC,EAAU,EACV,CAAS,EACT,SAAS,GAAG,KAAK,EACjB,OAAO,GAAG,EAAE,EACJ,EAAE;IACV,KAAK,MAAM,CAAC,IAAI,EAAE,CAAC,CAAC,IAAI;QACtB,CAAC,IAAI,EAAE,EAAE,CAAC;QACV,CAAC,GAAG,EAAE,CAAC,CAAC;QACR,CAAC,WAAW,EAAE,SAAS,CAAC;QACxB,CAAC,SAAS,EAAE,OAAO,CAAC;KACZ,EAAE,CAAC;QACX,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC;YACrB,MAAM,IAAI,SAAS,CAAC,GAAG,IAAI,uCAAuC,CAAC,EAAE,CAAC,CAAC;IAC3E,CAAC;IAED,IAAI,CAAC,IAAI,CAAC;QAAE,MAAM,IAAI,UAAU,CAAC,4BAA4B,CAAC,EAAE,CAAC,CAAC;IAElE,IAAI,SAAS,IAAI,CAAC;QAChB,MAAM,IAAI,UAAU,CAAC,oCAAoC,SAAS,EAAE,CAAC,CAAC;IAExE,IAAI,OAAO,GAAG,CAAC;QACb,MAAM,IAAI,UAAU,CAAC,mCAAmC,OAAO,EAAE,CAAC,CAAC;IAErE,MAAM,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC;IAEvB,8DAA8D;IAC9D,IAAI,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;IAElE,KAAK,IAAI,IAAI,GAAG,CAAC,EAAE,IAAI,GAAG,OAAO,EAAE,IAAI,EAAE,EAAE,CAAC;QAC1C,yEAAyE;QACzE,MAAM,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,UAAU,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;QAC9C,MAAM,IAAI,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;QAExC,CAAC,IAAI,IAAI,CAAC;QAEV,IAAI,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,IAAI,SAAS,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;YACxD,OAAO,IAAI,CAAC,IAAI,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC;IAC7B,CAAC;IAED,MAAM,IAAI,KAAK,CACb,6CAA6C,OAAO,wBAAwB,EAAE,OAAO,CAAC,EAAE,CACzF,CAAC;AACJ,CAAC,CAAC"}

package/dist/categories/manoeuvres/cw-two-impulse-transfer.d.ts

@@ -0,0 +1,58 @@
+import type { CwState, Vector3DTupleType } from '@interstellar-tools/types';
+export interface CwTwoImpulseTransferInterface {
+    /** First burn, at the start, in m/s (LVLH). */
+    dv1: Vector3DTupleType;
+    /** Second burn, on arrival, in m/s (LVLH). */
+    dv2: Vector3DTupleType;
+}
+/**
+ * Compute the **two burns of a Clohessy–Wiltshire rendezvous**: from a relative state to a target state in a given time.
+ *
+ * $$
+ * \mathbf v_0^+ = \Phi_{rv}^{-1}\left(\mathbf r_f - \Phi_{rr}\,\mathbf r_0\right), \quad
+ * \Delta\mathbf v_1 = \mathbf v_0^+ - \mathbf v_0, \quad
+ * \Delta\mathbf v_2 = \mathbf v_f - \left(\Phi_{vr}\,\mathbf r_0 + \Phi_{vv}\,\mathbf v_0^+\right)
+ * $$
+ *
+ * with the blocks of $ \Phi(t) $ from {@link cwStateTransition}.
+ *
+ * **Units**
+ * - Inputs: states in **m** and **m/s**, `n` in **rad/s**, `t` in **s**.
+ * - Output: both burns in **m/s**, as vectors in the LVLH frame.
+ *
+ * ::: info Notes
+ *
+ * - The first burn puts the deputy on a path that reaches the target position after `t`; the second matches the target velocity there. For a rendezvous with the chief, the target state is all zeros.
+ * - Some transfer times can't be targeted: in the orbital plane, every whole number of orbits ($ nt = 2\pi k $), and out of plane, every half orbit, where the cross-track motion returns to $ \pm z_0 $ whatever the burn. At a half orbit the cross-track burn is zero if the target is already on that path, and otherwise the function throws.
+ * - Valid near a circular reference orbit, for small separations.
+ *
+ * :::
+ *
+ * @param state0 - Deputy's relative state before the first burn `[x, y, z, xDot, yDot, zDot]` (finite).
+ * @param target - State to arrive at `[x, y, z, xDot, yDot, zDot]` (finite); all zeros to rendezvous with the chief.
+ * @param n - Mean motion of the circular reference orbit in **rad/s** (finite, > 0).
+ * @param t - Transfer time in **s** (finite, > 0).
+ * @returns The two burns `dv1` and `dv2` in **m/s**.
+ * @throws {TypeError} If any input is not a finite number.
+ * @throws {RangeError} If `n <= 0`, `t <= 0`, or the transfer time can't be targeted.
+ * @group Manoeuvres
+ *
+ * @example
+ * ```ts
+ * import type { CwState } from '@interstellar-tools/types';
+ *
+ * // 1 km behind the chief, at rest: rendezvous in a quarter of an orbit
+ * const n = 0.00113;
+ * const { dv1, dv2 } = cwTwoImpulseTransfer(
+ *   [0, -1000, 0, 0, 0, 0],
+ *   [0, 0, 0, 0, 0, 0],
+ *   n,
+ *   Math.PI / 2 / n
+ * );
+ * ```
+ *
+ * @see {@link cwPropagate}
+ * @see https://en.wikipedia.org/wiki/Clohessy%E2%80%93Wiltshire_equations Clohessy–Wiltshire equations
+ */
+export declare const cwTwoImpulseTransfer: (state0: CwState, target: CwState, n: number, t: number) => CwTwoImpulseTransferInterface;
+//# sourceMappingURL=cw-two-impulse-transfer.d.ts.map

package/dist/categories/manoeuvres/cw-two-impulse-transfer.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"cw-two-impulse-transfer.d.ts","sourceRoot":"","sources":["../../../src/categories/manoeuvres/cw-two-impulse-transfer.ts"],"names":[],"mappings":"AAAA,OAAO,KAAK,EAAE,OAAO,EAAE,iBAAiB,EAAE,MAAM,2BAA2B,CAAC;AAK5E,MAAM,WAAW,6BAA6B;IAC5C,+CAA+C;IAC/C,GAAG,EAAE,iBAAiB,CAAC;IACvB,8CAA8C;IAC9C,GAAG,EAAE,iBAAiB,CAAC;CACxB;AAED;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GAgDG;AACH,eAAO,MAAM,oBAAoB,GAC/B,QAAQ,OAAO,EACf,QAAQ,OAAO,EACf,GAAG,MAAM,EACT,GAAG,MAAM,KACR,6BA+DF,CAAC"}

package/dist/categories/manoeuvres/cw-two-impulse-transfer.js

@@ -0,0 +1,104 @@
+import { applyMatrix3 } from '../helpers/apply-matrix-3.js';
+import { cwStateTransition } from '../orbits/cw-state-transition.js';
+/**
+ * Compute the **two burns of a Clohessy–Wiltshire rendezvous**: from a relative state to a target state in a given time.
+ *
+ * $$
+ * \mathbf v_0^+ = \Phi_{rv}^{-1}\left(\mathbf r_f - \Phi_{rr}\,\mathbf r_0\right), \quad
+ * \Delta\mathbf v_1 = \mathbf v_0^+ - \mathbf v_0, \quad
+ * \Delta\mathbf v_2 = \mathbf v_f - \left(\Phi_{vr}\,\mathbf r_0 + \Phi_{vv}\,\mathbf v_0^+\right)
+ * $$
+ *
+ * with the blocks of $ \Phi(t) $ from {@link cwStateTransition}.
+ *
+ * **Units**
+ * - Inputs: states in **m** and **m/s**, `n` in **rad/s**, `t` in **s**.
+ * - Output: both burns in **m/s**, as vectors in the LVLH frame.
+ *
+ * ::: info Notes
+ *
+ * - The first burn puts the deputy on a path that reaches the target position after `t`; the second matches the target velocity there. For a rendezvous with the chief, the target state is all zeros.
+ * - Some transfer times can't be targeted: in the orbital plane, every whole number of orbits ($ nt = 2\pi k $), and out of plane, every half orbit, where the cross-track motion returns to $ \pm z_0 $ whatever the burn. At a half orbit the cross-track burn is zero if the target is already on that path, and otherwise the function throws.
+ * - Valid near a circular reference orbit, for small separations.
+ *
+ * :::
+ *
+ * @param state0 - Deputy's relative state before the first burn `[x, y, z, xDot, yDot, zDot]` (finite).
+ * @param target - State to arrive at `[x, y, z, xDot, yDot, zDot]` (finite); all zeros to rendezvous with the chief.
+ * @param n - Mean motion of the circular reference orbit in **rad/s** (finite, > 0).
+ * @param t - Transfer time in **s** (finite, > 0).
+ * @returns The two burns `dv1` and `dv2` in **m/s**.
+ * @throws {TypeError} If any input is not a finite number.
+ * @throws {RangeError} If `n <= 0`, `t <= 0`, or the transfer time can't be targeted.
+ * @group Manoeuvres
+ *
+ * @example
+ * ```ts
+ * import type { CwState } from '@interstellar-tools/types';
+ *
+ * // 1 km behind the chief, at rest: rendezvous in a quarter of an orbit
+ * const n = 0.00113;
+ * const { dv1, dv2 } = cwTwoImpulseTransfer(
+ *   [0, -1000, 0, 0, 0, 0],
+ *   [0, 0, 0, 0, 0, 0],
+ *   n,
+ *   Math.PI / 2 / n
+ * );
+ * ```
+ *
+ * @see {@link cwPropagate}
+ * @see https://en.wikipedia.org/wiki/Clohessy%E2%80%93Wiltshire_equations Clohessy–Wiltshire equations
+ */
+export const cwTwoImpulseTransfer = (state0, target, n, t) => {
+    for (const [name, state] of [
+        ['state0', state0],
+        ['target', target]
+    ]) {
+        for (const x of state)
+            if (!Number.isFinite(x))
+                throw new TypeError(`${name} must have finite components. Received: ${x}`);
+    }
+    if (!Number.isFinite(t))
+        throw new TypeError(`t must be a finite number. Received: ${t}`);
+    if (t <= 0)
+        throw new RangeError(`t must be > 0. Received: ${t}`);
+    const { rr, rv, vr, vv } = cwStateTransition(n, t);
+    const r0 = [state0[0], state0[1], state0[2]];
+    const v0 = [state0[3], state0[4], state0[5]];
+    const fromR0 = applyMatrix3(rr, r0);
+    // Position the first burn's velocity has to cover
+    const dx = target[0] - fromR0[0];
+    const dy = target[1] - fromR0[1];
+    const dz = target[2] - fromR0[2];
+    // In-plane: invert the 2×2 block of Φrv
+    const [[a, b], [c, d]] = rv;
+    const det = a * d - b * c;
+    if (Math.abs(det) * n * n <= 1e-10)
+        throw new RangeError(`The transfer time is a whole number of orbits (nt = 2πk), which can't be targeted in the orbital plane. Received: nt=${n * t}`);
+    const xDot = (d * dx - b * dy) / det;
+    const yDot = (-c * dx + a * dy) / det;
+    // Out of plane: z(t) = c z0 + (s / n) zDot0
+    const sOverN = rv[2][2];
+    let zDot;
+    if (Math.abs(sOverN) * n > 1e-10) {
+        zDot = dz / sOverN;
+    }
+    else if (Math.abs(dz) <= 1e-9 * Math.max(1, Math.abs(target[2]))) {
+        zDot = v0[2];
+    }
+    else {
+        throw new RangeError(`At a half orbit (nt = πk) the cross-track position can't be changed. Received: nt=${n * t}`);
+    }
+    const v0Plus = [xDot, yDot, zDot];
+    const fromR = applyMatrix3(vr, r0);
+    const fromV = applyMatrix3(vv, v0Plus);
+    return {
+        dv1: [xDot - v0[0], yDot - v0[1], zDot - v0[2]],
+        dv2: [
+            target[3] - (fromR[0] + fromV[0]),
+            target[4] - (fromR[1] + fromV[1]),
+            target[5] - (fromR[2] + fromV[2])
+        ]
+    };
+};
+//# sourceMappingURL=cw-two-impulse-transfer.js.map

package/dist/categories/manoeuvres/cw-two-impulse-transfer.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"cw-two-impulse-transfer.js","sourceRoot":"","sources":["../../../src/categories/manoeuvres/cw-two-impulse-transfer.ts"],"names":[],"mappings":"AAEA,OAAO,EAAE,YAAY,EAAE,MAAM,8BAA8B,CAAC;AAC5D,OAAO,EAAE,iBAAiB,EAAE,MAAM,kCAAkC,CAAC;AASrE;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GAgDG;AACH,MAAM,CAAC,MAAM,oBAAoB,GAAG,CAClC,MAAe,EACf,MAAe,EACf,CAAS,EACT,CAAS,EACsB,EAAE;IACjC,KAAK,MAAM,CAAC,IAAI,EAAE,KAAK,CAAC,IAAI;QAC1B,CAAC,QAAQ,EAAE,MAAM,CAAC;QAClB,CAAC,QAAQ,EAAE,MAAM,CAAC;KACV,EAAE,CAAC;QACX,KAAK,MAAM,CAAC,IAAI,KAAK;YACnB,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC;gBACrB,MAAM,IAAI,SAAS,CACjB,GAAG,IAAI,2CAA2C,CAAC,EAAE,CACtD,CAAC;IACR,CAAC;IAED,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC;QACrB,MAAM,IAAI,SAAS,CAAC,wCAAwC,CAAC,EAAE,CAAC,CAAC;IAEnE,IAAI,CAAC,IAAI,CAAC;QAAE,MAAM,IAAI,UAAU,CAAC,4BAA4B,CAAC,EAAE,CAAC,CAAC;IAElE,MAAM,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,GAAG,iBAAiB,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;IACnD,MAAM,EAAE,GAAsB,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC;IAChE,MAAM,EAAE,GAAsB,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC;IAChE,MAAM,MAAM,GAAG,YAAY,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC;IACpC,kDAAkD;IAClD,MAAM,EAAE,GAAG,MAAM,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC;IACjC,MAAM,EAAE,GAAG,MAAM,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC;IACjC,MAAM,EAAE,GAAG,MAAM,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC;IACjC,wCAAwC;IACxC,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC;IAC5B,MAAM,GAAG,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;IAE1B,IAAI,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,IAAI,KAAK;QAChC,MAAM,IAAI,UAAU,CAClB,wHAAwH,CAAC,GAAG,CAAC,EAAE,CAChI,CAAC;IAEJ,MAAM,IAAI,GAAG,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,GAAG,EAAE,CAAC,GAAG,GAAG,CAAC;IACrC,MAAM,IAAI,GAAG,CAAC,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,GAAG,EAAE,CAAC,GAAG,GAAG,CAAC;IACtC,4CAA4C;IAC5C,MAAM,MAAM,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IAExB,IAAI,IAAY,CAAC;IAEjB,IAAI,IAAI,CAAC,GAAG,CAAC,MAAM,CAAC,GAAG,CAAC,GAAG,KAAK,EAAE,CAAC;QACjC,IAAI,GAAG,EAAE,GAAG,MAAM,CAAC;IACrB,CAAC;SAAM,IAAI,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC;QACnE,IAAI,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC;IACf,CAAC;SAAM,CAAC;QACN,MAAM,IAAI,UAAU,CAClB,qFAAqF,CAAC,GAAG,CAAC,EAAE,CAC7F,CAAC;IACJ,CAAC;IAED,MAAM,MAAM,GAAsB,CAAC,IAAI,EAAE,IAAI,EAAE,IAAI,CAAC,CAAC;IACrD,MAAM,KAAK,GAAG,YAAY,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC;IACnC,MAAM,KAAK,GAAG,YAAY,CAAC,EAAE,EAAE,MAAM,CAAC,CAAC;IAEvC,OAAO;QACL,GAAG,EAAE,CAAC,IAAI,GAAG,EAAE,CAAC,CAAC,CAAC,EAAE,IAAI,GAAG,EAAE,CAAC,CAAC,CAAC,EAAE,IAAI,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC;QAC/C,GAAG,EAAE;YACH,MAAM,CAAC,CAAC,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC;YACjC,MAAM,CAAC,CAAC,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC;YACjC,MAAM,CAAC,CAAC,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC;SAClC;KACF,CAAC;AACJ,CAAC,CAAC"}