@interstellar-tools/equations 0.49.0 → 0.51.0

package/dist/categories/orbits/__tests__/peri-apoapsis-radii.spec.js

@@ -1,74 +0,0 @@
-import assert from 'node:assert/strict';
-import test, { describe } from 'node:test';
-import { periApoapsisRadii } from '../peri-apoapsis-radii';
-import { relClose } from 'packages/equations/src/__tests__/helpers';
-describe('periApoapsisRadii', () => {
-    test('ellipse: matches rp = a(1-e), ra = a(1+e)', () => {
-        const a = 7000e3; // m
-        const e = 0.01;
-        const { rp, ra } = periApoapsisRadii(a, e);
-        relClose(rp, a * (1 - e), 1e-15);
-        relClose(ra, a * (1 + e), 1e-15);
-        assert.ok(rp > 0 && ra > 0);
-        assert.ok(ra > rp);
-    });
-    test('circular: e = 0 → rp = ra = a', () => {
-        const a = 6778e3;
-        const e = 0;
-        const { rp, ra } = periApoapsisRadii(a, e);
-        relClose(rp, a, 1e-15);
-        relClose(ra, a, 1e-15);
-    });
-    test('ellipse with e close to 1 still yields positive rp and ra', () => {
-        const a = 10000000;
-        const e = 0.999999;
-        const { rp, ra } = periApoapsisRadii(a, e);
-        assert.ok(rp > 0 && ra > 0);
-        relClose(rp, a * (1 - e), 1e-12);
-        relClose(ra, a * (1 + e), 1e-12);
-    });
-    test('hyperbola: a < 0 and e > 1 → rp = |a|(e-1), ra = null', () => {
-        const a = -50000000; // m
-        const e = 1.2;
-        const { rp, ra } = periApoapsisRadii(a, e);
-        relClose(rp, Math.abs(a) * (e - 1), 1e-15);
-        assert.equal(ra, null);
-        assert.ok(rp > 0);
-    });
-    test('invalid combos: parabolic e=1 with a>0; and e>=1 for ellipse', () => {
-        const a = 7000e3;
-        assert.throws(() => periApoapsisRadii(a, 1), /Parabolic \(e=1\) is not defined with finite a/);
-        assert.throws(() => periApoapsisRadii(a, 1.1), /For a>0 \(elliptic\/circular\), eccentricity must satisfy 0 ≤ e < 1/);
-    });
-    test('invalid hyperbola: a < 0 with e ≤ 1 throws', () => {
-        const a = -20000000;
-        assert.throws(() => periApoapsisRadii(a, 1), /For a<0 \(hyperbolic\), eccentricity must satisfy e > 1/);
-        assert.throws(() => periApoapsisRadii(a, 0.5), /For a<0 \(hyperbolic\), eccentricity must satisfy e > 1/);
-    });
-    test('monotonicity (ellipse): increasing e decreases rp and increases ra', () => {
-        const a = 9000e3;
-        const e1 = 0.1;
-        const e2 = 0.3;
-        const r1 = periApoapsisRadii(a, e1);
-        const r2 = periApoapsisRadii(a, e2);
-        assert.ok(r2.rp < r1.rp); // rp decreases with e
-        assert.ok(r2.ra > r1.ra); // ra increases with e
-    });
-    test('input validation: a=0, non-finite a/e; negative e', () => {
-        const a = 7000e3;
-        const e = 0.1;
-        // a invalid
-        assert.throws(() => periApoapsisRadii(0, e), /non-zero/);
-        assert.throws(() => periApoapsisRadii(Number.NaN, e), /finite/);
-        // e invalid
-        assert.throws(() => periApoapsisRadii(a, -1e-3), /≥ 0/);
-        assert.throws(() => periApoapsisRadii(a, Number.NaN), /finite/);
-    });
-    test('numerical stability: very large magnitudes still finite', () => {
-        const { rp: rp1, ra: ra1 } = periApoapsisRadii(1e12, 0.2);
-        assert.ok(Number.isFinite(rp1) && Number.isFinite(ra1));
-        const { rp: rp2, ra: ra2 } = periApoapsisRadii(-1e12, 1.5);
-        assert.ok(Number.isFinite(rp2) && ra2 === null);
-    });
-});
-//# sourceMappingURL=peri-apoapsis-radii.spec.js.map

package/dist/categories/orbits/__tests__/peri-apoapsis-radii.spec.js.map

@@ -1 +0,0 @@
-{"version":3,"file":"peri-apoapsis-radii.spec.js","sourceRoot":"","sources":["../../../../src/categories/orbits/__tests__/peri-apoapsis-radii.spec.ts"],"names":[],"mappings":"AAAA,OAAO,MAAM,MAAM,oBAAoB,CAAC;AACxC,OAAO,IAAI,EAAE,EAAE,QAAQ,EAAE,MAAM,WAAW,CAAC;AAC3C,OAAO,EAAE,iBAAiB,EAAE,MAAM,wBAAwB,CAAC;AAC3D,OAAO,EAAE,QAAQ,EAAE,MAAM,0CAA0C,CAAC;AAEpE,QAAQ,CAAC,mBAAmB,EAAE,GAAG,EAAE;IACjC,IAAI,CAAC,2CAA2C,EAAE,GAAG,EAAE;QACrD,MAAM,CAAC,GAAG,MAAM,CAAC,CAAC,IAAI;QACtB,MAAM,CAAC,GAAG,IAAI,CAAC;QACf,MAAM,EAAE,EAAE,EAAE,EAAE,EAAE,GAAG,iBAAiB,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;QAC3C,QAAQ,CAAC,EAAE,EAAE,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC;QACjC,QAAQ,CAAC,EAAE,EAAE,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC;QACjC,MAAM,CAAC,EAAE,CAAC,EAAE,GAAG,CAAC,IAAI,EAAG,GAAG,CAAC,CAAC,CAAC;QAC7B,MAAM,CAAC,EAAE,CAAE,EAAa,GAAG,EAAE,CAAC,CAAC;IACjC,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,+BAA+B,EAAE,GAAG,EAAE;QACzC,MAAM,CAAC,GAAG,MAAM,CAAC;QACjB,MAAM,CAAC,GAAG,CAAC,CAAC;QACZ,MAAM,EAAE,EAAE,EAAE,EAAE,EAAE,GAAG,iBAAiB,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;QAC3C,QAAQ,CAAC,EAAE,EAAE,CAAC,EAAE,KAAK,CAAC,CAAC;QACvB,QAAQ,CAAC,EAAY,EAAE,CAAC,EAAE,KAAK,CAAC,CAAC;IACnC,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,2DAA2D,EAAE,GAAG,EAAE;QACrE,MAAM,CAAC,GAAG,QAAQ,CAAC;QACnB,MAAM,CAAC,GAAG,QAAQ,CAAC;QACnB,MAAM,EAAE,EAAE,EAAE,EAAE,EAAE,GAAG,iBAAiB,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;QAC3C,MAAM,CAAC,EAAE,CAAC,EAAE,GAAG,CAAC,IAAK,EAAa,GAAG,CAAC,CAAC,CAAC;QACxC,QAAQ,CAAC,EAAE,EAAE,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC;QACjC,QAAQ,CAAC,EAAY,EAAE,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC;IAC7C,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,uDAAuD,EAAE,GAAG,EAAE;QACjE,MAAM,CAAC,GAAG,CAAC,QAAQ,CAAC,CAAC,IAAI;QACzB,MAAM,CAAC,GAAG,GAAG,CAAC;QACd,MAAM,EAAE,EAAE,EAAE,EAAE,EAAE,GAAG,iBAAiB,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;QAC3C,QAAQ,CAAC,EAAE,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC;QAC3C,MAAM,CAAC,KAAK,CAAC,EAAE,EAAE,IAAI,CAAC,CAAC;QACvB,MAAM,CAAC,EAAE,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC;IACpB,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,8DAA8D,EAAE,GAAG,EAAE;QACxE,MAAM,CAAC,GAAG,MAAM,CAAC;QACjB,MAAM,CAAC,MAAM,CACX,GAAG,EAAE,CAAC,iBAAiB,CAAC,CAAC,EAAE,CAAC,CAAC,EAC7B,gDAAgD,CACjD,CAAC;QACF,MAAM,CAAC,MAAM,CACX,GAAG,EAAE,CAAC,iBAAiB,CAAC,CAAC,EAAE,GAAG,CAAC,EAC/B,qEAAqE,CACtE,CAAC;IACJ,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,4CAA4C,EAAE,GAAG,EAAE;QACtD,MAAM,CAAC,GAAG,CAAC,QAAQ,CAAC;QACpB,MAAM,CAAC,MAAM,CACX,GAAG,EAAE,CAAC,iBAAiB,CAAC,CAAC,EAAE,CAAC,CAAC,EAC7B,yDAAyD,CAC1D,CAAC;QACF,MAAM,CAAC,MAAM,CACX,GAAG,EAAE,CAAC,iBAAiB,CAAC,CAAC,EAAE,GAAG,CAAC,EAC/B,yDAAyD,CAC1D,CAAC;IACJ,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,oEAAoE,EAAE,GAAG,EAAE;QAC9E,MAAM,CAAC,GAAG,MAAM,CAAC;QACjB,MAAM,EAAE,GAAG,GAAG,CAAC;QACf,MAAM,EAAE,GAAG,GAAG,CAAC;QAEf,MAAM,EAAE,GAAG,iBAAiB,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC;QACpC,MAAM,EAAE,GAAG,iBAAiB,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC;QAEpC,MAAM,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,GAAG,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,sBAAsB;QAChD,MAAM,CAAC,EAAE,CAAE,EAAE,CAAC,EAAa,GAAI,EAAE,CAAC,EAAa,CAAC,CAAC,CAAC,sBAAsB;IAC1E,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,mDAAmD,EAAE,GAAG,EAAE;QAC7D,MAAM,CAAC,GAAG,MAAM,CAAC;QACjB,MAAM,CAAC,GAAG,GAAG,CAAC;QAEd,YAAY;QACZ,MAAM,CAAC,MAAM,CAAC,GAAG,EAAE,CAAC,iBAAiB,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,UAAU,CAAC,CAAC;QACzD,MAAM,CAAC,MAAM,CAAC,GAAG,EAAE,CAAC,iBAAiB,CAAC,MAAM,CAAC,GAAG,EAAE,CAAC,CAAC,EAAE,QAAQ,CAAC,CAAC;QAEhE,YAAY;QACZ,MAAM,CAAC,MAAM,CAAC,GAAG,EAAE,CAAC,iBAAiB,CAAC,CAAC,EAAE,CAAC,IAAI,CAAC,EAAE,KAAK,CAAC,CAAC;QACxD,MAAM,CAAC,MAAM,CAAC,GAAG,EAAE,CAAC,iBAAiB,CAAC,CAAC,EAAE,MAAM,CAAC,GAAG,CAAC,EAAE,QAAQ,CAAC,CAAC;IAClE,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,yDAAyD,EAAE,GAAG,EAAE;QACnE,MAAM,EAAE,EAAE,EAAE,GAAG,EAAE,EAAE,EAAE,GAAG,EAAE,GAAG,iBAAiB,CAAC,IAAI,EAAE,GAAG,CAAC,CAAC;QAC1D,MAAM,CAAC,EAAE,CAAC,MAAM,CAAC,QAAQ,CAAC,GAAG,CAAC,IAAI,MAAM,CAAC,QAAQ,CAAC,GAAa,CAAC,CAAC,CAAC;QAElE,MAAM,EAAE,EAAE,EAAE,GAAG,EAAE,EAAE,EAAE,GAAG,EAAE,GAAG,iBAAiB,CAAC,CAAC,IAAI,EAAE,GAAG,CAAC,CAAC;QAC3D,MAAM,CAAC,EAAE,CAAC,MAAM,CAAC,QAAQ,CAAC,GAAG,CAAC,IAAI,GAAG,KAAK,IAAI,CAAC,CAAC;IAClD,CAAC,CAAC,CAAC;AACL,CAAC,CAAC,CAAC"}

