@interstellar-tools/equations 0.1.2 → 0.2.0

package/dist/__tests__/gravitational-acceleration-on1-by2.int.spec.d.ts

@@ -0,0 +1,2 @@
+export {};
+//# sourceMappingURL=gravitational-acceleration-on1-by2.int.spec.d.ts.map

package/dist/__tests__/gravitational-acceleration-on1-by2.int.spec.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"gravitational-acceleration-on1-by2.int.spec.d.ts","sourceRoot":"","sources":["../../src/__tests__/gravitational-acceleration-on1-by2.int.spec.ts"],"names":[],"mappings":""}

package/dist/__tests__/gravitational-acceleration-on1-by2.int.spec.js

@@ -0,0 +1,86 @@
+import assert from 'node:assert/strict';
+import test, { describe } from 'node:test';
+import { gravitationalAccelerationOn1By2 } from '../law-of-gravitation';
+import { norm, relClose, scale, sub, vecRelClose } from './helpers';
+import { G_SI } from '@interstellar-tools/constants';
+describe('gravitationalAccelerationOn1By2', () => {
+    test('1D sanity (+x): ||a|| = G*m2/r^2 and points +x', () => {
+        const G = 1;
+        const m2 = 3;
+        const r1 = [0, 0, 0];
+        const r2 = [10, 0, 0];
+        const a = gravitationalAccelerationOn1By2(m2, r1, r2, G);
+        // Expected magnitude: 1*3/10^2 = 0.03, direction +x
+        relClose(norm(a), 0.03, 1e-14);
+        vecRelClose(a, [0.03, 0, 0], 1e-14);
+    });
+    test('1D sanity (-x): points toward body 2 at negative x', () => {
+        const G = 1;
+        const m2 = 3;
+        const r1 = [0, 0, 0];
+        const r2 = [-10, 0, 0];
+        const a = gravitationalAccelerationOn1By2(m2, r1, r2, G);
+        relClose(norm(a), 0.03, 1e-14);
+        vecRelClose(a, [-0.03, 0, 0], 1e-14);
+    });
+    test('direction equals r-hat for arbitrary 3D positions', () => {
+        const G = 2.5;
+        const m2 = 11;
+        const r1 = [0.2, -3, 5.5];
+        const r2 = [4.2, 1, -1];
+        const a = gravitationalAccelerationOn1By2(m2, r1, r2, G);
+        const dr = sub(r2, r1);
+        const r = Math.hypot(...dr);
+        const rhat = [dr[0] / r, dr[1] / r, dr[2] / r];
+        const ahat = scale(a, 1 / norm(a));
+        vecRelClose(ahat, rhat, 1e-12);
+    });
+    test('inverse-square scaling: doubling r → quarter acceleration', () => {
+        const G = 10;
+        const m2 = 5;
+        const r1 = [0, 0, 0];
+        const r2a = [2, 0, 0]; // r = 2
+        const r2b = [4, 0, 0]; // r = 4 (double → quarter a)
+        const aA = gravitationalAccelerationOn1By2(m2, r1, r2a, G);
+        const aB = gravitationalAccelerationOn1By2(m2, r1, r2b, G);
+        relClose(norm(aB), norm(aA) / 4, 1e-14);
+    });
+    test('linearity: scaling m2 or G scales a accordingly', () => {
+        const G = 3;
+        const m2 = 5;
+        const r1 = [1, 2, 3];
+        const r2 = [4, -1, 0];
+        const a = gravitationalAccelerationOn1By2(m2, r1, r2, G);
+        const a_m2x2 = gravitationalAccelerationOn1By2(2 * m2, r1, r2, G);
+        const a_Gx2 = gravitationalAccelerationOn1By2(m2, r1, r2, 2 * G);
+        vecRelClose(a_m2x2, scale(a, 2));
+        vecRelClose(a_Gx2, scale(a, 2));
+    });
+    test('real-world sanity: Earth due to Sun at 1 AU (~5.93e-3 m/s^2)', () => {
+        const mSun = 1.9885e30; // kg
+        const rEarth = [0, 0, 0];
+        const rSun = [1.495978707e11, 0, 0]; // m
+        const a = gravitationalAccelerationOn1By2(mSun, rEarth, rSun, G_SI);
+        // Closed-form expectation GM_sun / r^2
+        const dr = sub(rSun, rEarth);
+        const r = Math.hypot(...dr);
+        const expected = (G_SI * mSun) / (r * r);
+        relClose(norm(a), expected, 1e-12);
+        // Also within ~0.5% of commonly cited value
+        const ref = 5.93e-3;
+        assert.ok(Math.abs(norm(a) - ref) / ref < 5e-3, `|a| ${norm(a)} differs from ~${ref}`);
+    });
+    test('error propagation: invalid inputs cause throws', () => {
+        const r1 = [0, 0, 0];
+        const r2 = [1, 0, 0];
+        // negative / non-finite m2
+        assert.throws(() => gravitationalAccelerationOn1By2(-1, r1, r2, 1), /non-negative|finite/);
+        assert.throws(() => gravitationalAccelerationOn1By2(Number.NaN, r1, r2, 1), /finite/);
+        // coincident positions
+        assert.throws(() => gravitationalAccelerationOn1By2(1, r1, r1, 1), /singular \(r = 0\)/);
+        // invalid G
+        assert.throws(() => gravitationalAccelerationOn1By2(1, r1, r2, 0), /positive/);
+        assert.throws(() => gravitationalAccelerationOn1By2(1, r1, r2, Number.NaN), /finite/);
+    });
+});
+//# sourceMappingURL=gravitational-acceleration-on1-by2.int.spec.js.map

