@interstellar-tools/equations 0.1.1

package/README.md

@@ -0,0 +1,18 @@
+# @interstellar-tools/equations
+
+## Installation
+
+- [Node.js](https://nodejs.org/) version 22.9.0 or higher
+- npm version 11.5.1 or higher
+
+::: code-group
+
+```shell [npm]
+npm i --save @interstellar-tools/equations
+```
+
+```shell [yarn]
+yarn add @interstellar-tools/equations
+```
+
+:::

package/dist/__tests__/compute-angle.int.spec.d.ts

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

package/dist/__tests__/compute-angle.int.spec.d.ts.map

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

package/dist/__tests__/compute-angle.int.spec.js

@@ -0,0 +1,62 @@
+import assert from 'node:assert/strict';
+import test, { describe } from 'node:test';
+import { computeAngle } from '../compute-angle';
+import { wrapAngle } from '../wrap-angle';
+describe('computeAngle', () => {
+    // Test cases
+    const testCases = [
+        {
+            celestialBody: {
+                e: 0.0167,
+                angle: 0,
+                period: { value: 365.25, unit: 'day' }
+            },
+            timeStep: { value: 1, unit: 'day' }
+        },
+        {
+            celestialBody: {
+                e: 0.1,
+                angle: Math.PI / 3,
+                period: { value: 687, unit: 'day' }
+            },
+            timeStep: { value: 1, unit: 'day' }
+        },
+        // Retrograde motion test cases
+        {
+            celestialBody: {
+                e: 0.05,
+                angle: Math.PI / 2,
+                period: { value: -365.25, unit: 'day' } // Retrograde motion
+            },
+            timeStep: { value: 1, unit: 'day' }
+        },
+        {
+            celestialBody: {
+                e: 0.2,
+                angle: -Math.PI / 4,
+                period: { value: -4332.59, unit: 'day' } // Retrograde motion
+            },
+            timeStep: { value: 10, unit: 'day' }
+        }
+    ];
+    test('computes correct true anomaly', () => {
+        testCases.forEach(({ celestialBody, timeStep }) => {
+            const expected = wrapAngle(computeAngle(celestialBody, timeStep));
+            const result = wrapAngle(computeAngle(celestialBody, timeStep));
+            assert.ok(Math.abs(result - expected) < 1e-4, `Expected angle ≈ ${expected} for body=${JSON.stringify(celestialBody)}, timeStep=${timeStep.value}, but got ${result}`);
+        });
+    });
+    test('correctly adjusts for comets', () => {
+        const comet = {
+            e: 0.8,
+            q: 1.2,
+            angle: Math.PI / 4,
+            period: { value: 2000, unit: 'day' }
+        }; // Adding 'q' property to classify as comet
+        const timeStep = { value: 5, unit: 'day' };
+        const expected = wrapAngle(computeAngle(comet, timeStep));
+        const result = wrapAngle(computeAngle(comet, timeStep));
+        assert.ok(Math.abs(result - expected) < 1e-4, `Expected comet angle ≈ ${expected}, but got ${result}`);
+    });
+});
+//# sourceMappingURL=compute-angle.int.spec.js.map

package/dist/__tests__/compute-angle.int.spec.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"compute-angle.int.spec.js","sourceRoot":"","sources":["../../src/__tests__/compute-angle.int.spec.ts"],"names":[],"mappings":"AAAA,OAAO,MAAM,MAAM,oBAAoB,CAAC;AACxC,OAAO,IAAI,EAAE,EAAE,QAAQ,EAAE,MAAM,WAAW,CAAC;AAE3C,OAAO,EAAE,YAAY,EAAE,MAAM,kBAAkB,CAAC;AAChD,OAAO,EAAE,SAAS,EAAE,MAAM,eAAe,CAAC;AAM1C,QAAQ,CAAC,cAAc,EAAE,GAAG,EAAE;IAC5B,aAAa;IACb,MAAM,SAAS,GAAG;QAChB;YACE,aAAa,EAAE;gBACb,CAAC,EAAE,MAAM;gBACT,KAAK,EAAE,CAAC;gBACR,MAAM,EAAE,EAAE,KAAK,EAAE,MAAM,EAAE,IAAI,EAAE,KAAK,EAAE;aACvC;YACD,QAAQ,EAAE,EAAE,KAAK,EAAE,CAAC,EAAE,IAAI,EAAE,KAAK,EAAE;SACpC;QACD;YACE,aAAa,EAAE;gBACb,CAAC,EAAE,GAAG;gBACN,KAAK,EAAE,IAAI,CAAC,EAAE,GAAG,CAAC;gBAClB,MAAM,EAAE,EAAE,KAAK,EAAE,GAAG,EAAE,IAAI,EAAE,KAAK,EAAE;aACpC;YACD,QAAQ,EAAE,EAAE,KAAK,EAAE,CAAC,EAAE,IAAI,EAAE,KAAK,EAAE;SACpC;QACD,+BAA+B;QAC/B;YACE,aAAa,EAAE;gBACb,CAAC,EAAE,IAAI;gBACP,KAAK,EAAE,IAAI,CAAC,EAAE,GAAG,CAAC;gBAClB,MAAM,EAAE,EAAE,KAAK,EAAE,CAAC,MAAM,EAAE,IAAI,EAAE,KAAK,EAAE,CAAC,oBAAoB;aAC7D;YACD,QAAQ,EAAE,EAAE,KAAK,EAAE,CAAC,EAAE,IAAI,EAAE,KAAK,EAAE;SACpC;QACD;YACE,aAAa,EAAE;gBACb,CAAC,EAAE,GAAG;gBACN,KAAK,EAAE,CAAC,IAAI,CAAC,EAAE,GAAG,CAAC;gBACnB,MAAM,EAAE,EAAE,KAAK,EAAE,CAAC,OAAO,EAAE,IAAI,EAAE,KAAK,EAAE,CAAC,oBAAoB;aAC9D;YACD,QAAQ,EAAE,EAAE,KAAK,EAAE,EAAE,EAAE,IAAI,EAAE,KAAK,EAAE;SACrC;KACF,CAAC;IAEF,IAAI,CAAC,+BAA+B,EAAE,GAAG,EAAE;QACzC,SAAS,CAAC,OAAO,CAAC,CAAC,EAAE,aAAa,EAAE,QAAQ,EAAE,EAAE,EAAE;YAChD,MAAM,QAAQ,GAAG,SAAS,CACxB,YAAY,CACV,aAAkC,EAClC,QAA6B,CAC9B,CACF,CAAC;YACF,MAAM,MAAM,GAAG,SAAS,CACtB,YAAY,CACV,aAAkC,EAClC,QAA6B,CAC9B,CACF,CAAC;YAEF,MAAM,CAAC,EAAE,CACP,IAAI,CAAC,GAAG,CAAC,MAAM,GAAG,QAAQ,CAAC,GAAG,IAAI,EAClC,oBAAoB,QAAQ,aAAa,IAAI,CAAC,SAAS,CAAC,aAAa,CAAC,cAAc,QAAQ,CAAC,KAAK,aAAa,MAAM,EAAE,CACxH,CAAC;QACJ,CAAC,CAAC,CAAC;IACL,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,8BAA8B,EAAE,GAAG,EAAE;QACxC,MAAM,KAAK,GAAG;YACZ,CAAC,EAAE,GAAG;YACN,CAAC,EAAE,GAAG;YACN,KAAK,EAAE,IAAI,CAAC,EAAE,GAAG,CAAC;YAClB,MAAM,EAAE,EAAE,KAAK,EAAE,IAAI,EAAE,IAAI,EAAE,KAAK,EAAE;SACrC,CAAC,CAAC,2CAA2C;QAC9C,MAAM,QAAQ,GAAG,EAAE,KAAK,EAAE,CAAC,EAAE,IAAI,EAAE,KAAK,EAAE,CAAC;QAC3C,MAAM,QAAQ,GAAG,SAAS,CACxB,YAAY,CAAC,KAA0B,EAAE,QAA6B,CAAC,CACxE,CAAC;QACF,MAAM,MAAM,GAAG,SAAS,CACtB,YAAY,CAAC,KAA0B,EAAE,QAA6B,CAAC,CACxE,CAAC;QAEF,MAAM,CAAC,EAAE,CACP,IAAI,CAAC,GAAG,CAAC,MAAM,GAAG,QAAQ,CAAC,GAAG,IAAI,EAClC,0BAA0B,QAAQ,aAAa,MAAM,EAAE,CACxD,CAAC;IACJ,CAAC,CAAC,CAAC;AACL,CAAC,CAAC,CAAC"}

