@interstellar-tools/equations 0.49.0 → 0.50.0
This diff represents the content of publicly available package versions that have been released to one of the supported registries. The information contained in this diff is provided for informational purposes only and reflects changes between package versions as they appear in their respective public registries.
- package/dist/categories/angle/compute-angle.d.ts +1 -1
- package/dist/categories/angle/compute-angle.js +4 -4
- package/dist/categories/angle/compute-angle.js.map +1 -1
- package/dist/categories/anomalies/mean-to-eccentric-anomaly.d.ts +1 -1
- package/dist/categories/anomalies/mean-to-eccentric-anomaly.js +4 -4
- package/dist/categories/anomalies/mean-to-eccentric-anomaly.js.map +1 -1
- package/dist/categories/anomalies/true-to-mean-anomaly.d.ts +1 -1
- package/dist/categories/anomalies/true-to-mean-anomaly.js +2 -2
- package/dist/categories/anomalies/true-to-mean-anomaly.js.map +1 -1
- package/dist/categories/cartography/body-fixed-from-inertial-dcm-iau.js +3 -3
- package/dist/categories/cartography/body-fixed-from-inertial-dcm-iau.js.map +1 -1
- package/dist/categories/gravity/acceleration-on1-by2.js +1 -1
- package/dist/categories/gravity/acceleration-on1-by2.js.map +1 -1
- package/dist/categories/gravity/force-on1-by2.js +1 -1
- package/dist/categories/gravity/force-on1-by2.js.map +1 -1
- package/dist/categories/kepler/solve-kepler-high-eccentricity.d.ts +1 -1
- package/dist/categories/kepler/solve-kepler-high-eccentricity.js +2 -2
- package/dist/categories/kepler/solve-kepler-high-eccentricity.js.map +1 -1
- package/dist/categories/kepler/solve-kepler.d.ts +1 -1
- package/dist/categories/kepler/solve-kepler.js +5 -5
- package/dist/categories/kepler/solve-kepler.js.map +1 -1
- package/dist/index.d.ts +49 -49
- package/dist/index.d.ts.map +1 -1
- package/dist/index.js +49 -49
- package/dist/index.js.map +1 -1
- package/package.json +4 -4
- package/dist/__tests__/helpers/index.d.ts +0 -9
- package/dist/__tests__/helpers/index.d.ts.map +0 -1
- package/dist/__tests__/helpers/index.js +0 -45
- package/dist/__tests__/helpers/index.js.map +0 -1
- package/dist/categories/angle/__tests__/compute-angle.int.spec.d.ts +0 -2
- package/dist/categories/angle/__tests__/compute-angle.int.spec.d.ts.map +0 -1
- package/dist/categories/angle/__tests__/compute-angle.int.spec.js +0 -62
- package/dist/categories/angle/__tests__/compute-angle.int.spec.js.map +0 -1
- package/dist/categories/angle/__tests__/wrap-angle.spec.d.ts +0 -2
- package/dist/categories/angle/__tests__/wrap-angle.spec.d.ts.map +0 -1
- package/dist/categories/angle/__tests__/wrap-angle.spec.js +0 -60
- package/dist/categories/angle/__tests__/wrap-angle.spec.js.map +0 -1
- package/dist/categories/anomalies/__tests__/eccentric-to-true-anomaly.spec.d.ts +0 -2
- package/dist/categories/anomalies/__tests__/eccentric-to-true-anomaly.spec.d.ts.map +0 -1
- package/dist/categories/anomalies/__tests__/eccentric-to-true-anomaly.spec.js +0 -33
- package/dist/categories/anomalies/__tests__/eccentric-to-true-anomaly.spec.js.map +0 -1
- package/dist/categories/anomalies/__tests__/mean-to-eccentric-anomaly.spec.d.ts +0 -2
- package/dist/categories/anomalies/__tests__/mean-to-eccentric-anomaly.spec.d.ts.map +0 -1
- package/dist/categories/anomalies/__tests__/mean-to-eccentric-anomaly.spec.js +0 -197
- package/dist/categories/anomalies/__tests__/mean-to-eccentric-anomaly.spec.js.map +0 -1
- package/dist/categories/anomalies/__tests__/true-to-mean-anomaly.spec.d.ts +0 -2
- package/dist/categories/anomalies/__tests__/true-to-mean-anomaly.spec.d.ts.map +0 -1
- package/dist/categories/anomalies/__tests__/true-to-mean-anomaly.spec.js +0 -46
- package/dist/categories/anomalies/__tests__/true-to-mean-anomaly.spec.js.map +0 -1
- package/dist/categories/cartography/__tests__/body-fixed-from-inertial-dcm-iau.int.spec.d.ts +0 -2
- package/dist/categories/cartography/__tests__/body-fixed-from-inertial-dcm-iau.int.spec.d.ts.map +0 -1
- package/dist/categories/cartography/__tests__/body-fixed-from-inertial-dcm-iau.int.spec.js +0 -56
- package/dist/categories/cartography/__tests__/body-fixed-from-inertial-dcm-iau.int.spec.js.map +0 -1
- package/dist/categories/cartography/__tests__/eccentricity-squared-oblate-spheroid.spec.d.ts +0 -2
- package/dist/categories/cartography/__tests__/eccentricity-squared-oblate-spheroid.spec.d.ts.map +0 -1
- package/dist/categories/cartography/__tests__/eccentricity-squared-oblate-spheroid.spec.js +0 -60
- package/dist/categories/cartography/__tests__/eccentricity-squared-oblate-spheroid.spec.js.map +0 -1
- package/dist/categories/cartography/__tests__/flattening-oblate-spheroid.spec.d.ts +0 -2
- package/dist/categories/cartography/__tests__/flattening-oblate-spheroid.spec.d.ts.map +0 -1
- package/dist/categories/cartography/__tests__/flattening-oblate-spheroid.spec.js +0 -60
- package/dist/categories/cartography/__tests__/flattening-oblate-spheroid.spec.js.map +0 -1
- package/dist/categories/cartography/__tests__/is-on-triaxial-ellipsoid-surface.spec.d.ts +0 -2
- package/dist/categories/cartography/__tests__/is-on-triaxial-ellipsoid-surface.spec.d.ts.map +0 -1
- package/dist/categories/cartography/__tests__/is-on-triaxial-ellipsoid-surface.spec.js +0 -73
- package/dist/categories/cartography/__tests__/is-on-triaxial-ellipsoid-surface.spec.js.map +0 -1
- package/dist/categories/cartography/__tests__/planetocentric-latitude.spec.d.ts +0 -2
- package/dist/categories/cartography/__tests__/planetocentric-latitude.spec.d.ts.map +0 -1
- package/dist/categories/cartography/__tests__/planetocentric-latitude.spec.js +0 -69
- package/dist/categories/cartography/__tests__/planetocentric-latitude.spec.js.map +0 -1
- package/dist/categories/cartography/__tests__/planetographic-latitude-oblate.spec.d.ts +0 -2
- package/dist/categories/cartography/__tests__/planetographic-latitude-oblate.spec.d.ts.map +0 -1
- package/dist/categories/cartography/__tests__/planetographic-latitude-oblate.spec.js +0 -103
- package/dist/categories/cartography/__tests__/planetographic-latitude-oblate.spec.js.map +0 -1
- package/dist/categories/gravity/__tests__/gravitational-acceleration-on1-by2.int.spec.d.ts +0 -2
- package/dist/categories/gravity/__tests__/gravitational-acceleration-on1-by2.int.spec.d.ts.map +0 -1
- package/dist/categories/gravity/__tests__/gravitational-acceleration-on1-by2.int.spec.js +0 -87
- package/dist/categories/gravity/__tests__/gravitational-acceleration-on1-by2.int.spec.js.map +0 -1
- package/dist/categories/gravity/__tests__/gravitational-force-on1-by2.int.spec.d.ts +0 -2
- package/dist/categories/gravity/__tests__/gravitational-force-on1-by2.int.spec.d.ts.map +0 -1
- package/dist/categories/gravity/__tests__/gravitational-force-on1-by2.int.spec.js +0 -107
- package/dist/categories/gravity/__tests__/gravitational-force-on1-by2.int.spec.js.map +0 -1
- package/dist/categories/gravity/__tests__/gravitational-force.spec.d.ts +0 -2
- package/dist/categories/gravity/__tests__/gravitational-force.spec.d.ts.map +0 -1
- package/dist/categories/gravity/__tests__/gravitational-force.spec.js +0 -106
- package/dist/categories/gravity/__tests__/gravitational-force.spec.js.map +0 -1
- package/dist/categories/gravity/__tests__/gravitational-parameters.spec.d.ts +0 -2
- package/dist/categories/gravity/__tests__/gravitational-parameters.spec.d.ts.map +0 -1
- package/dist/categories/gravity/__tests__/gravitational-parameters.spec.js +0 -75
- package/dist/categories/gravity/__tests__/gravitational-parameters.spec.js.map +0 -1
- package/dist/categories/helpers/__tests__/apply-matrix-3.spec.d.ts +0 -2
- package/dist/categories/helpers/__tests__/apply-matrix-3.spec.d.ts.map +0 -1
- package/dist/categories/helpers/__tests__/apply-matrix-3.spec.js +0 -79
- package/dist/categories/helpers/__tests__/apply-matrix-3.spec.js.map +0 -1
- package/dist/categories/helpers/__tests__/det-3.spec.d.ts +0 -2
- package/dist/categories/helpers/__tests__/det-3.spec.d.ts.map +0 -1
- package/dist/categories/helpers/__tests__/det-3.spec.js +0 -73
- package/dist/categories/helpers/__tests__/det-3.spec.js.map +0 -1
- package/dist/categories/helpers/__tests__/mat-mul3.spec.d.ts +0 -2
- package/dist/categories/helpers/__tests__/mat-mul3.spec.d.ts.map +0 -1
- package/dist/categories/helpers/__tests__/mat-mul3.spec.js +0 -86
- package/dist/categories/helpers/__tests__/mat-mul3.spec.js.map +0 -1
- package/dist/categories/helpers/__tests__/misc.spec.d.ts +0 -2
- package/dist/categories/helpers/__tests__/misc.spec.d.ts.map +0 -1
- package/dist/categories/helpers/__tests__/misc.spec.js +0 -101
- package/dist/categories/helpers/__tests__/misc.spec.js.map +0 -1
- package/dist/categories/helpers/__tests__/radians.spec.d.ts +0 -2
- package/dist/categories/helpers/__tests__/radians.spec.d.ts.map +0 -1
- package/dist/categories/helpers/__tests__/radians.spec.js +0 -37
- package/dist/categories/helpers/__tests__/radians.spec.js.map +0 -1
- package/dist/categories/helpers/__tests__/rot-1.spec.d.ts +0 -2
- package/dist/categories/helpers/__tests__/rot-1.spec.d.ts.map +0 -1
- package/dist/categories/helpers/__tests__/rot-1.spec.js +0 -59
- package/dist/categories/helpers/__tests__/rot-1.spec.js.map +0 -1
- package/dist/categories/helpers/__tests__/rot-3.spec.d.ts +0 -2
- package/dist/categories/helpers/__tests__/rot-3.spec.d.ts.map +0 -1
- package/dist/categories/helpers/__tests__/rot-3.spec.js +0 -65
- package/dist/categories/helpers/__tests__/rot-3.spec.js.map +0 -1
- package/dist/categories/kepler/__tests__/solve-kepler-bisection.spec.d.ts +0 -2
- package/dist/categories/kepler/__tests__/solve-kepler-bisection.spec.d.ts.map +0 -1
- package/dist/categories/kepler/__tests__/solve-kepler-bisection.spec.js +0 -86
- package/dist/categories/kepler/__tests__/solve-kepler-bisection.spec.js.map +0 -1
- package/dist/categories/kepler/__tests__/solve-kepler-high-eccentricity.spec.d.ts +0 -2
- package/dist/categories/kepler/__tests__/solve-kepler-high-eccentricity.spec.d.ts.map +0 -1
- package/dist/categories/kepler/__tests__/solve-kepler-high-eccentricity.spec.js +0 -81
- package/dist/categories/kepler/__tests__/solve-kepler-high-eccentricity.spec.js.map +0 -1
- package/dist/categories/kepler/__tests__/solve-kepler-newton-raphson.spec.d.ts +0 -2
- package/dist/categories/kepler/__tests__/solve-kepler-newton-raphson.spec.d.ts.map +0 -1
- package/dist/categories/kepler/__tests__/solve-kepler-newton-raphson.spec.js +0 -50
- package/dist/categories/kepler/__tests__/solve-kepler-newton-raphson.spec.js.map +0 -1
- package/dist/categories/kepler/__tests__/solve-kepler.int.spec.d.ts +0 -2
- package/dist/categories/kepler/__tests__/solve-kepler.int.spec.d.ts.map +0 -1
- package/dist/categories/kepler/__tests__/solve-kepler.int.spec.js +0 -76
- package/dist/categories/kepler/__tests__/solve-kepler.int.spec.js.map +0 -1
- package/dist/categories/manoeuvres/__tests__/combine-burns-delta-v.spec.d.ts +0 -2
- package/dist/categories/manoeuvres/__tests__/combine-burns-delta-v.spec.d.ts.map +0 -1
- package/dist/categories/manoeuvres/__tests__/combine-burns-delta-v.spec.js +0 -103
- package/dist/categories/manoeuvres/__tests__/combine-burns-delta-v.spec.js.map +0 -1
- package/dist/categories/manoeuvres/__tests__/gravity-assist-turning-angle.spec.d.ts +0 -2
- package/dist/categories/manoeuvres/__tests__/gravity-assist-turning-angle.spec.d.ts.map +0 -1
- package/dist/categories/manoeuvres/__tests__/gravity-assist-turning-angle.spec.js +0 -89
- package/dist/categories/manoeuvres/__tests__/gravity-assist-turning-angle.spec.js.map +0 -1
- package/dist/categories/manoeuvres/__tests__/hohmann-transfer.spec.d.ts +0 -2
- package/dist/categories/manoeuvres/__tests__/hohmann-transfer.spec.d.ts.map +0 -1
- package/dist/categories/manoeuvres/__tests__/hohmann-transfer.spec.js +0 -101
- package/dist/categories/manoeuvres/__tests__/hohmann-transfer.spec.js.map +0 -1
- package/dist/categories/manoeuvres/__tests__/oberth-energy-gain.spec.d.ts +0 -2
- package/dist/categories/manoeuvres/__tests__/oberth-energy-gain.spec.d.ts.map +0 -1
- package/dist/categories/manoeuvres/__tests__/oberth-energy-gain.spec.js +0 -71
- package/dist/categories/manoeuvres/__tests__/oberth-energy-gain.spec.js.map +0 -1
- package/dist/categories/manoeuvres/__tests__/plane-change-delta-v.spec.d.ts +0 -2
- package/dist/categories/manoeuvres/__tests__/plane-change-delta-v.spec.d.ts.map +0 -1
- package/dist/categories/manoeuvres/__tests__/plane-change-delta-v.spec.js +0 -90
- package/dist/categories/manoeuvres/__tests__/plane-change-delta-v.spec.js.map +0 -1
- package/dist/categories/manoeuvres/__tests__/rocket-delta-v-from-isp.spec.d.ts +0 -2
- package/dist/categories/manoeuvres/__tests__/rocket-delta-v-from-isp.spec.d.ts.map +0 -1
- package/dist/categories/manoeuvres/__tests__/rocket-delta-v-from-isp.spec.js +0 -84
- package/dist/categories/manoeuvres/__tests__/rocket-delta-v-from-isp.spec.js.map +0 -1
- package/dist/categories/manoeuvres/__tests__/rocket-delta-v-from-ve.spec.d.ts +0 -2
- package/dist/categories/manoeuvres/__tests__/rocket-delta-v-from-ve.spec.d.ts.map +0 -1
- package/dist/categories/manoeuvres/__tests__/rocket-delta-v-from-ve.spec.js +0 -65
- package/dist/categories/manoeuvres/__tests__/rocket-delta-v-from-ve.spec.js.map +0 -1
- package/dist/categories/orbits/__tests__/atmospheric-drag-acceleration.spec.d.ts +0 -2
- package/dist/categories/orbits/__tests__/atmospheric-drag-acceleration.spec.d.ts.map +0 -1
- package/dist/categories/orbits/__tests__/atmospheric-drag-acceleration.spec.js +0 -71
- package/dist/categories/orbits/__tests__/atmospheric-drag-acceleration.spec.js.map +0 -1
- package/dist/categories/orbits/__tests__/characteristic-energy-c3.spec.d.ts +0 -2
- package/dist/categories/orbits/__tests__/characteristic-energy-c3.spec.d.ts.map +0 -1
- package/dist/categories/orbits/__tests__/characteristic-energy-c3.spec.js +0 -36
- package/dist/categories/orbits/__tests__/characteristic-energy-c3.spec.js.map +0 -1
- package/dist/categories/orbits/__tests__/circular-speed.spec.d.ts +0 -2
- package/dist/categories/orbits/__tests__/circular-speed.spec.d.ts.map +0 -1
- package/dist/categories/orbits/__tests__/circular-speed.spec.js +0 -69
- package/dist/categories/orbits/__tests__/circular-speed.spec.js.map +0 -1
- package/dist/categories/orbits/__tests__/cw-hill-derivatives.spec.d.ts +0 -2
- package/dist/categories/orbits/__tests__/cw-hill-derivatives.spec.d.ts.map +0 -1
- package/dist/categories/orbits/__tests__/cw-hill-derivatives.spec.js +0 -64
- package/dist/categories/orbits/__tests__/cw-hill-derivatives.spec.js.map +0 -1
- package/dist/categories/orbits/__tests__/escape-speed.spec.d.ts +0 -2
- package/dist/categories/orbits/__tests__/escape-speed.spec.d.ts.map +0 -1
- package/dist/categories/orbits/__tests__/escape-speed.spec.js +0 -79
- package/dist/categories/orbits/__tests__/escape-speed.spec.js.map +0 -1
- package/dist/categories/orbits/__tests__/flight-path-angle-from-true-anomaly.spec.d.ts +0 -2
- package/dist/categories/orbits/__tests__/flight-path-angle-from-true-anomaly.spec.d.ts.map +0 -1
- package/dist/categories/orbits/__tests__/flight-path-angle-from-true-anomaly.spec.js +0 -92
- package/dist/categories/orbits/__tests__/flight-path-angle-from-true-anomaly.spec.js.map +0 -1
- package/dist/categories/orbits/__tests__/hyperbolic-periapsis-speed.spec.d.ts +0 -2
- package/dist/categories/orbits/__tests__/hyperbolic-periapsis-speed.spec.d.ts.map +0 -1
- package/dist/categories/orbits/__tests__/hyperbolic-periapsis-speed.spec.js +0 -64
- package/dist/categories/orbits/__tests__/hyperbolic-periapsis-speed.spec.js.map +0 -1
- package/dist/categories/orbits/__tests__/j2-argument-of-perigee-rate.spec.d.ts +0 -2
- package/dist/categories/orbits/__tests__/j2-argument-of-perigee-rate.spec.d.ts.map +0 -1
- package/dist/categories/orbits/__tests__/j2-argument-of-perigee-rate.spec.js +0 -104
- package/dist/categories/orbits/__tests__/j2-argument-of-perigee-rate.spec.js.map +0 -1
- package/dist/categories/orbits/__tests__/j2-nodal-precession-rate.spec.d.ts +0 -2
- package/dist/categories/orbits/__tests__/j2-nodal-precession-rate.spec.d.ts.map +0 -1
- package/dist/categories/orbits/__tests__/j2-nodal-precession-rate.spec.js +0 -99
- package/dist/categories/orbits/__tests__/j2-nodal-precession-rate.spec.js.map +0 -1
- package/dist/categories/orbits/__tests__/kepler-period.spec.d.ts +0 -2
- package/dist/categories/orbits/__tests__/kepler-period.spec.d.ts.map +0 -1
- package/dist/categories/orbits/__tests__/kepler-period.spec.js +0 -72
- package/dist/categories/orbits/__tests__/kepler-period.spec.js.map +0 -1
- package/dist/categories/orbits/__tests__/mean-motion.spec.d.ts +0 -2
- package/dist/categories/orbits/__tests__/mean-motion.spec.d.ts.map +0 -1
- package/dist/categories/orbits/__tests__/mean-motion.spec.js +0 -48
- package/dist/categories/orbits/__tests__/mean-motion.spec.js.map +0 -1
- package/dist/categories/orbits/__tests__/peri-apoapsis-radii.spec.d.ts +0 -2
- package/dist/categories/orbits/__tests__/peri-apoapsis-radii.spec.d.ts.map +0 -1
- package/dist/categories/orbits/__tests__/peri-apoapsis-radii.spec.js +0 -74
- package/dist/categories/orbits/__tests__/peri-apoapsis-radii.spec.js.map +0 -1
- package/dist/categories/orbits/__tests__/specific-angular-momentum-from-elements.spec.d.ts +0 -2
- package/dist/categories/orbits/__tests__/specific-angular-momentum-from-elements.spec.d.ts.map +0 -1
- package/dist/categories/orbits/__tests__/specific-angular-momentum-from-elements.spec.js +0 -86
- package/dist/categories/orbits/__tests__/specific-angular-momentum-from-elements.spec.js.map +0 -1
- package/dist/categories/orbits/__tests__/specific-angular-momentum.spec.d.ts +0 -2
- package/dist/categories/orbits/__tests__/specific-angular-momentum.spec.d.ts.map +0 -1
- package/dist/categories/orbits/__tests__/specific-angular-momentum.spec.js +0 -85
- package/dist/categories/orbits/__tests__/specific-angular-momentum.spec.js.map +0 -1
- package/dist/categories/orbits/__tests__/specific-mechanical-energy.spec.d.ts +0 -2
- package/dist/categories/orbits/__tests__/specific-mechanical-energy.spec.d.ts.map +0 -1
- package/dist/categories/orbits/__tests__/specific-mechanical-energy.spec.js +0 -81
- package/dist/categories/orbits/__tests__/specific-mechanical-energy.spec.js.map +0 -1
- package/dist/categories/orbits/__tests__/sphere-of-influence-radius.spec.d.ts +0 -2
- package/dist/categories/orbits/__tests__/sphere-of-influence-radius.spec.d.ts.map +0 -1
- package/dist/categories/orbits/__tests__/sphere-of-influence-radius.spec.js +0 -62
- package/dist/categories/orbits/__tests__/sphere-of-influence-radius.spec.js.map +0 -1
- package/dist/categories/orbits/__tests__/vis-viva-speed.spec.d.ts +0 -2
- package/dist/categories/orbits/__tests__/vis-viva-speed.spec.d.ts.map +0 -1
- package/dist/categories/orbits/__tests__/vis-viva-speed.spec.js +0 -80
- package/dist/categories/orbits/__tests__/vis-viva-speed.spec.js.map +0 -1
|
@@ -1,62 +0,0 @@
|
|
|
1
|
-
import assert from 'node:assert/strict';
|
|
2
|
-
import test, { describe } from 'node:test';
|
|
3
|
-
import { computeAngle } from '../compute-angle';
|
|
4
|
-
import { wrapAngle } from '../wrap-angle';
|
|
5
|
-
describe('computeAngle', () => {
|
|
6
|
-
// Test cases
|
|
7
|
-
const testCases = [
|
|
8
|
-
{
|
|
9
|
-
celestialBody: {
|
|
10
|
-
e: 0.0167,
|
|
11
|
-
angle: 0,
|
|
12
|
-
period: { value: 365.25, unit: 'd' }
|
|
13
|
-
},
|
|
14
|
-
timeStep: { value: 1, unit: 'd' }
|
|
15
|
-
},
|
|
16
|
-
{
|
|
17
|
-
celestialBody: {
|
|
18
|
-
e: 0.1,
|
|
19
|
-
angle: Math.PI / 3,
|
|
20
|
-
period: { value: 687, unit: 'd' }
|
|
21
|
-
},
|
|
22
|
-
timeStep: { value: 1, unit: 'd' }
|
|
23
|
-
},
|
|
24
|
-
// Retrograde motion test cases
|
|
25
|
-
{
|
|
26
|
-
celestialBody: {
|
|
27
|
-
e: 0.05,
|
|
28
|
-
angle: Math.PI / 2,
|
|
29
|
-
period: { value: -365.25, unit: 'd' } // Retrograde motion
|
|
30
|
-
},
|
|
31
|
-
timeStep: { value: 1, unit: 'd' }
|
|
32
|
-
},
|
|
33
|
-
{
|
|
34
|
-
celestialBody: {
|
|
35
|
-
e: 0.2,
|
|
36
|
-
angle: -Math.PI / 4,
|
|
37
|
-
period: { value: -4332.59, unit: 'd' } // Retrograde motion
|
|
38
|
-
},
|
|
39
|
-
timeStep: { value: 10, unit: 'd' }
|
|
40
|
-
}
|
|
41
|
-
];
|
|
42
|
-
test('computes correct true anomaly', () => {
|
|
43
|
-
testCases.forEach(({ celestialBody, timeStep }) => {
|
|
44
|
-
const expected = wrapAngle(computeAngle(celestialBody, timeStep));
|
|
45
|
-
const result = wrapAngle(computeAngle(celestialBody, timeStep));
|
|
46
|
-
assert.ok(Math.abs(result - expected) < 1e-4, `Expected angle ≈ ${expected} for body=${JSON.stringify(celestialBody)}, timeStep=${timeStep.value}, but got ${result}`);
|
|
47
|
-
});
|
|
48
|
-
});
|
|
49
|
-
test('correctly adjusts for comets', () => {
|
|
50
|
-
const comet = {
|
|
51
|
-
e: 0.8,
|
|
52
|
-
q: 1.2,
|
|
53
|
-
angle: Math.PI / 4,
|
|
54
|
-
period: { value: 2000, unit: 'd' }
|
|
55
|
-
}; // Adding 'q' property to classify as comet
|
|
56
|
-
const timeStep = { value: 5, unit: 'd' };
|
|
57
|
-
const expected = wrapAngle(computeAngle(comet, timeStep));
|
|
58
|
-
const result = wrapAngle(computeAngle(comet, timeStep));
|
|
59
|
-
assert.ok(Math.abs(result - expected) < 1e-4, `Expected comet angle ≈ ${expected}, but got ${result}`);
|
|
60
|
-
});
|
|
61
|
-
});
|
|
62
|
-
//# sourceMappingURL=compute-angle.int.spec.js.map
|
|
@@ -1 +0,0 @@
|
|
|
1
|
-
{"version":3,"file":"compute-angle.int.spec.js","sourceRoot":"","sources":["../../../../src/categories/angle/__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,GAAG,EAAE;aACrC;YACD,QAAQ,EAAE,EAAE,KAAK,EAAE,CAAC,EAAE,IAAI,EAAE,GAAG,EAAE;SAClC;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,GAAG,EAAE;aAClC;YACD,QAAQ,EAAE,EAAE,KAAK,EAAE,CAAC,EAAE,IAAI,EAAE,GAAG,EAAE;SAClC;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,GAAG,EAAE,CAAC,oBAAoB;aAC3D;YACD,QAAQ,EAAE,EAAE,KAAK,EAAE,CAAC,EAAE,IAAI,EAAE,GAAG,EAAE;SAClC;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,GAAG,EAAE,CAAC,oBAAoB;aAC5D;YACD,QAAQ,EAAE,EAAE,KAAK,EAAE,EAAE,EAAE,IAAI,EAAE,GAAG,EAAE;SACnC;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,GAAG,EAAE;SACnC,CAAC,CAAC,2CAA2C;QAC9C,MAAM,QAAQ,GAAG,EAAE,KAAK,EAAE,CAAC,EAAE,IAAI,EAAE,GAAG,EAAE,CAAC;QACzC,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"}
|
|
@@ -1 +0,0 @@
|
|
|
1
|
-
{"version":3,"file":"wrap-angle.spec.d.ts","sourceRoot":"","sources":["../../../../src/categories/angle/__tests__/wrap-angle.spec.ts"],"names":[],"mappings":""}
|
|
@@ -1,60 +0,0 @@
|
|
|
1
|
-
import assert from 'node:assert';
|
|
2
|
-
import { describe, test } from 'node:test';
|
|
3
|
-
import { wrapAngle } from '../wrap-angle';
|
|
4
|
-
const EPSILON = 1e-10; // Small tolerance for floating-point comparisons
|
|
5
|
-
const assertApproxEqual = (actual, expected) => {
|
|
6
|
-
assert.ok(Math.abs(actual - expected) < EPSILON, `Expected ~${expected}, got ${actual}`);
|
|
7
|
-
};
|
|
8
|
-
describe('wrapAngle', () => {
|
|
9
|
-
test('angles already in range [-2π, 2π]', () => {
|
|
10
|
-
assert.equal(wrapAngle(0), 0);
|
|
11
|
-
assert.equal(wrapAngle(Math.PI), Math.PI);
|
|
12
|
-
assert.equal(wrapAngle(-Math.PI), -Math.PI);
|
|
13
|
-
assert.equal(wrapAngle(1.5 * Math.PI), 1.5 * Math.PI);
|
|
14
|
-
assert.equal(wrapAngle(-1.5 * Math.PI), -1.5 * Math.PI);
|
|
15
|
-
});
|
|
16
|
-
test('negative and retrograde angles', () => {
|
|
17
|
-
assertApproxEqual(wrapAngle(-2 * Math.PI), 0);
|
|
18
|
-
assertApproxEqual(wrapAngle(-3 * Math.PI), -Math.PI);
|
|
19
|
-
assertApproxEqual(wrapAngle(-4 * Math.PI), 0);
|
|
20
|
-
assertApproxEqual(wrapAngle(-5 * Math.PI), -Math.PI);
|
|
21
|
-
assertApproxEqual(wrapAngle(-6 * Math.PI), 0);
|
|
22
|
-
});
|
|
23
|
-
test('angles greater than 2π', () => {
|
|
24
|
-
assertApproxEqual(wrapAngle(2 * Math.PI), 0);
|
|
25
|
-
assertApproxEqual(wrapAngle(3 * Math.PI), Math.PI);
|
|
26
|
-
assertApproxEqual(wrapAngle(4 * Math.PI), 0);
|
|
27
|
-
assertApproxEqual(wrapAngle(5 * Math.PI), Math.PI);
|
|
28
|
-
assertApproxEqual(wrapAngle(6 * Math.PI), 0);
|
|
29
|
-
});
|
|
30
|
-
test('extreme values', () => {
|
|
31
|
-
assertApproxEqual(wrapAngle(100 * Math.PI), 0);
|
|
32
|
-
assertApproxEqual(wrapAngle(-100 * Math.PI), 0);
|
|
33
|
-
assertApproxEqual(wrapAngle(1000), wrapAngle(1000 % (2 * Math.PI)));
|
|
34
|
-
assertApproxEqual(wrapAngle(-1000), wrapAngle(-1000 % (2 * Math.PI)));
|
|
35
|
-
});
|
|
36
|
-
test('boundary values', () => {
|
|
37
|
-
assertApproxEqual(wrapAngle(2 * Math.PI), 0);
|
|
38
|
-
assertApproxEqual(wrapAngle(-2 * Math.PI), 0);
|
|
39
|
-
assertApproxEqual(wrapAngle(2 * Math.PI + EPSILON), EPSILON);
|
|
40
|
-
assertApproxEqual(wrapAngle(-2 * Math.PI - EPSILON), -EPSILON);
|
|
41
|
-
});
|
|
42
|
-
test('floating-point drift prevention', () => {
|
|
43
|
-
assertApproxEqual(wrapAngle(Math.PI + EPSILON / 2), Math.PI + EPSILON / 2);
|
|
44
|
-
assertApproxEqual(wrapAngle(-Math.PI - EPSILON / 2), -Math.PI - EPSILON / 2);
|
|
45
|
-
assertApproxEqual(wrapAngle(2 * Math.PI - EPSILON / 2), -EPSILON / 2);
|
|
46
|
-
assertApproxEqual(wrapAngle(-2 * Math.PI + EPSILON / 2), EPSILON / 2);
|
|
47
|
-
});
|
|
48
|
-
test('large positive and negative angles', () => {
|
|
49
|
-
assertApproxEqual(wrapAngle(1e6 * Math.PI), 0);
|
|
50
|
-
assertApproxEqual(wrapAngle(-1e6 * Math.PI), 0);
|
|
51
|
-
assertApproxEqual(wrapAngle(1e9), wrapAngle(1e9 % (2 * Math.PI)));
|
|
52
|
-
assertApproxEqual(wrapAngle(-1e9), wrapAngle(-1e9 % (2 * Math.PI)));
|
|
53
|
-
});
|
|
54
|
-
test('angles near zero', () => {
|
|
55
|
-
assertApproxEqual(wrapAngle(EPSILON), EPSILON);
|
|
56
|
-
assertApproxEqual(wrapAngle(-EPSILON), -EPSILON);
|
|
57
|
-
assertApproxEqual(wrapAngle(0), 0);
|
|
58
|
-
});
|
|
59
|
-
});
|
|
60
|
-
//# sourceMappingURL=wrap-angle.spec.js.map
|
|
@@ -1 +0,0 @@
|
|
|
1
|
-
{"version":3,"file":"wrap-angle.spec.js","sourceRoot":"","sources":["../../../../src/categories/angle/__tests__/wrap-angle.spec.ts"],"names":[],"mappings":"AAAA,OAAO,MAAM,MAAM,aAAa,CAAC;AACjC,OAAO,EAAE,QAAQ,EAAE,IAAI,EAAE,MAAM,WAAW,CAAC;AAE3C,OAAO,EAAE,SAAS,EAAE,MAAM,eAAe,CAAC;AAE1C,MAAM,OAAO,GAAG,KAAK,CAAC,CAAC,iDAAiD;AACxE,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,WAAW,EAAE,GAAG,EAAE;IACzB,IAAI,CAAC,mCAAmC,EAAE,GAAG,EAAE;QAC7C,MAAM,CAAC,KAAK,CAAC,SAAS,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;QAC9B,MAAM,CAAC,KAAK,CAAC,SAAS,CAAC,IAAI,CAAC,EAAE,CAAC,EAAE,IAAI,CAAC,EAAE,CAAC,CAAC;QAC1C,MAAM,CAAC,KAAK,CAAC,SAAS,CAAC,CAAC,IAAI,CAAC,EAAE,CAAC,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC;QAC5C,MAAM,CAAC,KAAK,CAAC,SAAS,CAAC,GAAG,GAAG,IAAI,CAAC,EAAE,CAAC,EAAE,GAAG,GAAG,IAAI,CAAC,EAAE,CAAC,CAAC;QACtD,MAAM,CAAC,KAAK,CAAC,SAAS,CAAC,CAAC,GAAG,GAAG,IAAI,CAAC,EAAE,CAAC,EAAE,CAAC,GAAG,GAAG,IAAI,CAAC,EAAE,CAAC,CAAC;IAC1D,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,gCAAgC,EAAE,GAAG,EAAE;QAC1C,iBAAiB,CAAC,SAAS,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;QAC9C,iBAAiB,CAAC,SAAS,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC;QACrD,iBAAiB,CAAC,SAAS,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;QAC9C,iBAAiB,CAAC,SAAS,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC;QACrD,iBAAiB,CAAC,SAAS,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;IAChD,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,wBAAwB,EAAE,GAAG,EAAE;QAClC,iBAAiB,CAAC,SAAS,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;QAC7C,iBAAiB,CAAC,SAAS,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,EAAE,IAAI,CAAC,EAAE,CAAC,CAAC;QACnD,iBAAiB,CAAC,SAAS,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;QAC7C,iBAAiB,CAAC,SAAS,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,EAAE,IAAI,CAAC,EAAE,CAAC,CAAC;QACnD,iBAAiB,CAAC,SAAS,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;IAC/C,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,gBAAgB,EAAE,GAAG,EAAE;QAC1B,iBAAiB,CAAC,SAAS,CAAC,GAAG,GAAG,IAAI,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;QAC/C,iBAAiB,CAAC,SAAS,CAAC,CAAC,GAAG,GAAG,IAAI,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;QAChD,iBAAiB,CAAC,SAAS,CAAC,IAAI,CAAC,EAAE,SAAS,CAAC,IAAI,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;QACpE,iBAAiB,CAAC,SAAS,CAAC,CAAC,IAAI,CAAC,EAAE,SAAS,CAAC,CAAC,IAAI,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;IACxE,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,iBAAiB,EAAE,GAAG,EAAE;QAC3B,iBAAiB,CAAC,SAAS,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;QAC7C,iBAAiB,CAAC,SAAS,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;QAC9C,iBAAiB,CAAC,SAAS,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,OAAO,CAAC,EAAE,OAAO,CAAC,CAAC;QAC7D,iBAAiB,CAAC,SAAS,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,OAAO,CAAC,EAAE,CAAC,OAAO,CAAC,CAAC;IACjE,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,iCAAiC,EAAE,GAAG,EAAE;QAC3C,iBAAiB,CAAC,SAAS,CAAC,IAAI,CAAC,EAAE,GAAG,OAAO,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,EAAE,GAAG,OAAO,GAAG,CAAC,CAAC,CAAC;QAC3E,iBAAiB,CACf,SAAS,CAAC,CAAC,IAAI,CAAC,EAAE,GAAG,OAAO,GAAG,CAAC,CAAC,EACjC,CAAC,IAAI,CAAC,EAAE,GAAG,OAAO,GAAG,CAAC,CACvB,CAAC;QACF,iBAAiB,CAAC,SAAS,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,OAAO,GAAG,CAAC,CAAC,EAAE,CAAC,OAAO,GAAG,CAAC,CAAC,CAAC;QACtE,iBAAiB,CAAC,SAAS,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,OAAO,GAAG,CAAC,CAAC,EAAE,OAAO,GAAG,CAAC,CAAC,CAAC;IACxE,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,oCAAoC,EAAE,GAAG,EAAE;QAC9C,iBAAiB,CAAC,SAAS,CAAC,GAAG,GAAG,IAAI,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;QAC/C,iBAAiB,CAAC,SAAS,CAAC,CAAC,GAAG,GAAG,IAAI,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;QAChD,iBAAiB,CAAC,SAAS,CAAC,GAAG,CAAC,EAAE,SAAS,CAAC,GAAG,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;QAClE,iBAAiB,CAAC,SAAS,CAAC,CAAC,GAAG,CAAC,EAAE,SAAS,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;IACtE,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,kBAAkB,EAAE,GAAG,EAAE;QAC5B,iBAAiB,CAAC,SAAS,CAAC,OAAO,CAAC,EAAE,OAAO,CAAC,CAAC;QAC/C,iBAAiB,CAAC,SAAS,CAAC,CAAC,OAAO,CAAC,EAAE,CAAC,OAAO,CAAC,CAAC;QACjD,iBAAiB,CAAC,SAAS,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;IACrC,CAAC,CAAC,CAAC;AACL,CAAC,CAAC,CAAC"}
|
|
@@ -1 +0,0 @@
|
|
|
1
|
-
{"version":3,"file":"eccentric-to-true-anomaly.spec.d.ts","sourceRoot":"","sources":["../../../../src/categories/anomalies/__tests__/eccentric-to-true-anomaly.spec.ts"],"names":[],"mappings":""}
|
|
@@ -1,33 +0,0 @@
|
|
|
1
|
-
import assert from 'node:assert/strict';
|
|
2
|
-
import test, { describe } from 'node:test';
|
|
3
|
-
import { eccentricToTrueAnomaly } from '../eccentric-to-true-anomaly';
|
|
4
|
-
describe('eccentricToTrueAnomaly', () => {
|
|
5
|
-
// Test cases for eccentricToTrueAnomaly
|
|
6
|
-
const testCases = [
|
|
7
|
-
{ E: 0, e: 0.5, expected: 0.0 }, // E = 0 should give V = 0
|
|
8
|
-
{ E: Math.PI / 4, e: 0.1, expected: 0.858858 }, // Approximate expected value
|
|
9
|
-
{ E: Math.PI / 2, e: 0.3, expected: 1.875489 }, // Approximate expected value
|
|
10
|
-
{ E: Math.PI, e: 0.2, expected: Math.PI }, // E = π should give V = π
|
|
11
|
-
{ E: (3 * Math.PI) / 2, e: 0.4, expected: 4.300872 } // Approximate expected value
|
|
12
|
-
];
|
|
13
|
-
test('computes correct true anomaly', () => {
|
|
14
|
-
testCases.forEach(({ E, e, expected }) => {
|
|
15
|
-
const result = eccentricToTrueAnomaly(E, e);
|
|
16
|
-
assert.ok(Math.abs(result - expected) < 1e-4, `Expected V ≈ ${expected} for E=${E}, e=${e}, but got ${result}`);
|
|
17
|
-
});
|
|
18
|
-
});
|
|
19
|
-
test('handles special cases', () => {
|
|
20
|
-
assert.strictEqual(eccentricToTrueAnomaly(0, 0), 0); // Circular orbit (e = 0) should return E
|
|
21
|
-
assert.strictEqual(eccentricToTrueAnomaly(Math.PI, 0), Math.PI); // Circular case with E = π
|
|
22
|
-
assert.strictEqual(eccentricToTrueAnomaly(0, 1), 0); // Parabolic orbit case
|
|
23
|
-
});
|
|
24
|
-
test('throws on invalid eccentricity', () => {
|
|
25
|
-
assert.throws(() => eccentricToTrueAnomaly((Math.PI / 3), -0.1), RangeError);
|
|
26
|
-
assert.throws(() => eccentricToTrueAnomaly((Math.PI / 3), 1.2), RangeError);
|
|
27
|
-
});
|
|
28
|
-
test('ensures numerical stability near E = π', () => {
|
|
29
|
-
assert.strictEqual(eccentricToTrueAnomaly((Math.PI + 1e-11), 0.5), Math.PI);
|
|
30
|
-
assert.strictEqual(eccentricToTrueAnomaly((Math.PI - 1e-11), 0.5), Math.PI);
|
|
31
|
-
});
|
|
32
|
-
});
|
|
33
|
-
//# sourceMappingURL=eccentric-to-true-anomaly.spec.js.map
|
|
@@ -1 +0,0 @@
|
|
|
1
|
-
{"version":3,"file":"eccentric-to-true-anomaly.spec.js","sourceRoot":"","sources":["../../../../src/categories/anomalies/__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;AAGtE,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,CAAY,EAAE,CAAC,CAAC,CAAC;YAEvD,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,CAAY,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,yCAAyC;QACzG,MAAM,CAAC,WAAW,CAAC,sBAAsB,CAAC,IAAI,CAAC,EAAa,EAAE,CAAC,CAAC,EAAE,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,2BAA2B;QACvG,MAAM,CAAC,WAAW,CAAC,sBAAsB,CAAC,CAAY,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,uBAAuB;IACzF,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,gCAAgC,EAAE,GAAG,EAAE;QAC1C,MAAM,CAAC,MAAM,CACX,GAAG,EAAE,CAAC,sBAAsB,CAAC,CAAC,IAAI,CAAC,EAAE,GAAG,CAAC,CAAY,EAAE,CAAC,GAAG,CAAC,EAC5D,UAAU,CACX,CAAC;QACF,MAAM,CAAC,MAAM,CACX,GAAG,EAAE,CAAC,sBAAsB,CAAC,CAAC,IAAI,CAAC,EAAE,GAAG,CAAC,CAAY,EAAE,GAAG,CAAC,EAC3D,UAAU,CACX,CAAC;IACJ,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,wCAAwC,EAAE,GAAG,EAAE;QAClD,MAAM,CAAC,WAAW,CAChB,sBAAsB,CAAC,CAAC,IAAI,CAAC,EAAE,GAAG,KAAK,CAAY,EAAE,GAAG,CAAC,EACzD,IAAI,CAAC,EAAE,CACR,CAAC;QACF,MAAM,CAAC,WAAW,CAChB,sBAAsB,CAAC,CAAC,IAAI,CAAC,EAAE,GAAG,KAAK,CAAY,EAAE,GAAG,CAAC,EACzD,IAAI,CAAC,EAAE,CACR,CAAC;IACJ,CAAC,CAAC,CAAC;AACL,CAAC,CAAC,CAAC"}
|
|
@@ -1 +0,0 @@
|
|
|
1
|
-
{"version":3,"file":"mean-to-eccentric-anomaly.spec.d.ts","sourceRoot":"","sources":["../../../../src/categories/anomalies/__tests__/mean-to-eccentric-anomaly.spec.ts"],"names":[],"mappings":""}
|
|
@@ -1,197 +0,0 @@
|
|
|
1
|
-
import assert from 'node:assert/strict';
|
|
2
|
-
import test, { describe } from 'node:test';
|
|
3
|
-
import { meanToEccentricAnomaly } from '../mean-to-eccentric-anomaly';
|
|
4
|
-
import { TWO_PI } from '@interstellar-tools/constants';
|
|
5
|
-
import { trueToMeanAnomaly } from '../true-to-mean-anomaly';
|
|
6
|
-
import { wrapAngle } from '../../angle/wrap-angle';
|
|
7
|
-
const EPSILON = 1e-10; // Increased tolerance for floating-point comparisons
|
|
8
|
-
/* node:coverage disable */
|
|
9
|
-
const assertApproxEqual = (actual, expected) => {
|
|
10
|
-
assert.ok(Math.abs(actual - expected) < EPSILON, `Expected ~${expected}, got ${actual}`);
|
|
11
|
-
};
|
|
12
|
-
/* node:coverage enable */
|
|
13
|
-
describe('meanToEccentricAnomaly', () => {
|
|
14
|
-
test('Circular Orbit (e=0)', () => {
|
|
15
|
-
const body = {
|
|
16
|
-
e: 0,
|
|
17
|
-
angle: Math.PI / 2,
|
|
18
|
-
period: { value: 365, unit: 'd' }
|
|
19
|
-
};
|
|
20
|
-
const timeStep = { value: 100, unit: 'd' };
|
|
21
|
-
const periodInDays = body.period.value;
|
|
22
|
-
const meanMotion = TWO_PI / Math.abs(periodInDays);
|
|
23
|
-
/* node:coverage disable */
|
|
24
|
-
const expectedM = wrapAngle(Math.PI / 2 + meanMotion * 100);
|
|
25
|
-
/* node:coverage enable */
|
|
26
|
-
const result = meanToEccentricAnomaly(body, timeStep);
|
|
27
|
-
assertApproxEqual(result, expectedM);
|
|
28
|
-
});
|
|
29
|
-
test('Elliptical Orbit (e > 0)', () => {
|
|
30
|
-
const body = {
|
|
31
|
-
e: 0.5,
|
|
32
|
-
angle: Math.PI / 2,
|
|
33
|
-
period: { value: 200, unit: 'd' }
|
|
34
|
-
};
|
|
35
|
-
const timeStep = { value: 50, unit: 'd' };
|
|
36
|
-
const periodInDays = body.period.value;
|
|
37
|
-
const meanMotion = TWO_PI / Math.abs(periodInDays);
|
|
38
|
-
const M0 = body.e === 0 ? body.angle : trueToMeanAnomaly(body.angle, body.e);
|
|
39
|
-
/* node:coverage disable */
|
|
40
|
-
const expectedM = wrapAngle(M0 + meanMotion * 50);
|
|
41
|
-
/* node:coverage enable */
|
|
42
|
-
const result = meanToEccentricAnomaly(body, timeStep);
|
|
43
|
-
assertApproxEqual(result, expectedM);
|
|
44
|
-
});
|
|
45
|
-
test('Large Time Step Clamping (Max 10 Orbits)', () => {
|
|
46
|
-
const body = {
|
|
47
|
-
e: 0.1,
|
|
48
|
-
angle: 1.0,
|
|
49
|
-
period: { value: 100, unit: 'd' }
|
|
50
|
-
};
|
|
51
|
-
const timeStep = { value: 2000, unit: 'd' };
|
|
52
|
-
const periodInDays = body.period.value;
|
|
53
|
-
const meanMotion = TWO_PI / Math.abs(periodInDays);
|
|
54
|
-
const maxAllowedTimeStep = 10 * periodInDays;
|
|
55
|
-
const M0 = trueToMeanAnomaly(body.angle, body.e);
|
|
56
|
-
/* node:coverage disable */
|
|
57
|
-
const expectedM = wrapAngle(M0 + meanMotion * maxAllowedTimeStep);
|
|
58
|
-
/* node:coverage enable */
|
|
59
|
-
const result = meanToEccentricAnomaly(body, timeStep);
|
|
60
|
-
assertApproxEqual(result, expectedM);
|
|
61
|
-
});
|
|
62
|
-
test('Eccentricity Out of Bounds Throws Error', () => {
|
|
63
|
-
const invalidBodies = [
|
|
64
|
-
{ e: -0.1, angle: 0, period: { value: 365, unit: 'd' } },
|
|
65
|
-
{ e: 1.1, angle: 0, period: { value: 365, unit: 'd' } }
|
|
66
|
-
];
|
|
67
|
-
const timeStep = { value: 10, unit: 'd' };
|
|
68
|
-
for (const body of invalidBodies) {
|
|
69
|
-
assert.throws(() => meanToEccentricAnomaly(body, timeStep), RangeError);
|
|
70
|
-
}
|
|
71
|
-
});
|
|
72
|
-
test('Zero Time Step Returns Initial Mean Anomaly', () => {
|
|
73
|
-
const body = {
|
|
74
|
-
e: 0.3,
|
|
75
|
-
angle: 2.0,
|
|
76
|
-
period: { value: 400, unit: 'd' }
|
|
77
|
-
};
|
|
78
|
-
const timeStep = { value: 0, unit: 'd' };
|
|
79
|
-
const M0 = body.e === 0 ? body.angle : trueToMeanAnomaly(body.angle, body.e);
|
|
80
|
-
const result = meanToEccentricAnomaly(body, timeStep);
|
|
81
|
-
assertApproxEqual(result, M0);
|
|
82
|
-
});
|
|
83
|
-
test('High Eccentricity Orbit (e ~ 1)', () => {
|
|
84
|
-
const body = {
|
|
85
|
-
e: 0.9951,
|
|
86
|
-
angle: -2.1196,
|
|
87
|
-
period: { value: 365.25, unit: 'd' }
|
|
88
|
-
};
|
|
89
|
-
const timeStep = { value: 0.25, unit: 'd' };
|
|
90
|
-
const periodInDays = body.period.value;
|
|
91
|
-
const meanMotion = TWO_PI / Math.abs(periodInDays);
|
|
92
|
-
const M0 = trueToMeanAnomaly(body.angle, body.e);
|
|
93
|
-
/* node:coverage disable */
|
|
94
|
-
const expectedM = wrapAngle(M0 + meanMotion * 0.25);
|
|
95
|
-
/* node:coverage enable */
|
|
96
|
-
const result = meanToEccentricAnomaly(body, timeStep);
|
|
97
|
-
assertApproxEqual(result, expectedM);
|
|
98
|
-
});
|
|
99
|
-
test('Angle Wrapping within [-2π, 2π]', () => {
|
|
100
|
-
const body = {
|
|
101
|
-
e: 0.5,
|
|
102
|
-
angle: 3 * Math.PI,
|
|
103
|
-
period: { value: 365, unit: 'd' }
|
|
104
|
-
};
|
|
105
|
-
const timeStep = { value: 10, unit: 'd' };
|
|
106
|
-
const periodInDays = body.period.value;
|
|
107
|
-
const meanMotion = TWO_PI / Math.abs(periodInDays);
|
|
108
|
-
const M0 = trueToMeanAnomaly(body.angle, body.e);
|
|
109
|
-
/* node:coverage disable */
|
|
110
|
-
const expectedM = wrapAngle(M0 + meanMotion * 10);
|
|
111
|
-
/* node:coverage enable */
|
|
112
|
-
const result = meanToEccentricAnomaly(body, timeStep);
|
|
113
|
-
assertApproxEqual(result, expectedM);
|
|
114
|
-
assert(result <= 2 * Math.PI && result >= -2 * Math.PI);
|
|
115
|
-
});
|
|
116
|
-
test('Minimum Time Step Clamping', () => {
|
|
117
|
-
const body = {
|
|
118
|
-
e: 0.5,
|
|
119
|
-
angle: Math.PI / 2,
|
|
120
|
-
period: { value: 365, unit: 'd' }
|
|
121
|
-
};
|
|
122
|
-
const timeStep = { value: 1e-6, unit: 'd' };
|
|
123
|
-
const periodInDays = body.period.value;
|
|
124
|
-
const meanMotion = TWO_PI / Math.abs(periodInDays);
|
|
125
|
-
const M0 = trueToMeanAnomaly(body.angle, body.e);
|
|
126
|
-
/* node:coverage disable */
|
|
127
|
-
const expectedM = wrapAngle(M0 + meanMotion * 1e-5); // Minimum time step
|
|
128
|
-
/* node:coverage enable */
|
|
129
|
-
const result = meanToEccentricAnomaly(body, timeStep);
|
|
130
|
-
assertApproxEqual(result, expectedM);
|
|
131
|
-
});
|
|
132
|
-
test("Halley's Comet", () => {
|
|
133
|
-
const body = {
|
|
134
|
-
e: 0.96714,
|
|
135
|
-
angle: (38.38 * Math.PI) / 180, // Convert degrees to radians
|
|
136
|
-
period: { value: 27576, unit: 'd' }
|
|
137
|
-
};
|
|
138
|
-
const timeStep = { value: 100, unit: 'd' };
|
|
139
|
-
const periodInDays = body.period.value;
|
|
140
|
-
const meanMotion = TWO_PI / Math.abs(periodInDays);
|
|
141
|
-
const M0 = trueToMeanAnomaly(body.angle, body.e);
|
|
142
|
-
/* node:coverage disable */
|
|
143
|
-
const expectedM = wrapAngle(M0 + meanMotion * 100);
|
|
144
|
-
/* node:coverage enable */
|
|
145
|
-
const result = meanToEccentricAnomaly(body, timeStep);
|
|
146
|
-
assertApproxEqual(result, expectedM);
|
|
147
|
-
});
|
|
148
|
-
test('Comet Hale-Bopp', () => {
|
|
149
|
-
const body = {
|
|
150
|
-
e: 0.9951,
|
|
151
|
-
angle: (180 * Math.PI) / 180, // Convert degrees to radians
|
|
152
|
-
period: { value: -253533, unit: 'd' } // Positive period for testing
|
|
153
|
-
};
|
|
154
|
-
const timeStep = { value: 0.25, unit: 'd' };
|
|
155
|
-
const periodInDays = body.period.value;
|
|
156
|
-
const meanMotion = TWO_PI / Math.abs(periodInDays);
|
|
157
|
-
const M0 = trueToMeanAnomaly(body.angle, body.e);
|
|
158
|
-
/* node:coverage disable */
|
|
159
|
-
const expectedM = wrapAngle(M0 + meanMotion * 0.25);
|
|
160
|
-
/* node:coverage enable */
|
|
161
|
-
const result = meanToEccentricAnomaly(body, timeStep);
|
|
162
|
-
assertApproxEqual(result, expectedM);
|
|
163
|
-
});
|
|
164
|
-
test('Comet 67P/Churyumov–Gerasimenko', () => {
|
|
165
|
-
const body = {
|
|
166
|
-
e: 0.641,
|
|
167
|
-
angle: 0, // Initial angle in radians
|
|
168
|
-
period: { value: -2484, unit: 'd' } // Positive period for testing
|
|
169
|
-
};
|
|
170
|
-
const timeStep = { value: 10, unit: 'd' };
|
|
171
|
-
const periodInDays = body.period.value;
|
|
172
|
-
const meanMotion = TWO_PI / Math.abs(periodInDays);
|
|
173
|
-
const M0 = trueToMeanAnomaly(body.angle, body.e);
|
|
174
|
-
/* node:coverage disable */
|
|
175
|
-
const expectedM = wrapAngle(M0 + meanMotion * 10);
|
|
176
|
-
/* node:coverage enable */
|
|
177
|
-
const result = meanToEccentricAnomaly(body, timeStep);
|
|
178
|
-
assertApproxEqual(result, expectedM);
|
|
179
|
-
});
|
|
180
|
-
test('Comet Encke', () => {
|
|
181
|
-
const body = {
|
|
182
|
-
e: 0.85,
|
|
183
|
-
angle: (160 * Math.PI) / 180, // Convert degrees to radians
|
|
184
|
-
period: { value: -1204, unit: 'd' } // Positive period for testing
|
|
185
|
-
};
|
|
186
|
-
const timeStep = { value: 5, unit: 'd' };
|
|
187
|
-
const periodInDays = body.period.value;
|
|
188
|
-
const meanMotion = TWO_PI / Math.abs(periodInDays);
|
|
189
|
-
const M0 = trueToMeanAnomaly(body.angle, body.e);
|
|
190
|
-
/* node:coverage disable */
|
|
191
|
-
const expectedM = wrapAngle(M0 + meanMotion * 5);
|
|
192
|
-
/* node:coverage enable */
|
|
193
|
-
const result = meanToEccentricAnomaly(body, timeStep);
|
|
194
|
-
assertApproxEqual(result, expectedM);
|
|
195
|
-
});
|
|
196
|
-
});
|
|
197
|
-
//# sourceMappingURL=mean-to-eccentric-anomaly.spec.js.map
|
|
@@ -1 +0,0 @@
|
|
|
1
|
-
{"version":3,"file":"mean-to-eccentric-anomaly.spec.js","sourceRoot":"","sources":["../../../../src/categories/anomalies/__tests__/mean-to-eccentric-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;AAMtE,OAAO,EAAE,MAAM,EAAE,MAAM,+BAA+B,CAAC;AACvD,OAAO,EAAE,iBAAiB,EAAE,MAAM,yBAAyB,CAAC;AAC5D,OAAO,EAAE,SAAS,EAAE,MAAM,wBAAwB,CAAC;AAEnD,MAAM,OAAO,GAAG,KAAK,CAAC,CAAC,qDAAqD;AAE5E,2BAA2B;AAC3B,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,0BAA0B;AAE1B,QAAQ,CAAC,wBAAwB,EAAE,GAAG,EAAE;IACtC,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,2BAA2B;QAC3B,MAAM,SAAS,GAAG,SAAS,CAAC,IAAI,CAAC,EAAE,GAAG,CAAC,GAAG,UAAU,GAAG,GAAG,CAAC,CAAC;QAC5D,0BAA0B;QAC1B,MAAM,MAAM,GAAG,sBAAsB,CAAC,IAAI,EAAE,QAAQ,CAAC,CAAC;QAEtD,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,iBAAiB,CAAC,IAAI,CAAC,KAAK,EAAE,IAAI,CAAC,CAAC,CAAC,CAAC;QACpE,2BAA2B;QAC3B,MAAM,SAAS,GAAG,SAAS,CAAC,EAAE,GAAG,UAAU,GAAG,EAAE,CAAC,CAAC;QAClD,0BAA0B;QAC1B,MAAM,MAAM,GAAG,sBAAsB,CAAC,IAAI,EAAE,QAAQ,CAAC,CAAC;QAEtD,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,iBAAiB,CAAC,IAAI,CAAC,KAAK,EAAE,IAAI,CAAC,CAAC,CAAC,CAAC;QACjD,2BAA2B;QAC3B,MAAM,SAAS,GAAG,SAAS,CAAC,EAAE,GAAG,UAAU,GAAG,kBAAkB,CAAC,CAAC;QAClE,0BAA0B;QAC1B,MAAM,MAAM,GAAG,sBAAsB,CAAC,IAAI,EAAE,QAAQ,CAAC,CAAC;QAEtD,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,sBAAsB,CAAC,IAAI,EAAE,QAAQ,CAAC,EAAE,UAAU,CAAC,CAAC;QAC1E,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,iBAAiB,CAAC,IAAI,CAAC,KAAK,EAAE,IAAI,CAAC,CAAC,CAAC,CAAC;QACpE,MAAM,MAAM,GAAG,sBAAsB,CAAC,IAAI,EAAE,QAAQ,CAAC,CAAC;QAEtD,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,iBAAiB,CAAC,IAAI,CAAC,KAAK,EAAE,IAAI,CAAC,CAAC,CAAC,CAAC;QACjD,2BAA2B;QAC3B,MAAM,SAAS,GAAG,SAAS,CAAC,EAAE,GAAG,UAAU,GAAG,IAAI,CAAC,CAAC;QACpD,0BAA0B;QAC1B,MAAM,MAAM,GAAG,sBAAsB,CAAC,IAAI,EAAE,QAAQ,CAAC,CAAC;QAEtD,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,iBAAiB,CAAC,IAAI,CAAC,KAAK,EAAE,IAAI,CAAC,CAAC,CAAC,CAAC;QACjD,2BAA2B;QAC3B,MAAM,SAAS,GAAG,SAAS,CAAC,EAAE,GAAG,UAAU,GAAG,EAAE,CAAC,CAAC;QAClD,0BAA0B;QAC1B,MAAM,MAAM,GAAG,sBAAsB,CAAC,IAAI,EAAE,QAAQ,CAAC,CAAC;QAEtD,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,iBAAiB,CAAC,IAAI,CAAC,KAAK,EAAE,IAAI,CAAC,CAAC,CAAC,CAAC;QACjD,2BAA2B;QAC3B,MAAM,SAAS,GAAG,SAAS,CAAC,EAAE,GAAG,UAAU,GAAG,IAAI,CAAC,CAAC,CAAC,oBAAoB;QACzE,0BAA0B;QAC1B,MAAM,MAAM,GAAG,sBAAsB,CAAC,IAAI,EAAE,QAAQ,CAAC,CAAC;QAEtD,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,iBAAiB,CAAC,IAAI,CAAC,KAAK,EAAE,IAAI,CAAC,CAAC,CAAC,CAAC;QACjD,2BAA2B;QAC3B,MAAM,SAAS,GAAG,SAAS,CAAC,EAAE,GAAG,UAAU,GAAG,GAAG,CAAC,CAAC;QACnD,0BAA0B;QAC1B,MAAM,MAAM,GAAG,sBAAsB,CAAC,IAAI,EAAE,QAAQ,CAAC,CAAC;QAEtD,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,iBAAiB,CAAC,IAAI,CAAC,KAAK,EAAE,IAAI,CAAC,CAAC,CAAC,CAAC;QACjD,2BAA2B;QAC3B,MAAM,SAAS,GAAG,SAAS,CAAC,EAAE,GAAG,UAAU,GAAG,IAAI,CAAC,CAAC;QACpD,0BAA0B;QAC1B,MAAM,MAAM,GAAG,sBAAsB,CAAC,IAAI,EAAE,QAAQ,CAAC,CAAC;QAEtD,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,iBAAiB,CAAC,IAAI,CAAC,KAAK,EAAE,IAAI,CAAC,CAAC,CAAC,CAAC;QACjD,2BAA2B;QAC3B,MAAM,SAAS,GAAG,SAAS,CAAC,EAAE,GAAG,UAAU,GAAG,EAAE,CAAC,CAAC;QAClD,0BAA0B;QAC1B,MAAM,MAAM,GAAG,sBAAsB,CAAC,IAAI,EAAE,QAAQ,CAAC,CAAC;QAEtD,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,iBAAiB,CAAC,IAAI,CAAC,KAAK,EAAE,IAAI,CAAC,CAAC,CAAC,CAAC;QACjD,2BAA2B;QAC3B,MAAM,SAAS,GAAG,SAAS,CAAC,EAAE,GAAG,UAAU,GAAG,CAAC,CAAC,CAAC;QACjD,0BAA0B;QAC1B,MAAM,MAAM,GAAG,sBAAsB,CAAC,IAAI,EAAE,QAAQ,CAAC,CAAC;QAEtD,iBAAiB,CAAC,MAAM,EAAE,SAAS,CAAC,CAAC;IACvC,CAAC,CAAC,CAAC;AACL,CAAC,CAAC,CAAC"}
|
|
@@ -1 +0,0 @@
|
|
|
1
|
-
{"version":3,"file":"true-to-mean-anomaly.spec.d.ts","sourceRoot":"","sources":["../../../../src/categories/anomalies/__tests__/true-to-mean-anomaly.spec.ts"],"names":[],"mappings":""}
|
|
@@ -1,46 +0,0 @@
|
|
|
1
|
-
import assert from 'node:assert/strict';
|
|
2
|
-
import test, { describe } from 'node:test';
|
|
3
|
-
import { trueToMeanAnomaly } from '../true-to-mean-anomaly';
|
|
4
|
-
import { wrapAngle } from '../../angle/wrap-angle';
|
|
5
|
-
describe('trueToMeanAnomaly', () => {
|
|
6
|
-
const EPSILON = 1e-10; // Floating-point tolerance
|
|
7
|
-
const assertApproxEqual = (actual, expected) => {
|
|
8
|
-
assert.ok(Math.abs(actual - expected) < EPSILON, `Expected ~${expected}, got ${actual}`);
|
|
9
|
-
};
|
|
10
|
-
test('Circular Orbit (e = 0)', () => {
|
|
11
|
-
const result = trueToMeanAnomaly((Math.PI / 3), 0); // e = 0 (should be identity)
|
|
12
|
-
assertApproxEqual(result, wrapAngle(Math.PI / 3));
|
|
13
|
-
});
|
|
14
|
-
test('Zero True Anomaly (V = 0)', () => {
|
|
15
|
-
const result = trueToMeanAnomaly(0, 0.5);
|
|
16
|
-
assertApproxEqual(result, 0);
|
|
17
|
-
});
|
|
18
|
-
test('True Anomaly at π/2 (90°)', () => {
|
|
19
|
-
const result = trueToMeanAnomaly((Math.PI / 2), 0.5);
|
|
20
|
-
const expectedM = 0.6141848493043778; // Corrected expected value
|
|
21
|
-
assertApproxEqual(result, expectedM);
|
|
22
|
-
});
|
|
23
|
-
test('True Anomaly at π (Apoapsis)', () => {
|
|
24
|
-
const result = trueToMeanAnomaly(Math.PI, 0.5);
|
|
25
|
-
const expectedM = wrapAngle(Math.PI); // Mean anomaly should be equal to π
|
|
26
|
-
assertApproxEqual(result, expectedM);
|
|
27
|
-
});
|
|
28
|
-
test('True Anomaly at 3π/2 (270°)', () => {
|
|
29
|
-
const result = trueToMeanAnomaly(((3 * Math.PI) / 2), 0.5);
|
|
30
|
-
const expectedM = -0.6141848493043782; // Corrected expected value
|
|
31
|
-
assertApproxEqual(result, expectedM);
|
|
32
|
-
});
|
|
33
|
-
test('True Anomaly at 2π (Periapsis)', () => {
|
|
34
|
-
const result = trueToMeanAnomaly((2 * Math.PI), 0.5);
|
|
35
|
-
assertApproxEqual(result, 0); // Should return to zero
|
|
36
|
-
});
|
|
37
|
-
test('Eccentricity Out of Bounds Throws Error', () => {
|
|
38
|
-
assert.throws(() => trueToMeanAnomaly((Math.PI / 3), -0.1), RangeError);
|
|
39
|
-
assert.throws(() => trueToMeanAnomaly((Math.PI / 3), 1.1), RangeError);
|
|
40
|
-
});
|
|
41
|
-
test('Large Eccentricity (e close to 1)', () => {
|
|
42
|
-
const result = trueToMeanAnomaly((Math.PI / 4), 0.99);
|
|
43
|
-
assert.ok(result >= 0 && result < 2 * Math.PI, `Expected result in range [0, 2π], got ${result}`);
|
|
44
|
-
});
|
|
45
|
-
});
|
|
46
|
-
//# sourceMappingURL=true-to-mean-anomaly.spec.js.map
|
|
@@ -1 +0,0 @@
|
|
|
1
|
-
{"version":3,"file":"true-to-mean-anomaly.spec.js","sourceRoot":"","sources":["../../../../src/categories/anomalies/__tests__/true-to-mean-anomaly.spec.ts"],"names":[],"mappings":"AAAA,OAAO,MAAM,MAAM,oBAAoB,CAAC;AACxC,OAAO,IAAI,EAAE,EAAE,QAAQ,EAAE,MAAM,WAAW,CAAC;AAE3C,OAAO,EAAE,iBAAiB,EAAE,MAAM,yBAAyB,CAAC;AAG5D,OAAO,EAAE,SAAS,EAAE,MAAM,wBAAwB,CAAC;AAEnD,QAAQ,CAAC,mBAAmB,EAAE,GAAG,EAAE;IACjC,MAAM,OAAO,GAAG,KAAK,CAAC,CAAC,2BAA2B;IAClD,MAAM,iBAAiB,GAAG,CAAC,MAAc,EAAE,QAAgB,EAAE,EAAE;QAC7D,MAAM,CAAC,EAAE,CACP,IAAI,CAAC,GAAG,CAAC,MAAM,GAAG,QAAQ,CAAC,GAAG,OAAO,EACrC,aAAa,QAAQ,SAAS,MAAM,EAAE,CACvC,CAAC;IACJ,CAAC,CAAC;IAEF,IAAI,CAAC,wBAAwB,EAAE,GAAG,EAAE;QAClC,MAAM,MAAM,GAAG,iBAAiB,CAAC,CAAC,IAAI,CAAC,EAAE,GAAG,CAAC,CAAY,EAAE,CAAC,CAAC,CAAC,CAAC,6BAA6B;QAE5F,iBAAiB,CAAC,MAAM,EAAE,SAAS,CAAC,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC;IACpD,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,2BAA2B,EAAE,GAAG,EAAE;QACrC,MAAM,MAAM,GAAG,iBAAiB,CAAC,CAAY,EAAE,GAAG,CAAC,CAAC;QAEpD,iBAAiB,CAAC,MAAM,EAAE,CAAC,CAAC,CAAC;IAC/B,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,2BAA2B,EAAE,GAAG,EAAE;QACrC,MAAM,MAAM,GAAG,iBAAiB,CAAC,CAAC,IAAI,CAAC,EAAE,GAAG,CAAC,CAAY,EAAE,GAAG,CAAC,CAAC;QAChE,MAAM,SAAS,GAAG,kBAAkB,CAAC,CAAC,2BAA2B;QAEjE,iBAAiB,CAAC,MAAM,EAAE,SAAS,CAAC,CAAC;IACvC,CAAC,CAAC,CAAC;IACH,IAAI,CAAC,8BAA8B,EAAE,GAAG,EAAE;QACxC,MAAM,MAAM,GAAG,iBAAiB,CAAC,IAAI,CAAC,EAAa,EAAE,GAAG,CAAC,CAAC;QAC1D,MAAM,SAAS,GAAG,SAAS,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,oCAAoC;QAE1E,iBAAiB,CAAC,MAAM,EAAE,SAAS,CAAC,CAAC;IACvC,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,6BAA6B,EAAE,GAAG,EAAE;QACvC,MAAM,MAAM,GAAG,iBAAiB,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,GAAG,CAAC,CAAY,EAAE,GAAG,CAAC,CAAC;QACtE,MAAM,SAAS,GAAG,CAAC,kBAAkB,CAAC,CAAC,2BAA2B;QAElE,iBAAiB,CAAC,MAAM,EAAE,SAAS,CAAC,CAAC;IACvC,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,gCAAgC,EAAE,GAAG,EAAE;QAC1C,MAAM,MAAM,GAAG,iBAAiB,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,CAAY,EAAE,GAAG,CAAC,CAAC;QAEhE,iBAAiB,CAAC,MAAM,EAAE,CAAC,CAAC,CAAC,CAAC,wBAAwB;IACxD,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,yCAAyC,EAAE,GAAG,EAAE;QACnD,MAAM,CAAC,MAAM,CACX,GAAG,EAAE,CAAC,iBAAiB,CAAC,CAAC,IAAI,CAAC,EAAE,GAAG,CAAC,CAAY,EAAE,CAAC,GAAG,CAAC,EACvD,UAAU,CACX,CAAC;QACF,MAAM,CAAC,MAAM,CACX,GAAG,EAAE,CAAC,iBAAiB,CAAC,CAAC,IAAI,CAAC,EAAE,GAAG,CAAC,CAAY,EAAE,GAAG,CAAC,EACtD,UAAU,CACX,CAAC;IACJ,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,mCAAmC,EAAE,GAAG,EAAE;QAC7C,MAAM,MAAM,GAAG,iBAAiB,CAAC,CAAC,IAAI,CAAC,EAAE,GAAG,CAAC,CAAY,EAAE,IAAI,CAAC,CAAC;QAEjE,MAAM,CAAC,EAAE,CACP,MAAM,IAAI,CAAC,IAAI,MAAM,GAAG,CAAC,GAAG,IAAI,CAAC,EAAE,EACnC,yCAAyC,MAAM,EAAE,CAClD,CAAC;IACJ,CAAC,CAAC,CAAC;AACL,CAAC,CAAC,CAAC"}
|
package/dist/categories/cartography/__tests__/body-fixed-from-inertial-dcm-iau.int.spec.d.ts.map
DELETED
|
@@ -1 +0,0 @@
|
|
|
1
|
-
{"version":3,"file":"body-fixed-from-inertial-dcm-iau.int.spec.d.ts","sourceRoot":"","sources":["../../../../src/categories/cartography/__tests__/body-fixed-from-inertial-dcm-iau.int.spec.ts"],"names":[],"mappings":""}
|
|
@@ -1,56 +0,0 @@
|
|
|
1
|
-
import assert from 'node:assert/strict';
|
|
2
|
-
import test, { describe } from 'node:test';
|
|
3
|
-
import { absClose, matAbsClose } from 'packages/equations/src/__tests__/helpers';
|
|
4
|
-
import { bodyFixedFromInertialDcmIAU } from '../body-fixed-from-inertial-dcm-iau';
|
|
5
|
-
import { transpose3 } from '../../helpers/transpose-3';
|
|
6
|
-
import { matMul3 } from '../../helpers/mat-mul3';
|
|
7
|
-
import { det3 } from '../../helpers/det-3';
|
|
8
|
-
const I3 = [
|
|
9
|
-
[1, 0, 0],
|
|
10
|
-
[0, 1, 0],
|
|
11
|
-
[0, 0, 1]
|
|
12
|
-
];
|
|
13
|
-
describe('bodyFixedFromInertialDcmIAU', () => {
|
|
14
|
-
test('returns identity for a configuration that collapses to I (alphaP=-π/2, deltaP=π/2, W=0)', () => {
|
|
15
|
-
const alphaP = (-Math.PI / 2);
|
|
16
|
-
const deltaP = (Math.PI / 2);
|
|
17
|
-
const W = 0;
|
|
18
|
-
const R = bodyFixedFromInertialDcmIAU(alphaP, deltaP, W);
|
|
19
|
-
matAbsClose(R, I3, 1e-12);
|
|
20
|
-
});
|
|
21
|
-
test('produces an orthonormal DCM: RᵀR ≈ I', () => {
|
|
22
|
-
const alphaP = 1.234;
|
|
23
|
-
const deltaP = 0.456;
|
|
24
|
-
const W = 2.789;
|
|
25
|
-
const R = bodyFixedFromInertialDcmIAU(alphaP, deltaP, W);
|
|
26
|
-
const Rt = transpose3(R);
|
|
27
|
-
const RtR = matMul3(Rt, R);
|
|
28
|
-
matAbsClose(RtR, I3, 1e-12);
|
|
29
|
-
});
|
|
30
|
-
test('has determinant +1 (proper rotation)', () => {
|
|
31
|
-
const alphaP = 0.12;
|
|
32
|
-
const deltaP = -0.34;
|
|
33
|
-
const W = 0.56;
|
|
34
|
-
const R = bodyFixedFromInertialDcmIAU(alphaP, deltaP, W);
|
|
35
|
-
const d = det3(R);
|
|
36
|
-
absClose(d, 1, 1e-12, `det(R) expected ~1, got ${d}`);
|
|
37
|
-
});
|
|
38
|
-
test('is 2π-periodic in W: R(W + 2π) ≈ R(W)', () => {
|
|
39
|
-
const alphaP = 0.9;
|
|
40
|
-
const deltaP = 0.3;
|
|
41
|
-
const W = 1.1;
|
|
42
|
-
const R1 = bodyFixedFromInertialDcmIAU(alphaP, deltaP, W);
|
|
43
|
-
const R2 = bodyFixedFromInertialDcmIAU(alphaP, deltaP, (W + 2 * Math.PI));
|
|
44
|
-
matAbsClose(R1, R2, 1e-12);
|
|
45
|
-
});
|
|
46
|
-
test('throws TypeError for non-finite inputs', () => {
|
|
47
|
-
const ok = 0;
|
|
48
|
-
assert.throws(() => bodyFixedFromInertialDcmIAU(Number.NaN, ok, ok), (err) => err instanceof TypeError &&
|
|
49
|
-
err.message === 'alphaP must be finite. Received: NaN');
|
|
50
|
-
assert.throws(() => bodyFixedFromInertialDcmIAU(ok, Number.POSITIVE_INFINITY, ok), (err) => err instanceof TypeError &&
|
|
51
|
-
err.message === 'deltaP must be finite. Received: Infinity');
|
|
52
|
-
assert.throws(() => bodyFixedFromInertialDcmIAU(ok, ok, Number.NEGATIVE_INFINITY), (err) => err instanceof TypeError &&
|
|
53
|
-
err.message === 'W must be finite. Received: -Infinity');
|
|
54
|
-
});
|
|
55
|
-
});
|
|
56
|
-
//# sourceMappingURL=body-fixed-from-inertial-dcm-iau.int.spec.js.map
|
package/dist/categories/cartography/__tests__/body-fixed-from-inertial-dcm-iau.int.spec.js.map
DELETED
|
@@ -1 +0,0 @@
|
|
|
1
|
-
{"version":3,"file":"body-fixed-from-inertial-dcm-iau.int.spec.js","sourceRoot":"","sources":["../../../../src/categories/cartography/__tests__/body-fixed-from-inertial-dcm-iau.int.spec.ts"],"names":[],"mappings":"AACA,OAAO,MAAM,MAAM,oBAAoB,CAAC;AACxC,OAAO,IAAI,EAAE,EAAE,QAAQ,EAAE,MAAM,WAAW,CAAC;AAC3C,OAAO,EACL,QAAQ,EACR,WAAW,EACZ,MAAM,0CAA0C,CAAC;AAClD,OAAO,EAAE,2BAA2B,EAAE,MAAM,qCAAqC,CAAC;AAClF,OAAO,EAAE,UAAU,EAAE,MAAM,2BAA2B,CAAC;AACvD,OAAO,EAAE,OAAO,EAAE,MAAM,wBAAwB,CAAC;AACjD,OAAO,EAAE,IAAI,EAAE,MAAM,qBAAqB,CAAC;AAE3C,MAAM,EAAE,GAAkB;IACxB,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;IACT,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;IACT,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;CACD,CAAC;AAEX,QAAQ,CAAC,6BAA6B,EAAE,GAAG,EAAE;IAC3C,IAAI,CAAC,yFAAyF,EAAE,GAAG,EAAE;QACnG,MAAM,MAAM,GAAG,CAAC,CAAC,IAAI,CAAC,EAAE,GAAG,CAAC,CAAY,CAAC;QACzC,MAAM,MAAM,GAAG,CAAC,IAAI,CAAC,EAAE,GAAG,CAAC,CAAY,CAAC;QACxC,MAAM,CAAC,GAAG,CAAY,CAAC;QAEvB,MAAM,CAAC,GAAG,2BAA2B,CAAC,MAAM,EAAE,MAAM,EAAE,CAAC,CAAC,CAAC;QACzD,WAAW,CAAC,CAAC,EAAE,EAAE,EAAE,KAAK,CAAC,CAAC;IAC5B,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,sCAAsC,EAAE,GAAG,EAAE;QAChD,MAAM,MAAM,GAAG,KAAgB,CAAC;QAChC,MAAM,MAAM,GAAG,KAAgB,CAAC;QAChC,MAAM,CAAC,GAAG,KAAgB,CAAC;QAE3B,MAAM,CAAC,GAAG,2BAA2B,CAAC,MAAM,EAAE,MAAM,EAAE,CAAC,CAAC,CAAC;QAEzD,MAAM,EAAE,GAAG,UAAU,CAAC,CAAC,CAAC,CAAC;QACzB,MAAM,GAAG,GAAG,OAAO,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC;QAE3B,WAAW,CAAC,GAAG,EAAE,EAAE,EAAE,KAAK,CAAC,CAAC;IAC9B,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,sCAAsC,EAAE,GAAG,EAAE;QAChD,MAAM,MAAM,GAAG,IAAe,CAAC;QAC/B,MAAM,MAAM,GAAG,CAAC,IAAe,CAAC;QAChC,MAAM,CAAC,GAAG,IAAe,CAAC;QAE1B,MAAM,CAAC,GAAG,2BAA2B,CAAC,MAAM,EAAE,MAAM,EAAE,CAAC,CAAC,CAAC;QACzD,MAAM,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC;QAElB,QAAQ,CAAC,CAAC,EAAE,CAAC,EAAE,KAAK,EAAE,2BAA2B,CAAC,EAAE,CAAC,CAAC;IACxD,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,uCAAuC,EAAE,GAAG,EAAE;QACjD,MAAM,MAAM,GAAG,GAAc,CAAC;QAC9B,MAAM,MAAM,GAAG,GAAc,CAAC;QAC9B,MAAM,CAAC,GAAG,GAAc,CAAC;QAEzB,MAAM,EAAE,GAAG,2BAA2B,CAAC,MAAM,EAAE,MAAM,EAAE,CAAC,CAAC,CAAC;QAC1D,MAAM,EAAE,GAAG,2BAA2B,CACpC,MAAM,EACN,MAAM,EACN,CAAC,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,EAAE,CAAY,CAC7B,CAAC;QAEF,WAAW,CAAC,EAAE,EAAE,EAAE,EAAE,KAAK,CAAC,CAAC;IAC7B,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,wCAAwC,EAAE,GAAG,EAAE;QAClD,MAAM,EAAE,GAAG,CAAY,CAAC;QAExB,MAAM,CAAC,MAAM,CACX,GAAG,EAAE,CAAC,2BAA2B,CAAC,MAAM,CAAC,GAAc,EAAE,EAAE,EAAE,EAAE,CAAC,EAChE,CAAC,GAAG,EAAE,EAAE,CACN,GAAG,YAAY,SAAS;YACxB,GAAG,CAAC,OAAO,KAAK,sCAAsC,CACzD,CAAC;QAEF,MAAM,CAAC,MAAM,CACX,GAAG,EAAE,CACH,2BAA2B,CACzB,EAAE,EACF,MAAM,CAAC,iBAA4B,EACnC,EAAE,CACH,EACH,CAAC,GAAG,EAAE,EAAE,CACN,GAAG,YAAY,SAAS;YACxB,GAAG,CAAC,OAAO,KAAK,2CAA2C,CAC9D,CAAC;QAEF,MAAM,CAAC,MAAM,CACX,GAAG,EAAE,CACH,2BAA2B,CACzB,EAAE,EACF,EAAE,EACF,MAAM,CAAC,iBAA4B,CACpC,EACH,CAAC,GAAG,EAAE,EAAE,CACN,GAAG,YAAY,SAAS;YACxB,GAAG,CAAC,OAAO,KAAK,uCAAuC,CAC1D,CAAC;IACJ,CAAC,CAAC,CAAC;AACL,CAAC,CAAC,CAAC"}
|
package/dist/categories/cartography/__tests__/eccentricity-squared-oblate-spheroid.spec.d.ts.map
DELETED
|
@@ -1 +0,0 @@
|
|
|
1
|
-
{"version":3,"file":"eccentricity-squared-oblate-spheroid.spec.d.ts","sourceRoot":"","sources":["../../../../src/categories/cartography/__tests__/eccentricity-squared-oblate-spheroid.spec.ts"],"names":[],"mappings":""}
|
|
@@ -1,60 +0,0 @@
|
|
|
1
|
-
import assert from 'node:assert/strict';
|
|
2
|
-
import test, { describe } from 'node:test';
|
|
3
|
-
import { eccentricitySquaredOblateSpheroid } from '../eccentricity-squared-oblate-spheroid';
|
|
4
|
-
import { relClose } from 'packages/equations/src/__tests__/helpers';
|
|
5
|
-
describe('eccentricitySquaredOblateSpheroid', () => {
|
|
6
|
-
test('computes e² = 1 - (c²/a²) for typical Earth-like radii (WGS84-ish)', () => {
|
|
7
|
-
const a = 6378137.0; // m (equatorial)
|
|
8
|
-
const c = 6356752.314245; // m (polar)
|
|
9
|
-
const e2 = eccentricitySquaredOblateSpheroid(a, c);
|
|
10
|
-
const expected = 1 - (c * c) / (a * a);
|
|
11
|
-
relClose(e2, expected);
|
|
12
|
-
// sanity: Earth e² is ~0.00669438 (ballpark)
|
|
13
|
-
assert.ok(e2 > 0 && e2 < 0.01, `unexpected e2 magnitude: ${e2}`);
|
|
14
|
-
});
|
|
15
|
-
test('returns 0 for a perfect sphere (a === c)', () => {
|
|
16
|
-
const a = 10;
|
|
17
|
-
const c = 10;
|
|
18
|
-
const e2 = eccentricitySquaredOblateSpheroid(a, c);
|
|
19
|
-
// This should be exactly 0 for identical inputs, but keep it tolerant anyway.
|
|
20
|
-
relClose(e2, 0);
|
|
21
|
-
});
|
|
22
|
-
test('increases as c decreases for fixed a (more flattened => larger e²)', () => {
|
|
23
|
-
const a = 10;
|
|
24
|
-
const e2LessFlat = eccentricitySquaredOblateSpheroid(a, 9.9);
|
|
25
|
-
const e2MoreFlat = eccentricitySquaredOblateSpheroid(a, 9.0);
|
|
26
|
-
assert.ok(e2MoreFlat > e2LessFlat, `expected e2 to increase with flattening: ${e2LessFlat} -> ${e2MoreFlat}`);
|
|
27
|
-
});
|
|
28
|
-
test('is dimensionless and within [0, 1) for physically valid oblate spheroids (a >= c > 0)', () => {
|
|
29
|
-
const a = 10;
|
|
30
|
-
const c = 7;
|
|
31
|
-
const e2 = eccentricitySquaredOblateSpheroid(a, c);
|
|
32
|
-
assert.ok(e2 >= 0, `expected e2 >= 0, got ${e2}`);
|
|
33
|
-
assert.ok(e2 < 1, `expected e2 < 1, got ${e2}`);
|
|
34
|
-
});
|
|
35
|
-
test('throws TypeError when a is not finite (NaN)', () => {
|
|
36
|
-
assert.throws(() => eccentricitySquaredOblateSpheroid(Number.NaN, 1), (err) => err instanceof TypeError &&
|
|
37
|
-
err.message === 'a must be finite. Received: NaN');
|
|
38
|
-
});
|
|
39
|
-
test('throws TypeError when c is not finite (Infinity)', () => {
|
|
40
|
-
assert.throws(() => eccentricitySquaredOblateSpheroid(1, Number.POSITIVE_INFINITY), (err) => err instanceof TypeError &&
|
|
41
|
-
err.message === 'c must be finite. Received: Infinity');
|
|
42
|
-
});
|
|
43
|
-
test('throws RangeError when a <= 0', () => {
|
|
44
|
-
assert.throws(() => eccentricitySquaredOblateSpheroid(0, 1), (err) => err instanceof RangeError &&
|
|
45
|
-
err.message === 'a must be > 0. Received: 0');
|
|
46
|
-
});
|
|
47
|
-
test('throws RangeError when c <= 0', () => {
|
|
48
|
-
assert.throws(() => eccentricitySquaredOblateSpheroid(1, 0), (err) => err instanceof RangeError &&
|
|
49
|
-
err.message === 'c must be > 0. Received: 0');
|
|
50
|
-
});
|
|
51
|
-
test('for prolate inputs (c > a), e² becomes negative (still computed)', () => {
|
|
52
|
-
// The function intentionally doesn't enforce a >= c; it only enforces positivity.
|
|
53
|
-
// This test documents that behavior.
|
|
54
|
-
const a = 10;
|
|
55
|
-
const c = 11;
|
|
56
|
-
const e2 = eccentricitySquaredOblateSpheroid(a, c);
|
|
57
|
-
assert.ok(e2 < 0, `expected e2 < 0 for c > a, got ${e2}`);
|
|
58
|
-
});
|
|
59
|
-
});
|
|
60
|
-
//# sourceMappingURL=eccentricity-squared-oblate-spheroid.spec.js.map
|