package/dist/__tests__/gravitational-acceleration-on1-by2.int.spec.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"gravitational-acceleration-on1-by2.int.spec.js","sourceRoot":"","sources":["../../src/__tests__/gravitational-acceleration-on1-by2.int.spec.ts"],"names":[],"mappings":"AAAA,OAAO,MAAM,MAAM,oBAAoB,CAAC;AACxC,OAAO,IAAI,EAAE,EAAE,QAAQ,EAAE,MAAM,WAAW,CAAC;AAC3C,OAAO,EAAE,+BAA+B,EAAE,MAAM,uBAAuB,CAAC;AACxE,OAAO,EAAE,IAAI,EAAE,QAAQ,EAAE,KAAK,EAAE,GAAG,EAAE,WAAW,EAAE,MAAM,WAAW,CAAC;AACpE,OAAO,EAAE,IAAI,EAAE,MAAM,+BAA+B,CAAC;AAGrD,QAAQ,CAAC,iCAAiC,EAAE,GAAG,EAAE;IAC/C,IAAI,CAAC,gDAAgD,EAAE,GAAG,EAAE;QAC1D,MAAM,CAAC,GAAG,CAAC,CAAC;QACZ,MAAM,EAAE,GAAG,CAAC,CAAC;QACb,MAAM,EAAE,GAAsB,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;QACxC,MAAM,EAAE,GAAsB,CAAC,EAAE,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;QAEzC,MAAM,CAAC,GAAG,+BAA+B,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC,CAAC;QAEzD,oDAAoD;QACpD,QAAQ,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,IAAI,EAAE,KAAK,CAAC,CAAC;QAC/B,WAAW,CAAC,CAAC,EAAE,CAAC,IAAI,EAAE,CAAC,EAAE,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC;IACtC,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,oDAAoD,EAAE,GAAG,EAAE;QAC9D,MAAM,CAAC,GAAG,CAAC,CAAC;QACZ,MAAM,EAAE,GAAG,CAAC,CAAC;QACb,MAAM,EAAE,GAAsB,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;QACxC,MAAM,EAAE,GAAsB,CAAC,CAAC,EAAE,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;QAE1C,MAAM,CAAC,GAAG,+BAA+B,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC,CAAC;QAEzD,QAAQ,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,IAAI,EAAE,KAAK,CAAC,CAAC;QAC/B,WAAW,CAAC,CAAC,EAAE,CAAC,CAAC,IAAI,EAAE,CAAC,EAAE,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC;IACvC,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,mDAAmD,EAAE,GAAG,EAAE;QAC7D,MAAM,CAAC,GAAG,GAAG,CAAC;QACd,MAAM,EAAE,GAAG,EAAE,CAAC;QACd,MAAM,EAAE,GAAsB,CAAC,GAAG,EAAE,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC;QAC7C,MAAM,EAAE,GAAsB,CAAC,GAAG,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;QAE3C,MAAM,CAAC,GAAG,+BAA+B,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC,CAAC;QACzD,MAAM,EAAE,GAAG,GAAG,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC;QACvB,MAAM,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,GAAG,EAAE,CAAC,CAAC;QAC5B,MAAM,IAAI,GAAsB,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;QAElE,MAAM,IAAI,GAAG,KAAK,CAAC,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,CAAsB,CAAC;QACxD,WAAW,CAAC,IAAI,EAAE,IAAI,EAAE,KAAK,CAAC,CAAC;IACjC,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,2DAA2D,EAAE,GAAG,EAAE;QACrE,MAAM,CAAC,GAAG,EAAE,CAAC;QACb,MAAM,EAAE,GAAG,CAAC,CAAC;QACb,MAAM,EAAE,GAAsB,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;QAExC,MAAM,GAAG,GAAsB,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,QAAQ;QAClD,MAAM,GAAG,GAAsB,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,6BAA6B;QAEvE,MAAM,EAAE,GAAG,+BAA+B,CAAC,EAAE,EAAE,EAAE,EAAE,GAAG,EAAE,CAAC,CAAC,CAAC;QAC3D,MAAM,EAAE,GAAG,+BAA+B,CAAC,EAAE,EAAE,EAAE,EAAE,GAAG,EAAE,CAAC,CAAC,CAAC;QAE3D,QAAQ,CAAC,IAAI,CAAC,EAAE,CAAC,EAAE,IAAI,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,KAAK,CAAC,CAAC;IAC1C,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,iDAAiD,EAAE,GAAG,EAAE;QAC3D,MAAM,CAAC,GAAG,CAAC,CAAC;QACZ,MAAM,EAAE,GAAG,CAAC,CAAC;QACb,MAAM,EAAE,GAAsB,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;QACxC,MAAM,EAAE,GAAsB,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;QAEzC,MAAM,CAAC,GAAG,+BAA+B,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC,CAAC;QACzD,MAAM,MAAM,GAAG,+BAA+B,CAAC,CAAC,GAAG,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC,CAAC;QAClE,MAAM,KAAK,GAAG,+BAA+B,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC;QAEjE,WAAW,CAAC,MAAM,EAAE,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;QACjC,WAAW,CAAC,KAAK,EAAE,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;IAClC,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,8DAA8D,EAAE,GAAG,EAAE;QACxE,MAAM,IAAI,GAAG,SAAS,CAAC,CAAC,KAAK;QAC7B,MAAM,MAAM,GAAsB,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;QAC5C,MAAM,IAAI,GAAsB,CAAC,cAAc,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,IAAI;QAE5D,MAAM,CAAC,GAAG,+BAA+B,CAAC,IAAI,EAAE,MAAM,EAAE,IAAI,EAAE,IAAI,CAAC,CAAC;QAEpE,uCAAuC;QACvC,MAAM,EAAE,GAAG,GAAG,CAAC,IAAI,EAAE,MAAM,CAAC,CAAC;QAC7B,MAAM,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,GAAG,EAAE,CAAC,CAAC;QAC5B,MAAM,QAAQ,GAAG,CAAC,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;QAEzC,QAAQ,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,QAAQ,EAAE,KAAK,CAAC,CAAC;QAEnC,4CAA4C;QAC5C,MAAM,GAAG,GAAG,OAAO,CAAC;QACpB,MAAM,CAAC,EAAE,CACP,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,GAAG,GAAG,GAAG,IAAI,EACpC,OAAO,IAAI,CAAC,CAAC,CAAC,kBAAkB,GAAG,EAAE,CACtC,CAAC;IACJ,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,gDAAgD,EAAE,GAAG,EAAE;QAC1D,MAAM,EAAE,GAAsB,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;QACxC,MAAM,EAAE,GAAsB,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;QAExC,2BAA2B;QAC3B,MAAM,CAAC,MAAM,CACX,GAAG,EAAE,CAAC,+BAA+B,CAAC,CAAC,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC,EACpD,qBAAqB,CACtB,CAAC;QACF,MAAM,CAAC,MAAM,CACX,GAAG,EAAE,CAAC,+BAA+B,CAAC,MAAM,CAAC,GAAG,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC,EAC5D,QAAQ,CACT,CAAC;QAEF,uBAAuB;QACvB,MAAM,CAAC,MAAM,CACX,GAAG,EAAE,CAAC,+BAA+B,CAAC,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC,EACnD,oBAAoB,CACrB,CAAC;QAEF,YAAY;QACZ,MAAM,CAAC,MAAM,CACX,GAAG,EAAE,CAAC,+BAA+B,CAAC,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC,EACnD,UAAU,CACX,CAAC;QACF,MAAM,CAAC,MAAM,CACX,GAAG,EAAE,CAAC,+BAA+B,CAAC,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,MAAM,CAAC,GAAG,CAAC,EAC5D,QAAQ,CACT,CAAC;IACJ,CAAC,CAAC,CAAC;AACL,CAAC,CAAC,CAAC"}

package/dist/__tests__/gravitational-force-on1-by2.int.spec.d.ts

@@ -0,0 +1,2 @@
+export {};
+//# sourceMappingURL=gravitational-force-on1-by2.int.spec.d.ts.map

package/dist/__tests__/gravitational-force-on1-by2.int.spec.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"gravitational-force-on1-by2.int.spec.d.ts","sourceRoot":"","sources":["../../src/__tests__/gravitational-force-on1-by2.int.spec.ts"],"names":[],"mappings":""}

package/dist/__tests__/gravitational-force-on1-by2.int.spec.js

@@ -0,0 +1,106 @@
+import assert from 'node:assert/strict';
+import test, { describe } from 'node:test';
+import { gravitationalForceOn1By2 } from '../law-of-gravitation';
+import { norm, relClose, scale, sub, vecRelClose } from './helpers';
+import { G_SI } from '@interstellar-tools/constants';
+describe('gravitationalForceOn1By2', () => {
+    test('1D sanity (+x): |F| = G*m1*m2/r^2 and points +x', () => {
+        const G = 1;
+        const m1 = 2;
+        const m2 = 3;
+        const r1 = [0, 0, 0];
+        const r2 = [10, 0, 0];
+        const F = gravitationalForceOn1By2(m1, m2, r1, r2, G);
+        // Expected magnitude and direction
+        const expectedMag = (G * m1 * m2) / 10 ** 2; // 0.06
+        relClose(norm(F), expectedMag, 1e-14);
+        vecRelClose(F, [expectedMag, 0, 0], 1e-14);
+    });
+    test('1D sanity (-x): points toward body 2 (at -x)', () => {
+        const G = 1;
+        const m1 = 2;
+        const m2 = 3;
+        const r1 = [0, 0, 0];
+        const r2 = [-10, 0, 0];
+        const F = gravitationalForceOn1By2(m1, m2, r1, r2, G);
+        relClose(norm(F), 0.06, 1e-14);
+        vecRelClose(F, [-0.06, 0, 0], 1e-14);
+    });
+    test('direction equals r-hat for arbitrary 3D positions', () => {
+        const G = 2.5;
+        const m1 = 7;
+        const m2 = 11;
+        const r1 = [0.2, -3, 5.5];
+        const r2 = [4.2, 1, -1];
+        const F = gravitationalForceOn1By2(m1, m2, r1, r2, G);
+        const dr = sub(r2, r1);
+        const r = Math.hypot(...dr);
+        const rhat = [dr[0] / r, dr[1] / r, dr[2] / r];
+        // Normalize F and compare to rhat
+        const Fhat = scale(F, 1 / norm(F));
+        vecRelClose(Fhat, rhat, 1e-12);
+    });
+    test('Newton’s 3rd law: F_1<-2 = - F_2<-1', () => {
+        const G = 9.87;
+        const m1 = 7;
+        const m2 = 11;
+        const r1 = [0, -2, 5];
+        const r2 = [3, 1, -1];
+        const F12 = gravitationalForceOn1By2(m1, m2, r1, r2, G);
+        const F21 = gravitationalForceOn1By2(m2, m1, r2, r1, G);
+        vecRelClose(F12, [-F21[0], -F21[1], -F21[2]]);
+    });
+    test('inverse-square scaling: doubling r → quarter force', () => {
+        const G = 10;
+        const m1 = 2;
+        const m2 = 5;
+        const r1 = [0, 0, 0];
+        const r2a = [2, 0, 0]; // r = 2
+        const r2b = [4, 0, 0]; // r = 4 (double → quarter force)
+        const Fa = gravitationalForceOn1By2(m1, m2, r1, r2a, G);
+        const Fb = gravitationalForceOn1By2(m1, m2, r1, r2b, G);
+        relClose(norm(Fb), norm(Fa) / 4, 1e-14);
+    });
+    test('linearity: scaling m1, m2, or G scales F accordingly', () => {
+        const G = 3;
+        const m1 = 2;
+        const m2 = 5;
+        const r1 = [1, 2, 3];
+        const r2 = [4, -1, 0];
+        const F = gravitationalForceOn1By2(m1, m2, r1, r2, G);
+        const F_m1x2 = gravitationalForceOn1By2(2 * m1, m2, r1, r2, G);
+        const F_m2x2 = gravitationalForceOn1By2(m1, 2 * m2, r1, r2, G);
+        const F_Gx2 = gravitationalForceOn1By2(m1, m2, r1, r2, 2 * G);
+        // Doubling any of m1, m2, or G doubles the force vector
+        vecRelClose(F_m1x2, scale(F, 2));
+        vecRelClose(F_m2x2, scale(F, 2));
+        vecRelClose(F_Gx2, scale(F, 2));
+    });
+    test('real-world sanity: Earth–Sun at ~1 AU (~3.54e22 N)', () => {
+        const mEarth = 5.972e24; // kg
+        const mSun = 1.9885e30; // kg
+        const rEarth = [0, 0, 0];
+        const rSun = [1.495978707e11, 0, 0]; // m
+        const F = gravitationalForceOn1By2(mEarth, mSun, rEarth, rSun, G_SI);
+        const r = Math.hypot(rSun[0] - rEarth[0], rSun[1] - rEarth[1], rSun[2] - rEarth[2]);
+        const expectedMag = (G_SI * mEarth * mSun) / (r * r);
+        relClose(norm(F), expectedMag, 1e-12);
+        // also within ~0.5% of the commonly quoted magnitude
+        const ref = 3.542e22;
+        assert.ok(Math.abs(norm(F) - ref) / ref < 5e-3, `|F| ${norm(F)} differs from ~${ref}`);
+    });
+    test('error propagation: invalid inputs cause throws', () => {
+        const r1 = [0, 0, 0];
+        const r2 = [1, 0, 0];
+        // negative / non-finite masses
+        assert.throws(() => gravitationalForceOn1By2(-1, 1, r1, r2, 1), /non-negative|finite/);
+        assert.throws(() => gravitationalForceOn1By2(1, -1, r1, r2, 1), /non-negative|finite/);
+        assert.throws(() => gravitationalForceOn1By2(Number.NaN, 1, r1, r2, 1), /finite/);
+        // coincident positions
+        assert.throws(() => gravitationalForceOn1By2(1, 1, r1, r1, 1), /singular \(r = 0\)/);
+        // invalid G
+        assert.throws(() => gravitationalForceOn1By2(1, 1, r1, r2, 0), /positive/);
+        assert.throws(() => gravitationalForceOn1By2(1, 1, r1, r2, Number.NaN), /finite/);
+    });
+});
+//# sourceMappingURL=gravitational-force-on1-by2.int.spec.js.map