package/dist/__tests__/compute-mean-anomaly.int.spec.d.ts

@@ -0,0 +1,2 @@
+export {};
+//# sourceMappingURL=compute-mean-anomaly.int.spec.d.ts.map

package/dist/__tests__/compute-mean-anomaly.int.spec.d.ts.map

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

package/dist/__tests__/compute-mean-anomaly.int.spec.js

@@ -0,0 +1,175 @@
+import assert from 'node:assert/strict';
+import test, { describe } from 'node:test';
+import { computeMeanAnomaly } from '../compute-mean-anomaly';
+import { trueAnomalyToMeanAnomaly } from '../true-anomaly-to-mean-anomaly';
+import { wrapAngle } from '../wrap-angle';
+import { TWO_PI } from '@interstellar-tools/constants';
+const EPSILON = 1e-10; // Increased tolerance for floating-point comparisons
+const assertApproxEqual = (actual, expected) => {
+    assert.ok(Math.abs(actual - expected) < EPSILON, `Expected ~${expected}, got ${actual}`);
+};
+describe('computeMeanAnomaly', () => {
+    test('Circular Orbit (e=0)', () => {
+        const body = {
+            e: 0,
+            angle: Math.PI / 2,
+            period: { value: 365, unit: 'd' }
+        };
+        const timeStep = { value: 100, unit: 'd' };
+        const periodInDays = body.period.value;
+        const meanMotion = TWO_PI / Math.abs(periodInDays);
+        const expectedM = wrapAngle(Math.PI / 2 + meanMotion * 100);
+        const result = computeMeanAnomaly(body, timeStep);
+        assertApproxEqual(result, expectedM);
+    });
+    test('Elliptical Orbit (e > 0)', () => {
+        const body = {
+            e: 0.5,
+            angle: Math.PI / 2,
+            period: { value: 200, unit: 'd' }
+        };
+        const timeStep = { value: 50, unit: 'd' };
+        const periodInDays = body.period.value;
+        const meanMotion = TWO_PI / Math.abs(periodInDays);
+        const M0 = body.e === 0 ? body.angle : trueAnomalyToMeanAnomaly(body.angle, body.e);
+        const expectedM = wrapAngle(M0 + meanMotion * 50);
+        const result = computeMeanAnomaly(body, timeStep);
+        assertApproxEqual(result, expectedM);
+    });
+    test('Large Time Step Clamping (Max 10 Orbits)', () => {
+        const body = {
+            e: 0.1,
+            angle: 1.0,
+            period: { value: 100, unit: 'd' }
+        };
+        const timeStep = { value: 2000, unit: 'd' };
+        const periodInDays = body.period.value;
+        const meanMotion = TWO_PI / Math.abs(periodInDays);
+        const maxAllowedTimeStep = 10 * periodInDays;
+        const M0 = trueAnomalyToMeanAnomaly(body.angle, body.e);
+        const expectedM = wrapAngle(M0 + meanMotion * maxAllowedTimeStep);
+        const result = computeMeanAnomaly(body, timeStep);
+        assertApproxEqual(result, expectedM);
+    });
+    test('Eccentricity Out of Bounds Throws Error', () => {
+        const invalidBodies = [
+            { e: -0.1, angle: 0, period: { value: 365, unit: 'd' } },
+            { e: 1.1, angle: 0, period: { value: 365, unit: 'd' } }
+        ];
+        const timeStep = { value: 10, unit: 'd' };
+        for (const body of invalidBodies) {
+            assert.throws(() => computeMeanAnomaly(body, timeStep), RangeError);
+        }
+    });
+    test('Zero Time Step Returns Initial Mean Anomaly', () => {
+        const body = {
+            e: 0.3,
+            angle: 2.0,
+            period: { value: 400, unit: 'd' }
+        };
+        const timeStep = { value: 0, unit: 'd' };
+        const M0 = body.e === 0 ? body.angle : trueAnomalyToMeanAnomaly(body.angle, body.e);
+        const result = computeMeanAnomaly(body, timeStep);
+        assertApproxEqual(result, M0);
+    });
+    test('High Eccentricity Orbit (e ~ 1)', () => {
+        const body = {
+            e: 0.9951,
+            angle: -2.1196,
+            period: { value: 365.25, unit: 'd' }
+        };
+        const timeStep = { value: 0.25, unit: 'd' };
+        const periodInDays = body.period.value;
+        const meanMotion = TWO_PI / Math.abs(periodInDays);
+        const M0 = trueAnomalyToMeanAnomaly(body.angle, body.e);
+        const expectedM = wrapAngle(M0 + meanMotion * 0.25);
+        const result = computeMeanAnomaly(body, timeStep);
+        assertApproxEqual(result, expectedM);
+    });
+    test('Angle Wrapping within [-2π, 2π]', () => {
+        const body = {
+            e: 0.5,
+            angle: 3 * Math.PI,
+            period: { value: 365, unit: 'd' }
+        };
+        const timeStep = { value: 10, unit: 'd' };
+        const periodInDays = body.period.value;
+        const meanMotion = TWO_PI / Math.abs(periodInDays);
+        const M0 = trueAnomalyToMeanAnomaly(body.angle, body.e);
+        const expectedM = wrapAngle(M0 + meanMotion * 10);
+        const result = computeMeanAnomaly(body, timeStep);
+        assertApproxEqual(result, expectedM);
+        assert(result <= 2 * Math.PI && result >= -2 * Math.PI);
+    });
+    test('Minimum Time Step Clamping', () => {
+        const body = {
+            e: 0.5,
+            angle: Math.PI / 2,
+            period: { value: 365, unit: 'd' }
+        };
+        const timeStep = { value: 1e-6, unit: 'd' };
+        const periodInDays = body.period.value;
+        const meanMotion = TWO_PI / Math.abs(periodInDays);
+        const M0 = trueAnomalyToMeanAnomaly(body.angle, body.e);
+        const expectedM = wrapAngle(M0 + meanMotion * 1e-5); // Minimum time step
+        const result = computeMeanAnomaly(body, timeStep);
+        assertApproxEqual(result, expectedM);
+    });
+    test("Halley's Comet", () => {
+        const body = {
+            e: 0.96714,
+            angle: (38.38 * Math.PI) / 180, // Convert degrees to radians
+            period: { value: 27576, unit: 'd' }
+        };
+        const timeStep = { value: 100, unit: 'd' };
+        const periodInDays = body.period.value;
+        const meanMotion = TWO_PI / Math.abs(periodInDays);
+        const M0 = trueAnomalyToMeanAnomaly(body.angle, body.e);
+        const expectedM = wrapAngle(M0 + meanMotion * 100);
+        const result = computeMeanAnomaly(body, timeStep);
+        assertApproxEqual(result, expectedM);
+    });
+    test('Comet Hale-Bopp', () => {
+        const body = {
+            e: 0.9951,
+            angle: (180 * Math.PI) / 180, // Convert degrees to radians
+            period: { value: -253533, unit: 'd' } // Positive period for testing
+        };
+        const timeStep = { value: 0.25, unit: 'd' };
+        const periodInDays = body.period.value;
+        const meanMotion = TWO_PI / Math.abs(periodInDays);
+        const M0 = trueAnomalyToMeanAnomaly(body.angle, body.e);
+        const expectedM = wrapAngle(M0 + meanMotion * 0.25);
+        const result = computeMeanAnomaly(body, timeStep);
+        assertApproxEqual(result, expectedM);
+    });
+    test('Comet 67P/Churyumov–Gerasimenko', () => {
+        const body = {
+            e: 0.641,
+            angle: 0, // Initial angle in radians
+            period: { value: -2484, unit: 'd' } // Positive period for testing
+        };
+        const timeStep = { value: 10, unit: 'd' };
+        const periodInDays = body.period.value;
+        const meanMotion = TWO_PI / Math.abs(periodInDays);
+        const M0 = trueAnomalyToMeanAnomaly(body.angle, body.e);
+        const expectedM = wrapAngle(M0 + meanMotion * 10);
+        const result = computeMeanAnomaly(body, timeStep);
+        assertApproxEqual(result, expectedM);
+    });
+    test('Comet Encke', () => {
+        const body = {
+            e: 0.85,
+            angle: (160 * Math.PI) / 180, // Convert degrees to radians
+            period: { value: -1204, unit: 'd' } // Positive period for testing
+        };
+        const timeStep = { value: 5, unit: 'd' };
+        const periodInDays = body.period.value;
+        const meanMotion = TWO_PI / Math.abs(periodInDays);
+        const M0 = trueAnomalyToMeanAnomaly(body.angle, body.e);
+        const expectedM = wrapAngle(M0 + meanMotion * 5);
+        const result = computeMeanAnomaly(body, timeStep);
+        assertApproxEqual(result, expectedM);
+    });
+});
+//# sourceMappingURL=compute-mean-anomaly.int.spec.js.map