package/dist/categories/orbits/__tests__/specific-angular-momentum-from-elements.spec.d.ts

@@ -1,2 +0,0 @@
-export {};
-//# sourceMappingURL=specific-angular-momentum-from-elements.spec.d.ts.map

package/dist/categories/orbits/__tests__/specific-angular-momentum-from-elements.spec.d.ts.map

@@ -1 +0,0 @@
-{"version":3,"file":"specific-angular-momentum-from-elements.spec.d.ts","sourceRoot":"","sources":["../../../../src/categories/orbits/__tests__/specific-angular-momentum-from-elements.spec.ts"],"names":[],"mappings":""}

package/dist/categories/orbits/__tests__/specific-angular-momentum-from-elements.spec.js

@@ -1,86 +0,0 @@
-import assert from 'node:assert/strict';
-import test, { describe } from 'node:test';
-import { absClose, relClose } from 'packages/equations/src/__tests__/helpers';
-import { specificAngularMomentumFromElements } from '../specific-angular-momentum-from-elements';
-describe('specificAngularMomentumFromElements', () => {
-    test('definition (elliptic): h = sqrt(μ a (1-e²))', () => {
-        const mu = 3.986004418e14; // m³/s² (Earth GM)
-        const a = 7000e3; // m
-        const e = 0.1; // -
-        const h = specificAngularMomentumFromElements(a, e, mu);
-        const expected = Math.sqrt(mu * a * (1 - e * e));
-        relClose(h, expected, 1e-15);
-    });
-    test('circular limit (e=0): h = sqrt(μ a)', () => {
-        const mu = 3.986004418e14;
-        const a = 6778e3;
-        const e = 0;
-        const h = specificAngularMomentumFromElements(a, e, mu);
-        const expected = Math.sqrt(mu * a);
-        relClose(h, expected, 1e-15);
-    });
-    test('hyperbolic case (a<0, e>1): radicand positive → valid h', () => {
-        const mu = 3.986004418e14;
-        const a = -50000000; // m
-        const e = 1.2; // -
-        const h = specificAngularMomentumFromElements(a, e, mu);
-        const expected = Math.sqrt(mu * a * (1 - e * e)); // positive since a<0 and (1-e²)<0
-        relClose(h, expected, 1e-15);
-    });
-    test('tolerance clamp: tiny negative radicand → returns ~0', () => {
-        const mu = 3.986004418e14;
-        const a = 7000e3; // m (positive)
-        // Function tolerance: tol = |1e-14 * mu * |a||
-        const tol = Math.abs(1e-14 * mu * Math.abs(a));
-        // Target radicand = -tol/2: μ*a*(1-e²) = -tol/2  =>  1-e² = -(tol/2)/(μ*a)
-        const oneMinusE2 = -(tol / 2) / (mu * a);
-        const e = Math.sqrt(1 - oneMinusE2); // slightly > 1 (physically inconsistent with a>0), used only to hit tolerance path
-        const h = specificAngularMomentumFromElements(a, e, mu);
-        absClose(h, 0, 1e-12);
-    });
-    test('invalid: clearly negative radicand throws (a>0, e>1 far from 1)', () => {
-        const mu = 3.986004418e14;
-        const a = 7000e3; // m
-        const e = 1.1; // makes 1 - e² < 0 with a>0 → negative radicand
-        assert.throws(() => specificAngularMomentumFromElements(a, e, mu), /Invalid \(a,e\).*μ·a·\(1-e²\) < 0/);
-    });
-    test('μ = 0 → h = 0 (degenerate but mathematically consistent)', () => {
-        const mu = 0;
-        const a = 7000e3;
-        const e = 0.3;
-        const h = specificAngularMomentumFromElements(a, e, mu);
-        absClose(h, 0, 0);
-    });
-    test('scaling: for fixed e, h ∝ sqrt(|a|)', () => {
-        const mu = 3.986004418e14;
-        const e = 0.2;
-        const a1 = 8000e3;
-        const a2 = 2 * a1;
-        const h1 = specificAngularMomentumFromElements(a1, e, mu);
-        const h2 = specificAngularMomentumFromElements(a2, e, mu);
-        const expectedRatio = Math.sqrt(Math.abs(a2) / Math.abs(a1));
-        relClose(h2 / h1, expectedRatio, 1e-12);
-    });
-    test('input validation: a=0, non-finite a/e/μ, negative e, negative μ', () => {
-        const a = 7000e3;
-        const e = 0.1;
-        const mu = 3.986004418e14;
-        // a invalid
-        assert.throws(() => specificAngularMomentumFromElements(0, e, mu), /non-zero/);
-        assert.throws(() => specificAngularMomentumFromElements(Number.NaN, e, mu), /finite/);
-        // e invalid
-        assert.throws(() => specificAngularMomentumFromElements(a, -1e-3, mu), /≥ 0/);
-        assert.throws(() => specificAngularMomentumFromElements(a, Number.NaN, mu), /finite/);
-        // μ invalid
-        assert.throws(() => specificAngularMomentumFromElements(a, e, -1), /≥ 0/);
-        assert.throws(() => specificAngularMomentumFromElements(a, e, Number.NaN), /finite/);
-    });
-    test('numerical sanity: very large magnitudes still finite', () => {
-        const mu = 1e20;
-        const a = 1e9;
-        const e = 0.5;
-        const h = specificAngularMomentumFromElements(a, e, mu);
-        assert.ok(Number.isFinite(h) && h > 0);
-    });
-});
-//# sourceMappingURL=specific-angular-momentum-from-elements.spec.js.map

package/dist/categories/orbits/__tests__/specific-angular-momentum-from-elements.spec.js.map

