@interstellar-tools/equations 0.1.2 → 0.3.0

package/dist/categories/gravity/gravitational-force.d.ts

@@ -0,0 +1,94 @@
+import type { ForceResultInterface, Vector3DTupleType } from '@interstellar-tools/types';
+/**
+ * Compute the Newtonian **gravitational force** on body 1 due to body 2.
+ *
+ * **Definitions**
+ *
+ * $$
+ * \mathbf r=\mathbf r_2-\mathbf r_1,\quad
+ * r=\lVert\mathbf r\rVert,\quad
+ * \hat{\mathbf r}=\frac{\mathbf r}{r}
+ * $$
+ *
+ * **Law & Decomposition**
+ *
+ * $$
+ * \mathbf F_{1\leftarrow2}
+ * = G\,\frac{m_1 m_2}{r^2}\,\hat{\mathbf r}
+ * $$
+ *
+ * $$
+ * \lVert\mathbf F\rVert
+ * = G\,\frac{m_1 m_2}{r^2},\qquad
+ * \mathbf F=\langle F_x,F_y,F_z\rangle
+ * = \lVert\mathbf F\rVert\,\hat{\mathbf r}
+ * $$
+ *
+ * $$
+ * F_x=\lVert\mathbf F\rVert\,\hat r_x,\quad
+ * F_y=\lVert\mathbf F\rVert\,\hat r_y,\quad
+ * F_z=\lVert\mathbf F\rVert\,\hat r_z
+ * $$
+ *
+ * **Newton's 3rd law, for context**
+ *
+ * $$
+ * \mathbf F_{1\leftarrow2}=-\,\mathbf F_{2\leftarrow1}
+ * $$
+ *
+ * **Conventions & units**
+ * - Input masses in **kilograms (kg)**; positions in **meters (m)**.
+ * - Output force in **newtons (N)**.
+ * - Returned `vector` and `direction` point **from body 1 toward body 2**; i.e., `direction = \hat{\mathbf r}`.
+ *
+ * **Invariants (within numerical tolerance)**
+ * - `vector = magnitude * direction`
+ * - `magnitude = sqrt(Fx² + Fy² + Fz²)`
+ * - `||direction|| = 1`
+ *
+ *
+ * @param {number} m1 Mass of body 1 (kg).
+ * @param {number} m2 Mass of body 2 (kg).
+ * @param {Vector3DTupleType} r1 Position of body 1 as a 3-tuple `[x, y, z]` in meters.
+ * @param {Vector3DTupleType} r2 Position of body 2 as a 3-tuple `[x, y, z]` in meters.
+ * @param {number} [G]  Gravitational constant (defaults to `G_SI`, m³·kg⁻¹·s⁻²).
+ * @returns {ForceResultInterface} `{ vector, magnitude, direction }`.
+ * @throws {@link https://developer.mozilla.org/en-US/docs/Web/JavaScript/Reference/Global_Objects/Error | Error} If `m1` or `m2` are non-finite or negative.
+ * @throws {@link https://developer.mozilla.org/en-US/docs/Web/JavaScript/Reference/Global_Objects/Error | Error} If positions coincide (`r = 0`).
+ * @throws {@link https://developer.mozilla.org/en-US/docs/Web/JavaScript/Reference/Global_Objects/Error | Error} If `G` is non-finite or not positive.
+ *
+ * @example
+ * ```ts
+ * import { gravitationalForce } from "@interstellar-tools/equations";
+ * import { G_SI } from "@interstellar-tools/constants";
+ *
+ * // Earth–Sun at ~1 AU
+ * const mEarth = 5.972e24;     // kg
+ * const mSun   = 1.9885e30;    // kg
+ * const rEarth: [number, number, number] = [0, 0, 0];
+ * const rSun  : [number, number, number] = [1.495978707e11, 0, 0]; // m
+ *
+ * const F = gravitationalForce(mEarth, mSun, rEarth, rSun, G_SI);
+ * console.log(F.magnitude.toExponential(3)); // ≈ "3.542e+22" N
+ * console.log(F.vector); // Points from Earth toward Sun (−x here)
+ * ```
+ *
+ * **Get F directly (recommended)**
+ *
+ * ```ts
+ * import { gravitationalForce } from "@interstellar-tools/equations";
+ * import type { Vector3DTupleType, ForceResultInterface } from "@interstellar-tools/types";
+ *
+ * const m1 = 5.972e24; // kg
+ * const m2 = 1.9885e30; // kg
+ * const r1: Vector3DTupleType = [0, 0, 0];                 // m
+ * const r2: Vector3DTupleType = [1.495978707e11, 0, 0];    // m
+ *
+ * const { vector: F, magnitude, direction }: ForceResultInterface = gravitationalForce(m1, m2, r1, r2);
+ * ```
+ *
+ * @see https://en.wikipedia.org/wiki/Newton%27s_law_of_universal_gravitation
+ * @group Gravity
+ */
+export declare const gravitationalForce: (m1: number, m2: number, r1: Vector3DTupleType, r2: Vector3DTupleType, G?: number) => ForceResultInterface;
+//# sourceMappingURL=gravitational-force.d.ts.map

package/dist/categories/gravity/gravitational-force.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"gravitational-force.d.ts","sourceRoot":"","sources":["../../../src/categories/gravity/gravitational-force.ts"],"names":[],"mappings":"AACA,OAAO,KAAK,EACV,oBAAoB,EACpB,iBAAiB,EAClB,MAAM,2BAA2B,CAAC;AAEnC;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GA0FG;AACH,eAAO,MAAM,kBAAkB,GAC7B,IAAI,MAAM,EACV,IAAI,MAAM,EACV,IAAI,iBAAiB,EACrB,IAAI,iBAAiB,EACrB,IAAG,MAAa,KACf,oBAgCF,CAAC"}

package/dist/categories/gravity/gravitational-force.js