package/dist/__tests__/compute-mean-anomaly.int.spec.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"compute-mean-anomaly.int.spec.js","sourceRoot":"","sources":["../../src/__tests__/compute-mean-anomaly.int.spec.ts"],"names":[],"mappings":"AAAA,OAAO,MAAM,MAAM,oBAAoB,CAAC;AACxC,OAAO,IAAI,EAAE,EAAE,QAAQ,EAAE,MAAM,WAAW,CAAC;AAE3C,OAAO,EAAE,kBAAkB,EAAE,MAAM,yBAAyB,CAAC;AAC7D,OAAO,EAAE,wBAAwB,EAAE,MAAM,iCAAiC,CAAC;AAC3E,OAAO,EAAE,SAAS,EAAE,MAAM,eAAe,CAAC;AAK1C,OAAO,EAAE,MAAM,EAAE,MAAM,+BAA+B,CAAC;AAEvD,MAAM,OAAO,GAAG,KAAK,CAAC,CAAC,qDAAqD;AAC5E,MAAM,iBAAiB,GAAG,CAAC,MAAc,EAAE,QAAgB,EAAE,EAAE;IAC7D,MAAM,CAAC,EAAE,CACP,IAAI,CAAC,GAAG,CAAC,MAAM,GAAG,QAAQ,CAAC,GAAG,OAAO,EACrC,aAAa,QAAQ,SAAS,MAAM,EAAE,CACvC,CAAC;AACJ,CAAC,CAAC;AAEF,QAAQ,CAAC,oBAAoB,EAAE,GAAG,EAAE;IAClC,IAAI,CAAC,sBAAsB,EAAE,GAAG,EAAE;QAChC,MAAM,IAAI,GAAsB;YAC9B,CAAC,EAAE,CAAC;YACJ,KAAK,EAAE,IAAI,CAAC,EAAE,GAAG,CAAC;YAClB,MAAM,EAAE,EAAE,KAAK,EAAE,GAAG,EAAE,IAAI,EAAE,GAAG,EAAE;SACb,CAAC;QACvB,MAAM,QAAQ,GAAsB,EAAE,KAAK,EAAE,GAAG,EAAE,IAAI,EAAE,GAAG,EAAE,CAAC;QAC9D,MAAM,YAAY,GAAG,IAAI,CAAC,MAAM,CAAC,KAAK,CAAC;QACvC,MAAM,UAAU,GAAG,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,YAAY,CAAC,CAAC;QACnD,MAAM,SAAS,GAAG,SAAS,CAAC,IAAI,CAAC,EAAE,GAAG,CAAC,GAAG,UAAU,GAAG,GAAG,CAAC,CAAC;QAC5D,MAAM,MAAM,GAAG,kBAAkB,CAAC,IAAI,EAAE,QAAQ,CAAC,CAAC;QAElD,iBAAiB,CAAC,MAAM,EAAE,SAAS,CAAC,CAAC;IACvC,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,0BAA0B,EAAE,GAAG,EAAE;QACpC,MAAM,IAAI,GAAsB;YAC9B,CAAC,EAAE,GAAG;YACN,KAAK,EAAE,IAAI,CAAC,EAAE,GAAG,CAAC;YAClB,MAAM,EAAE,EAAE,KAAK,EAAE,GAAG,EAAE,IAAI,EAAE,GAAG,EAAE;SACb,CAAC;QACvB,MAAM,QAAQ,GAAsB,EAAE,KAAK,EAAE,EAAE,EAAE,IAAI,EAAE,GAAG,EAAE,CAAC;QAC7D,MAAM,YAAY,GAAG,IAAI,CAAC,MAAM,CAAC,KAAK,CAAC;QACvC,MAAM,UAAU,GAAG,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,YAAY,CAAC,CAAC;QACnD,MAAM,EAAE,GACN,IAAI,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,wBAAwB,CAAC,IAAI,CAAC,KAAK,EAAE,IAAI,CAAC,CAAC,CAAC,CAAC;QAC3E,MAAM,SAAS,GAAG,SAAS,CAAC,EAAE,GAAG,UAAU,GAAG,EAAE,CAAC,CAAC;QAClD,MAAM,MAAM,GAAG,kBAAkB,CAAC,IAAI,EAAE,QAAQ,CAAC,CAAC;QAElD,iBAAiB,CAAC,MAAM,EAAE,SAAS,CAAC,CAAC;IACvC,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,0CAA0C,EAAE,GAAG,EAAE;QACpD,MAAM,IAAI,GAAsB;YAC9B,CAAC,EAAE,GAAG;YACN,KAAK,EAAE,GAAG;YACV,MAAM,EAAE,EAAE,KAAK,EAAE,GAAG,EAAE,IAAI,EAAE,GAAG,EAAE;SACb,CAAC;QACvB,MAAM,QAAQ,GAAsB,EAAE,KAAK,EAAE,IAAI,EAAE,IAAI,EAAE,GAAG,EAAE,CAAC;QAC/D,MAAM,YAAY,GAAG,IAAI,CAAC,MAAM,CAAC,KAAK,CAAC;QACvC,MAAM,UAAU,GAAG,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,YAAY,CAAC,CAAC;QACnD,MAAM,kBAAkB,GAAG,EAAE,GAAG,YAAY,CAAC;QAC7C,MAAM,EAAE,GAAG,wBAAwB,CAAC,IAAI,CAAC,KAAK,EAAE,IAAI,CAAC,CAAC,CAAC,CAAC;QACxD,MAAM,SAAS,GAAG,SAAS,CAAC,EAAE,GAAG,UAAU,GAAG,kBAAkB,CAAC,CAAC;QAClE,MAAM,MAAM,GAAG,kBAAkB,CAAC,IAAI,EAAE,QAAQ,CAAC,CAAC;QAElD,iBAAiB,CAAC,MAAM,EAAE,SAAS,CAAC,CAAC;IACvC,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,yCAAyC,EAAE,GAAG,EAAE;QACnD,MAAM,aAAa,GAAwB;YACzC,EAAE,CAAC,EAAE,CAAC,GAAG,EAAE,KAAK,EAAE,CAAC,EAAE,MAAM,EAAE,EAAE,KAAK,EAAE,GAAG,EAAE,IAAI,EAAE,GAAG,EAAE,EAAE;YACxD,EAAE,CAAC,EAAE,GAAG,EAAE,KAAK,EAAE,CAAC,EAAE,MAAM,EAAE,EAAE,KAAK,EAAE,GAAG,EAAE,IAAI,EAAE,GAAG,EAAE,EAAE;SACjC,CAAC;QACzB,MAAM,QAAQ,GAAsB,EAAE,KAAK,EAAE,EAAE,EAAE,IAAI,EAAE,GAAG,EAAE,CAAC;QAE7D,KAAK,MAAM,IAAI,IAAI,aAAa,EAAE,CAAC;YACjC,MAAM,CAAC,MAAM,CAAC,GAAG,EAAE,CAAC,kBAAkB,CAAC,IAAI,EAAE,QAAQ,CAAC,EAAE,UAAU,CAAC,CAAC;QACtE,CAAC;IACH,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,6CAA6C,EAAE,GAAG,EAAE;QACvD,MAAM,IAAI,GAAsB;YAC9B,CAAC,EAAE,GAAG;YACN,KAAK,EAAE,GAAG;YACV,MAAM,EAAE,EAAE,KAAK,EAAE,GAAG,EAAE,IAAI,EAAE,GAAG,EAAE;SACb,CAAC;QACvB,MAAM,QAAQ,GAAsB,EAAE,KAAK,EAAE,CAAC,EAAE,IAAI,EAAE,GAAG,EAAE,CAAC;QAC5D,MAAM,EAAE,GACN,IAAI,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,wBAAwB,CAAC,IAAI,CAAC,KAAK,EAAE,IAAI,CAAC,CAAC,CAAC,CAAC;QAC3E,MAAM,MAAM,GAAG,kBAAkB,CAAC,IAAI,EAAE,QAAQ,CAAC,CAAC;QAElD,iBAAiB,CAAC,MAAM,EAAE,EAAE,CAAC,CAAC;IAChC,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,iCAAiC,EAAE,GAAG,EAAE;QAC3C,MAAM,IAAI,GAAsB;YAC9B,CAAC,EAAE,MAAM;YACT,KAAK,EAAE,CAAC,MAAM;YACd,MAAM,EAAE,EAAE,KAAK,EAAE,MAAM,EAAE,IAAI,EAAE,GAAG,EAAE;SAChB,CAAC;QACvB,MAAM,QAAQ,GAAsB,EAAE,KAAK,EAAE,IAAI,EAAE,IAAI,EAAE,GAAG,EAAE,CAAC;QAC/D,MAAM,YAAY,GAAG,IAAI,CAAC,MAAM,CAAC,KAAK,CAAC;QACvC,MAAM,UAAU,GAAG,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,YAAY,CAAC,CAAC;QACnD,MAAM,EAAE,GAAG,wBAAwB,CAAC,IAAI,CAAC,KAAK,EAAE,IAAI,CAAC,CAAC,CAAC,CAAC;QACxD,MAAM,SAAS,GAAG,SAAS,CAAC,EAAE,GAAG,UAAU,GAAG,IAAI,CAAC,CAAC;QACpD,MAAM,MAAM,GAAG,kBAAkB,CAAC,IAAI,EAAE,QAAQ,CAAC,CAAC;QAElD,iBAAiB,CAAC,MAAM,EAAE,SAAS,CAAC,CAAC;IACvC,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,iCAAiC,EAAE,GAAG,EAAE;QAC3C,MAAM,IAAI,GAAsB;YAC9B,CAAC,EAAE,GAAG;YACN,KAAK,EAAE,CAAC,GAAG,IAAI,CAAC,EAAE;YAClB,MAAM,EAAE,EAAE,KAAK,EAAE,GAAG,EAAE,IAAI,EAAE,GAAG,EAAE;SACb,CAAC;QACvB,MAAM,QAAQ,GAAsB,EAAE,KAAK,EAAE,EAAE,EAAE,IAAI,EAAE,GAAG,EAAE,CAAC;QAC7D,MAAM,YAAY,GAAG,IAAI,CAAC,MAAM,CAAC,KAAK,CAAC;QACvC,MAAM,UAAU,GAAG,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,YAAY,CAAC,CAAC;QACnD,MAAM,EAAE,GAAG,wBAAwB,CAAC,IAAI,CAAC,KAAK,EAAE,IAAI,CAAC,CAAC,CAAC,CAAC;QACxD,MAAM,SAAS,GAAG,SAAS,CAAC,EAAE,GAAG,UAAU,GAAG,EAAE,CAAC,CAAC;QAClD,MAAM,MAAM,GAAG,kBAAkB,CAAC,IAAI,EAAE,QAAQ,CAAC,CAAC;QAElD,iBAAiB,CAAC,MAAM,EAAE,SAAS,CAAC,CAAC;QACrC,MAAM,CAAC,MAAM,IAAI,CAAC,GAAG,IAAI,CAAC,EAAE,IAAI,MAAM,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,CAAC;IAC1D,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,4BAA4B,EAAE,GAAG,EAAE;QACtC,MAAM,IAAI,GAAsB;YAC9B,CAAC,EAAE,GAAG;YACN,KAAK,EAAE,IAAI,CAAC,EAAE,GAAG,CAAC;YAClB,MAAM,EAAE,EAAE,KAAK,EAAE,GAAG,EAAE,IAAI,EAAE,GAAG,EAAE;SACb,CAAC;QACvB,MAAM,QAAQ,GAAsB,EAAE,KAAK,EAAE,IAAI,EAAE,IAAI,EAAE,GAAG,EAAE,CAAC;QAC/D,MAAM,YAAY,GAAG,IAAI,CAAC,MAAM,CAAC,KAAK,CAAC;QACvC,MAAM,UAAU,GAAG,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,YAAY,CAAC,CAAC;QACnD,MAAM,EAAE,GAAG,wBAAwB,CAAC,IAAI,CAAC,KAAK,EAAE,IAAI,CAAC,CAAC,CAAC,CAAC;QACxD,MAAM,SAAS,GAAG,SAAS,CAAC,EAAE,GAAG,UAAU,GAAG,IAAI,CAAC,CAAC,CAAC,oBAAoB;QACzE,MAAM,MAAM,GAAG,kBAAkB,CAAC,IAAI,EAAE,QAAQ,CAAC,CAAC;QAElD,iBAAiB,CAAC,MAAM,EAAE,SAAS,CAAC,CAAC;IACvC,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,gBAAgB,EAAE,GAAG,EAAE;QAC1B,MAAM,IAAI,GAAsB;YAC9B,CAAC,EAAE,OAAO;YACV,KAAK,EAAE,CAAC,KAAK,GAAG,IAAI,CAAC,EAAE,CAAC,GAAG,GAAG,EAAE,6BAA6B;YAC7D,MAAM,EAAE,EAAE,KAAK,EAAE,KAAK,EAAE,IAAI,EAAE,GAAG,EAAE;SACf,CAAC;QACvB,MAAM,QAAQ,GAAsB,EAAE,KAAK,EAAE,GAAG,EAAE,IAAI,EAAE,GAAG,EAAE,CAAC;QAC9D,MAAM,YAAY,GAAG,IAAI,CAAC,MAAM,CAAC,KAAK,CAAC;QACvC,MAAM,UAAU,GAAG,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,YAAY,CAAC,CAAC;QACnD,MAAM,EAAE,GAAG,wBAAwB,CAAC,IAAI,CAAC,KAAK,EAAE,IAAI,CAAC,CAAC,CAAC,CAAC;QACxD,MAAM,SAAS,GAAG,SAAS,CAAC,EAAE,GAAG,UAAU,GAAG,GAAG,CAAC,CAAC;QACnD,MAAM,MAAM,GAAG,kBAAkB,CAAC,IAAI,EAAE,QAAQ,CAAC,CAAC;QAElD,iBAAiB,CAAC,MAAM,EAAE,SAAS,CAAC,CAAC;IACvC,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,iBAAiB,EAAE,GAAG,EAAE;QAC3B,MAAM,IAAI,GAAsB;YAC9B,CAAC,EAAE,MAAM;YACT,KAAK,EAAE,CAAC,GAAG,GAAG,IAAI,CAAC,EAAE,CAAC,GAAG,GAAG,EAAE,6BAA6B;YAC3D,MAAM,EAAE,EAAE,KAAK,EAAE,CAAC,MAAM,EAAE,IAAI,EAAE,GAAG,EAAE,CAAC,8BAA8B;SAChD,CAAC;QACvB,MAAM,QAAQ,GAAsB,EAAE,KAAK,EAAE,IAAI,EAAE,IAAI,EAAE,GAAG,EAAE,CAAC;QAC/D,MAAM,YAAY,GAAG,IAAI,CAAC,MAAM,CAAC,KAAK,CAAC;QACvC,MAAM,UAAU,GAAG,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,YAAY,CAAC,CAAC;QACnD,MAAM,EAAE,GAAG,wBAAwB,CAAC,IAAI,CAAC,KAAK,EAAE,IAAI,CAAC,CAAC,CAAC,CAAC;QACxD,MAAM,SAAS,GAAG,SAAS,CAAC,EAAE,GAAG,UAAU,GAAG,IAAI,CAAC,CAAC;QACpD,MAAM,MAAM,GAAG,kBAAkB,CAAC,IAAI,EAAE,QAAQ,CAAC,CAAC;QAElD,iBAAiB,CAAC,MAAM,EAAE,SAAS,CAAC,CAAC;IACvC,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,iCAAiC,EAAE,GAAG,EAAE;QAC3C,MAAM,IAAI,GAAsB;YAC9B,CAAC,EAAE,KAAK;YACR,KAAK,EAAE,CAAC,EAAE,2BAA2B;YACrC,MAAM,EAAE,EAAE,KAAK,EAAE,CAAC,IAAI,EAAE,IAAI,EAAE,GAAG,EAAE,CAAC,8BAA8B;SAC9C,CAAC;QACvB,MAAM,QAAQ,GAAsB,EAAE,KAAK,EAAE,EAAE,EAAE,IAAI,EAAE,GAAG,EAAE,CAAC;QAC7D,MAAM,YAAY,GAAG,IAAI,CAAC,MAAM,CAAC,KAAK,CAAC;QACvC,MAAM,UAAU,GAAG,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,YAAY,CAAC,CAAC;QACnD,MAAM,EAAE,GAAG,wBAAwB,CAAC,IAAI,CAAC,KAAK,EAAE,IAAI,CAAC,CAAC,CAAC,CAAC;QACxD,MAAM,SAAS,GAAG,SAAS,CAAC,EAAE,GAAG,UAAU,GAAG,EAAE,CAAC,CAAC;QAClD,MAAM,MAAM,GAAG,kBAAkB,CAAC,IAAI,EAAE,QAAQ,CAAC,CAAC;QAElD,iBAAiB,CAAC,MAAM,EAAE,SAAS,CAAC,CAAC;IACvC,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,aAAa,EAAE,GAAG,EAAE;QACvB,MAAM,IAAI,GAAsB;YAC9B,CAAC,EAAE,IAAI;YACP,KAAK,EAAE,CAAC,GAAG,GAAG,IAAI,CAAC,EAAE,CAAC,GAAG,GAAG,EAAE,6BAA6B;YAC3D,MAAM,EAAE,EAAE,KAAK,EAAE,CAAC,IAAI,EAAE,IAAI,EAAE,GAAG,EAAE,CAAC,8BAA8B;SAC9C,CAAC;QACvB,MAAM,QAAQ,GAAsB,EAAE,KAAK,EAAE,CAAC,EAAE,IAAI,EAAE,GAAG,EAAE,CAAC;QAC5D,MAAM,YAAY,GAAG,IAAI,CAAC,MAAM,CAAC,KAAK,CAAC;QACvC,MAAM,UAAU,GAAG,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,YAAY,CAAC,CAAC;QACnD,MAAM,EAAE,GAAG,wBAAwB,CAAC,IAAI,CAAC,KAAK,EAAE,IAAI,CAAC,CAAC,CAAC,CAAC;QACxD,MAAM,SAAS,GAAG,SAAS,CAAC,EAAE,GAAG,UAAU,GAAG,CAAC,CAAC,CAAC;QACjD,MAAM,MAAM,GAAG,kBAAkB,CAAC,IAAI,EAAE,QAAQ,CAAC,CAAC;QAElD,iBAAiB,CAAC,MAAM,EAAE,SAAS,CAAC,CAAC;IACvC,CAAC,CAAC,CAAC;AACL,CAAC,CAAC,CAAC"}