@@ -1 +0,0 @@
-{"version":3,"file":"specific-angular-momentum-from-elements.spec.js","sourceRoot":"","sources":["../../../../src/categories/orbits/__tests__/specific-angular-momentum-from-elements.spec.ts"],"names":[],"mappings":"AAAA,OAAO,MAAM,MAAM,oBAAoB,CAAC;AACxC,OAAO,IAAI,EAAE,EAAE,QAAQ,EAAE,MAAM,WAAW,CAAC;AAC3C,OAAO,EAAE,QAAQ,EAAE,QAAQ,EAAE,MAAM,0CAA0C,CAAC;AAE9E,OAAO,EAAE,mCAAmC,EAAE,MAAM,4CAA4C,CAAC;AAEjG,QAAQ,CAAC,qCAAqC,EAAE,GAAG,EAAE;IACnD,IAAI,CAAC,6CAA6C,EAAE,GAAG,EAAE;QACvD,MAAM,EAAE,GAAG,cAAc,CAAC,CAAC,mBAAmB;QAC9C,MAAM,CAAC,GAAG,MAAM,CAAC,CAAC,IAAI;QACtB,MAAM,CAAC,GAAG,GAAG,CAAC,CAAC,IAAI;QACnB,MAAM,CAAC,GAAG,mCAAmC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC;QACxD,MAAM,QAAQ,GAAG,IAAI,CAAC,IAAI,CAAC,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;QACjD,QAAQ,CAAC,CAAC,EAAE,QAAQ,EAAE,KAAK,CAAC,CAAC;IAC/B,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,qCAAqC,EAAE,GAAG,EAAE;QAC/C,MAAM,EAAE,GAAG,cAAc,CAAC;QAC1B,MAAM,CAAC,GAAG,MAAM,CAAC;QACjB,MAAM,CAAC,GAAG,CAAC,CAAC;QACZ,MAAM,CAAC,GAAG,mCAAmC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC;QACxD,MAAM,QAAQ,GAAG,IAAI,CAAC,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC;QACnC,QAAQ,CAAC,CAAC,EAAE,QAAQ,EAAE,KAAK,CAAC,CAAC;IAC/B,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,yDAAyD,EAAE,GAAG,EAAE;QACnE,MAAM,EAAE,GAAG,cAAc,CAAC;QAC1B,MAAM,CAAC,GAAG,CAAC,QAAQ,CAAC,CAAC,IAAI;QACzB,MAAM,CAAC,GAAG,GAAG,CAAC,CAAC,IAAI;QACnB,MAAM,CAAC,GAAG,mCAAmC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC;QACxD,MAAM,QAAQ,GAAG,IAAI,CAAC,IAAI,CAAC,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,kCAAkC;QACpF,QAAQ,CAAC,CAAC,EAAE,QAAQ,EAAE,KAAK,CAAC,CAAC;IAC/B,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,sDAAsD,EAAE,GAAG,EAAE;QAChE,MAAM,EAAE,GAAG,cAAc,CAAC;QAC1B,MAAM,CAAC,GAAG,MAAM,CAAC,CAAC,eAAe;QACjC,+CAA+C;QAC/C,MAAM,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,GAAG,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;QAC/C,2EAA2E;QAC3E,MAAM,UAAU,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,GAAG,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC;QACzC,MAAM,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,UAAU,CAAC,CAAC,CAAC,mFAAmF;QAExH,MAAM,CAAC,GAAG,mCAAmC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC;QACxD,QAAQ,CAAC,CAAC,EAAE,CAAC,EAAE,KAAK,CAAC,CAAC;IACxB,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,iEAAiE,EAAE,GAAG,EAAE;QAC3E,MAAM,EAAE,GAAG,cAAc,CAAC;QAC1B,MAAM,CAAC,GAAG,MAAM,CAAC,CAAC,IAAI;QACtB,MAAM,CAAC,GAAG,GAAG,CAAC,CAAC,gDAAgD;QAC/D,MAAM,CAAC,MAAM,CACX,GAAG,EAAE,CAAC,mCAAmC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,EACnD,mCAAmC,CACpC,CAAC;IACJ,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,0DAA0D,EAAE,GAAG,EAAE;QACpE,MAAM,EAAE,GAAG,CAAC,CAAC;QACb,MAAM,CAAC,GAAG,MAAM,CAAC;QACjB,MAAM,CAAC,GAAG,GAAG,CAAC;QACd,MAAM,CAAC,GAAG,mCAAmC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC;QACxD,QAAQ,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;IACpB,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,qCAAqC,EAAE,GAAG,EAAE;QAC/C,MAAM,EAAE,GAAG,cAAc,CAAC;QAC1B,MAAM,CAAC,GAAG,GAAG,CAAC;QACd,MAAM,EAAE,GAAG,MAAM,CAAC;QAClB,MAAM,EAAE,GAAG,CAAC,GAAG,EAAE,CAAC;QAElB,MAAM,EAAE,GAAG,mCAAmC,CAAC,EAAE,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC;QAC1D,MAAM,EAAE,GAAG,mCAAmC,CAAC,EAAE,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC;QAE1D,MAAM,aAAa,GAAG,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC;QAC7D,QAAQ,CAAC,EAAE,GAAG,EAAE,EAAE,aAAa,EAAE,KAAK,CAAC,CAAC;IAC1C,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,iEAAiE,EAAE,GAAG,EAAE;QAC3E,MAAM,CAAC,GAAG,MAAM,CAAC;QACjB,MAAM,CAAC,GAAG,GAAG,CAAC;QACd,MAAM,EAAE,GAAG,cAAc,CAAC;QAE1B,YAAY;QACZ,MAAM,CAAC,MAAM,CACX,GAAG,EAAE,CAAC,mCAAmC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,EACnD,UAAU,CACX,CAAC;QACF,MAAM,CAAC,MAAM,CACX,GAAG,EAAE,CAAC,mCAAmC,CAAC,MAAM,CAAC,GAAG,EAAE,CAAC,EAAE,EAAE,CAAC,EAC5D,QAAQ,CACT,CAAC;QAEF,YAAY;QACZ,MAAM,CAAC,MAAM,CACX,GAAG,EAAE,CAAC,mCAAmC,CAAC,CAAC,EAAE,CAAC,IAAI,EAAE,EAAE,CAAC,EACvD,KAAK,CACN,CAAC;QACF,MAAM,CAAC,MAAM,CACX,GAAG,EAAE,CAAC,mCAAmC,CAAC,CAAC,EAAE,MAAM,CAAC,GAAG,EAAE,EAAE,CAAC,EAC5D,QAAQ,CACT,CAAC;QAEF,YAAY;QACZ,MAAM,CAAC,MAAM,CAAC,GAAG,EAAE,CAAC,mCAAmC,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC;QAC1E,MAAM,CAAC,MAAM,CACX,GAAG,EAAE,CAAC,mCAAmC,CAAC,CAAC,EAAE,CAAC,EAAE,MAAM,CAAC,GAAG,CAAC,EAC3D,QAAQ,CACT,CAAC;IACJ,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,sDAAsD,EAAE,GAAG,EAAE;QAChE,MAAM,EAAE,GAAG,IAAI,CAAC;QAChB,MAAM,CAAC,GAAG,GAAG,CAAC;QACd,MAAM,CAAC,GAAG,GAAG,CAAC;QACd,MAAM,CAAC,GAAG,mCAAmC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC;QACxD,MAAM,CAAC,EAAE,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC;IACzC,CAAC,CAAC,CAAC;AACL,CAAC,CAAC,CAAC"}

package/dist/categories/orbits/__tests__/specific-angular-momentum.spec.d.ts

@@ -1,2 +0,0 @@
-export {};
-//# sourceMappingURL=specific-angular-momentum.spec.d.ts.map

package/dist/categories/orbits/__tests__/specific-angular-momentum.spec.d.ts.map

@@ -1 +0,0 @@
-{"version":3,"file":"specific-angular-momentum.spec.d.ts","sourceRoot":"","sources":["../../../../src/categories/orbits/__tests__/specific-angular-momentum.spec.ts"],"names":[],"mappings":""}

package/dist/categories/orbits/__tests__/specific-angular-momentum.spec.js