package/dist/__tests__/gravitational-force-on1-by2.int.spec.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"gravitational-force-on1-by2.int.spec.js","sourceRoot":"","sources":["../../src/__tests__/gravitational-force-on1-by2.int.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,uBAAuB,CAAC;AACjE,OAAO,EAAE,IAAI,EAAE,QAAQ,EAAE,KAAK,EAAE,GAAG,EAAE,WAAW,EAAE,MAAM,WAAW,CAAC;AACpE,OAAO,EAAE,IAAI,EAAE,MAAM,+BAA+B,CAAC;AAGrD,QAAQ,CAAC,0BAA0B,EAAE,GAAG,EAAE;IACxC,IAAI,CAAC,iDAAiD,EAAE,GAAG,EAAE;QAC3D,MAAM,CAAC,GAAG,CAAC,CAAC;QACZ,MAAM,EAAE,GAAG,CAAC,CAAC;QACb,MAAM,EAAE,GAAG,CAAC,CAAC;QACb,MAAM,EAAE,GAAsB,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;QACxC,MAAM,EAAE,GAAsB,CAAC,EAAE,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;QAEzC,MAAM,CAAC,GAAG,wBAAwB,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC,CAAC;QAEtD,mCAAmC;QACnC,MAAM,WAAW,GAAG,CAAC,CAAC,GAAG,EAAE,GAAG,EAAE,CAAC,GAAG,EAAE,IAAI,CAAC,CAAC,CAAC,OAAO;QACpD,QAAQ,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,WAAW,EAAE,KAAK,CAAC,CAAC;QACtC,WAAW,CAAC,CAAC,EAAE,CAAC,WAAW,EAAE,CAAC,EAAE,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC;IAC7C,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,8CAA8C,EAAE,GAAG,EAAE;QACxD,MAAM,CAAC,GAAG,CAAC,CAAC;QACZ,MAAM,EAAE,GAAG,CAAC,CAAC;QACb,MAAM,EAAE,GAAG,CAAC,CAAC;QACb,MAAM,EAAE,GAAsB,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;QACxC,MAAM,EAAE,GAAsB,CAAC,CAAC,EAAE,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;QAE1C,MAAM,CAAC,GAAG,wBAAwB,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC,CAAC;QACtD,QAAQ,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,IAAI,EAAE,KAAK,CAAC,CAAC;QAC/B,WAAW,CAAC,CAAC,EAAE,CAAC,CAAC,IAAI,EAAE,CAAC,EAAE,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC;IACvC,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,mDAAmD,EAAE,GAAG,EAAE;QAC7D,MAAM,CAAC,GAAG,GAAG,CAAC;QACd,MAAM,EAAE,GAAG,CAAC,CAAC;QACb,MAAM,EAAE,GAAG,EAAE,CAAC;QACd,MAAM,EAAE,GAAsB,CAAC,GAAG,EAAE,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC;QAC7C,MAAM,EAAE,GAAsB,CAAC,GAAG,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;QAE3C,MAAM,CAAC,GAAG,wBAAwB,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC,CAAC;QACtD,MAAM,EAAE,GAAG,GAAG,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC;QACvB,MAAM,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,GAAG,EAAE,CAAC,CAAC;QAC5B,MAAM,IAAI,GAAsB,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;QAElE,kCAAkC;QAClC,MAAM,IAAI,GAAG,KAAK,CAAC,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,CAAsB,CAAC;QACxD,WAAW,CAAC,IAAI,EAAE,IAAI,EAAE,KAAK,CAAC,CAAC;IACjC,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,qCAAqC,EAAE,GAAG,EAAE;QAC/C,MAAM,CAAC,GAAG,IAAI,CAAC;QACf,MAAM,EAAE,GAAG,CAAC,CAAC;QACb,MAAM,EAAE,GAAG,EAAE,CAAC;QACd,MAAM,EAAE,GAAsB,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;QACzC,MAAM,EAAE,GAAsB,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;QAEzC,MAAM,GAAG,GAAG,wBAAwB,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC,CAAC;QACxD,MAAM,GAAG,GAAG,wBAAwB,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC,CAAC;QAExD,WAAW,CAAC,GAAG,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IAChD,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,oDAAoD,EAAE,GAAG,EAAE;QAC9D,MAAM,CAAC,GAAG,EAAE,CAAC;QACb,MAAM,EAAE,GAAG,CAAC,CAAC;QACb,MAAM,EAAE,GAAG,CAAC,CAAC;QACb,MAAM,EAAE,GAAsB,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;QAExC,MAAM,GAAG,GAAsB,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,QAAQ;QAClD,MAAM,GAAG,GAAsB,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,iCAAiC;QAE3E,MAAM,EAAE,GAAG,wBAAwB,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,GAAG,EAAE,CAAC,CAAC,CAAC;QACxD,MAAM,EAAE,GAAG,wBAAwB,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,GAAG,EAAE,CAAC,CAAC,CAAC;QAExD,QAAQ,CAAC,IAAI,CAAC,EAAE,CAAC,EAAE,IAAI,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,KAAK,CAAC,CAAC;IAC1C,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,sDAAsD,EAAE,GAAG,EAAE;QAChE,MAAM,CAAC,GAAG,CAAC,CAAC;QACZ,MAAM,EAAE,GAAG,CAAC,CAAC;QACb,MAAM,EAAE,GAAG,CAAC,CAAC;QACb,MAAM,EAAE,GAAsB,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;QACxC,MAAM,EAAE,GAAsB,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;QAEzC,MAAM,CAAC,GAAG,wBAAwB,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC,CAAC;QACtD,MAAM,MAAM,GAAG,wBAAwB,CAAC,CAAC,GAAG,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC,CAAC;QAC/D,MAAM,MAAM,GAAG,wBAAwB,CAAC,EAAE,EAAE,CAAC,GAAG,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC,CAAC;QAC/D,MAAM,KAAK,GAAG,wBAAwB,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC;QAE9D,wDAAwD;QACxD,WAAW,CAAC,MAAM,EAAE,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;QACjC,WAAW,CAAC,MAAM,EAAE,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;QACjC,WAAW,CAAC,KAAK,EAAE,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;IAClC,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,oDAAoD,EAAE,GAAG,EAAE;QAC9D,MAAM,MAAM,GAAG,QAAQ,CAAC,CAAC,KAAK;QAC9B,MAAM,IAAI,GAAG,SAAS,CAAC,CAAC,KAAK;QAC7B,MAAM,MAAM,GAAsB,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;QAC5C,MAAM,IAAI,GAAsB,CAAC,cAAc,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,IAAI;QAE5D,MAAM,CAAC,GAAG,wBAAwB,CAAC,MAAM,EAAE,IAAI,EAAE,MAAM,EAAE,IAAI,EAAE,IAAI,CAAC,CAAC;QACrE,MAAM,CAAC,GAAG,IAAI,CAAC,KAAK,CAClB,IAAI,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,EACnB,IAAI,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,EACnB,IAAI,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,CACpB,CAAC;QACF,MAAM,WAAW,GAAG,CAAC,IAAI,GAAG,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;QAErD,QAAQ,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,WAAW,EAAE,KAAK,CAAC,CAAC;QAEtC,qDAAqD;QACrD,MAAM,GAAG,GAAG,QAAQ,CAAC;QACrB,MAAM,CAAC,EAAE,CACP,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,GAAG,GAAG,GAAG,IAAI,EACpC,OAAO,IAAI,CAAC,CAAC,CAAC,kBAAkB,GAAG,EAAE,CACtC,CAAC;IACJ,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,gDAAgD,EAAE,GAAG,EAAE;QAC1D,MAAM,EAAE,GAAsB,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;QACxC,MAAM,EAAE,GAAsB,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;QAExC,+BAA+B;QAC/B,MAAM,CAAC,MAAM,CACX,GAAG,EAAE,CAAC,wBAAwB,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC,EAChD,qBAAqB,CACtB,CAAC;QACF,MAAM,CAAC,MAAM,CACX,GAAG,EAAE,CAAC,wBAAwB,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC,EAChD,qBAAqB,CACtB,CAAC;QACF,MAAM,CAAC,MAAM,CACX,GAAG,EAAE,CAAC,wBAAwB,CAAC,MAAM,CAAC,GAAG,EAAE,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC,EACxD,QAAQ,CACT,CAAC;QAEF,uBAAuB;QACvB,MAAM,CAAC,MAAM,CACX,GAAG,EAAE,CAAC,wBAAwB,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC,EAC/C,oBAAoB,CACrB,CAAC;QAEF,YAAY;QACZ,MAAM,CAAC,MAAM,CAAC,GAAG,EAAE,CAAC,wBAAwB,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC,EAAE,UAAU,CAAC,CAAC;QAC3E,MAAM,CAAC,MAAM,CACX,GAAG,EAAE,CAAC,wBAAwB,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,MAAM,CAAC,GAAG,CAAC,EACxD,QAAQ,CACT,CAAC;IACJ,CAAC,CAAC,CAAC;AACL,CAAC,CAAC,CAAC"}