package/dist/__tests__/eccentric-to-true-anomaly.spec.d.ts

@@ -0,0 +1,2 @@
+export {};
+//# sourceMappingURL=eccentric-to-true-anomaly.spec.d.ts.map

package/dist/__tests__/eccentric-to-true-anomaly.spec.d.ts.map

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

package/dist/__tests__/eccentric-to-true-anomaly.spec.js

@@ -0,0 +1,33 @@
+import assert from 'node:assert/strict';
+import test, { describe } from 'node:test';
+import { eccentricToTrueAnomaly } from '../eccentric-to-true-anomaly';
+describe('eccentricToTrueAnomaly', () => {
+    // Test cases for eccentricToTrueAnomaly
+    const testCases = [
+        { E: 0, e: 0.5, expected: 0.0 }, // E = 0 should give V = 0
+        { E: Math.PI / 4, e: 0.1, expected: 0.858858 }, // Approximate expected value
+        { E: Math.PI / 2, e: 0.3, expected: 1.875489 }, // Approximate expected value
+        { E: Math.PI, e: 0.2, expected: Math.PI }, // E = π should give V = π
+        { E: (3 * Math.PI) / 2, e: 0.4, expected: 4.300872 } // Approximate expected value
+    ];
+    test('computes correct true anomaly', () => {
+        testCases.forEach(({ E, e, expected }) => {
+            const result = eccentricToTrueAnomaly(E, e);
+            assert.ok(Math.abs(result - expected) < 1e-4, `Expected V ≈ ${expected} for E=${E}, e=${e}, but got ${result}`);
+        });
+    });
+    test('handles special cases', () => {
+        assert.strictEqual(eccentricToTrueAnomaly(0, 0), 0); // Circular orbit (e = 0) should return E
+        assert.strictEqual(eccentricToTrueAnomaly(Math.PI, 0), Math.PI); // Circular case with E = π
+        assert.strictEqual(eccentricToTrueAnomaly(0, 1), 0); // Parabolic orbit case
+    });
+    test('throws on invalid eccentricity', () => {
+        assert.throws(() => eccentricToTrueAnomaly(Math.PI / 3, -0.1), RangeError);
+        assert.throws(() => eccentricToTrueAnomaly(Math.PI / 3, 1.2), RangeError);
+    });
+    test('ensures numerical stability near E = π', () => {
+        assert.strictEqual(eccentricToTrueAnomaly(Math.PI + 1e-11, 0.5), Math.PI);
+        assert.strictEqual(eccentricToTrueAnomaly(Math.PI - 1e-11, 0.5), Math.PI);
+    });
+});
+//# sourceMappingURL=eccentric-to-true-anomaly.spec.js.map