@@ -1,85 +0,0 @@
-import assert from 'node:assert/strict';
-import test, { describe } from 'node:test';
-import { absClose, relClose } from 'packages/equations/src/__tests__/helpers';
-import { specificAngularMomentum } from '../specific-angular-momentum';
-import { dot, norm } from '../../helpers/misc';
-describe('specificAngularMomentum', () => {
-    test('orthogonal state: r ⟂ v → h = r * v', () => {
-        const r = [6778e3, 0, 0]; // m
-        const v = [0, 7.67e3, 0]; // m/s
-        const h = specificAngularMomentum(r, v);
-        const expected = norm(r) * norm(v); // r*v_perp (since perp)
-        relClose(h, expected, 1e-14);
-    });
-    test('manual cross vs function: ||r×v|| matches', () => {
-        const r = [1.2e7, -3.4e6, 5.6e6];
-        const v = [-7.8e3, 9.1e3, 2.3e3];
-        const hx = r[1] * v[2] - r[2] * v[1];
-        const hy = r[2] * v[0] - r[0] * v[2];
-        const hz = r[0] * v[1] - r[1] * v[0];
-        const expected = Math.hypot(hx, hy, hz);
-        const h = specificAngularMomentum(r, v);
-        relClose(h, expected, 1e-15);
-    });
-    test('scalar identity: h = |r| |v| sin(theta)', () => {
-        const r = [7e6, 2e6, -1e6];
-        const v = [1e3, -7e3, 4e3];
-        const rr = norm(r);
-        const vv = norm(v);
-        const cosT = dot(r, v) / (rr * vv);
-        const sinT = Math.sqrt(Math.max(0, 1 - cosT * cosT));
-        const expected = rr * vv * sinT;
-        const h = specificAngularMomentum(r, v);
-        relClose(h, expected, 1e-12);
-    });
-    test('parallel vectors (purely radial) → h = 0', () => {
-        const r = [8e6, -2e6, 4e6];
-        const k = 3.5e-3; // arbitrary scale to mimic m/s
-        const v = [k * r[0], k * r[1], k * r[2]];
-        const h = specificAngularMomentum(r, v);
-        absClose(h, 0, 1e-12);
-    });
-    test('zero velocity → h = 0', () => {
-        const r = [7e6, 0, 0];
-        const v = [0, 0, 0];
-        const h = specificAngularMomentum(r, v);
-        absClose(h, 0, 1e-12);
-    });
-    test('rotation invariance of magnitude (rotate 90° about z)', () => {
-        const r = [5e6, 2e6, 0];
-        const v = [1e3, 3e3, 0];
-        const h1 = specificAngularMomentum(r, v);
-        // 90° rotation about z: (x,y) -> (-y, x)
-        const r2 = [-r[1], r[0], 0];
-        const v2 = [-v[1], v[0], 0];
-        const h2 = specificAngularMomentum(r2, v2);
-        relClose(h2, h1, 1e-15);
-    });
-    test('units sanity (circular identity): h ≈ sqrt(mu * r) with v_c = sqrt(mu/r)', () => {
-        const mu = 3.986004418e14; // Earth GM
-        const r = [6778e3, 0, 0];
-        const v = [0, Math.sqrt(mu / r[0]), 0];
-        const h = specificAngularMomentum(r, v);
-        const expected = Math.sqrt(mu * r[0]);
-        relClose(h, expected, 1e-12);
-    });
-    test('input validation: non-finite components throw', () => {
-        const rOK = [7e6, 0, 0];
-        const vOK = [0, 7.5e3, 0];
-        assert.throws(() => specificAngularMomentum([Number.NaN, 0, 0], vOK), /finite/);
-        assert.throws(() => specificAngularMomentum([Infinity, 0, 0], vOK), /finite/);
-        assert.throws(() => specificAngularMomentum(rOK, [0, Number.NaN, 0]), /finite/);
-    });
-    test('input validation: zero-length position vector throws', () => {
-        const rZero = [0, 0, 0];
-        const v = [1e3, 2e3, 3e3];
-        assert.throws(() => specificAngularMomentum(rZero, v), /non-zero/);
-    });
-    test('numerical stability: very large magnitudes still finite', () => {
-        const r = [1e9, 1e9, 1e9];
-        const v = [1e5, -2e5, 3e5];
-        const h = specificAngularMomentum(r, v);
-        assert.ok(Number.isFinite(h) && h >= 0);
-    });
-});
-//# sourceMappingURL=specific-angular-momentum.spec.js.map

package/dist/categories/orbits/__tests__/specific-angular-momentum.spec.js.map

@@ -1 +0,0 @@
-{"version":3,"file":"specific-angular-momentum.spec.js","sourceRoot":"","sources":["../../../../src/categories/orbits/__tests__/specific-angular-momentum.spec.ts"],"names":[],"mappings":"AAAA,OAAO,MAAM,MAAM,oBAAoB,CAAC;AACxC,OAAO,IAAI,EAAE,EAAE,QAAQ,EAAE,MAAM,WAAW,CAAC;AAC3C,OAAO,EAAE,QAAQ,EAAE,QAAQ,EAAE,MAAM,0CAA0C,CAAC;AAC9E,OAAO,EAAE,uBAAuB,EAAE,MAAM,8BAA8B,CAAC;AAEvE,OAAO,EAAE,GAAG,EAAE,IAAI,EAAE,MAAM,oBAAoB,CAAC;AAE/C,QAAQ,CAAC,yBAAyB,EAAE,GAAG,EAAE;IACvC,IAAI,CAAC,qCAAqC,EAAE,GAAG,EAAE;QAC/C,MAAM,CAAC,GAAsB,CAAC,MAAM,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,IAAI;QACjD,MAAM,CAAC,GAAsB,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,CAAC,CAAC,CAAC,MAAM;QACnD,MAAM,CAAC,GAAG,uBAAuB,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;QACxC,MAAM,QAAQ,GAAG,IAAI,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,wBAAwB;QAC5D,QAAQ,CAAC,CAAC,EAAE,QAAQ,EAAE,KAAK,CAAC,CAAC;IAC/B,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,2CAA2C,EAAE,GAAG,EAAE;QACrD,MAAM,CAAC,GAAsB,CAAC,KAAK,EAAE,CAAC,KAAK,EAAE,KAAK,CAAC,CAAC;QACpD,MAAM,CAAC,GAAsB,CAAC,CAAC,KAAK,EAAE,KAAK,EAAE,KAAK,CAAC,CAAC;QAEpD,MAAM,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;QACrC,MAAM,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;QACrC,MAAM,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;QACrC,MAAM,QAAQ,GAAG,IAAI,CAAC,KAAK,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC;QAExC,MAAM,CAAC,GAAG,uBAAuB,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;QACxC,QAAQ,CAAC,CAAC,EAAE,QAAQ,EAAE,KAAK,CAAC,CAAC;IAC/B,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,yCAAyC,EAAE,GAAG,EAAE;QACnD,MAAM,CAAC,GAAsB,CAAC,GAAG,EAAE,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC;QAC9C,MAAM,CAAC,GAAsB,CAAC,GAAG,EAAE,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC;QAC9C,MAAM,EAAE,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC;QACnB,MAAM,EAAE,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC;QACnB,MAAM,IAAI,GAAG,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,EAAE,GAAG,EAAE,CAAC,CAAC;QACnC,MAAM,IAAI,GAAG,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,GAAG,IAAI,GAAG,IAAI,CAAC,CAAC,CAAC;QACrD,MAAM,QAAQ,GAAG,EAAE,GAAG,EAAE,GAAG,IAAI,CAAC;QAEhC,MAAM,CAAC,GAAG,uBAAuB,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;QACxC,QAAQ,CAAC,CAAC,EAAE,QAAQ,EAAE,KAAK,CAAC,CAAC;IAC/B,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,0CAA0C,EAAE,GAAG,EAAE;QACpD,MAAM,CAAC,GAAsB,CAAC,GAAG,EAAE,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC;QAC9C,MAAM,CAAC,GAAG,MAAM,CAAC,CAAC,+BAA+B;QACjD,MAAM,CAAC,GAAsB,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QAC5D,MAAM,CAAC,GAAG,uBAAuB,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;QACxC,QAAQ,CAAC,CAAC,EAAE,CAAC,EAAE,KAAK,CAAC,CAAC;IACxB,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,uBAAuB,EAAE,GAAG,EAAE;QACjC,MAAM,CAAC,GAAsB,CAAC,GAAG,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;QACzC,MAAM,CAAC,GAAsB,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;QACvC,MAAM,CAAC,GAAG,uBAAuB,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;QACxC,QAAQ,CAAC,CAAC,EAAE,CAAC,EAAE,KAAK,CAAC,CAAC;IACxB,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,uDAAuD,EAAE,GAAG,EAAE;QACjE,MAAM,CAAC,GAAsB,CAAC,GAAG,EAAE,GAAG,EAAE,CAAC,CAAC,CAAC;QAC3C,MAAM,CAAC,GAAsB,CAAC,GAAG,EAAE,GAAG,EAAE,CAAC,CAAC,CAAC;QAE3C,MAAM,EAAE,GAAG,uBAAuB,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;QAEzC,yCAAyC;QACzC,MAAM,EAAE,GAAsB,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;QAC/C,MAAM,EAAE,GAAsB,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;QAE/C,MAAM,EAAE,GAAG,uBAAuB,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC;QAC3C,QAAQ,CAAC,EAAE,EAAE,EAAE,EAAE,KAAK,CAAC,CAAC;IAC1B,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,0EAA0E,EAAE,GAAG,EAAE;QACpF,MAAM,EAAE,GAAG,cAAc,CAAC,CAAC,WAAW;QACtC,MAAM,CAAC,GAAsB,CAAC,MAAM,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;QAC5C,MAAM,CAAC,GAAsB,CAAC,CAAC,EAAE,IAAI,CAAC,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;QAC1D,MAAM,CAAC,GAAG,uBAAuB,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;QACxC,MAAM,QAAQ,GAAG,IAAI,CAAC,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QACtC,QAAQ,CAAC,CAAC,EAAE,QAAQ,EAAE,KAAK,CAAC,CAAC;IAC/B,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,+CAA+C,EAAE,GAAG,EAAE;QACzD,MAAM,GAAG,GAAsB,CAAC,GAAG,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;QAC3C,MAAM,GAAG,GAAsB,CAAC,CAAC,EAAE,KAAK,EAAE,CAAC,CAAC,CAAC;QAE7C,MAAM,CAAC,MAAM,CACX,GAAG,EAAE,CACH,uBAAuB,CAAC,CAAC,MAAM,CAAC,GAAG,EAAE,CAAC,EAAE,CAAC,CAAsB,EAAE,GAAG,CAAC,EACvE,QAAQ,CACT,CAAC;QACF,MAAM,CAAC,MAAM,CACX,GAAG,EAAE,CACH,uBAAuB,CACrB,CAAC,QAAQ,EAAE,CAAC,EAAE,CAAC,CAAiC,EAChD,GAAG,CACJ,EACH,QAAQ,CACT,CAAC;QACF,MAAM,CAAC,MAAM,CACX,GAAG,EAAE,CACH,uBAAuB,CAAC,GAAG,EAAE,CAAC,CAAC,EAAE,MAAM,CAAC,GAAG,EAAE,CAAC,CAAsB,CAAC,EACvE,QAAQ,CACT,CAAC;IACJ,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,sDAAsD,EAAE,GAAG,EAAE;QAChE,MAAM,KAAK,GAAsB,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;QAC3C,MAAM,CAAC,GAAsB,CAAC,GAAG,EAAE,GAAG,EAAE,GAAG,CAAC,CAAC;QAC7C,MAAM,CAAC,MAAM,CAAC,GAAG,EAAE,CAAC,uBAAuB,CAAC,KAAK,EAAE,CAAC,CAAC,EAAE,UAAU,CAAC,CAAC;IACrE,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,yDAAyD,EAAE,GAAG,EAAE;QACnE,MAAM,CAAC,GAAsB,CAAC,GAAG,EAAE,GAAG,EAAE,GAAG,CAAC,CAAC;QAC7C,MAAM,CAAC,GAAsB,CAAC,GAAG,EAAE,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC;QAC9C,MAAM,CAAC,GAAG,uBAAuB,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;QACxC,MAAM,CAAC,EAAE,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC;IAC1C,CAAC,CAAC,CAAC;AACL,CAAC,CAAC,CAAC"}