package/dist/__tests__/gravitational-force.spec.d.ts

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

package/dist/__tests__/gravitational-force.spec.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"gravitational-force.spec.d.ts","sourceRoot":"","sources":["../../src/__tests__/gravitational-force.spec.ts"],"names":[],"mappings":""}

package/dist/__tests__/gravitational-force.spec.js

@@ -0,0 +1,105 @@
+import assert from 'node:assert/strict';
+import test, { describe } from 'node:test';
+import { gravitationalForce } from '../law-of-gravitation';
+import { G_SI } from '@interstellar-tools/constants';
+import { norm, relClose, vecRelClose } from './helpers';
+describe('gravitationalForce', () => {
+    test('computes simple 1D case (+x) with custom G', () => {
+        const G = 1; // simple numbers to check algebra
+        const m1 = 2;
+        const m2 = 3;
+        const r1 = [0, 0, 0];
+        const r2 = [10, 0, 0];
+        const res = gravitationalForce(m1, m2, r1, r2, G);
+        // Expected: |F| = G m1 m2 / r^2 = 6 / 100 = 0.06; direction = +x
+        relClose(res.magnitude, 0.06);
+        vecRelClose(res.vector, [0.06, 0, 0]);
+        // invariants
+        relClose(norm(res.vector), res.magnitude);
+        relClose(norm(res.direction), 1);
+        vecRelClose(res.vector, [
+            res.magnitude * res.direction[0],
+            res.magnitude * res.direction[1],
+            res.magnitude * res.direction[2]
+        ]);
+    });
+    test('direction points from body 1 toward body 2 (arbitrary 3D)', () => {
+        const G = 1.23;
+        const m1 = 1;
+        const m2 = 4;
+        const r1 = [1, 2, 3];
+        const r2 = [-2, 6, 3]; // displacement = (-3, 4, 0)
+        const res = gravitationalForce(m1, m2, r1, r2, G);
+        // direction should align with normalized (r2 - r1)
+        const dr = [
+            r2[0] - r1[0],
+            r2[1] - r1[1],
+            r2[2] - r1[2]
+        ];
+        const r = Math.hypot(...dr);
+        const rhat = [dr[0] / r, dr[1] / r, dr[2] / r];
+        vecRelClose(res.direction, rhat, 1e-12);
+        // vector == magnitude * rhat
+        vecRelClose(res.vector, [
+            res.magnitude * rhat[0],
+            res.magnitude * rhat[1],
+            res.magnitude * rhat[2]
+        ]);
+    });
+    test("Newton's 3rd law: F_1<-2 = - F_2<-1", () => {
+        const G = 9.87;
+        const m1 = 7;
+        const m2 = 11;
+        const r1 = [0, -2, 5];
+        const r2 = [3, 1, -1];
+        const F12 = gravitationalForce(m1, m2, r1, r2, G).vector;
+        const F21 = gravitationalForce(m2, m1, r2, r1, G).vector;
+        vecRelClose(F12, [-F21[0], -F21[1], -F21[2]]);
+    });
+    test('real-world sanity: Earth-Sun at ~1 AU (matches known ~3.54e22 N)', () => {
+        const mEarth = 5.972e24; // kg
+        const mSun = 1.9885e30; // kg
+        const r1 = [0, 0, 0];
+        const r2 = [1.495978707e11, 0, 0]; // m
+        const res = gravitationalForce(mEarth, mSun, r1, r2, G_SI);
+        // Compare to closed-form expectation (not using the function's internal dir)
+        const r = Math.hypot(r2[0] - r1[0], r2[1] - r1[1], r2[2] - r1[2]);
+        const expectedMag = (G_SI * mEarth * mSun) / (r * r);
+        // exact equivalence (up to FP)
+        relClose(res.magnitude, expectedMag, 1e-15);
+        // also compare to the commonly quoted magnitude ~3.542e22 N within 0.5%
+        const ref = 3.542e22;
+        assert.ok(Math.abs(res.magnitude - ref) / ref < 5e-3, `|F| ${res.magnitude} differs from ~${ref}`);
+    });
+    test('error: negative or non-finite masses', () => {
+        const r1 = [0, 0, 0];
+        const r2 = [1, 0, 0];
+        assert.throws(() => gravitationalForce(-1, 1, r1, r2, 1), /non-negative/);
+        assert.throws(() => gravitationalForce(1, -1, r1, r2, 1), /non-negative/);
+        assert.throws(() => gravitationalForce(Number.NaN, 1, r1, r2, 1), /finite/);
+        assert.throws(() => gravitationalForce(1, Number.POSITIVE_INFINITY, r1, r2, 1), /finite/);
+    });
+    test('error: coincident positions (r = 0)', () => {
+        const r = [3, -4, 5];
+        assert.throws(() => gravitationalForce(1, 1, r, r, 1), /singular \(r = 0\)/);
+    });
+    test('error: invalid G (non-finite or <= 0)', () => {
+        const r1 = [0, 0, 0];
+        const r2 = [1, 0, 0];
+        assert.throws(() => gravitationalForce(1, 1, r1, r2, 0), /positive/);
+        assert.throws(() => gravitationalForce(1, 1, r1, r2, Number.NaN), /finite/);
+        assert.throws(() => gravitationalForce(1, 1, r1, r2, -1), /positive/);
+    });
+    test('scales with 1/r^2', () => {
+        const G = 10;
+        const m1 = 2;
+        const m2 = 5;
+        const r1 = [0, 0, 0];
+        const r2a = [2, 0, 0]; // r = 2
+        const r2b = [4, 0, 0]; // r = 4 (double -> quarter force)
+        const Fa = gravitationalForce(m1, m2, r1, r2a, G).magnitude;
+        const Fb = gravitationalForce(m1, m2, r1, r2b, G).magnitude;
+        relClose(Fb, Fa / 4);
+    });
+});
+//# sourceMappingURL=gravitational-force.spec.js.map