@@ -0,0 +1,119 @@
+import { G_SI } from '@interstellar-tools/constants';
+/**
+ * Compute the Newtonian **gravitational force** on body 1 due to body 2.
+ *
+ * **Definitions**
+ *
+ * $$
+ * \mathbf r=\mathbf r_2-\mathbf r_1,\quad
+ * r=\lVert\mathbf r\rVert,\quad
+ * \hat{\mathbf r}=\frac{\mathbf r}{r}
+ * $$
+ *
+ * **Law & Decomposition**
+ *
+ * $$
+ * \mathbf F_{1\leftarrow2}
+ * = G\,\frac{m_1 m_2}{r^2}\,\hat{\mathbf r}
+ * $$
+ *
+ * $$
+ * \lVert\mathbf F\rVert
+ * = G\,\frac{m_1 m_2}{r^2},\qquad
+ * \mathbf F=\langle F_x,F_y,F_z\rangle
+ * = \lVert\mathbf F\rVert\,\hat{\mathbf r}
+ * $$
+ *
+ * $$
+ * F_x=\lVert\mathbf F\rVert\,\hat r_x,\quad
+ * F_y=\lVert\mathbf F\rVert\,\hat r_y,\quad
+ * F_z=\lVert\mathbf F\rVert\,\hat r_z
+ * $$
+ *
+ * **Newton's 3rd law, for context**
+ *
+ * $$
+ * \mathbf F_{1\leftarrow2}=-\,\mathbf F_{2\leftarrow1}
+ * $$
+ *
+ * **Conventions & units**
+ * - Input masses in **kilograms (kg)**; positions in **meters (m)**.
+ * - Output force in **newtons (N)**.
+ * - Returned `vector` and `direction` point **from body 1 toward body 2**; i.e., `direction = \hat{\mathbf r}`.
+ *
+ * **Invariants (within numerical tolerance)**
+ * - `vector = magnitude * direction`
+ * - `magnitude = sqrt(Fx² + Fy² + Fz²)`
+ * - `||direction|| = 1`
+ *
+ *
+ * @param {number} m1 Mass of body 1 (kg).
+ * @param {number} m2 Mass of body 2 (kg).
+ * @param {Vector3DTupleType} r1 Position of body 1 as a 3-tuple `[x, y, z]` in meters.
+ * @param {Vector3DTupleType} r2 Position of body 2 as a 3-tuple `[x, y, z]` in meters.
+ * @param {number} [G]  Gravitational constant (defaults to `G_SI`, m³·kg⁻¹·s⁻²).
+ * @returns {ForceResultInterface} `{ vector, magnitude, direction }`.
+ * @throws {@link https://developer.mozilla.org/en-US/docs/Web/JavaScript/Reference/Global_Objects/Error | Error} If `m1` or `m2` are non-finite or negative.
+ * @throws {@link https://developer.mozilla.org/en-US/docs/Web/JavaScript/Reference/Global_Objects/Error | Error} If positions coincide (`r = 0`).
+ * @throws {@link https://developer.mozilla.org/en-US/docs/Web/JavaScript/Reference/Global_Objects/Error | Error} If `G` is non-finite or not positive.
+ *
+ * @example
+ * ```ts
+ * import { gravitationalForce } from "@interstellar-tools/equations";
+ * import { G_SI } from "@interstellar-tools/constants";
+ *
+ * // Earth–Sun at ~1 AU
+ * const mEarth = 5.972e24;     // kg
+ * const mSun   = 1.9885e30;    // kg
+ * const rEarth: [number, number, number] = [0, 0, 0];
+ * const rSun  : [number, number, number] = [1.495978707e11, 0, 0]; // m
+ *
+ * const F = gravitationalForce(mEarth, mSun, rEarth, rSun, G_SI);
+ * console.log(F.magnitude.toExponential(3)); // ≈ "3.542e+22" N
+ * console.log(F.vector); // Points from Earth toward Sun (−x here)
+ * ```
+ *
+ * **Get F directly (recommended)**
+ *
+ * ```ts
+ * import { gravitationalForce } from "@interstellar-tools/equations";
+ * import type { Vector3DTupleType, ForceResultInterface } from "@interstellar-tools/types";
+ *
+ * const m1 = 5.972e24; // kg
+ * const m2 = 1.9885e30; // kg
+ * const r1: Vector3DTupleType = [0, 0, 0];                 // m
+ * const r2: Vector3DTupleType = [1.495978707e11, 0, 0];    // m
+ *
+ * const { vector: F, magnitude, direction }: ForceResultInterface = gravitationalForce(m1, m2, r1, r2);
+ * ```
+ *
+ * @see https://en.wikipedia.org/wiki/Newton%27s_law_of_universal_gravitation
+ * @group Gravity
+ */
+export const gravitationalForce = (m1, m2, r1, r2, G = G_SI) => {
+    if (!Number.isFinite(m1) || !Number.isFinite(m2) || m1 < 0 || m2 < 0) {
+        throw new Error('m1 and m2 must be finite, non-negative numbers (kg).');
+    }
+    if (!Number.isFinite(G) || G <= 0) {
+        throw new Error('G must be a finite, positive number (m^3·kg^-1·s^-2).');
+    }
+    const rx = r2[0] - r1[0];
+    const ry = r2[1] - r1[1];
+    const rz = r2[2] - r1[2];
+    const r2mag = rx * rx + ry * ry + rz * rz;
+    if (r2mag === 0) {
+        throw new Error('Bodies are at the same position: gravitational force is singular (r = 0).');
+    }
+    const invR = 1 / Math.sqrt(r2mag);
+    const dir = [rx * invR, ry * invR, rz * invR];
+    // |F| = G m1 m2 / r^2
+    const magnitude = (G * m1 * m2) / r2mag;
+    // Vector points from body 1 toward body 2.
+    const Fx = dir[0] * magnitude;
+    const Fy = dir[1] * magnitude;
+    const Fz = dir[2] * magnitude;
+    // Ensure the tuple type is preserved for callers.
+    const vector = [Fx, Fy, Fz];
+    return { vector, magnitude, direction: dir };
+};
+//# sourceMappingURL=gravitational-force.js.map