package/dist/categories/orbits/__tests__/specific-mechanical-energy.spec.d.ts

@@ -1,2 +0,0 @@
-export {};
-//# sourceMappingURL=specific-mechanical-energy.spec.d.ts.map

package/dist/categories/orbits/__tests__/specific-mechanical-energy.spec.d.ts.map

@@ -1 +0,0 @@
-{"version":3,"file":"specific-mechanical-energy.spec.d.ts","sourceRoot":"","sources":["../../../../src/categories/orbits/__tests__/specific-mechanical-energy.spec.ts"],"names":[],"mappings":""}

package/dist/categories/orbits/__tests__/specific-mechanical-energy.spec.js

@@ -1,81 +0,0 @@
-import assert from 'node:assert/strict';
-import test, { describe } from 'node:test';
-import { specificMechanicalEnergy } from '../specific-mechanical-energy';
-import { absClose, relClose } from 'packages/equations/src/__tests__/helpers';
-describe('specificMechanicalEnergy', () => {
-    test('circular orbit: ε = -μ/(2r) when v = sqrt(μ/r)', () => {
-        const mu = 3.986004418e14; // m³/s² (Earth GM)
-        const r = 6378e3 + 400e3; // m (LEO-ish)
-        const v = Math.sqrt(mu / r); // m/s
-        const eps = specificMechanicalEnergy(v, r, mu);
-        const expected = -mu / (2 * r);
-        relClose(eps, expected, 1e-12);
-    });
-    test('parabolic case: v = sqrt(2μ/r) → ε ≈ 0', () => {
-        const mu = 3.986004418e14; // m³/s²
-        const r = 7000e3; // m
-        const v = Math.sqrt((2 * mu) / r);
-        const eps = specificMechanicalEnergy(v, r, mu);
-        // Use absolute tolerance because the expected value is 0
-        absClose(eps, 0, 1e-9);
-    });
-    test('hyperbolic case: v slightly above escape → ε > 0', () => {
-        const mu = 3.986004418e14;
-        const r = 7000e3;
-        const v_escape = Math.sqrt((2 * mu) / r);
-        const v = 1.01 * v_escape;
-        const eps = specificMechanicalEnergy(v, r, mu);
-        assert.ok(eps > 0, 'ε should be positive for hyperbolic cases');
-    });
-    test('elliptic case: v below circular → ε < 0', () => {
-        const mu = 3.986004418e14;
-        const r = 7000e3;
-        const v_circ = Math.sqrt(mu / r);
-        const v = 0.95 * v_circ;
-        const eps = specificMechanicalEnergy(v, r, mu);
-        assert.ok(eps < 0, 'ε should be negative for bound (elliptic) cases');
-    });
-    test('semi-major axis identity: ε = -μ/(2a) (using vis-viva for v)', () => {
-        const mu = 3.986004418e14;
-        const a = 10000e3; // m (semi-major axis)
-        const r = 8000e3; // m (some point on the ellipse)
-        const v = Math.sqrt(mu * (2 / r - 1 / a)); // vis-viva
-        const eps = specificMechanicalEnergy(v, r, mu);
-        const expected = -mu / (2 * a);
-        relClose(eps, expected, 1e-12);
-    });
-    test('monotonic in speed: increasing v increases ε (same r, μ)', () => {
-        const mu = 3.986004418e14;
-        const r = 8000e3;
-        const v1 = 7000;
-        const v2 = 7500;
-        const eps1 = specificMechanicalEnergy(v1, r, mu);
-        const eps2 = specificMechanicalEnergy(v2, r, mu);
-        assert.ok(eps2 > eps1);
-    });
-    test('increases with radius for fixed v: larger r → larger ε (less negative)', () => {
-        const mu = 3.986004418e14;
-        const v = 7500;
-        const r1 = 7000e3;
-        const r2 = 8000e3;
-        const eps1 = specificMechanicalEnergy(v, r1, mu);
-        const eps2 = specificMechanicalEnergy(v, r2, mu);
-        assert.ok(eps2 > eps1);
-    });
-    test('throws on invalid inputs', () => {
-        const mu = 3.986004418e14;
-        const r = 7000e3;
-        const v = 7500;
-        // v invalid
-        assert.throws(() => specificMechanicalEnergy(Number.NaN, r, mu), /non-negative/);
-        assert.throws(() => specificMechanicalEnergy(-1, r, mu), /non-negative/);
-        // r invalid
-        assert.throws(() => specificMechanicalEnergy(v, Number.NaN, mu), /positive/);
-        assert.throws(() => specificMechanicalEnergy(v, 0, mu), /positive/);
-        assert.throws(() => specificMechanicalEnergy(v, -1, mu), /positive/);
-        // mu invalid
-        assert.throws(() => specificMechanicalEnergy(v, r, Number.NaN), /non-negative/);
-        assert.throws(() => specificMechanicalEnergy(v, r, -1), /non-negative/);
-    });
-});
-//# sourceMappingURL=specific-mechanical-energy.spec.js.map

package/dist/categories/orbits/__tests__/specific-mechanical-energy.spec.js.map