package/dist/__tests__/gravitational-force.spec.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"gravitational-force.spec.js","sourceRoot":"","sources":["../../src/__tests__/gravitational-force.spec.ts"],"names":[],"mappings":"AAAA,OAAO,MAAM,MAAM,oBAAoB,CAAC;AACxC,OAAO,IAAI,EAAE,EAAE,QAAQ,EAAE,MAAM,WAAW,CAAC;AAC3C,OAAO,EAAE,kBAAkB,EAAE,MAAM,uBAAuB,CAAC;AAC3D,OAAO,EAAE,IAAI,EAAE,MAAM,+BAA+B,CAAC;AACrD,OAAO,EAAE,IAAI,EAAE,QAAQ,EAAE,WAAW,EAAE,MAAM,WAAW,CAAC;AAExD,QAAQ,CAAC,oBAAoB,EAAE,GAAG,EAAE;IAClC,IAAI,CAAC,4CAA4C,EAAE,GAAG,EAAE;QACtD,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,kCAAkC;QAC/C,MAAM,EAAE,GAAG,CAAC,CAAC;QACb,MAAM,EAAE,GAAG,CAAC,CAAC;QACb,MAAM,EAAE,GAA6B,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;QAC/C,MAAM,EAAE,GAA6B,CAAC,EAAE,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;QAEhD,MAAM,GAAG,GAAG,kBAAkB,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC,CAAC;QAElD,iEAAiE;QACjE,QAAQ,CAAC,GAAG,CAAC,SAAS,EAAE,IAAI,CAAC,CAAC;QAC9B,WAAW,CAAC,GAAG,CAAC,MAAM,EAAE,CAAC,IAAI,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;QACtC,aAAa;QACb,QAAQ,CAAC,IAAI,CAAC,GAAG,CAAC,MAAM,CAAC,EAAE,GAAG,CAAC,SAAS,CAAC,CAAC;QAC1C,QAAQ,CAAC,IAAI,CAAC,GAAG,CAAC,SAAS,CAAC,EAAE,CAAC,CAAC,CAAC;QACjC,WAAW,CAAC,GAAG,CAAC,MAAM,EAAE;YACtB,GAAG,CAAC,SAAS,GAAG,GAAG,CAAC,SAAS,CAAC,CAAC,CAAC;YAChC,GAAG,CAAC,SAAS,GAAG,GAAG,CAAC,SAAS,CAAC,CAAC,CAAC;YAChC,GAAG,CAAC,SAAS,GAAG,GAAG,CAAC,SAAS,CAAC,CAAC,CAAC;SACjC,CAAC,CAAC;IACL,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,2DAA2D,EAAE,GAAG,EAAE;QACrE,MAAM,CAAC,GAAG,IAAI,CAAC;QACf,MAAM,EAAE,GAAG,CAAC,CAAC;QACb,MAAM,EAAE,GAAG,CAAC,CAAC;QACb,MAAM,EAAE,GAA6B,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;QAC/C,MAAM,EAAE,GAA6B,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,4BAA4B;QAC7E,MAAM,GAAG,GAAG,kBAAkB,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC,CAAC;QAElD,mDAAmD;QACnD,MAAM,EAAE,GAA6B;YACnC,EAAE,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC;YACb,EAAE,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC;YACb,EAAE,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC;SACd,CAAC;QACF,MAAM,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,GAAG,EAAE,CAAC,CAAC;QAC5B,MAAM,IAAI,GAA6B,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;QAEzE,WAAW,CAAC,GAAG,CAAC,SAAS,EAAE,IAAI,EAAE,KAAK,CAAC,CAAC;QACxC,6BAA6B;QAC7B,WAAW,CAAC,GAAG,CAAC,MAAM,EAAE;YACtB,GAAG,CAAC,SAAS,GAAG,IAAI,CAAC,CAAC,CAAC;YACvB,GAAG,CAAC,SAAS,GAAG,IAAI,CAAC,CAAC,CAAC;YACvB,GAAG,CAAC,SAAS,GAAG,IAAI,CAAC,CAAC,CAAC;SACxB,CAAC,CAAC;IACL,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,qCAAqC,EAAE,GAAG,EAAE;QAC/C,MAAM,CAAC,GAAG,IAAI,CAAC;QACf,MAAM,EAAE,GAAG,CAAC,CAAC;QACb,MAAM,EAAE,GAAG,EAAE,CAAC;QACd,MAAM,EAAE,GAA6B,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;QAChD,MAAM,EAAE,GAA6B,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;QAEhD,MAAM,GAAG,GAAG,kBAAkB,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC,CAAC,MAAM,CAAC;QACzD,MAAM,GAAG,GAAG,kBAAkB,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC,CAAC,MAAM,CAAC;QAEzD,WAAW,CAAC,GAAG,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IAChD,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,kEAAkE,EAAE,GAAG,EAAE;QAC5E,MAAM,MAAM,GAAG,QAAQ,CAAC,CAAC,KAAK;QAC9B,MAAM,IAAI,GAAG,SAAS,CAAC,CAAC,KAAK;QAC7B,MAAM,EAAE,GAA6B,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;QAC/C,MAAM,EAAE,GAA6B,CAAC,cAAc,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,IAAI;QAEjE,MAAM,GAAG,GAAG,kBAAkB,CAAC,MAAM,EAAE,IAAI,EAAE,EAAE,EAAE,EAAE,EAAE,IAAI,CAAC,CAAC;QAE3D,6EAA6E;QAC7E,MAAM,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;QAClE,MAAM,WAAW,GAAG,CAAC,IAAI,GAAG,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;QAErD,+BAA+B;QAC/B,QAAQ,CAAC,GAAG,CAAC,SAAS,EAAE,WAAW,EAAE,KAAK,CAAC,CAAC;QAE5C,wEAAwE;QACxE,MAAM,GAAG,GAAG,QAAQ,CAAC;QACrB,MAAM,CAAC,EAAE,CACP,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,SAAS,GAAG,GAAG,CAAC,GAAG,GAAG,GAAG,IAAI,EAC1C,OAAO,GAAG,CAAC,SAAS,kBAAkB,GAAG,EAAE,CAC5C,CAAC;IACJ,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,sCAAsC,EAAE,GAAG,EAAE;QAChD,MAAM,EAAE,GAA6B,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;QAC/C,MAAM,EAAE,GAA6B,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;QAE/C,MAAM,CAAC,MAAM,CAAC,GAAG,EAAE,CAAC,kBAAkB,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC,EAAE,cAAc,CAAC,CAAC;QAC1E,MAAM,CAAC,MAAM,CAAC,GAAG,EAAE,CAAC,kBAAkB,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC,EAAE,cAAc,CAAC,CAAC;QAC1E,MAAM,CAAC,MAAM,CAAC,GAAG,EAAE,CAAC,kBAAkB,CAAC,MAAM,CAAC,GAAG,EAAE,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC,EAAE,QAAQ,CAAC,CAAC;QAC5E,MAAM,CAAC,MAAM,CACX,GAAG,EAAE,CAAC,kBAAkB,CAAC,CAAC,EAAE,MAAM,CAAC,iBAAiB,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC,EAChE,QAAQ,CACT,CAAC;IACJ,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,qCAAqC,EAAE,GAAG,EAAE;QAC/C,MAAM,CAAC,GAA6B,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;QAC/C,MAAM,CAAC,MAAM,CACX,GAAG,EAAE,CAAC,kBAAkB,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,EACvC,oBAAoB,CACrB,CAAC;IACJ,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,uCAAuC,EAAE,GAAG,EAAE;QACjD,MAAM,EAAE,GAA6B,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;QAC/C,MAAM,EAAE,GAA6B,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;QAE/C,MAAM,CAAC,MAAM,CAAC,GAAG,EAAE,CAAC,kBAAkB,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC,EAAE,UAAU,CAAC,CAAC;QACrE,MAAM,CAAC,MAAM,CAAC,GAAG,EAAE,CAAC,kBAAkB,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,MAAM,CAAC,GAAG,CAAC,EAAE,QAAQ,CAAC,CAAC;QAC5E,MAAM,CAAC,MAAM,CAAC,GAAG,EAAE,CAAC,kBAAkB,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC,CAAC,EAAE,UAAU,CAAC,CAAC;IACxE,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,mBAAmB,EAAE,GAAG,EAAE;QAC7B,MAAM,CAAC,GAAG,EAAE,CAAC;QACb,MAAM,EAAE,GAAG,CAAC,CAAC;QACb,MAAM,EAAE,GAAG,CAAC,CAAC;QACb,MAAM,EAAE,GAA6B,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;QAE/C,MAAM,GAAG,GAA6B,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,QAAQ;QACzD,MAAM,GAAG,GAA6B,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,kCAAkC;QAEnF,MAAM,EAAE,GAAG,kBAAkB,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,GAAG,EAAE,CAAC,CAAC,CAAC,SAAS,CAAC;QAC5D,MAAM,EAAE,GAAG,kBAAkB,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,GAAG,EAAE,CAAC,CAAC,CAAC,SAAS,CAAC;QAE5D,QAAQ,CAAC,EAAE,EAAE,EAAE,GAAG,CAAC,CAAC,CAAC;IACvB,CAAC,CAAC,CAAC;AACL,CAAC,CAAC,CAAC"}

package/dist/__tests__/gravitational-parameters.spec.d.ts

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

package/dist/__tests__/gravitational-parameters.spec.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"gravitational-parameters.spec.d.ts","sourceRoot":"","sources":["../../src/__tests__/gravitational-parameters.spec.ts"],"names":[],"mappings":""}

package/dist/__tests__/gravitational-parameters.spec.js