package/dist/categories/gravity/gravitational-force.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"gravitational-force.js","sourceRoot":"","sources":["../../../src/categories/gravity/gravitational-force.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,IAAI,EAAE,MAAM,+BAA+B,CAAC;AAMrD;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GA0FG;AACH,MAAM,CAAC,MAAM,kBAAkB,GAAG,CAChC,EAAU,EACV,EAAU,EACV,EAAqB,EACrB,EAAqB,EACrB,IAAY,IAAI,EACM,EAAE;IACxB,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,EAAE,CAAC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,EAAE,CAAC,IAAI,EAAE,GAAG,CAAC,IAAI,EAAE,GAAG,CAAC,EAAE,CAAC;QACrE,MAAM,IAAI,KAAK,CAAC,sDAAsD,CAAC,CAAC;IAC1E,CAAC;IAED,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,IAAI,CAAC,IAAI,CAAC,EAAE,CAAC;QAClC,MAAM,IAAI,KAAK,CAAC,uDAAuD,CAAC,CAAC;IAC3E,CAAC;IAED,MAAM,EAAE,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC;IACzB,MAAM,EAAE,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC;IACzB,MAAM,EAAE,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC;IACzB,MAAM,KAAK,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC;IAE1C,IAAI,KAAK,KAAK,CAAC,EAAE,CAAC;QAChB,MAAM,IAAI,KAAK,CACb,2EAA2E,CAC5E,CAAC;IACJ,CAAC;IAED,MAAM,IAAI,GAAG,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC;IAClC,MAAM,GAAG,GAAsB,CAAC,EAAE,GAAG,IAAI,EAAE,EAAE,GAAG,IAAI,EAAE,EAAE,GAAG,IAAI,CAAC,CAAC;IACjE,sBAAsB;IACtB,MAAM,SAAS,GAAG,CAAC,CAAC,GAAG,EAAE,GAAG,EAAE,CAAC,GAAG,KAAK,CAAC;IACxC,2CAA2C;IAC3C,MAAM,EAAE,GAAG,GAAG,CAAC,CAAC,CAAC,GAAG,SAAS,CAAC;IAC9B,MAAM,EAAE,GAAG,GAAG,CAAC,CAAC,CAAC,GAAG,SAAS,CAAC;IAC9B,MAAM,EAAE,GAAG,GAAG,CAAC,CAAC,CAAC,GAAG,SAAS,CAAC;IAC9B,kDAAkD;IAClD,MAAM,MAAM,GAAsB,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC;IAE/C,OAAO,EAAE,MAAM,EAAE,SAAS,EAAE,SAAS,EAAE,GAAG,EAAE,CAAC;AAC/C,CAAC,CAAC"}

package/dist/categories/gravity/gravitational-parameter.d.ts

@@ -0,0 +1,49 @@
+/**
+ * Two-body **gravitational parameter** ($ \mu $).
+ *
+ * **Definition**
+ *
+ * $$
+ * \mu = G\,(M + m)
+ * $$
+ *
+ * **Common approximation, when** ($ M \gg m $)
+ *
+ * $$
+ * \mu \approx G\,M
+ * $$
+ *
+ * **Units**
+ * - Inputs: masses in **kg**, ( $ G $ ) in **m³·kg⁻¹·s⁻²** (defaults to `G_SI`).
+ * - Output: ($ \mu $) in **m³/s²**.
+ *
+ * ::: tip
+ *
+ * - Use `m = 0` (default) for the standard published ($ \mu $) of a central body.
+ * - When modelling a binary with comparable masses, pass both `M` and `m`.
+ *
+ * :::
+ *
+ * @param {number} M Mass of the primary body (kg).
+ * @param {number} m Mass of the secondary body (kg). Defaults to `0`.
+ * @param {number} [G] Gravitational constant. Defaults to `G_SI`.
+ * @returns {number} Gravitational parameter ($ \mu $) in m³/s².
+ * @throws {@link https://developer.mozilla.org/en-US/docs/Web/JavaScript/Reference/Global_Objects/Error | Error} If `M` is non-finite or negative.
+ * @throws {@link https://developer.mozilla.org/en-US/docs/Web/JavaScript/Reference/Global_Objects/Error | Error} If `m` is non-finite or negative.
+ * @throws {@link https://developer.mozilla.org/en-US/docs/Web/JavaScript/Reference/Global_Objects/Error | Error} If `G` is non-finite or negative.
+ *
+ * @example
+ * ```ts
+ * // Standard gravitational parameter for Earth (approximate: m ~ 0)
+ * const muEarth = gravitationalParameter(5.97219e24); // ≈ 3.986004e14 m^3/s^2
+ *
+ * // Earth + 1000 kg satellite: practically identical to GM
+ * const muExact = gravitationalParameter(5.97219e24, 1000);
+ *
+ * // Earth–Moon system (two-body μ)
+ * const muEarthMoon = gravitationalParameter(5.97219e24, 7.342e22);
+ * ```
+ * @group Gravity
+ */
+export declare const gravitationalParameter: (M: number, m?: number, G?: number) => number;
+//# sourceMappingURL=gravitational-parameter.d.ts.map

package/dist/categories/gravity/gravitational-parameter.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"gravitational-parameter.d.ts","sourceRoot":"","sources":["../../../src/categories/gravity/gravitational-parameter.ts"],"names":[],"mappings":"AAEA;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GA8CG;AACH,eAAO,MAAM,sBAAsB,GACjC,GAAG,MAAM,EACT,IAAG,MAAU,EACb,IAAG,MAAa,KACf,MAcF,CAAC"}

package/dist/categories/gravity/gravitational-parameter.js