@@ -1 +0,0 @@
-{"version":3,"file":"specific-mechanical-energy.spec.js","sourceRoot":"","sources":["../../../../src/categories/orbits/__tests__/specific-mechanical-energy.spec.ts"],"names":[],"mappings":"AAAA,OAAO,MAAM,MAAM,oBAAoB,CAAC;AACxC,OAAO,IAAI,EAAE,EAAE,QAAQ,EAAE,MAAM,WAAW,CAAC;AAC3C,OAAO,EAAE,wBAAwB,EAAE,MAAM,+BAA+B,CAAC;AACzE,OAAO,EAAE,QAAQ,EAAE,QAAQ,EAAE,MAAM,0CAA0C,CAAC;AAE9E,QAAQ,CAAC,0BAA0B,EAAE,GAAG,EAAE;IACxC,IAAI,CAAC,gDAAgD,EAAE,GAAG,EAAE;QAC1D,MAAM,EAAE,GAAG,cAAc,CAAC,CAAC,mBAAmB;QAC9C,MAAM,CAAC,GAAG,MAAM,GAAG,KAAK,CAAC,CAAC,cAAc;QACxC,MAAM,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,MAAM;QAEnC,MAAM,GAAG,GAAG,wBAAwB,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC;QAC/C,MAAM,QAAQ,GAAG,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;QAE/B,QAAQ,CAAC,GAAG,EAAE,QAAQ,EAAE,KAAK,CAAC,CAAC;IACjC,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,wCAAwC,EAAE,GAAG,EAAE;QAClD,MAAM,EAAE,GAAG,cAAc,CAAC,CAAC,QAAQ;QACnC,MAAM,CAAC,GAAG,MAAM,CAAC,CAAC,IAAI;QACtB,MAAM,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC;QAElC,MAAM,GAAG,GAAG,wBAAwB,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC;QAE/C,yDAAyD;QACzD,QAAQ,CAAC,GAAG,EAAE,CAAC,EAAE,IAAI,CAAC,CAAC;IACzB,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,kDAAkD,EAAE,GAAG,EAAE;QAC5D,MAAM,EAAE,GAAG,cAAc,CAAC;QAC1B,MAAM,CAAC,GAAG,MAAM,CAAC;QACjB,MAAM,QAAQ,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC;QACzC,MAAM,CAAC,GAAG,IAAI,GAAG,QAAQ,CAAC;QAE1B,MAAM,GAAG,GAAG,wBAAwB,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC;QAC/C,MAAM,CAAC,EAAE,CAAC,GAAG,GAAG,CAAC,EAAE,2CAA2C,CAAC,CAAC;IAClE,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,yCAAyC,EAAE,GAAG,EAAE;QACnD,MAAM,EAAE,GAAG,cAAc,CAAC;QAC1B,MAAM,CAAC,GAAG,MAAM,CAAC;QACjB,MAAM,MAAM,GAAG,IAAI,CAAC,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC;QACjC,MAAM,CAAC,GAAG,IAAI,GAAG,MAAM,CAAC;QAExB,MAAM,GAAG,GAAG,wBAAwB,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC;QAC/C,MAAM,CAAC,EAAE,CAAC,GAAG,GAAG,CAAC,EAAE,iDAAiD,CAAC,CAAC;IACxE,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,8DAA8D,EAAE,GAAG,EAAE;QACxE,MAAM,EAAE,GAAG,cAAc,CAAC;QAC1B,MAAM,CAAC,GAAG,OAAO,CAAC,CAAC,sBAAsB;QACzC,MAAM,CAAC,GAAG,MAAM,CAAC,CAAC,gCAAgC;QAClD,MAAM,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,WAAW;QAEtD,MAAM,GAAG,GAAG,wBAAwB,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC;QAC/C,MAAM,QAAQ,GAAG,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;QAE/B,QAAQ,CAAC,GAAG,EAAE,QAAQ,EAAE,KAAK,CAAC,CAAC;IACjC,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,0DAA0D,EAAE,GAAG,EAAE;QACpE,MAAM,EAAE,GAAG,cAAc,CAAC;QAC1B,MAAM,CAAC,GAAG,MAAM,CAAC;QACjB,MAAM,EAAE,GAAG,IAAI,CAAC;QAChB,MAAM,EAAE,GAAG,IAAI,CAAC;QAEhB,MAAM,IAAI,GAAG,wBAAwB,CAAC,EAAE,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC;QACjD,MAAM,IAAI,GAAG,wBAAwB,CAAC,EAAE,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC;QAEjD,MAAM,CAAC,EAAE,CAAC,IAAI,GAAG,IAAI,CAAC,CAAC;IACzB,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,wEAAwE,EAAE,GAAG,EAAE;QAClF,MAAM,EAAE,GAAG,cAAc,CAAC;QAC1B,MAAM,CAAC,GAAG,IAAI,CAAC;QACf,MAAM,EAAE,GAAG,MAAM,CAAC;QAClB,MAAM,EAAE,GAAG,MAAM,CAAC;QAElB,MAAM,IAAI,GAAG,wBAAwB,CAAC,CAAC,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC;QACjD,MAAM,IAAI,GAAG,wBAAwB,CAAC,CAAC,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC;QAEjD,MAAM,CAAC,EAAE,CAAC,IAAI,GAAG,IAAI,CAAC,CAAC;IACzB,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,0BAA0B,EAAE,GAAG,EAAE;QACpC,MAAM,EAAE,GAAG,cAAc,CAAC;QAC1B,MAAM,CAAC,GAAG,MAAM,CAAC;QACjB,MAAM,CAAC,GAAG,IAAI,CAAC;QAEf,YAAY;QACZ,MAAM,CAAC,MAAM,CACX,GAAG,EAAE,CAAC,wBAAwB,CAAC,MAAM,CAAC,GAAG,EAAE,CAAC,EAAE,EAAE,CAAC,EACjD,cAAc,CACf,CAAC;QACF,MAAM,CAAC,MAAM,CAAC,GAAG,EAAE,CAAC,wBAAwB,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,EAAE,cAAc,CAAC,CAAC;QAEzE,YAAY;QACZ,MAAM,CAAC,MAAM,CACX,GAAG,EAAE,CAAC,wBAAwB,CAAC,CAAC,EAAE,MAAM,CAAC,GAAG,EAAE,EAAE,CAAC,EACjD,UAAU,CACX,CAAC;QACF,MAAM,CAAC,MAAM,CAAC,GAAG,EAAE,CAAC,wBAAwB,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,EAAE,UAAU,CAAC,CAAC;QACpE,MAAM,CAAC,MAAM,CAAC,GAAG,EAAE,CAAC,wBAAwB,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,EAAE,CAAC,EAAE,UAAU,CAAC,CAAC;QAErE,aAAa;QACb,MAAM,CAAC,MAAM,CACX,GAAG,EAAE,CAAC,wBAAwB,CAAC,CAAC,EAAE,CAAC,EAAE,MAAM,CAAC,GAAG,CAAC,EAChD,cAAc,CACf,CAAC;QACF,MAAM,CAAC,MAAM,CAAC,GAAG,EAAE,CAAC,wBAAwB,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,cAAc,CAAC,CAAC;IAC1E,CAAC,CAAC,CAAC;AACL,CAAC,CAAC,CAAC"}

package/dist/categories/orbits/__tests__/sphere-of-influence-radius.spec.d.ts

@@ -1,2 +0,0 @@
-export {};
-//# sourceMappingURL=sphere-of-influence-radius.spec.d.ts.map

package/dist/categories/orbits/__tests__/sphere-of-influence-radius.spec.d.ts.map

@@ -1 +0,0 @@
-{"version":3,"file":"sphere-of-influence-radius.spec.d.ts","sourceRoot":"","sources":["../../../../src/categories/orbits/__tests__/sphere-of-influence-radius.spec.ts"],"names":[],"mappings":""}

package/dist/categories/orbits/__tests__/sphere-of-influence-radius.spec.js

@@ -1,62 +0,0 @@
-import assert from 'node:assert/strict';
-import test, { describe } from 'node:test';
-import { sphereOfInfluenceRadius } from '../sphere-of-influence-radius';
-import { relClose } from 'packages/equations/src/__tests__/helpers';
-describe('sphereOfInfluenceRadius', () => {
-    test('computes rSOI = a * (m/M)^(2/5) for typical inputs (Earth around Sun, km-based)', () => {
-        const aKm = 149597870.7; // km
-        const mEarth = 5.972e24; // kg
-        const mSun = 1.9885e30; // kg
-        const rSoi = sphereOfInfluenceRadius(aKm, mEarth, mSun);
-        const expected = aKm * Math.pow(mEarth / mSun, 2 / 5);
-        relClose(rSoi, expected);
-        assert.ok(rSoi > 0);
-    });
-    test('scales linearly with a (doubling a doubles rSOI)', () => {
-        const a = 10;
-        const m = 2;
-        const M = 100;
-        const r1 = sphereOfInfluenceRadius(a, m, M);
-        const r2 = sphereOfInfluenceRadius(2 * a, m, M);
-        relClose(r2, 2 * r1);
-    });
-    test('increases with secondary mass m (all else equal)', () => {
-        const a = 10;
-        const M = 100;
-        const rSmall = sphereOfInfluenceRadius(a, 1, M);
-        const rLarge = sphereOfInfluenceRadius(a, 10, M);
-        assert.ok(rLarge > rSmall, `expected rSOI to increase with m: ${rSmall} -> ${rLarge}`);
-    });
-    test('decreases with primary mass M (all else equal)', () => {
-        const a = 10;
-        const m = 10;
-        const rSmallM = sphereOfInfluenceRadius(a, m, 100);
-        const rLargeM = sphereOfInfluenceRadius(a, m, 1000);
-        assert.ok(rLargeM < rSmallM, `expected rSOI to decrease with M: ${rSmallM} -> ${rLargeM}`);
-    });
-    test('throws TypeError when a is not finite (NaN)', () => {
-        assert.throws(() => sphereOfInfluenceRadius(Number.NaN, 1, 1), (err) => err instanceof TypeError &&
-            err.message === 'a must be a finite number. Received: NaN');
-    });
-    test('throws TypeError when m is not finite (Infinity)', () => {
-        assert.throws(() => sphereOfInfluenceRadius(1, Number.POSITIVE_INFINITY, 1), (err) => err instanceof TypeError &&
-            err.message === 'm must be a finite number. Received: Infinity');
-    });
-    test('throws TypeError when M is not finite (-Infinity)', () => {
-        assert.throws(() => sphereOfInfluenceRadius(1, 1, Number.NEGATIVE_INFINITY), (err) => err instanceof TypeError &&
-            err.message === 'M must be a finite number. Received: -Infinity');
-    });
-    test('throws RangeError when a <= 0', () => {
-        assert.throws(() => sphereOfInfluenceRadius(0, 1, 1), (err) => err instanceof RangeError &&
-            err.message === 'a must be > 0. Received: 0');
-    });
-    test('throws RangeError when m <= 0', () => {
-        assert.throws(() => sphereOfInfluenceRadius(1, 0, 1), (err) => err instanceof RangeError &&
-            err.message === 'm must be > 0. Received: 0');
-    });
-    test('throws RangeError when M <= 0', () => {
-        assert.throws(() => sphereOfInfluenceRadius(1, 1, 0), (err) => err instanceof RangeError &&
-            err.message === 'M must be > 0. Received: 0');
-    });
-});
-//# sourceMappingURL=sphere-of-influence-radius.spec.js.map

package/dist/categories/orbits/__tests__/sphere-of-influence-radius.spec.js.map