@@ -0,0 +1,75 @@
+import assert from 'node:assert/strict';
+import test, { describe } from 'node:test';
+import { relClose } from './helpers';
+import { gravitationalParameter } from '../gravitational-parameter';
+describe('gravitationalParameter', () => {
+    test('computes μ = G(M + m) for simple numbers', () => {
+        const G = 10;
+        const M = 2;
+        const m = 3;
+        const mu = gravitationalParameter(M, m, G);
+        assert.equal(mu, 10 * (2 + 3)); // 50
+    });
+    test('matches canonical GM_Earth when M = GM/G', () => {
+        const GM_canonical = 3.986004418e14; // m^3/s^2
+        const M_from_GM = GM_canonical / 6.6743e-11;
+        const mu = gravitationalParameter(M_from_GM); // m=0, G=G_SI
+        relClose(mu, GM_canonical, 1e-12);
+    });
+    test('M = 0 returns μ = G * m', () => {
+        const G = 6.7;
+        const M = 0;
+        const m = 123;
+        assert.equal(gravitationalParameter(M, m, G), G * m);
+    });
+    test('commutative in masses: μ(M,m) = μ(m,M)', () => {
+        const G = 3.21;
+        const M = 7;
+        const m = 11;
+        const mu1 = gravitationalParameter(M, m, G);
+        const mu2 = gravitationalParameter(m, M, G);
+        assert.equal(mu1, mu2);
+    });
+    test('approximation: when m ≪ M, μ ≈ G M within ~m/M relative error', () => {
+        const G = 2;
+        const M = 1e8;
+        const m = 1e2; // m/M = 1e-6
+        const exact = gravitationalParameter(M, m, G);
+        const approx = G * M;
+        // Relative error should be about m/M
+        const relErr = Math.abs(exact - approx) / approx;
+        assert.ok(relErr < 2e-6, `relative error ${relErr} not < 2e-6`);
+    });
+    test('scales linearly with G and with (M + m)', () => {
+        const G = 5;
+        const M = 10;
+        const m = 4;
+        const mu = gravitationalParameter(M, m, G);
+        const mu_doubleG = gravitationalParameter(M, m, 2 * G);
+        const mu_doubleMasses = gravitationalParameter(2 * M, 2 * m, G);
+        // Doubling G doubles μ
+        relClose(mu_doubleG, 2 * mu, 1e-15, 'doubling G should double μ');
+        // Doubling both masses doubles (M + m) and thus doubles μ
+        relClose(mu_doubleMasses, 2 * mu, 1e-15, 'doubling masses should double μ');
+    });
+    test('throws on invalid M (negative or non-finite)', () => {
+        const G = 6.7;
+        const m = 1;
+        assert.throws(() => gravitationalParameter(-1, m, G), /finite|positive/);
+        assert.throws(() => gravitationalParameter(Number.NaN, m, G), /finite/);
+    });
+    test('throws on invalid m (negative or non-finite)', () => {
+        const G = 6.7;
+        const M = 1;
+        assert.throws(() => gravitationalParameter(M, -1, G), /finite|positive/);
+        assert.throws(() => gravitationalParameter(M, Number.NaN, G), /finite/);
+    });
+    test('throws on invalid G (non-finite or ≤ 0)', () => {
+        const M = 1;
+        const m = 1;
+        assert.throws(() => gravitationalParameter(M, m, 0), /positive/);
+        assert.throws(() => gravitationalParameter(M, m, -1), /positive/);
+        assert.throws(() => gravitationalParameter(M, m, Number.NaN), /finite/);
+    });
+});
+//# sourceMappingURL=gravitational-parameters.spec.js.map

package/dist/__tests__/gravitational-parameters.spec.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"gravitational-parameters.spec.js","sourceRoot":"","sources":["../../src/__tests__/gravitational-parameters.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,MAAM,WAAW,CAAC;AACrC,OAAO,EAAE,sBAAsB,EAAE,MAAM,4BAA4B,CAAC;AAEpE,QAAQ,CAAC,wBAAwB,EAAE,GAAG,EAAE;IACtC,IAAI,CAAC,0CAA0C,EAAE,GAAG,EAAE;QACpD,MAAM,CAAC,GAAG,EAAE,CAAC;QACb,MAAM,CAAC,GAAG,CAAC,CAAC;QACZ,MAAM,CAAC,GAAG,CAAC,CAAC;QACZ,MAAM,EAAE,GAAG,sBAAsB,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;QAC3C,MAAM,CAAC,KAAK,CAAC,EAAE,EAAE,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,KAAK;IACvC,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,0CAA0C,EAAE,GAAG,EAAE;QACpD,MAAM,YAAY,GAAG,cAAc,CAAC,CAAC,UAAU;QAC/C,MAAM,SAAS,GAAG,YAAY,GAAG,UAAU,CAAC;QAC5C,MAAM,EAAE,GAAG,sBAAsB,CAAC,SAAS,CAAC,CAAC,CAAC,cAAc;QAC5D,QAAQ,CAAC,EAAE,EAAE,YAAY,EAAE,KAAK,CAAC,CAAC;IACpC,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,yBAAyB,EAAE,GAAG,EAAE;QACnC,MAAM,CAAC,GAAG,GAAG,CAAC;QACd,MAAM,CAAC,GAAG,CAAC,CAAC;QACZ,MAAM,CAAC,GAAG,GAAG,CAAC;QACd,MAAM,CAAC,KAAK,CAAC,sBAAsB,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC;IACvD,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,wCAAwC,EAAE,GAAG,EAAE;QAClD,MAAM,CAAC,GAAG,IAAI,CAAC;QACf,MAAM,CAAC,GAAG,CAAC,CAAC;QACZ,MAAM,CAAC,GAAG,EAAE,CAAC;QACb,MAAM,GAAG,GAAG,sBAAsB,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;QAC5C,MAAM,GAAG,GAAG,sBAAsB,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;QAC5C,MAAM,CAAC,KAAK,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC;IACzB,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,+DAA+D,EAAE,GAAG,EAAE;QACzE,MAAM,CAAC,GAAG,CAAC,CAAC;QACZ,MAAM,CAAC,GAAG,GAAG,CAAC;QACd,MAAM,CAAC,GAAG,GAAG,CAAC,CAAC,aAAa;QAC5B,MAAM,KAAK,GAAG,sBAAsB,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;QAC9C,MAAM,MAAM,GAAG,CAAC,GAAG,CAAC,CAAC;QACrB,qCAAqC;QACrC,MAAM,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,GAAG,MAAM,CAAC,GAAG,MAAM,CAAC;QACjD,MAAM,CAAC,EAAE,CAAC,MAAM,GAAG,IAAI,EAAE,kBAAkB,MAAM,aAAa,CAAC,CAAC;IAClE,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,yCAAyC,EAAE,GAAG,EAAE;QACnD,MAAM,CAAC,GAAG,CAAC,CAAC;QACZ,MAAM,CAAC,GAAG,EAAE,CAAC;QACb,MAAM,CAAC,GAAG,CAAC,CAAC;QAEZ,MAAM,EAAE,GAAG,sBAAsB,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;QAC3C,MAAM,UAAU,GAAG,sBAAsB,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC;QACvD,MAAM,eAAe,GAAG,sBAAsB,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC;QAEhE,uBAAuB;QACvB,QAAQ,CAAC,UAAU,EAAE,CAAC,GAAG,EAAE,EAAE,KAAK,EAAE,4BAA4B,CAAC,CAAC;QAClE,0DAA0D;QAC1D,QAAQ,CAAC,eAAe,EAAE,CAAC,GAAG,EAAE,EAAE,KAAK,EAAE,iCAAiC,CAAC,CAAC;IAC9E,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,8CAA8C,EAAE,GAAG,EAAE;QACxD,MAAM,CAAC,GAAG,GAAG,CAAC;QACd,MAAM,CAAC,GAAG,CAAC,CAAC;QACZ,MAAM,CAAC,MAAM,CAAC,GAAG,EAAE,CAAC,sBAAsB,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,EAAE,iBAAiB,CAAC,CAAC;QACzE,MAAM,CAAC,MAAM,CAAC,GAAG,EAAE,CAAC,sBAAsB,CAAC,MAAM,CAAC,GAAG,EAAE,CAAC,EAAE,CAAC,CAAC,EAAE,QAAQ,CAAC,CAAC;IAC1E,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,8CAA8C,EAAE,GAAG,EAAE;QACxD,MAAM,CAAC,GAAG,GAAG,CAAC;QACd,MAAM,CAAC,GAAG,CAAC,CAAC;QACZ,MAAM,CAAC,MAAM,CAAC,GAAG,EAAE,CAAC,sBAAsB,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,iBAAiB,CAAC,CAAC;QACzE,MAAM,CAAC,MAAM,CAAC,GAAG,EAAE,CAAC,sBAAsB,CAAC,CAAC,EAAE,MAAM,CAAC,GAAG,EAAE,CAAC,CAAC,EAAE,QAAQ,CAAC,CAAC;IAC1E,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,yCAAyC,EAAE,GAAG,EAAE;QACnD,MAAM,CAAC,GAAG,CAAC,CAAC;QACZ,MAAM,CAAC,GAAG,CAAC,CAAC;QACZ,MAAM,CAAC,MAAM,CAAC,GAAG,EAAE,CAAC,sBAAsB,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,EAAE,UAAU,CAAC,CAAC;QACjE,MAAM,CAAC,MAAM,CAAC,GAAG,EAAE,CAAC,sBAAsB,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,UAAU,CAAC,CAAC;QAClE,MAAM,CAAC,MAAM,CAAC,GAAG,EAAE,CAAC,sBAAsB,CAAC,CAAC,EAAE,CAAC,EAAE,MAAM,CAAC,GAAG,CAAC,EAAE,QAAQ,CAAC,CAAC;IAC1E,CAAC,CAAC,CAAC;AACL,CAAC,CAAC,CAAC"}