package/dist/__tests__/eccentric-to-true-anomaly.spec.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"eccentric-to-true-anomaly.spec.js","sourceRoot":"","sources":["../../src/__tests__/eccentric-to-true-anomaly.spec.ts"],"names":[],"mappings":"AAAA,OAAO,MAAM,MAAM,oBAAoB,CAAC;AACxC,OAAO,IAAI,EAAE,EAAE,QAAQ,EAAE,MAAM,WAAW,CAAC;AAE3C,OAAO,EAAE,sBAAsB,EAAE,MAAM,8BAA8B,CAAC;AAEtE,QAAQ,CAAC,wBAAwB,EAAE,GAAG,EAAE;IACtC,wCAAwC;IACxC,MAAM,SAAS,GAAG;QAChB,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,GAAG,EAAE,QAAQ,EAAE,GAAG,EAAE,EAAE,0BAA0B;QAC3D,EAAE,CAAC,EAAE,IAAI,CAAC,EAAE,GAAG,CAAC,EAAE,CAAC,EAAE,GAAG,EAAE,QAAQ,EAAE,QAAQ,EAAE,EAAE,6BAA6B;QAC7E,EAAE,CAAC,EAAE,IAAI,CAAC,EAAE,GAAG,CAAC,EAAE,CAAC,EAAE,GAAG,EAAE,QAAQ,EAAE,QAAQ,EAAE,EAAE,6BAA6B;QAC7E,EAAE,CAAC,EAAE,IAAI,CAAC,EAAE,EAAE,CAAC,EAAE,GAAG,EAAE,QAAQ,EAAE,IAAI,CAAC,EAAE,EAAE,EAAE,0BAA0B;QACrE,EAAE,CAAC,EAAE,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,GAAG,EAAE,QAAQ,EAAE,QAAQ,EAAE,CAAC,6BAA6B;KACnF,CAAC;IAEF,IAAI,CAAC,+BAA+B,EAAE,GAAG,EAAE;QACzC,SAAS,CAAC,OAAO,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,QAAQ,EAAE,EAAE,EAAE;YACvC,MAAM,MAAM,GAAG,sBAAsB,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;YAE5C,MAAM,CAAC,EAAE,CACP,IAAI,CAAC,GAAG,CAAC,MAAM,GAAG,QAAQ,CAAC,GAAG,IAAI,EAClC,gBAAgB,QAAQ,UAAU,CAAC,OAAO,CAAC,aAAa,MAAM,EAAE,CACjE,CAAC;QACJ,CAAC,CAAC,CAAC;IACL,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,uBAAuB,EAAE,GAAG,EAAE;QACjC,MAAM,CAAC,WAAW,CAAC,sBAAsB,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,yCAAyC;QAC9F,MAAM,CAAC,WAAW,CAAC,sBAAsB,CAAC,IAAI,CAAC,EAAE,EAAE,CAAC,CAAC,EAAE,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,2BAA2B;QAC5F,MAAM,CAAC,WAAW,CAAC,sBAAsB,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,uBAAuB;IAC9E,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,gCAAgC,EAAE,GAAG,EAAE;QAC1C,MAAM,CAAC,MAAM,CAAC,GAAG,EAAE,CAAC,sBAAsB,CAAC,IAAI,CAAC,EAAE,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,UAAU,CAAC,CAAC;QAC3E,MAAM,CAAC,MAAM,CAAC,GAAG,EAAE,CAAC,sBAAsB,CAAC,IAAI,CAAC,EAAE,GAAG,CAAC,EAAE,GAAG,CAAC,EAAE,UAAU,CAAC,CAAC;IAC5E,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,wCAAwC,EAAE,GAAG,EAAE;QAClD,MAAM,CAAC,WAAW,CAAC,sBAAsB,CAAC,IAAI,CAAC,EAAE,GAAG,KAAK,EAAE,GAAG,CAAC,EAAE,IAAI,CAAC,EAAE,CAAC,CAAC;QAC1E,MAAM,CAAC,WAAW,CAAC,sBAAsB,CAAC,IAAI,CAAC,EAAE,GAAG,KAAK,EAAE,GAAG,CAAC,EAAE,IAAI,CAAC,EAAE,CAAC,CAAC;IAC5E,CAAC,CAAC,CAAC;AACL,CAAC,CAAC,CAAC"}