@@ -1 +0,0 @@
-{"version":3,"file":"sphere-of-influence-radius.spec.js","sourceRoot":"","sources":["../../../../src/categories/orbits/__tests__/sphere-of-influence-radius.spec.ts"],"names":[],"mappings":"AAAA,OAAO,MAAM,MAAM,oBAAoB,CAAC;AACxC,OAAO,IAAI,EAAE,EAAE,QAAQ,EAAE,MAAM,WAAW,CAAC;AAC3C,OAAO,EAAE,uBAAuB,EAAE,MAAM,+BAA+B,CAAC;AACxE,OAAO,EAAE,QAAQ,EAAE,MAAM,0CAA0C,CAAC;AAEpE,QAAQ,CAAC,yBAAyB,EAAE,GAAG,EAAE;IACvC,IAAI,CAAC,iFAAiF,EAAE,GAAG,EAAE;QAC3F,MAAM,GAAG,GAAG,WAAa,CAAC,CAAC,KAAK;QAChC,MAAM,MAAM,GAAG,QAAQ,CAAC,CAAC,KAAK;QAC9B,MAAM,IAAI,GAAG,SAAS,CAAC,CAAC,KAAK;QAE7B,MAAM,IAAI,GAAG,uBAAuB,CAAC,GAAG,EAAE,MAAM,EAAE,IAAI,CAAC,CAAC;QAExD,MAAM,QAAQ,GAAG,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,MAAM,GAAG,IAAI,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC;QACtD,QAAQ,CAAC,IAAI,EAAE,QAAQ,CAAC,CAAC;QACzB,MAAM,CAAC,EAAE,CAAC,IAAI,GAAG,CAAC,CAAC,CAAC;IACtB,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,kDAAkD,EAAE,GAAG,EAAE;QAC5D,MAAM,CAAC,GAAG,EAAE,CAAC;QACb,MAAM,CAAC,GAAG,CAAC,CAAC;QACZ,MAAM,CAAC,GAAG,GAAG,CAAC;QAEd,MAAM,EAAE,GAAG,uBAAuB,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;QAC5C,MAAM,EAAE,GAAG,uBAAuB,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;QAEhD,QAAQ,CAAC,EAAE,EAAE,CAAC,GAAG,EAAE,CAAC,CAAC;IACvB,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,kDAAkD,EAAE,GAAG,EAAE;QAC5D,MAAM,CAAC,GAAG,EAAE,CAAC;QACb,MAAM,CAAC,GAAG,GAAG,CAAC;QAEd,MAAM,MAAM,GAAG,uBAAuB,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;QAChD,MAAM,MAAM,GAAG,uBAAuB,CAAC,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC,CAAC;QAEjD,MAAM,CAAC,EAAE,CACP,MAAM,GAAG,MAAM,EACf,qCAAqC,MAAM,OAAO,MAAM,EAAE,CAC3D,CAAC;IACJ,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,gDAAgD,EAAE,GAAG,EAAE;QAC1D,MAAM,CAAC,GAAG,EAAE,CAAC;QACb,MAAM,CAAC,GAAG,EAAE,CAAC;QAEb,MAAM,OAAO,GAAG,uBAAuB,CAAC,CAAC,EAAE,CAAC,EAAE,GAAG,CAAC,CAAC;QACnD,MAAM,OAAO,GAAG,uBAAuB,CAAC,CAAC,EAAE,CAAC,EAAE,IAAI,CAAC,CAAC;QAEpD,MAAM,CAAC,EAAE,CACP,OAAO,GAAG,OAAO,EACjB,qCAAqC,OAAO,OAAO,OAAO,EAAE,CAC7D,CAAC;IACJ,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,6CAA6C,EAAE,GAAG,EAAE;QACvD,MAAM,CAAC,MAAM,CACX,GAAG,EAAE,CAAC,uBAAuB,CAAC,MAAM,CAAC,GAAG,EAAE,CAAC,EAAE,CAAC,CAAC,EAC/C,CAAC,GAAG,EAAE,EAAE,CACN,GAAG,YAAY,SAAS;YACxB,GAAG,CAAC,OAAO,KAAK,0CAA0C,CAC7D,CAAC;IACJ,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,kDAAkD,EAAE,GAAG,EAAE;QAC5D,MAAM,CAAC,MAAM,CACX,GAAG,EAAE,CAAC,uBAAuB,CAAC,CAAC,EAAE,MAAM,CAAC,iBAAiB,EAAE,CAAC,CAAC,EAC7D,CAAC,GAAG,EAAE,EAAE,CACN,GAAG,YAAY,SAAS;YACxB,GAAG,CAAC,OAAO,KAAK,+CAA+C,CAClE,CAAC;IACJ,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,mDAAmD,EAAE,GAAG,EAAE;QAC7D,MAAM,CAAC,MAAM,CACX,GAAG,EAAE,CAAC,uBAAuB,CAAC,CAAC,EAAE,CAAC,EAAE,MAAM,CAAC,iBAAiB,CAAC,EAC7D,CAAC,GAAG,EAAE,EAAE,CACN,GAAG,YAAY,SAAS;YACxB,GAAG,CAAC,OAAO,KAAK,gDAAgD,CACnE,CAAC;IACJ,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,+BAA+B,EAAE,GAAG,EAAE;QACzC,MAAM,CAAC,MAAM,CACX,GAAG,EAAE,CAAC,uBAAuB,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,EACtC,CAAC,GAAG,EAAE,EAAE,CACN,GAAG,YAAY,UAAU;YACzB,GAAG,CAAC,OAAO,KAAK,4BAA4B,CAC/C,CAAC;IACJ,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,+BAA+B,EAAE,GAAG,EAAE;QACzC,MAAM,CAAC,MAAM,CACX,GAAG,EAAE,CAAC,uBAAuB,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,EACtC,CAAC,GAAG,EAAE,EAAE,CACN,GAAG,YAAY,UAAU;YACzB,GAAG,CAAC,OAAO,KAAK,4BAA4B,CAC/C,CAAC;IACJ,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,+BAA+B,EAAE,GAAG,EAAE;QACzC,MAAM,CAAC,MAAM,CACX,GAAG,EAAE,CAAC,uBAAuB,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,EACtC,CAAC,GAAG,EAAE,EAAE,CACN,GAAG,YAAY,UAAU;YACzB,GAAG,CAAC,OAAO,KAAK,4BAA4B,CAC/C,CAAC;IACJ,CAAC,CAAC,CAAC;AACL,CAAC,CAAC,CAAC"}

package/dist/categories/orbits/__tests__/vis-viva-speed.spec.d.ts

@@ -1,2 +0,0 @@
-export {};
-//# sourceMappingURL=vis-viva-speed.spec.d.ts.map

package/dist/categories/orbits/__tests__/vis-viva-speed.spec.d.ts.map

@@ -1 +0,0 @@
-{"version":3,"file":"vis-viva-speed.spec.d.ts","sourceRoot":"","sources":["../../../../src/categories/orbits/__tests__/vis-viva-speed.spec.ts"],"names":[],"mappings":""}

package/dist/categories/orbits/__tests__/vis-viva-speed.spec.js

@@ -1,80 +0,0 @@
-import assert from 'node:assert/strict';
-import test, { describe } from 'node:test';
-import { absClose, relClose } from 'packages/equations/src/__tests__/helpers';
-import { visVivaSpeed } from '../vis-viva-speed';
-describe('visVivaSpeed', () => {
-    test('circular orbit: a = r → v = sqrt(μ/r)', () => {
-        const mu = 3.986004418e14; // m³/s² (Earth GM)
-        const r = 6778e3; // m (LEO-ish)
-        const a = r;
-        const v = visVivaSpeed(r, a, mu);
-        const expected = Math.sqrt(mu / r);
-        relClose(v, expected, 1e-12);
-    });
-    test('parabolic trajectory: a = Infinity → v = sqrt(2μ/r)', () => {
-        const mu = 3.986004418e14;
-        const r = 7000e3;
-        const v = visVivaSpeed(r, Infinity, mu);
-        const expected = Math.sqrt((2 * mu) / r);
-        relClose(v, expected, 1e-12);
-    });
-    test('hyperbolic case: a < 0 uses same formula and yields valid speed', () => {
-        const mu = 3.986004418e14;
-        const r = 10000e3;
-        const a = -15000e3; // negative semi-major axis
-        const v = visVivaSpeed(r, a, mu);
-        const expected = Math.sqrt(mu * (2 / r - 1 / a));
-        relClose(v, expected, 1e-12);
-    });
-    test('radicand zero boundary: choose a such that 2/r - 1/a = 0 → v = 0', () => {
-        const mu = 3.986004418e14;
-        const r = 8000e3;
-        const a = r / 2; // 1/a = 2/r
-        const v = visVivaSpeed(r, a, mu);
-        absClose(v, 0, 1e-12);
-    });
-    test('tiny negative radicand (within tolerance) is clamped to v ≈ 0', () => {
-        const mu = 3.986004418e14;
-        const r = 8000e3;
-        // Function's tolerance: tol = |1e-14 * mu / r|
-        const tol = Math.abs((1e-14 * mu) / r);
-        const delta = tol / (2 * mu); // makes radicand = μ*(2/r - 1/a) = -tol/2
-        // Choose 1/a = 2/r + delta  → a = 1 / (2/r + delta)
-        const a = 1 / (2 / r + delta);
-        const v = visVivaSpeed(r, a, mu);
-        absClose(v, 0, 1e-12);
-    });
-    test('significantly negative radicand throws (physically impossible state)', () => {
-        const mu = 3.986004418e14;
-        const r = 8000e3;
-        const a = r / 4; // 1/a = 4/r -> 2/r - 1/a = -2/r (negative)
-        assert.throws(() => visVivaSpeed(r, a, mu), /Invalid state: μ\*\(2\/r - 1\/a\) < 0/);
-    });
-    test('monotonic with radius for fixed a, μ: larger r → smaller v (ellipse domain)', () => {
-        const mu = 3.986004418e14;
-        const a = 10000e3; // ellipse: r must be < 2a
-        const r1 = 8000e3;
-        const r2 = 15000e3; // < 2a = 20000 km
-        const v1 = visVivaSpeed(r1, a, mu);
-        const v2 = visVivaSpeed(r2, a, mu);
-        assert.ok(v1 > v2, 'speed should decrease as radius increases (for fixed a)');
-    });
-    test('input validation: r<=0, mu<0/NaN, a non-finite (except +Infinity), a=0', () => {
-        const mu = 3.986004418e14;
-        const r = 7000e3;
-        const a = r;
-        // r invalid
-        assert.throws(() => visVivaSpeed(0, a, mu), /positive/);
-        assert.throws(() => visVivaSpeed(-1, a, mu), /positive/);
-        assert.throws(() => visVivaSpeed(Number.NaN, a, mu), /finite/);
-        // mu invalid
-        assert.throws(() => visVivaSpeed(r, a, -1), /non-negative/);
-        assert.throws(() => visVivaSpeed(r, a, Number.NaN), /finite/);
-        // a invalid (non-finite other than +Infinity)
-        assert.throws(() => visVivaSpeed(r, Number.NaN, mu), /finite/);
-        assert.throws(() => visVivaSpeed(r, Number.NEGATIVE_INFINITY, mu), /finite/);
-        // a = 0 invalid
-        assert.throws(() => visVivaSpeed(r, 0, mu), /non-zero/);
-    });
-});
-//# sourceMappingURL=vis-viva-speed.spec.js.map