package/dist/__tests__/helpers/index.d.ts

@@ -0,0 +1,8 @@
+import { Vector3DTupleType } from '@interstellar-tools/types';
+export declare const relClose: (a: number, b: number, eps?: number, msg?: string) => void;
+export declare const vecRelClose: (a: Vector3DTupleType, b: Vector3DTupleType, eps?: number, msg?: string) => void;
+export declare const norm: (v: Vector3DTupleType) => number;
+export declare const scale: (v: Vector3DTupleType, s: number) => Vector3DTupleType;
+export declare const sub: (a: Vector3DTupleType, b: Vector3DTupleType) => Vector3DTupleType;
+export declare const absClose: (a: number, b: number, eps?: number, msg?: string) => void;
+//# sourceMappingURL=index.d.ts.map

package/dist/__tests__/helpers/index.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"index.d.ts","sourceRoot":"","sources":["../../../src/__tests__/helpers/index.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,iBAAiB,EAAE,MAAM,2BAA2B,CAAC;AAG9D,eAAO,MAAM,QAAQ,GAAI,GAAG,MAAM,EAAE,GAAG,MAAM,EAAE,YAAW,EAAE,MAAM,MAAM,SAMvE,CAAC;AAEF,eAAO,MAAM,WAAW,GACtB,GAAG,iBAAiB,EACpB,GAAG,iBAAiB,EACpB,YAAW,EACX,MAAM,MAAM,SAOb,CAAC;AAEF,eAAO,MAAM,IAAI,GAAI,GAAG,iBAAiB,WAExC,CAAC;AAEF,eAAO,MAAM,KAAK,GAAI,GAAG,iBAAiB,EAAE,GAAG,MAAM,KAAG,iBAEvD,CAAC;AAEF,eAAO,MAAM,GAAG,GACd,GAAG,iBAAiB,EACpB,GAAG,iBAAiB,KACnB,iBAEF,CAAC;AAEF,eAAO,MAAM,QAAQ,GAAI,GAAG,MAAM,EAAE,GAAG,MAAM,EAAE,YAAW,EAAE,MAAM,MAAM,SAEvE,CAAC"}

package/dist/__tests__/helpers/index.js

@@ -0,0 +1,25 @@
+import assert from 'node:assert/strict';
+export const relClose = (a, b, eps = 1e-12, msg) => {
+    const scale = Math.max(1, Math.abs(a), Math.abs(b));
+    assert.ok(Math.abs(a - b) <= eps * scale, msg ?? `|${a} - ${b}| > ${eps} * ${scale}`);
+};
+export const vecRelClose = (a, b, eps = 1e-12, msg) => {
+    assert.equal(a.length, 3);
+    assert.equal(b.length, 3);
+    for (let i = 0; i < 3; i++) {
+        relClose(a[i], b[i], eps, msg);
+    }
+};
+export const norm = (v) => {
+    return Math.hypot(v[0], v[1], v[2]);
+};
+export const scale = (v, s) => {
+    return [v[0] * s, v[1] * s, v[2] * s];
+};
+export const sub = (a, b) => {
+    return [a[0] - b[0], a[1] - b[1], a[2] - b[2]];
+};
+export const absClose = (a, b, eps = 1e-12, msg) => {
+    assert.ok(Math.abs(a - b) <= eps, msg ?? `|${a} - ${b}| > ${eps}`);
+};
+//# sourceMappingURL=index.js.map

package/dist/__tests__/helpers/index.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"index.js","sourceRoot":"","sources":["../../../src/__tests__/helpers/index.ts"],"names":[],"mappings":"AACA,OAAO,MAAM,MAAM,oBAAoB,CAAC;AAExC,MAAM,CAAC,MAAM,QAAQ,GAAG,CAAC,CAAS,EAAE,CAAS,EAAE,GAAG,GAAG,KAAK,EAAE,GAAY,EAAE,EAAE;IAC1E,MAAM,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;IACpD,MAAM,CAAC,EAAE,CACP,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,IAAI,GAAG,GAAG,KAAK,EAC9B,GAAG,IAAI,IAAI,CAAC,MAAM,CAAC,OAAO,GAAG,MAAM,KAAK,EAAE,CAC3C,CAAC;AACJ,CAAC,CAAC;AAEF,MAAM,CAAC,MAAM,WAAW,GAAG,CACzB,CAAoB,EACpB,CAAoB,EACpB,GAAG,GAAG,KAAK,EACX,GAAY,EACZ,EAAE;IACF,MAAM,CAAC,KAAK,CAAC,CAAC,CAAC,MAAM,EAAE,CAAC,CAAC,CAAC;IAC1B,MAAM,CAAC,KAAK,CAAC,CAAC,CAAC,MAAM,EAAE,CAAC,CAAC,CAAC;IAC1B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;QAC3B,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,GAAG,EAAE,GAAG,CAAC,CAAC;IACjC,CAAC;AACH,CAAC,CAAC;AAEF,MAAM,CAAC,MAAM,IAAI,GAAG,CAAC,CAAoB,EAAE,EAAE;IAC3C,OAAO,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;AACtC,CAAC,CAAC;AAEF,MAAM,CAAC,MAAM,KAAK,GAAG,CAAC,CAAoB,EAAE,CAAS,EAAqB,EAAE;IAC1E,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;AACxC,CAAC,CAAC;AAEF,MAAM,CAAC,MAAM,GAAG,GAAG,CACjB,CAAoB,EACpB,CAAoB,EACD,EAAE;IACrB,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;AACjD,CAAC,CAAC;AAEF,MAAM,CAAC,MAAM,QAAQ,GAAG,CAAC,CAAS,EAAE,CAAS,EAAE,GAAG,GAAG,KAAK,EAAE,GAAY,EAAE,EAAE;IAC1E,MAAM,CAAC,EAAE,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,IAAI,GAAG,EAAE,GAAG,IAAI,IAAI,CAAC,MAAM,CAAC,OAAO,GAAG,EAAE,CAAC,CAAC;AACrE,CAAC,CAAC"}

package/dist/__tests__/kepler-period.spec.d.ts

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

package/dist/__tests__/kepler-period.spec.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"kepler-period.spec.d.ts","sourceRoot":"","sources":["../../src/__tests__/kepler-period.spec.ts"],"names":[],"mappings":""}

package/dist/__tests__/kepler-period.spec.js