@@ -0,0 +1,61 @@
+import { G_SI } from '@interstellar-tools/constants';
+/**
+ * Two-body **gravitational parameter** ($ \mu $).
+ *
+ * **Definition**
+ *
+ * $$
+ * \mu = G\,(M + m)
+ * $$
+ *
+ * **Common approximation, when** ($ M \gg m $)
+ *
+ * $$
+ * \mu \approx G\,M
+ * $$
+ *
+ * **Units**
+ * - Inputs: masses in **kg**, ( $ G $ ) in **m³·kg⁻¹·s⁻²** (defaults to `G_SI`).
+ * - Output: ($ \mu $) in **m³/s²**.
+ *
+ * ::: tip
+ *
+ * - Use `m = 0` (default) for the standard published ($ \mu $) of a central body.
+ * - When modelling a binary with comparable masses, pass both `M` and `m`.
+ *
+ * :::
+ *
+ * @param {number} M Mass of the primary body (kg).
+ * @param {number} m Mass of the secondary body (kg). Defaults to `0`.
+ * @param {number} [G] Gravitational constant. Defaults to `G_SI`.
+ * @returns {number} Gravitational parameter ($ \mu $) in m³/s².
+ * @throws {@link https://developer.mozilla.org/en-US/docs/Web/JavaScript/Reference/Global_Objects/Error | Error} If `M` is non-finite or negative.
+ * @throws {@link https://developer.mozilla.org/en-US/docs/Web/JavaScript/Reference/Global_Objects/Error | Error} If `m` is non-finite or negative.
+ * @throws {@link https://developer.mozilla.org/en-US/docs/Web/JavaScript/Reference/Global_Objects/Error | Error} If `G` is non-finite or negative.
+ *
+ * @example
+ * ```ts
+ * // Standard gravitational parameter for Earth (approximate: m ~ 0)
+ * const muEarth = gravitationalParameter(5.97219e24); // ≈ 3.986004e14 m^3/s^2
+ *
+ * // Earth + 1000 kg satellite: practically identical to GM
+ * const muExact = gravitationalParameter(5.97219e24, 1000);
+ *
+ * // Earth–Moon system (two-body μ)
+ * const muEarthMoon = gravitationalParameter(5.97219e24, 7.342e22);
+ * ```
+ * @group Gravity
+ */
+export const gravitationalParameter = (M, m = 0, G = G_SI) => {
+    if (!Number.isFinite(M) || M < 0) {
+        throw new Error('M must be a finite, positive number (kg).');
+    }
+    if (!Number.isFinite(m) || m < 0) {
+        throw new Error('m must be a finite, positive number (kg).');
+    }
+    if (!Number.isFinite(G) || G <= 0) {
+        throw new Error('G must be a finite, positive number (m^3·kg^-1·s^-2).');
+    }
+    return G * (M + m);
+};
+//# sourceMappingURL=gravitational-parameter.js.map

package/dist/categories/gravity/gravitational-parameter.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"gravitational-parameter.js","sourceRoot":"","sources":["../../../src/categories/gravity/gravitational-parameter.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,IAAI,EAAE,MAAM,+BAA+B,CAAC;AAErD;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GA8CG;AACH,MAAM,CAAC,MAAM,sBAAsB,GAAG,CACpC,CAAS,EACT,IAAY,CAAC,EACb,IAAY,IAAI,EACR,EAAE;IACV,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC;QACjC,MAAM,IAAI,KAAK,CAAC,2CAA2C,CAAC,CAAC;IAC/D,CAAC;IAED,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC;QACjC,MAAM,IAAI,KAAK,CAAC,2CAA2C,CAAC,CAAC;IAC/D,CAAC;IAED,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,IAAI,CAAC,IAAI,CAAC,EAAE,CAAC;QAClC,MAAM,IAAI,KAAK,CAAC,uDAAuD,CAAC,CAAC;IAC3E,CAAC;IAED,OAAO,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;AACrB,CAAC,CAAC"}

package/dist/{solve-kepler-bisection.d.ts → categories/kepler/solve-kepler-bisection.d.ts}