package/dist/__tests__/solve-kepler-bisection.spec.d.ts

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

package/dist/__tests__/solve-kepler-bisection.spec.d.ts.map

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

package/dist/__tests__/solve-kepler-bisection.spec.js

@@ -0,0 +1,41 @@
+import assert from 'node:assert/strict';
+import test, { describe } from 'node:test';
+import { solveKeplerBisection } from '../solve-kepler-bisection';
+describe('solveKeplerBisection', () => {
+    test("Solves Kepler's equation for small eccentricities", () => {
+        const M = Math.PI / 4; // 45 degrees in radians
+        const e = 0.1;
+        const expected = solveKeplerBisection(M, e, 50, 1e-9);
+        const result = solveKeplerBisection(M, e, 50, 1e-9);
+        assert.ok(Math.abs(result - expected) < 1e-9, `Expected ${expected}, got ${result}`);
+    });
+    test("Solves Kepler's equation for moderate eccentricities", () => {
+        const M = Math.PI / 2; // 90 degrees in radians
+        const e = 0.5;
+        const expected = solveKeplerBisection(M, e, 50, 1e-9);
+        const result = solveKeplerBisection(M, e, 50, 1e-9);
+        assert.ok(Math.abs(result - expected) < 1e-9, `Expected ${expected}, got ${result}`);
+    });
+    test("Solves Kepler's equation for high eccentricities", () => {
+        const M = Math.PI / 3; // 60 degrees in radians
+        const e = 0.9;
+        const expected = solveKeplerBisection(M, e, 50, 1e-9);
+        const result = solveKeplerBisection(M, e, 50, 1e-9);
+        assert.ok(Math.abs(result - expected) < 1e-9, `Expected ${expected}, got ${result}`);
+    });
+    test('Ensures solution falls within expected bounds', () => {
+        const M = Math.PI / 6;
+        const e = 0.4;
+        const result = solveKeplerBisection(M, e, 50, 1e-9);
+        assert.ok(result >= 0 && result <= Math.PI, `Eccentric anomaly ${result} out of bounds`);
+    });
+    test('Handles convergence within the given tolerance', () => {
+        const M = Math.PI / 2;
+        const e = 0.7;
+        const tolerance = 1e-12;
+        const result = solveKeplerBisection(M, e, 50, tolerance);
+        const F = result - e * Math.sin(result) - M;
+        assert.ok(Math.abs(F) < tolerance, `Convergence issue: Residual ${F}`);
+    });
+});
+//# sourceMappingURL=solve-kepler-bisection.spec.js.map

package/dist/__tests__/solve-kepler-bisection.spec.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"solve-kepler-bisection.spec.js","sourceRoot":"","sources":["../../src/__tests__/solve-kepler-bisection.spec.ts"],"names":[],"mappings":"AAAA,OAAO,MAAM,MAAM,oBAAoB,CAAC;AACxC,OAAO,IAAI,EAAE,EAAE,QAAQ,EAAE,MAAM,WAAW,CAAC;AAE3C,OAAO,EAAE,oBAAoB,EAAE,MAAM,2BAA2B,CAAC;AAEjE,QAAQ,CAAC,sBAAsB,EAAE,GAAG,EAAE;IACpC,IAAI,CAAC,mDAAmD,EAAE,GAAG,EAAE;QAC7D,MAAM,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,wBAAwB;QAC/C,MAAM,CAAC,GAAG,GAAG,CAAC;QACd,MAAM,QAAQ,GAAG,oBAAoB,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,EAAE,IAAI,CAAC,CAAC;QACtD,MAAM,MAAM,GAAG,oBAAoB,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,EAAE,IAAI,CAAC,CAAC;QAEpD,MAAM,CAAC,EAAE,CACP,IAAI,CAAC,GAAG,CAAC,MAAM,GAAG,QAAQ,CAAC,GAAG,IAAI,EAClC,YAAY,QAAQ,SAAS,MAAM,EAAE,CACtC,CAAC;IACJ,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,sDAAsD,EAAE,GAAG,EAAE;QAChE,MAAM,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,wBAAwB;QAC/C,MAAM,CAAC,GAAG,GAAG,CAAC;QACd,MAAM,QAAQ,GAAG,oBAAoB,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,EAAE,IAAI,CAAC,CAAC;QACtD,MAAM,MAAM,GAAG,oBAAoB,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,EAAE,IAAI,CAAC,CAAC;QAEpD,MAAM,CAAC,EAAE,CACP,IAAI,CAAC,GAAG,CAAC,MAAM,GAAG,QAAQ,CAAC,GAAG,IAAI,EAClC,YAAY,QAAQ,SAAS,MAAM,EAAE,CACtC,CAAC;IACJ,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,kDAAkD,EAAE,GAAG,EAAE;QAC5D,MAAM,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,wBAAwB;QAC/C,MAAM,CAAC,GAAG,GAAG,CAAC;QACd,MAAM,QAAQ,GAAG,oBAAoB,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,EAAE,IAAI,CAAC,CAAC;QACtD,MAAM,MAAM,GAAG,oBAAoB,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,EAAE,IAAI,CAAC,CAAC;QAEpD,MAAM,CAAC,EAAE,CACP,IAAI,CAAC,GAAG,CAAC,MAAM,GAAG,QAAQ,CAAC,GAAG,IAAI,EAClC,YAAY,QAAQ,SAAS,MAAM,EAAE,CACtC,CAAC;IACJ,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,+CAA+C,EAAE,GAAG,EAAE;QACzD,MAAM,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC;QACtB,MAAM,CAAC,GAAG,GAAG,CAAC;QACd,MAAM,MAAM,GAAG,oBAAoB,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,EAAE,IAAI,CAAC,CAAC;QAEpD,MAAM,CAAC,EAAE,CACP,MAAM,IAAI,CAAC,IAAI,MAAM,IAAI,IAAI,CAAC,EAAE,EAChC,qBAAqB,MAAM,gBAAgB,CAC5C,CAAC;IACJ,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,gDAAgD,EAAE,GAAG,EAAE;QAC1D,MAAM,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC;QACtB,MAAM,CAAC,GAAG,GAAG,CAAC;QACd,MAAM,SAAS,GAAG,KAAK,CAAC;QACxB,MAAM,MAAM,GAAG,oBAAoB,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,EAAE,SAAS,CAAC,CAAC;QACzD,MAAM,CAAC,GAAG,MAAM,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,MAAM,CAAC,GAAG,CAAC,CAAC;QAE5C,MAAM,CAAC,EAAE,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,SAAS,EAAE,+BAA+B,CAAC,EAAE,CAAC,CAAC;IACzE,CAAC,CAAC,CAAC;AACL,CAAC,CAAC,CAAC"}

package/dist/__tests__/solve-kepler-high-eccentricity.spec.d.ts

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

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

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

package/dist/__tests__/solve-kepler-high-eccentricity.spec.js