package/dist/categories/orbits/__tests__/vis-viva-speed.spec.js.map

@@ -1 +0,0 @@
-{"version":3,"file":"vis-viva-speed.spec.js","sourceRoot":"","sources":["../../../../src/categories/orbits/__tests__/vis-viva-speed.spec.ts"],"names":[],"mappings":"AAAA,OAAO,MAAM,MAAM,oBAAoB,CAAC;AACxC,OAAO,IAAI,EAAE,EAAE,QAAQ,EAAE,MAAM,WAAW,CAAC;AAC3C,OAAO,EAAE,QAAQ,EAAE,QAAQ,EAAE,MAAM,0CAA0C,CAAC;AAC9E,OAAO,EAAE,YAAY,EAAE,MAAM,mBAAmB,CAAC;AAEjD,QAAQ,CAAC,cAAc,EAAE,GAAG,EAAE;IAC5B,IAAI,CAAC,uCAAuC,EAAE,GAAG,EAAE;QACjD,MAAM,EAAE,GAAG,cAAc,CAAC,CAAC,mBAAmB;QAC9C,MAAM,CAAC,GAAG,MAAM,CAAC,CAAC,cAAc;QAChC,MAAM,CAAC,GAAG,CAAC,CAAC;QAEZ,MAAM,CAAC,GAAG,YAAY,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC;QACjC,MAAM,QAAQ,GAAG,IAAI,CAAC,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC;QAEnC,QAAQ,CAAC,CAAC,EAAE,QAAQ,EAAE,KAAK,CAAC,CAAC;IAC/B,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,qDAAqD,EAAE,GAAG,EAAE;QAC/D,MAAM,EAAE,GAAG,cAAc,CAAC;QAC1B,MAAM,CAAC,GAAG,MAAM,CAAC;QAEjB,MAAM,CAAC,GAAG,YAAY,CAAC,CAAC,EAAE,QAAQ,EAAE,EAAE,CAAC,CAAC;QACxC,MAAM,QAAQ,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC;QAEzC,QAAQ,CAAC,CAAC,EAAE,QAAQ,EAAE,KAAK,CAAC,CAAC;IAC/B,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,iEAAiE,EAAE,GAAG,EAAE;QAC3E,MAAM,EAAE,GAAG,cAAc,CAAC;QAC1B,MAAM,CAAC,GAAG,OAAO,CAAC;QAClB,MAAM,CAAC,GAAG,CAAC,OAAO,CAAC,CAAC,2BAA2B;QAE/C,MAAM,CAAC,GAAG,YAAY,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC;QACjC,MAAM,QAAQ,GAAG,IAAI,CAAC,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;QAEjD,QAAQ,CAAC,CAAC,EAAE,QAAQ,EAAE,KAAK,CAAC,CAAC;IAC/B,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,kEAAkE,EAAE,GAAG,EAAE;QAC5E,MAAM,EAAE,GAAG,cAAc,CAAC;QAC1B,MAAM,CAAC,GAAG,MAAM,CAAC;QACjB,MAAM,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,YAAY;QAE7B,MAAM,CAAC,GAAG,YAAY,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC;QACjC,QAAQ,CAAC,CAAC,EAAE,CAAC,EAAE,KAAK,CAAC,CAAC;IACxB,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,+DAA+D,EAAE,GAAG,EAAE;QACzE,MAAM,EAAE,GAAG,cAAc,CAAC;QAC1B,MAAM,CAAC,GAAG,MAAM,CAAC;QAEjB,+CAA+C;QAC/C,MAAM,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,KAAK,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC;QACvC,MAAM,KAAK,GAAG,GAAG,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,0CAA0C;QAExE,oDAAoD;QACpD,MAAM,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,KAAK,CAAC,CAAC;QAE9B,MAAM,CAAC,GAAG,YAAY,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC;QACjC,QAAQ,CAAC,CAAC,EAAE,CAAC,EAAE,KAAK,CAAC,CAAC;IACxB,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,sEAAsE,EAAE,GAAG,EAAE;QAChF,MAAM,EAAE,GAAG,cAAc,CAAC;QAC1B,MAAM,CAAC,GAAG,MAAM,CAAC;QACjB,MAAM,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,2CAA2C;QAE5D,MAAM,CAAC,MAAM,CACX,GAAG,EAAE,CAAC,YAAY,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,EAC5B,uCAAuC,CACxC,CAAC;IACJ,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,6EAA6E,EAAE,GAAG,EAAE;QACvF,MAAM,EAAE,GAAG,cAAc,CAAC;QAC1B,MAAM,CAAC,GAAG,OAAO,CAAC,CAAC,0BAA0B;QAC7C,MAAM,EAAE,GAAG,MAAM,CAAC;QAClB,MAAM,EAAE,GAAG,OAAO,CAAC,CAAC,kBAAkB;QAEtC,MAAM,EAAE,GAAG,YAAY,CAAC,EAAE,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC;QACnC,MAAM,EAAE,GAAG,YAAY,CAAC,EAAE,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC;QAEnC,MAAM,CAAC,EAAE,CACP,EAAE,GAAG,EAAE,EACP,yDAAyD,CAC1D,CAAC;IACJ,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,wEAAwE,EAAE,GAAG,EAAE;QAClF,MAAM,EAAE,GAAG,cAAc,CAAC;QAC1B,MAAM,CAAC,GAAG,MAAM,CAAC;QACjB,MAAM,CAAC,GAAG,CAAC,CAAC;QAEZ,YAAY;QACZ,MAAM,CAAC,MAAM,CAAC,GAAG,EAAE,CAAC,YAAY,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,EAAE,UAAU,CAAC,CAAC;QACxD,MAAM,CAAC,MAAM,CAAC,GAAG,EAAE,CAAC,YAAY,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,EAAE,UAAU,CAAC,CAAC;QACzD,MAAM,CAAC,MAAM,CAAC,GAAG,EAAE,CAAC,YAAY,CAAC,MAAM,CAAC,GAAG,EAAE,CAAC,EAAE,EAAE,CAAC,EAAE,QAAQ,CAAC,CAAC;QAE/D,aAAa;QACb,MAAM,CAAC,MAAM,CAAC,GAAG,EAAE,CAAC,YAAY,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,cAAc,CAAC,CAAC;QAC5D,MAAM,CAAC,MAAM,CAAC,GAAG,EAAE,CAAC,YAAY,CAAC,CAAC,EAAE,CAAC,EAAE,MAAM,CAAC,GAAG,CAAC,EAAE,QAAQ,CAAC,CAAC;QAE9D,8CAA8C;QAC9C,MAAM,CAAC,MAAM,CAAC,GAAG,EAAE,CAAC,YAAY,CAAC,CAAC,EAAE,MAAM,CAAC,GAAG,EAAE,EAAE,CAAC,EAAE,QAAQ,CAAC,CAAC;QAC/D,MAAM,CAAC,MAAM,CACX,GAAG,EAAE,CAAC,YAAY,CAAC,CAAC,EAAE,MAAM,CAAC,iBAAiB,EAAE,EAAE,CAAC,EACnD,QAAQ,CACT,CAAC;QAEF,gBAAgB;QAChB,MAAM,CAAC,MAAM,CAAC,GAAG,EAAE,CAAC,YAAY,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,EAAE,UAAU,CAAC,CAAC;IAC1D,CAAC,CAAC,CAAC;AACL,CAAC,CAAC,CAAC"}