@@ -69,7 +69,7 @@ import { Radians } from '@interstellar-tools/types';
  *
  * @see [Kepler's Equation (Wikipedia)](https://en.wikipedia.org/wiki/Kepler%27s_equation)
  * @see [Bisection Method (Wikipedia)](https://en.wikipedia.org/wiki/Bisection_method)
- * @category Solve for Kepler
+ * @group Kepler Solvers
  */
 export declare const solveKeplerBisection: (M: Radians, e: number, maxIter: number, tolerance: number) => Radians;
 //# sourceMappingURL=solve-kepler-bisection.d.ts.map

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

@@ -0,0 +1 @@
+{"version":3,"file":"solve-kepler-bisection.d.ts","sourceRoot":"","sources":["../../../src/categories/kepler/solve-kepler-bisection.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,OAAO,EAAE,MAAM,2BAA2B,CAAC;AAEpD;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GAuEG;AACH,eAAO,MAAM,oBAAoB,GAC/B,GAAG,OAAO,EACV,GAAG,MAAM,EACT,SAAS,MAAM,EACf,WAAW,MAAM,KAChB,OAwBF,CAAC"}

package/dist/{solve-kepler-bisection.js → categories/kepler/solve-kepler-bisection.js}

@@ -68,7 +68,7 @@
  *
  * @see [Kepler's Equation (Wikipedia)](https://en.wikipedia.org/wiki/Kepler%27s_equation)
  * @see [Bisection Method (Wikipedia)](https://en.wikipedia.org/wiki/Bisection_method)
- * @category Solve for Kepler
+ * @group Kepler Solvers
  */
 export const solveKeplerBisection = (M, e, maxIter, tolerance) => {
     let E_low = 0;

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

@@ -0,0 +1 @@
+{"version":3,"file":"solve-kepler-bisection.js","sourceRoot":"","sources":["../../../src/categories/kepler/solve-kepler-bisection.ts"],"names":[],"mappings":"AAEA;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GAuEG;AACH,MAAM,CAAC,MAAM,oBAAoB,GAAG,CAClC,CAAU,EACV,CAAS,EACT,OAAe,EACf,SAAiB,EACR,EAAE;IACX,IAAI,KAAK,GAAG,CAAC,CAAC;IACd,IAAI,MAAM,GAAG,IAAI,CAAC,EAAE,CAAC;IACrB,IAAI,CAAC,GAAG,CAAC,KAAK,GAAG,MAAM,CAAC,GAAG,CAAC,CAAC;IAC7B,IAAI,IAAI,GAAG,CAAC,CAAC;IAEb,OAAO,IAAI,GAAG,OAAO,EAAE,CAAC;QACtB,MAAM,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;QAElC,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,SAAS,EAAE,CAAC;YAC5B,OAAO,CAAC,CAAC,CAAC,YAAY;QACxB,CAAC;QAED,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC;YACV,MAAM,GAAG,CAAC,CAAC;QACb,CAAC;aAAM,CAAC;YACN,KAAK,GAAG,CAAC,CAAC;QACZ,CAAC;QAED,CAAC,GAAG,CAAC,KAAK,GAAG,MAAM,CAAC,GAAG,CAAC,CAAC;QACzB,IAAI,EAAE,CAAC;IACT,CAAC;IAED,OAAO,CAAC,CAAC,CAAC,qBAAqB;AACjC,CAAC,CAAC"}

package/dist/{solve-kepler-high-eccentricity.d.ts → categories/kepler/solve-kepler-high-eccentricity.d.ts}

@@ -3,8 +3,6 @@ import { Radians } from '@interstellar-tools/types';
  * Solves **Kepler's Equation** for the **Eccentric Anomaly** ($E$) in highly eccentric orbits
  * ($e > 0.9$) using an iterative numerical approach.
  *
- * ---
- *
  * **Mathematical Explanation:**
  *
  * Kepler's equation for eccentric anomaly ($E$) is:
@@ -55,13 +53,11 @@ import { Radians } from '@interstellar-tools/types';
  * 3. **Angle Wrapping (Elliptical Only):**
  *    - The result is wrapped using `wrapAngle()` for consistency.
  *
- * ---
  *
  * **Performance Considerations:**
  * - **Typically converges in fewer than 10 iterations for $e > 0.9$.**
  * - **Time complexity:** $O(1)$ for standard convergence.
  *
- * ---
  *
  * @param {Radians} M - Mean anomaly ($M$) in **radians**.
  * @param {number} e - Orbital eccentricity ($e > 0.9$ for high-eccentricity orbits).
@@ -69,7 +65,6 @@ import { Radians } from '@interstellar-tools/types';
  * @param {number} [tolerance=1e-9] - Convergence criterion for stopping the iteration.
  * @returns {Radians} The **eccentric anomaly** ($E$) in **radians** (wrapped to $[-\pi, \pi]$ for elliptical orbits).
  *
- * ---
  *
  * @example
  * ```ts
@@ -89,11 +84,10 @@ import { Radians } from '@interstellar-tools/types';
  * console.log(solveKeplerHighEccentricity(M_hyperbolic, e_hyperbolic));
  * ```
  *
- * ---
- *
  * @see [Kepler's Equation (Wikipedia)](https://en.wikipedia.org/wiki/Kepler%27s_equation)
  * @see [Hyperbolic Anomaly (Wikipedia)](https://en.wikipedia.org/wiki/Hyperbolic_trajectory#Hyperbolic_anomaly)
- * @category Solve for Kepler
+ * @group Kepler Solvers
+ *
  */
 export declare const solveKeplerHighEccentricity: (M: Radians, e: number, maxIter?: number, tolerance?: number) => Radians;
 //# sourceMappingURL=solve-kepler-high-eccentricity.d.ts.map

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

@@ -0,0 +1 @@
+{"version":3,"file":"solve-kepler-high-eccentricity.d.ts","sourceRoot":"","sources":["../../../src/categories/kepler/solve-kepler-high-eccentricity.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,OAAO,EAAE,MAAM,2BAA2B,CAAC;AAKpD;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GAyFG;AAEH,eAAO,MAAM,2BAA2B,GACtC,GAAG,OAAO,EACV,GAAG,MAAM,EACT,gBAAsE,EACtE,kBAAgB,KACf,OAgEF,CAAC"}

package/dist/{solve-kepler-high-eccentricity.js → categories/kepler/solve-kepler-high-eccentricity.js}

@@ -1,10 +1,9 @@
-import { wrapAngle } from './wrap-angle';
+import { wrapAngle } from '../angle/wrap-angle';
+/* node:coverage disable */
 /**
  * Solves **Kepler's Equation** for the **Eccentric Anomaly** ($E$) in highly eccentric orbits
  * ($e > 0.9$) using an iterative numerical approach.
  *
- * ---
- *
  * **Mathematical Explanation:**
  *
  * Kepler's equation for eccentric anomaly ($E$) is:
@@ -55,13 +54,11 @@ import { wrapAngle } from './wrap-angle';
  * 3. **Angle Wrapping (Elliptical Only):**
  *    - The result is wrapped using `wrapAngle()` for consistency.
  *
- * ---
  *
  * **Performance Considerations:**
  * - **Typically converges in fewer than 10 iterations for $e > 0.9$.**
  * - **Time complexity:** $O(1)$ for standard convergence.
  *
- * ---
  *
  * @param {Radians} M - Mean anomaly ($M$) in **radians**.
  * @param {number} e - Orbital eccentricity ($e > 0.9$ for high-eccentricity orbits).
@@ -69,7 +66,6 @@ import { wrapAngle } from './wrap-angle';
  * @param {number} [tolerance=1e-9] - Convergence criterion for stopping the iteration.
  * @returns {Radians} The **eccentric anomaly** ($E$) in **radians** (wrapped to $[-\pi, \pi]$ for elliptical orbits).
  *
- * ---
  *
  * @example
  * ```ts
@@ -89,12 +85,12 @@ import { wrapAngle } from './wrap-angle';
  * console.log(solveKeplerHighEccentricity(M_hyperbolic, e_hyperbolic));
  * ```
  *
- * ---
- *
  * @see [Kepler's Equation (Wikipedia)](https://en.wikipedia.org/wiki/Kepler%27s_equation)
  * @see [Hyperbolic Anomaly (Wikipedia)](https://en.wikipedia.org/wiki/Hyperbolic_trajectory#Hyperbolic_anomaly)
- * @category Solve for Kepler
+ * @group Kepler Solvers
+ *
  */
+/* node:coverage enable */
 export const solveKeplerHighEccentricity = (M, e, maxIter = Math.max(300, Math.floor(5 + 3 * Math.log(1 + Math.abs(M)))), tolerance = 1e-9) => {
     if (e < 0) {
         throw new RangeError(`Invalid eccentricity: ${e}. Must be e ≥ 0.`);

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

@@ -0,0 +1 @@
+{"version":3,"file":"solve-kepler-high-eccentricity.js","sourceRoot":"","sources":["../../../src/categories/kepler/solve-kepler-high-eccentricity.ts"],"names":[],"mappings":"AAEA,OAAO,EAAE,SAAS,EAAE,MAAM,qBAAqB,CAAC;AAEhD,2BAA2B;AAC3B;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GAyFG;AACH,0BAA0B;AAC1B,MAAM,CAAC,MAAM,2BAA2B,GAAG,CACzC,CAAU,EACV,CAAS,EACT,OAAO,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,EAAE,IAAI,CAAC,KAAK,CAAC,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EACtE,SAAS,GAAG,IAAI,EACP,EAAE;IACX,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC;QACV,MAAM,IAAI,UAAU,CAAC,yBAAyB,CAAC,kBAAkB,CAAC,CAAC;IACrE,CAAC;IAED,IAAI,CAAS,CAAC;IAEd,gEAAgE;IAChE,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC;QACV,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;IACpD,CAAC;SAAM,CAAC;QACN,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;IACpD,CAAC;IAED,IAAI,KAAK,GAAG,CAAC,CAAC;IACd,IAAI,IAAI,GAAG,CAAC,CAAC;IACb,IAAI,UAAU,GAAG,CAAC,CAAC;IAEnB,OAAO,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,GAAG,SAAS,IAAI,IAAI,GAAG,OAAO,EAAE,CAAC;QACrD,IAAI,EAAU,CAAC;QACf,IAAI,GAAW,CAAC;QAEhB,qCAAqC;QACrC,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC;YACV,EAAE,GAAG,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;YAC7B,GAAG,GAAG,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;QAC5B,CAAC;QACD,+BAA+B;aAC1B,CAAC;YACJ,MAAM,KAAK,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;YAC3B,MAAM,KAAK,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;YAE3B,EAAE,GAAG,CAAC,GAAG,KAAK,GAAG,CAAC,GAAG,CAAC,CAAC;YACvB,GAAG,GAAG,CAAC,GAAG,KAAK,GAAG,CAAC,CAAC;QACtB,CAAC;QAED,oCAAoC;QACpC,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,IAAI,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC;QAErD,KAAK,GAAG,EAAE,GAAG,GAAG,CAAC;QAEjB,yCAAyC;QACzC,KAAK,GAAG,IAAI,CAAC,IAAI,CAAC,KAAK,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,EAAE,GAAG,CAAC,CAAC;QAE1D,CAAC,IAAI,KAAK,CAAC;QACX,IAAI,EAAE,CAAC;QAEP,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,IAAI,QAAQ,CAAC,CAAC,CAAC,EAAE,CAAC;YAC7B,UAAU,GAAG,CAAC,CAAC;QACjB,CAAC;QAED,iDAAiD;QACjD,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,GAAG,EAAE,CAAC;YACtB,CAAC,GAAG,GAAG,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;YACvB,MAAM;QACR,CAAC;IACH,CAAC;IAED,6DAA6D;IAC7D,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC;QACV,UAAU,GAAG,SAAS,CAAC,UAAU,CAAC,CAAC;IACrC,CAAC;IAED,OAAO,UAAU,CAAC;AACpB,CAAC,CAAC"}

package/dist/{solve-kepler-newton-raphson.d.ts → categories/kepler/solve-kepler-newton-raphson.d.ts}

@@ -81,7 +81,7 @@ import { Radians } from '@interstellar-tools/types';
  * @see [Kepler's Equation (Wikipedia)](https://en.wikipedia.org/wiki/Kepler%27s_equation)
  * @see [Newton-Raphson Method (Wikipedia)](https://en.wikipedia.org/wiki/Newton%27s_method)
  * @see [Eccentric Anomaly (Wikipedia)](https://en.wikipedia.org/wiki/Mean_anomaly#Eccentric_anomaly)
- * @category Solve for Kepler
+ * @group Kepler Solvers
  */
 export declare const solveKeplerNewtonRaphson: (M: Radians, e: number, maxIter: number, tolerance: number) => Radians;
 //# sourceMappingURL=solve-kepler-newton-raphson.d.ts.map

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

@@ -0,0 +1 @@
+{"version":3,"file":"solve-kepler-newton-raphson.d.ts","sourceRoot":"","sources":["../../../src/categories/kepler/solve-kepler-newton-raphson.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,OAAO,EAAE,MAAM,2BAA2B,CAAC;AAEpD;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GAmFG;AACH,eAAO,MAAM,wBAAwB,GACnC,GAAG,OAAO,EACV,GAAG,MAAM,EACT,SAAS,MAAM,EACf,WAAW,MAAM,KAChB,OAwCF,CAAC"}

package/dist/{solve-kepler-newton-raphson.js → categories/kepler/solve-kepler-newton-raphson.js}

@@ -80,7 +80,7 @@
  * @see [Kepler's Equation (Wikipedia)](https://en.wikipedia.org/wiki/Kepler%27s_equation)
  * @see [Newton-Raphson Method (Wikipedia)](https://en.wikipedia.org/wiki/Newton%27s_method)
  * @see [Eccentric Anomaly (Wikipedia)](https://en.wikipedia.org/wiki/Mean_anomaly#Eccentric_anomaly)
- * @category Solve for Kepler
+ * @group Kepler Solvers
  */
 export const solveKeplerNewtonRaphson = (M, e, maxIter, tolerance) => {
     if (e < 0 || e >= 1) {

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

@@ -0,0 +1 @@
+{"version":3,"file":"solve-kepler-newton-raphson.js","sourceRoot":"","sources":["../../../src/categories/kepler/solve-kepler-newton-raphson.ts"],"names":[],"mappings":"AAEA;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GAmFG;AACH,MAAM,CAAC,MAAM,wBAAwB,GAAG,CACtC,CAAU,EACV,CAAS,EACT,OAAe,EACf,SAAiB,EACR,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,IAAI,CAAS,CAAC;IAEd,sBAAsB;IACtB,IAAI,CAAC,GAAG,GAAG,EAAE,CAAC;QACZ,CAAC,GAAG,CAAC,CAAC;IACR,CAAC;SAAM,IAAI,CAAC,GAAG,IAAI,EAAE,CAAC;QACpB,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;IAClD,CAAC;SAAM,CAAC;QACN,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,wBAAwB;IAC3C,CAAC;IAED,IAAI,IAAI,GAAG,CAAC,CAAC;IAEb,OAAO,IAAI,GAAG,OAAO,EAAE,CAAC;QACtB,MAAM,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;QAClC,MAAM,EAAE,GAAG,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;QAC/B,MAAM,GAAG,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;QAC5B,MAAM,GAAG,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;QAC5B,MAAM,MAAM,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,gBAAgB;QAExC,oDAAoD;QACpD,IAAI,UAAU,GAAG,MAAM,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,GAAG,MAAM,GAAG,GAAG,CAAC,GAAG,EAAE,CAAC,CAAC;QAE1D,UAAU,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,UAAU,GAAG,UAAU,GAAG,GAAG,CAAC,GAAG,EAAE,CAAC;QAEjE,CAAC,IAAI,UAAU,CAAC;QAEhB,IAAI,IAAI,CAAC,GAAG,CAAC,UAAU,CAAC,GAAG,SAAS,EAAE,CAAC;YACrC,OAAO,CAAC,CAAC,CAAC,YAAY;QACxB,CAAC;QAED,IAAI,EAAE,CAAC;IACT,CAAC;IAED,OAAO,GAAG,CAAC,CAAC,mBAAmB;AACjC,CAAC,CAAC"}

package/dist/{solve-kepler.d.ts → categories/kepler/solve-kepler.d.ts}

@@ -76,7 +76,7 @@ import { Radians } from '@interstellar-tools/types';
  * @see [Newton-Raphson Method (Wikipedia)](https://en.wikipedia.org/wiki/Newton%27s_method)
  * @see [Eccentric Anomaly (Wikipedia)](https://en.wikipedia.org/wiki/Mean_anomaly#Eccentric_anomaly)
  * @see [RangeError](https://developer.mozilla.org/en-US/docs/Web/JavaScript/Reference/Global_Objects/RangeError)
- * @category Solve for Kepler
+ * @group Kepler Solvers
  */
 export declare const solveKepler: (M: Radians, e: number, maxIter?: number, tolerance?: number) => Radians;
 //# sourceMappingURL=solve-kepler.d.ts.map

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

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

package/dist/{solve-kepler.js → categories/kepler/solve-kepler.js}

@@ -1,7 +1,7 @@
+import { wrapAngle } from '../angle/wrap-angle';
 import { solveKeplerBisection } from './solve-kepler-bisection';
 import { solveKeplerHighEccentricity } from './solve-kepler-high-eccentricity';
 import { solveKeplerNewtonRaphson } from './solve-kepler-newton-raphson';
-import { wrapAngle } from './wrap-angle';
 /**
  * Solves **Kepler's Equation** for the **Eccentric Anomaly** ($E$) using an adaptive approach:
  * - **Newton-Raphson method** for fast convergence.
@@ -79,7 +79,7 @@ import { wrapAngle } from './wrap-angle';
  * @see [Newton-Raphson Method (Wikipedia)](https://en.wikipedia.org/wiki/Newton%27s_method)
  * @see [Eccentric Anomaly (Wikipedia)](https://en.wikipedia.org/wiki/Mean_anomaly#Eccentric_anomaly)
  * @see [RangeError](https://developer.mozilla.org/en-US/docs/Web/JavaScript/Reference/Global_Objects/RangeError)
- * @category Solve for Kepler
+ * @group Kepler Solvers
  */
 export const solveKepler = (M, e, maxIter = 50, tolerance = 1e-9) => {
     if (e < 0 || e >= 1) {

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

@@ -0,0 +1 @@
+{"version":3,"file":"solve-kepler.js","sourceRoot":"","sources":["../../../src/categories/kepler/solve-kepler.ts"],"names":[],"mappings":"AAEA,OAAO,EAAE,SAAS,EAAE,MAAM,qBAAqB,CAAC;AAEhD,OAAO,EAAE,oBAAoB,EAAE,MAAM,0BAA0B,CAAC;AAChE,OAAO,EAAE,2BAA2B,EAAE,MAAM,kCAAkC,CAAC;AAC/E,OAAO,EAAE,wBAAwB,EAAE,MAAM,+BAA+B,CAAC;AAEzE;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GA8EG;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,wDAAwD;IACxD,IAAI,CAAC,GAAG,GAAG,EAAE,CAAC;QACZ,OAAO,2BAA2B,CAAC,CAAC,EAAE,CAAC,EAAE,OAAO,EAAE,SAAS,CAAC,CAAC;IAC/D,CAAC;IAED,+DAA+D;IAC/D,IAAI,CAAC,GAAG,wBAAwB,CAAC,CAAC,EAAE,CAAC,EAAE,OAAO,EAAE,SAAS,CAAC,CAAC;IAE3D,IAAI,KAAK,CAAC,CAAC,CAAC,EAAE,CAAC;QACb,CAAC,GAAG,oBAAoB,CAAC,CAAC,EAAE,CAAC,EAAE,OAAO,EAAE,SAAS,CAAC,CAAC;IACrD,CAAC;IAED,OAAO,SAAS,CAAC,CAAC,CAAC,CAAC;AACtB,CAAC,CAAC"}

package/dist/categories/orbits/circular-speed.d.ts

@@ -0,0 +1,42 @@
+/**
+ * **Circular speed** ($ v_c $) for a body in a circular orbit at radius ($ r $).
+ *
+ * **Definition**
+ *
+ * $$
+ * v_c=\sqrt{\frac{\mu}{r}}
+ * $$
+ *
+ * where ($ \mu $) is the standard gravitational parameter and ($ r $) is the distance
+ * from the central body's center to the orbiting body.
+ *
+ * **Units**
+ * - Inputs: ($ r $) in **meters (m)**, ($ \mu $) in **m³/s²**.
+ * - Output: ($ v_c $) in **m/s**.
+ *
+ * ::: info
+ *
+ * - This is the special case of the vis-viva equation for circular orbits (($ a=r $)).
+ * - If you need a non-circular orbit, use `visVivaSpeed(r, a, mu)`.
+ *
+ * :::
+ *
+ * @param {number} r  Orbital radius ($ r $) (m), must be finite and **> 0**.
+ * @param {number} mu Gravitational parameter ($ \mu $) (m³/s²), must be finite and **≥ 0**.
+ * @returns {number} Circular speed ($ v_c $) (m/s).
+ * @throws {@link https://developer.mozilla.org/en-US/docs/Web/JavaScript/Reference/Global_Objects/Error | Error} If inputs are non-finite or out of domain.
+ *
+ * @see {@link visVivaSpeed}
+ *
+ * @example
+ * ```ts
+ * // LEO ~400 km altitude around Earth
+ * const muEarth = 3.986004418e14; // m^3/s^2
+ * const r = 6378e3 + 400e3;       // m
+ * const v = circularSpeed(r, muEarth); // ≈ 7670 m/s
+ * ```
+ *
+ * @group Orbits
+ */
+export declare const circularSpeed: (r: number, mu: number) => number;
+//# sourceMappingURL=circular-speed.d.ts.map

package/dist/categories/orbits/circular-speed.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"circular-speed.d.ts","sourceRoot":"","sources":["../../../src/categories/orbits/circular-speed.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GAuCG;AACH,eAAO,MAAM,aAAa,GAAI,GAAG,MAAM,EAAE,IAAI,MAAM,KAAG,MAUrD,CAAC"}

package/dist/categories/orbits/circular-speed.js

@@ -0,0 +1,50 @@
+/**
+ * **Circular speed** ($ v_c $) for a body in a circular orbit at radius ($ r $).
+ *
+ * **Definition**
+ *
+ * $$
+ * v_c=\sqrt{\frac{\mu}{r}}
+ * $$
+ *
+ * where ($ \mu $) is the standard gravitational parameter and ($ r $) is the distance
+ * from the central body's center to the orbiting body.
+ *
+ * **Units**
+ * - Inputs: ($ r $) in **meters (m)**, ($ \mu $) in **m³/s²**.
+ * - Output: ($ v_c $) in **m/s**.
+ *
+ * ::: info
+ *
+ * - This is the special case of the vis-viva equation for circular orbits (($ a=r $)).
+ * - If you need a non-circular orbit, use `visVivaSpeed(r, a, mu)`.
+ *
+ * :::
+ *
+ * @param {number} r  Orbital radius ($ r $) (m), must be finite and **> 0**.
+ * @param {number} mu Gravitational parameter ($ \mu $) (m³/s²), must be finite and **≥ 0**.
+ * @returns {number} Circular speed ($ v_c $) (m/s).
+ * @throws {@link https://developer.mozilla.org/en-US/docs/Web/JavaScript/Reference/Global_Objects/Error | Error} If inputs are non-finite or out of domain.
+ *
+ * @see {@link visVivaSpeed}
+ *
+ * @example
+ * ```ts
+ * // LEO ~400 km altitude around Earth
+ * const muEarth = 3.986004418e14; // m^3/s^2
+ * const r = 6378e3 + 400e3;       // m
+ * const v = circularSpeed(r, muEarth); // ≈ 7670 m/s
+ * ```
+ *
+ * @group Orbits
+ */
+export const circularSpeed = (r, mu) => {
+    if (!Number.isFinite(r) || r <= 0) {
+        throw new Error('r must be a finite, positive number (m).');
+    }
+    if (!Number.isFinite(mu) || mu < 0) {
+        throw new Error('mu must be a finite, non-negative number (m^3/s^2).');
+    }
+    return Math.sqrt(mu / r);
+};
+//# sourceMappingURL=circular-speed.js.map

package/dist/categories/orbits/circular-speed.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"circular-speed.js","sourceRoot":"","sources":["../../../src/categories/orbits/circular-speed.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GAuCG;AACH,MAAM,CAAC,MAAM,aAAa,GAAG,CAAC,CAAS,EAAE,EAAU,EAAU,EAAE;IAC7D,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,IAAI,CAAC,IAAI,CAAC,EAAE,CAAC;QAClC,MAAM,IAAI,KAAK,CAAC,0CAA0C,CAAC,CAAC;IAC9D,CAAC;IAED,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,EAAE,CAAC,IAAI,EAAE,GAAG,CAAC,EAAE,CAAC;QACnC,MAAM,IAAI,KAAK,CAAC,qDAAqD,CAAC,CAAC;IACzE,CAAC;IAED,OAAO,IAAI,CAAC,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC;AAC3B,CAAC,CAAC"}