@@ -0,0 +1,81 @@
+import assert from 'node:assert/strict';
+import test, { describe } from 'node:test';
+import { solveKeplerHighEccentricity } from '../solve-kepler-high-eccentricity';
+const EPSILON = 1e-9;
+const assertApproxEqual = (actual, expected) => {
+    assert.ok(Math.abs(actual - expected) < EPSILON, `Expected ~${expected}, got ${actual}`);
+};
+// **Define test cases based on Halley's Comet & Comet Hale-Bopp**
+const testCases = [
+    {
+        name: "Halley's Comet (Typical Case)",
+        M: (38.38 * Math.PI) / 180, // Mean anomaly converted to radians
+        e: 0.96714,
+        expected: 1.63500444611954 // Placeholder expected value (needs verification)
+    },
+    {
+        name: 'Comet Hale-Bopp (Extreme Eccentricity)',
+        M: (180 * Math.PI) / 180,
+        e: 0.9951,
+        expected: Math.PI // Placeholder expected value (Kepler's equation should return π for M=π)
+    },
+    {
+        name: "Halley's Comet (Small M)",
+        M: 0.01, // Small anomaly case
+        e: 0.96714,
+        expected: 0.238206497655823 // Placeholder expected value (needs verification)
+    },
+    {
+        name: 'Comet Hale-Bopp (Large M)',
+        M: 5, // Large anomaly case
+        e: 0.9951,
+        expected: 4.1553401654182 // Placeholder expected value (needs verification)
+    }
+];
+describe('solveKeplerHighEccentricity (Real Data Tests)', () => {
+    for (const { name, M, e, expected } of testCases) {
+        test(name, () => {
+            const result = solveKeplerHighEccentricity(M, e);
+            assertApproxEqual(result, expected);
+        });
+    }
+});
+describe('solveKeplerHighEccentricity (General Behavior Tests)', () => {
+    test('converges for typical values', () => {
+        const M = Math.PI / 4; // 45 degrees in radians
+        const e = 0.95; // High eccentricity
+        const result = solveKeplerHighEccentricity(M, e);
+        assertApproxEqual(result, 1.7242367109653982);
+    });
+    test('works for small M values', () => {
+        const M = 0.01;
+        const e = 0.99;
+        const result = solveKeplerHighEccentricity(M, e);
+        assertApproxEqual(result, 0.3422703164917758);
+    });
+    test('handles M = 0 case correctly', () => {
+        const M = 0;
+        const e = 0.98;
+        const result = solveKeplerHighEccentricity(M, e);
+        assertApproxEqual(result, 0); // Eccentric anomaly should also be zero
+    });
+    test('handles large M values', () => {
+        const M = 5; // Large mean anomaly
+        const e = 0.95;
+        const result = solveKeplerHighEccentricity(M, e);
+        assertApproxEqual(result, 4.181000769932799);
+    });
+    test('ensures correct output for extremely large M', () => {
+        const M = 10 * Math.PI; // Large anomaly
+        const e = 0.97;
+        const result = solveKeplerHighEccentricity(M, e, Math.max(50, Math.floor(M)));
+        assertApproxEqual(result, 0);
+    });
+    test('does not exceed max iterations', () => {
+        const M = 3;
+        const e = 0.999; // Extremely high eccentricity
+        const result = solveKeplerHighEccentricity(M, e, Math.max(50, Math.floor(M)), 1e-9);
+        assert.ok(result !== null); // Ensure it still produces a valid number
+    });
+});
+//# sourceMappingURL=solve-kepler-high-eccentricity.spec.js.map

package/dist/__tests__/solve-kepler-high-eccentricity.spec.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"solve-kepler-high-eccentricity.spec.js","sourceRoot":"","sources":["../../src/__tests__/solve-kepler-high-eccentricity.spec.ts"],"names":[],"mappings":"AAAA,OAAO,MAAM,MAAM,oBAAoB,CAAC;AACxC,OAAO,IAAI,EAAE,EAAE,QAAQ,EAAE,MAAM,WAAW,CAAC;AAE3C,OAAO,EAAE,2BAA2B,EAAE,MAAM,mCAAmC,CAAC;AAEhF,MAAM,OAAO,GAAG,IAAI,CAAC;AACrB,MAAM,iBAAiB,GAAG,CAAC,MAAc,EAAE,QAAgB,EAAE,EAAE;IAC7D,MAAM,CAAC,EAAE,CACP,IAAI,CAAC,GAAG,CAAC,MAAM,GAAG,QAAQ,CAAC,GAAG,OAAO,EACrC,aAAa,QAAQ,SAAS,MAAM,EAAE,CACvC,CAAC;AACJ,CAAC,CAAC;AACF,kEAAkE;AAClE,MAAM,SAAS,GAAG;IAChB;QACE,IAAI,EAAE,+BAA+B;QACrC,CAAC,EAAE,CAAC,KAAK,GAAG,IAAI,CAAC,EAAE,CAAC,GAAG,GAAG,EAAE,oCAAoC;QAChE,CAAC,EAAE,OAAO;QACV,QAAQ,EAAE,gBAAgB,CAAC,kDAAkD;KAC9E;IACD;QACE,IAAI,EAAE,wCAAwC;QAC9C,CAAC,EAAE,CAAC,GAAG,GAAG,IAAI,CAAC,EAAE,CAAC,GAAG,GAAG;QACxB,CAAC,EAAE,MAAM;QACT,QAAQ,EAAE,IAAI,CAAC,EAAE,CAAC,yEAAyE;KAC5F;IACD;QACE,IAAI,EAAE,0BAA0B;QAChC,CAAC,EAAE,IAAI,EAAE,qBAAqB;QAC9B,CAAC,EAAE,OAAO;QACV,QAAQ,EAAE,iBAAiB,CAAC,kDAAkD;KAC/E;IACD;QACE,IAAI,EAAE,2BAA2B;QACjC,CAAC,EAAE,CAAC,EAAE,qBAAqB;QAC3B,CAAC,EAAE,MAAM;QACT,QAAQ,EAAE,eAAe,CAAC,kDAAkD;KAC7E;CACF,CAAC;AAEF,QAAQ,CAAC,+CAA+C,EAAE,GAAG,EAAE;IAC7D,KAAK,MAAM,EAAE,IAAI,EAAE,CAAC,EAAE,CAAC,EAAE,QAAQ,EAAE,IAAI,SAAS,EAAE,CAAC;QACjD,IAAI,CAAC,IAAI,EAAE,GAAG,EAAE;YACd,MAAM,MAAM,GAAG,2BAA2B,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;YAEjD,iBAAiB,CAAC,MAAM,EAAE,QAAQ,CAAC,CAAC;QACtC,CAAC,CAAC,CAAC;IACL,CAAC;AACH,CAAC,CAAC,CAAC;AAEH,QAAQ,CAAC,sDAAsD,EAAE,GAAG,EAAE;IACpE,IAAI,CAAC,8BAA8B,EAAE,GAAG,EAAE;QACxC,MAAM,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,wBAAwB;QAC/C,MAAM,CAAC,GAAG,IAAI,CAAC,CAAC,oBAAoB;QACpC,MAAM,MAAM,GAAG,2BAA2B,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;QAEjD,iBAAiB,CAAC,MAAM,EAAE,kBAAkB,CAAC,CAAC;IAChD,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,0BAA0B,EAAE,GAAG,EAAE;QACpC,MAAM,CAAC,GAAG,IAAI,CAAC;QACf,MAAM,CAAC,GAAG,IAAI,CAAC;QACf,MAAM,MAAM,GAAG,2BAA2B,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;QAEjD,iBAAiB,CAAC,MAAM,EAAE,kBAAkB,CAAC,CAAC;IAChD,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,8BAA8B,EAAE,GAAG,EAAE;QACxC,MAAM,CAAC,GAAG,CAAC,CAAC;QACZ,MAAM,CAAC,GAAG,IAAI,CAAC;QACf,MAAM,MAAM,GAAG,2BAA2B,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;QAEjD,iBAAiB,CAAC,MAAM,EAAE,CAAC,CAAC,CAAC,CAAC,wCAAwC;IACxE,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,wBAAwB,EAAE,GAAG,EAAE;QAClC,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,qBAAqB;QAClC,MAAM,CAAC,GAAG,IAAI,CAAC;QACf,MAAM,MAAM,GAAG,2BAA2B,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;QAEjD,iBAAiB,CAAC,MAAM,EAAE,iBAAiB,CAAC,CAAC;IAC/C,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,8CAA8C,EAAE,GAAG,EAAE;QACxD,MAAM,CAAC,GAAG,EAAE,GAAG,IAAI,CAAC,EAAE,CAAC,CAAC,gBAAgB;QACxC,MAAM,CAAC,GAAG,IAAI,CAAC;QACf,MAAM,MAAM,GAAG,2BAA2B,CACxC,CAAC,EACD,CAAC,EACD,IAAI,CAAC,GAAG,CAAC,EAAE,EAAE,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAC5B,CAAC;QAEF,iBAAiB,CAAC,MAAM,EAAE,CAAC,CAAC,CAAC;IAC/B,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,gCAAgC,EAAE,GAAG,EAAE;QAC1C,MAAM,CAAC,GAAG,CAAC,CAAC;QACZ,MAAM,CAAC,GAAG,KAAK,CAAC,CAAC,8BAA8B;QAC/C,MAAM,MAAM,GAAG,2BAA2B,CACxC,CAAC,EACD,CAAC,EACD,IAAI,CAAC,GAAG,CAAC,EAAE,EAAE,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,EAC3B,IAAI,CACL,CAAC;QAEF,MAAM,CAAC,EAAE,CAAC,MAAM,KAAK,IAAI,CAAC,CAAC,CAAC,0CAA0C;IACxE,CAAC,CAAC,CAAC;AACL,CAAC,CAAC,CAAC"}