@@ -0,0 +1,72 @@
+import assert from 'node:assert/strict';
+import test, { describe } from 'node:test';
+import { relClose } from './helpers';
+import { keplerPeriod } from '../kepler-period';
+describe('keplerPeriod', () => {
+    test('GEO sanity: a ≈ 42,164 km around Earth → T ≈ 86,164 s', () => {
+        const a = 42164e3; // m
+        const muEarth = 3.986004418e14; // m^3/s^2
+        const T = keplerPeriod(a, muEarth);
+        const ref = 86164.0905; // sidereal day (s)
+        // Keep tolerance modest: orbit isn't exactly Keplerian in reality; our formula is ideal
+        relClose(T, ref, 1e-3); // 0.1% relative tolerance
+    });
+    test('LEO circular (~400 km): a = Re + 400 km → T ~ 5550 s', () => {
+        const a = 6378e3 + 400e3; // m
+        const muEarth = 3.986004418e14; // m^3/s^2
+        const T = keplerPeriod(a, muEarth);
+        const ref = 5550; // ~92.5 min
+        relClose(T, ref, 5e-2); // 5% slack; depends on precise constants
+    });
+    test('Scaling with semi-major axis: T ∝ a^(3/2) (same μ)', () => {
+        const mu = 3.986004418e14;
+        const a1 = 7000e3;
+        const a2 = 2 * a1;
+        const T1 = keplerPeriod(a1, mu);
+        const T2 = keplerPeriod(a2, mu);
+        const expectedRatio = Math.pow(2, 1.5); // 2^(3/2)
+        relClose(T2 / T1, expectedRatio, 1e-12);
+    });
+    test('Scaling with μ: doubling μ reduces T by 1/√2 (same a)', () => {
+        const a = 10000e3;
+        const mu = 3.986004418e14;
+        const T1 = keplerPeriod(a, mu);
+        const T2 = keplerPeriod(a, 2 * mu);
+        const expectedRatio = 1 / Math.SQRT2;
+        relClose(T2 / T1, expectedRatio, 1e-12);
+    });
+    test('Heliocentric: a = 1 AU around Sun → T ≈ 365.256 days', () => {
+        const AU = 149597870700; // m (IAU 2012)
+        const muSun = 1.32712440018e20; // m^3/s^2
+        const T = keplerPeriod(AU, muSun);
+        const siderealYear = 365.256363004 * 86400; // s
+        relClose(T, siderealYear, 1e-3); // 0.1% tolerance
+    });
+    test('Input validation: a <= 0, μ <= 0, or non-finite throws', () => {
+        const a = 7000e3;
+        const mu = 3.986004418e14;
+        // a invalid
+        assert.throws(() => keplerPeriod(0, mu), /positive/);
+        assert.throws(() => keplerPeriod(-1, mu), /positive/);
+        assert.throws(() => keplerPeriod(Number.NaN, mu), /finite/);
+        // μ invalid
+        assert.throws(() => keplerPeriod(a, 0), /positive/);
+        assert.throws(() => keplerPeriod(a, -1), /positive/);
+        assert.throws(() => keplerPeriod(a, Number.NaN), /finite/);
+    });
+    test('Numerical stability: very large a still yields finite result', () => {
+        const a = 1e9; // 1,000,000 km
+        const muEarth = 3.986004418e14; // m^3/s^2
+        const T = keplerPeriod(a, muEarth);
+        assert.ok(Number.isFinite(T) && T > 0);
+    });
+    // Optional tight identity test against the formula directly
+    test('Definition check: T = 2π * sqrt(a^3/μ)', () => {
+        const a = 12345678; // m
+        const mu = 7.654321e13; // m^3/s^2
+        const T = keplerPeriod(a, mu);
+        const expected = 2 * Math.PI * Math.sqrt((a * a * a) / mu);
+        relClose(T, expected, 1e-15);
+    });
+});
+//# sourceMappingURL=kepler-period.spec.js.map

package/dist/__tests__/kepler-period.spec.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"kepler-period.spec.js","sourceRoot":"","sources":["../../src/__tests__/kepler-period.spec.ts"],"names":[],"mappings":"AAAA,OAAO,MAAM,MAAM,oBAAoB,CAAC;AACxC,OAAO,IAAI,EAAE,EAAE,QAAQ,EAAE,MAAM,WAAW,CAAC;AAE3C,OAAO,EAAE,QAAQ,EAAE,MAAM,WAAW,CAAC;AACrC,OAAO,EAAE,YAAY,EAAE,MAAM,kBAAkB,CAAC;AAEhD,QAAQ,CAAC,cAAc,EAAE,GAAG,EAAE;IAC5B,IAAI,CAAC,uDAAuD,EAAE,GAAG,EAAE;QACjE,MAAM,CAAC,GAAG,OAAO,CAAC,CAAC,IAAI;QACvB,MAAM,OAAO,GAAG,cAAc,CAAC,CAAC,UAAU;QAC1C,MAAM,CAAC,GAAG,YAAY,CAAC,CAAC,EAAE,OAAO,CAAC,CAAC;QAEnC,MAAM,GAAG,GAAG,UAAU,CAAC,CAAC,mBAAmB;QAC3C,wFAAwF;QACxF,QAAQ,CAAC,CAAC,EAAE,GAAG,EAAE,IAAI,CAAC,CAAC,CAAC,0BAA0B;IACpD,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,sDAAsD,EAAE,GAAG,EAAE;QAChE,MAAM,CAAC,GAAG,MAAM,GAAG,KAAK,CAAC,CAAC,IAAI;QAC9B,MAAM,OAAO,GAAG,cAAc,CAAC,CAAC,UAAU;QAC1C,MAAM,CAAC,GAAG,YAAY,CAAC,CAAC,EAAE,OAAO,CAAC,CAAC;QAEnC,MAAM,GAAG,GAAG,IAAI,CAAC,CAAC,YAAY;QAC9B,QAAQ,CAAC,CAAC,EAAE,GAAG,EAAE,IAAI,CAAC,CAAC,CAAC,yCAAyC;IACnE,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,oDAAoD,EAAE,GAAG,EAAE;QAC9D,MAAM,EAAE,GAAG,cAAc,CAAC;QAC1B,MAAM,EAAE,GAAG,MAAM,CAAC;QAClB,MAAM,EAAE,GAAG,CAAC,GAAG,EAAE,CAAC;QAElB,MAAM,EAAE,GAAG,YAAY,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC;QAChC,MAAM,EAAE,GAAG,YAAY,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC;QAEhC,MAAM,aAAa,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,UAAU;QAClD,QAAQ,CAAC,EAAE,GAAG,EAAE,EAAE,aAAa,EAAE,KAAK,CAAC,CAAC;IAC1C,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,uDAAuD,EAAE,GAAG,EAAE;QACjE,MAAM,CAAC,GAAG,OAAO,CAAC;QAClB,MAAM,EAAE,GAAG,cAAc,CAAC;QAE1B,MAAM,EAAE,GAAG,YAAY,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC;QAC/B,MAAM,EAAE,GAAG,YAAY,CAAC,CAAC,EAAE,CAAC,GAAG,EAAE,CAAC,CAAC;QAEnC,MAAM,aAAa,GAAG,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC;QACrC,QAAQ,CAAC,EAAE,GAAG,EAAE,EAAE,aAAa,EAAE,KAAK,CAAC,CAAC;IAC1C,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,sDAAsD,EAAE,GAAG,EAAE;QAChE,MAAM,EAAE,GAAG,YAAe,CAAC,CAAC,eAAe;QAC3C,MAAM,KAAK,GAAG,gBAAgB,CAAC,CAAC,UAAU;QAC1C,MAAM,CAAC,GAAG,YAAY,CAAC,EAAE,EAAE,KAAK,CAAC,CAAC;QAElC,MAAM,YAAY,GAAG,aAAa,GAAG,KAAK,CAAC,CAAC,IAAI;QAChD,QAAQ,CAAC,CAAC,EAAE,YAAY,EAAE,IAAI,CAAC,CAAC,CAAC,iBAAiB;IACpD,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,wDAAwD,EAAE,GAAG,EAAE;QAClE,MAAM,CAAC,GAAG,MAAM,CAAC;QACjB,MAAM,EAAE,GAAG,cAAc,CAAC;QAE1B,YAAY;QACZ,MAAM,CAAC,MAAM,CAAC,GAAG,EAAE,CAAC,YAAY,CAAC,CAAC,EAAE,EAAE,CAAC,EAAE,UAAU,CAAC,CAAC;QACrD,MAAM,CAAC,MAAM,CAAC,GAAG,EAAE,CAAC,YAAY,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,EAAE,UAAU,CAAC,CAAC;QACtD,MAAM,CAAC,MAAM,CAAC,GAAG,EAAE,CAAC,YAAY,CAAC,MAAM,CAAC,GAAG,EAAE,EAAE,CAAC,EAAE,QAAQ,CAAC,CAAC;QAE5D,YAAY;QACZ,MAAM,CAAC,MAAM,CAAC,GAAG,EAAE,CAAC,YAAY,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,UAAU,CAAC,CAAC;QACpD,MAAM,CAAC,MAAM,CAAC,GAAG,EAAE,CAAC,YAAY,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,UAAU,CAAC,CAAC;QACrD,MAAM,CAAC,MAAM,CAAC,GAAG,EAAE,CAAC,YAAY,CAAC,CAAC,EAAE,MAAM,CAAC,GAAG,CAAC,EAAE,QAAQ,CAAC,CAAC;IAC7D,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,8DAA8D,EAAE,GAAG,EAAE;QACxE,MAAM,CAAC,GAAG,GAAG,CAAC,CAAC,eAAe;QAC9B,MAAM,OAAO,GAAG,cAAc,CAAC,CAAC,UAAU;QAC1C,MAAM,CAAC,GAAG,YAAY,CAAC,CAAC,EAAE,OAAO,CAAC,CAAC;QAEnC,MAAM,CAAC,EAAE,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC;IACzC,CAAC,CAAC,CAAC;IAEH,4DAA4D;IAC5D,IAAI,CAAC,wCAAwC,EAAE,GAAG,EAAE;QAClD,MAAM,CAAC,GAAG,QAAU,CAAC,CAAC,IAAI;QAC1B,MAAM,EAAE,GAAG,WAAW,CAAC,CAAC,UAAU;QAElC,MAAM,CAAC,GAAG,YAAY,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC;QAC9B,MAAM,QAAQ,GAAG,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC;QAE3D,QAAQ,CAAC,CAAC,EAAE,QAAQ,EAAE,KAAK,CAAC,CAAC;IAC/B,CAAC,CAAC,CAAC;AACL,CAAC,CAAC,CAAC"}

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

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

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

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