package/dist/__tests__/solve-kepler-newton-raphson.spec.d.ts

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

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

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

package/dist/__tests__/solve-kepler-newton-raphson.spec.js

@@ -0,0 +1,50 @@
+import assert from 'node:assert/strict';
+import test, { describe } from 'node:test';
+import { solveKeplerNewtonRaphson } from '../solve-kepler-newton-raphson';
+const EPSILON = 1e-9;
+const assertApproxEqual = (actual, expected) => {
+    assert.ok(Math.abs(actual - expected) < EPSILON, `Expected ~${expected}, got ${actual}`);
+};
+describe('solveKeplerNewtonRaphson', () => {
+    test('converges for typical values', () => {
+        const M = Math.PI / 4; // 45 degrees in radians
+        const e = 0.1; // Low eccentricity
+        const result = solveKeplerNewtonRaphson(M, e, 50, 1e-9);
+        assertApproxEqual(result, 0.8612648849); // Expected numerical result
+    });
+    test('works for small M values', () => {
+        const M = 0.01;
+        const e = 0.2;
+        const result = solveKeplerNewtonRaphson(M, e, 50, 1e-9);
+        assertApproxEqual(result, 0.0124999186);
+    });
+    test('handles M = 0 case correctly', () => {
+        const M = 0;
+        const e = 0.5;
+        const result = solveKeplerNewtonRaphson(M, e, 50, 1e-9);
+        assertApproxEqual(result, 0); // Eccentric anomaly should also be zero
+    });
+    test('handles high eccentricity values correctly', () => {
+        const M = Math.PI / 3;
+        const e = 0.95;
+        const result = solveKeplerNewtonRaphson(M, e, 50, 1e-9);
+        assertApproxEqual(result, 1.9349139832);
+    });
+    test('ensures correct output for nearly parabolic orbits', () => {
+        const M = Math.PI / 6;
+        const e = 0.97; // Nearly parabolic
+        const result = solveKeplerNewtonRaphson(M, e, 50, 1e-9);
+        assertApproxEqual(result, 1.4904711747);
+    });
+    test('throws RangeError for invalid eccentricities', () => {
+        assert.throws(() => solveKeplerNewtonRaphson(Math.PI / 4, -0.1, 50, 1e-9), RangeError);
+        assert.throws(() => solveKeplerNewtonRaphson(Math.PI / 4, 1.2, 50, 1e-9), RangeError);
+    });
+    test('returns NaN if the method does not converge', () => {
+        const M = 3;
+        const e = 0.9999; // Extremely high eccentricity
+        const result = solveKeplerNewtonRaphson(M, e, 5, 1e-9); // Very low maxIter
+        assert.ok(Number.isNaN(result));
+    });
+});
+//# sourceMappingURL=solve-kepler-newton-raphson.spec.js.map

package/dist/__tests__/solve-kepler-newton-raphson.spec.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"solve-kepler-newton-raphson.spec.js","sourceRoot":"","sources":["../../src/__tests__/solve-kepler-newton-raphson.spec.ts"],"names":[],"mappings":"AAAA,OAAO,MAAM,MAAM,oBAAoB,CAAC;AACxC,OAAO,IAAI,EAAE,EAAE,QAAQ,EAAE,MAAM,WAAW,CAAC;AAE3C,OAAO,EAAE,wBAAwB,EAAE,MAAM,gCAAgC,CAAC;AAE1E,MAAM,OAAO,GAAG,IAAI,CAAC;AACrB,MAAM,iBAAiB,GAAG,CAAC,MAAc,EAAE,QAAgB,EAAE,EAAE;IAC7D,MAAM,CAAC,EAAE,CACP,IAAI,CAAC,GAAG,CAAC,MAAM,GAAG,QAAQ,CAAC,GAAG,OAAO,EACrC,aAAa,QAAQ,SAAS,MAAM,EAAE,CACvC,CAAC;AACJ,CAAC,CAAC;AAEF,QAAQ,CAAC,0BAA0B,EAAE,GAAG,EAAE;IACxC,IAAI,CAAC,8BAA8B,EAAE,GAAG,EAAE;QACxC,MAAM,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,wBAAwB;QAC/C,MAAM,CAAC,GAAG,GAAG,CAAC,CAAC,mBAAmB;QAClC,MAAM,MAAM,GAAG,wBAAwB,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,EAAE,IAAI,CAAC,CAAC;QAExD,iBAAiB,CAAC,MAAM,EAAE,YAAY,CAAC,CAAC,CAAC,4BAA4B;IACvE,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,0BAA0B,EAAE,GAAG,EAAE;QACpC,MAAM,CAAC,GAAG,IAAI,CAAC;QACf,MAAM,CAAC,GAAG,GAAG,CAAC;QACd,MAAM,MAAM,GAAG,wBAAwB,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,EAAE,IAAI,CAAC,CAAC;QAExD,iBAAiB,CAAC,MAAM,EAAE,YAAY,CAAC,CAAC;IAC1C,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,8BAA8B,EAAE,GAAG,EAAE;QACxC,MAAM,CAAC,GAAG,CAAC,CAAC;QACZ,MAAM,CAAC,GAAG,GAAG,CAAC;QACd,MAAM,MAAM,GAAG,wBAAwB,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,EAAE,IAAI,CAAC,CAAC;QAExD,iBAAiB,CAAC,MAAM,EAAE,CAAC,CAAC,CAAC,CAAC,wCAAwC;IACxE,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,4CAA4C,EAAE,GAAG,EAAE;QACtD,MAAM,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC;QACtB,MAAM,CAAC,GAAG,IAAI,CAAC;QACf,MAAM,MAAM,GAAG,wBAAwB,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,EAAE,IAAI,CAAC,CAAC;QAExD,iBAAiB,CAAC,MAAM,EAAE,YAAY,CAAC,CAAC;IAC1C,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,oDAAoD,EAAE,GAAG,EAAE;QAC9D,MAAM,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC;QACtB,MAAM,CAAC,GAAG,IAAI,CAAC,CAAC,mBAAmB;QACnC,MAAM,MAAM,GAAG,wBAAwB,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,EAAE,IAAI,CAAC,CAAC;QAExD,iBAAiB,CAAC,MAAM,EAAE,YAAY,CAAC,CAAC;IAC1C,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,8CAA8C,EAAE,GAAG,EAAE;QACxD,MAAM,CAAC,MAAM,CACX,GAAG,EAAE,CAAC,wBAAwB,CAAC,IAAI,CAAC,EAAE,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,EAAE,EAAE,IAAI,CAAC,EAC3D,UAAU,CACX,CAAC;QACF,MAAM,CAAC,MAAM,CACX,GAAG,EAAE,CAAC,wBAAwB,CAAC,IAAI,CAAC,EAAE,GAAG,CAAC,EAAE,GAAG,EAAE,EAAE,EAAE,IAAI,CAAC,EAC1D,UAAU,CACX,CAAC;IACJ,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,6CAA6C,EAAE,GAAG,EAAE;QACvD,MAAM,CAAC,GAAG,CAAC,CAAC;QACZ,MAAM,CAAC,GAAG,MAAM,CAAC,CAAC,8BAA8B;QAChD,MAAM,MAAM,GAAG,wBAAwB,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,IAAI,CAAC,CAAC,CAAC,mBAAmB;QAE3E,MAAM,CAAC,EAAE,CAAC,MAAM,CAAC,KAAK,CAAC,MAAM,CAAC,CAAC,CAAC;IAClC,CAAC,CAAC,CAAC;AACL,CAAC,CAAC,CAAC"}

package/dist/__tests__/solve-kepler.int.spec.d.ts

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

package/dist/__tests__/solve-kepler.int.spec.d.ts.map

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