@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.
Files changed (230) hide show
  1. package/dist/categories/angle/compute-angle.d.ts +1 -1
  2. package/dist/categories/angle/compute-angle.js +4 -4
  3. package/dist/categories/angle/compute-angle.js.map +1 -1
  4. package/dist/categories/anomalies/mean-to-eccentric-anomaly.d.ts +1 -1
  5. package/dist/categories/anomalies/mean-to-eccentric-anomaly.js +4 -4
  6. package/dist/categories/anomalies/mean-to-eccentric-anomaly.js.map +1 -1
  7. package/dist/categories/anomalies/true-to-mean-anomaly.d.ts +1 -1
  8. package/dist/categories/anomalies/true-to-mean-anomaly.js +2 -2
  9. package/dist/categories/anomalies/true-to-mean-anomaly.js.map +1 -1
  10. package/dist/categories/cartography/body-fixed-from-inertial-dcm-iau.js +3 -3
  11. package/dist/categories/cartography/body-fixed-from-inertial-dcm-iau.js.map +1 -1
  12. package/dist/categories/gravity/acceleration-on1-by2.js +1 -1
  13. package/dist/categories/gravity/acceleration-on1-by2.js.map +1 -1
  14. package/dist/categories/gravity/force-on1-by2.js +1 -1
  15. package/dist/categories/gravity/force-on1-by2.js.map +1 -1
  16. package/dist/categories/kepler/solve-kepler-high-eccentricity.d.ts +1 -1
  17. package/dist/categories/kepler/solve-kepler-high-eccentricity.js +2 -2
  18. package/dist/categories/kepler/solve-kepler-high-eccentricity.js.map +1 -1
  19. package/dist/categories/kepler/solve-kepler.d.ts +1 -1
  20. package/dist/categories/kepler/solve-kepler.js +5 -5
  21. package/dist/categories/kepler/solve-kepler.js.map +1 -1
  22. package/dist/index.d.ts +49 -49
  23. package/dist/index.d.ts.map +1 -1
  24. package/dist/index.js +49 -49
  25. package/dist/index.js.map +1 -1
  26. package/package.json +4 -4
  27. package/dist/__tests__/helpers/index.d.ts +0 -9
  28. package/dist/__tests__/helpers/index.d.ts.map +0 -1
  29. package/dist/__tests__/helpers/index.js +0 -45
  30. package/dist/__tests__/helpers/index.js.map +0 -1
  31. package/dist/categories/angle/__tests__/compute-angle.int.spec.d.ts +0 -2
  32. package/dist/categories/angle/__tests__/compute-angle.int.spec.d.ts.map +0 -1
  33. package/dist/categories/angle/__tests__/compute-angle.int.spec.js +0 -62
  34. package/dist/categories/angle/__tests__/compute-angle.int.spec.js.map +0 -1
  35. package/dist/categories/angle/__tests__/wrap-angle.spec.d.ts +0 -2
  36. package/dist/categories/angle/__tests__/wrap-angle.spec.d.ts.map +0 -1
  37. package/dist/categories/angle/__tests__/wrap-angle.spec.js +0 -60
  38. package/dist/categories/angle/__tests__/wrap-angle.spec.js.map +0 -1
  39. package/dist/categories/anomalies/__tests__/eccentric-to-true-anomaly.spec.d.ts +0 -2
  40. package/dist/categories/anomalies/__tests__/eccentric-to-true-anomaly.spec.d.ts.map +0 -1
  41. package/dist/categories/anomalies/__tests__/eccentric-to-true-anomaly.spec.js +0 -33
  42. package/dist/categories/anomalies/__tests__/eccentric-to-true-anomaly.spec.js.map +0 -1
  43. package/dist/categories/anomalies/__tests__/mean-to-eccentric-anomaly.spec.d.ts +0 -2
  44. package/dist/categories/anomalies/__tests__/mean-to-eccentric-anomaly.spec.d.ts.map +0 -1
  45. package/dist/categories/anomalies/__tests__/mean-to-eccentric-anomaly.spec.js +0 -197
  46. package/dist/categories/anomalies/__tests__/mean-to-eccentric-anomaly.spec.js.map +0 -1
  47. package/dist/categories/anomalies/__tests__/true-to-mean-anomaly.spec.d.ts +0 -2
  48. package/dist/categories/anomalies/__tests__/true-to-mean-anomaly.spec.d.ts.map +0 -1
  49. package/dist/categories/anomalies/__tests__/true-to-mean-anomaly.spec.js +0 -46
  50. package/dist/categories/anomalies/__tests__/true-to-mean-anomaly.spec.js.map +0 -1
  51. package/dist/categories/cartography/__tests__/body-fixed-from-inertial-dcm-iau.int.spec.d.ts +0 -2
  52. package/dist/categories/cartography/__tests__/body-fixed-from-inertial-dcm-iau.int.spec.d.ts.map +0 -1
  53. package/dist/categories/cartography/__tests__/body-fixed-from-inertial-dcm-iau.int.spec.js +0 -56
  54. package/dist/categories/cartography/__tests__/body-fixed-from-inertial-dcm-iau.int.spec.js.map +0 -1
  55. package/dist/categories/cartography/__tests__/eccentricity-squared-oblate-spheroid.spec.d.ts +0 -2
  56. package/dist/categories/cartography/__tests__/eccentricity-squared-oblate-spheroid.spec.d.ts.map +0 -1
  57. package/dist/categories/cartography/__tests__/eccentricity-squared-oblate-spheroid.spec.js +0 -60
  58. package/dist/categories/cartography/__tests__/eccentricity-squared-oblate-spheroid.spec.js.map +0 -1
  59. package/dist/categories/cartography/__tests__/flattening-oblate-spheroid.spec.d.ts +0 -2
  60. package/dist/categories/cartography/__tests__/flattening-oblate-spheroid.spec.d.ts.map +0 -1
  61. package/dist/categories/cartography/__tests__/flattening-oblate-spheroid.spec.js +0 -60
  62. package/dist/categories/cartography/__tests__/flattening-oblate-spheroid.spec.js.map +0 -1
  63. package/dist/categories/cartography/__tests__/is-on-triaxial-ellipsoid-surface.spec.d.ts +0 -2
  64. package/dist/categories/cartography/__tests__/is-on-triaxial-ellipsoid-surface.spec.d.ts.map +0 -1
  65. package/dist/categories/cartography/__tests__/is-on-triaxial-ellipsoid-surface.spec.js +0 -73
  66. package/dist/categories/cartography/__tests__/is-on-triaxial-ellipsoid-surface.spec.js.map +0 -1
  67. package/dist/categories/cartography/__tests__/planetocentric-latitude.spec.d.ts +0 -2
  68. package/dist/categories/cartography/__tests__/planetocentric-latitude.spec.d.ts.map +0 -1
  69. package/dist/categories/cartography/__tests__/planetocentric-latitude.spec.js +0 -69
  70. package/dist/categories/cartography/__tests__/planetocentric-latitude.spec.js.map +0 -1
  71. package/dist/categories/cartography/__tests__/planetographic-latitude-oblate.spec.d.ts +0 -2
  72. package/dist/categories/cartography/__tests__/planetographic-latitude-oblate.spec.d.ts.map +0 -1
  73. package/dist/categories/cartography/__tests__/planetographic-latitude-oblate.spec.js +0 -103
  74. package/dist/categories/cartography/__tests__/planetographic-latitude-oblate.spec.js.map +0 -1
  75. package/dist/categories/gravity/__tests__/gravitational-acceleration-on1-by2.int.spec.d.ts +0 -2
  76. package/dist/categories/gravity/__tests__/gravitational-acceleration-on1-by2.int.spec.d.ts.map +0 -1
  77. package/dist/categories/gravity/__tests__/gravitational-acceleration-on1-by2.int.spec.js +0 -87
  78. package/dist/categories/gravity/__tests__/gravitational-acceleration-on1-by2.int.spec.js.map +0 -1
  79. package/dist/categories/gravity/__tests__/gravitational-force-on1-by2.int.spec.d.ts +0 -2
  80. package/dist/categories/gravity/__tests__/gravitational-force-on1-by2.int.spec.d.ts.map +0 -1
  81. package/dist/categories/gravity/__tests__/gravitational-force-on1-by2.int.spec.js +0 -107
  82. package/dist/categories/gravity/__tests__/gravitational-force-on1-by2.int.spec.js.map +0 -1
  83. package/dist/categories/gravity/__tests__/gravitational-force.spec.d.ts +0 -2
  84. package/dist/categories/gravity/__tests__/gravitational-force.spec.d.ts.map +0 -1
  85. package/dist/categories/gravity/__tests__/gravitational-force.spec.js +0 -106
  86. package/dist/categories/gravity/__tests__/gravitational-force.spec.js.map +0 -1
  87. package/dist/categories/gravity/__tests__/gravitational-parameters.spec.d.ts +0 -2
  88. package/dist/categories/gravity/__tests__/gravitational-parameters.spec.d.ts.map +0 -1
  89. package/dist/categories/gravity/__tests__/gravitational-parameters.spec.js +0 -75
  90. package/dist/categories/gravity/__tests__/gravitational-parameters.spec.js.map +0 -1
  91. package/dist/categories/helpers/__tests__/apply-matrix-3.spec.d.ts +0 -2
  92. package/dist/categories/helpers/__tests__/apply-matrix-3.spec.d.ts.map +0 -1
  93. package/dist/categories/helpers/__tests__/apply-matrix-3.spec.js +0 -79
  94. package/dist/categories/helpers/__tests__/apply-matrix-3.spec.js.map +0 -1
  95. package/dist/categories/helpers/__tests__/det-3.spec.d.ts +0 -2
  96. package/dist/categories/helpers/__tests__/det-3.spec.d.ts.map +0 -1
  97. package/dist/categories/helpers/__tests__/det-3.spec.js +0 -73
  98. package/dist/categories/helpers/__tests__/det-3.spec.js.map +0 -1
  99. package/dist/categories/helpers/__tests__/mat-mul3.spec.d.ts +0 -2
  100. package/dist/categories/helpers/__tests__/mat-mul3.spec.d.ts.map +0 -1
  101. package/dist/categories/helpers/__tests__/mat-mul3.spec.js +0 -86
  102. package/dist/categories/helpers/__tests__/mat-mul3.spec.js.map +0 -1
  103. package/dist/categories/helpers/__tests__/misc.spec.d.ts +0 -2
  104. package/dist/categories/helpers/__tests__/misc.spec.d.ts.map +0 -1
  105. package/dist/categories/helpers/__tests__/misc.spec.js +0 -101
  106. package/dist/categories/helpers/__tests__/misc.spec.js.map +0 -1
  107. package/dist/categories/helpers/__tests__/radians.spec.d.ts +0 -2
  108. package/dist/categories/helpers/__tests__/radians.spec.d.ts.map +0 -1
  109. package/dist/categories/helpers/__tests__/radians.spec.js +0 -37
  110. package/dist/categories/helpers/__tests__/radians.spec.js.map +0 -1
  111. package/dist/categories/helpers/__tests__/rot-1.spec.d.ts +0 -2
  112. package/dist/categories/helpers/__tests__/rot-1.spec.d.ts.map +0 -1
  113. package/dist/categories/helpers/__tests__/rot-1.spec.js +0 -59
  114. package/dist/categories/helpers/__tests__/rot-1.spec.js.map +0 -1
  115. package/dist/categories/helpers/__tests__/rot-3.spec.d.ts +0 -2
  116. package/dist/categories/helpers/__tests__/rot-3.spec.d.ts.map +0 -1
  117. package/dist/categories/helpers/__tests__/rot-3.spec.js +0 -65
  118. package/dist/categories/helpers/__tests__/rot-3.spec.js.map +0 -1
  119. package/dist/categories/kepler/__tests__/solve-kepler-bisection.spec.d.ts +0 -2
  120. package/dist/categories/kepler/__tests__/solve-kepler-bisection.spec.d.ts.map +0 -1
  121. package/dist/categories/kepler/__tests__/solve-kepler-bisection.spec.js +0 -86
  122. package/dist/categories/kepler/__tests__/solve-kepler-bisection.spec.js.map +0 -1
  123. package/dist/categories/kepler/__tests__/solve-kepler-high-eccentricity.spec.d.ts +0 -2
  124. package/dist/categories/kepler/__tests__/solve-kepler-high-eccentricity.spec.d.ts.map +0 -1
  125. package/dist/categories/kepler/__tests__/solve-kepler-high-eccentricity.spec.js +0 -81
  126. package/dist/categories/kepler/__tests__/solve-kepler-high-eccentricity.spec.js.map +0 -1
  127. package/dist/categories/kepler/__tests__/solve-kepler-newton-raphson.spec.d.ts +0 -2
  128. package/dist/categories/kepler/__tests__/solve-kepler-newton-raphson.spec.d.ts.map +0 -1
  129. package/dist/categories/kepler/__tests__/solve-kepler-newton-raphson.spec.js +0 -50
  130. package/dist/categories/kepler/__tests__/solve-kepler-newton-raphson.spec.js.map +0 -1
  131. package/dist/categories/kepler/__tests__/solve-kepler.int.spec.d.ts +0 -2
  132. package/dist/categories/kepler/__tests__/solve-kepler.int.spec.d.ts.map +0 -1
  133. package/dist/categories/kepler/__tests__/solve-kepler.int.spec.js +0 -76
  134. package/dist/categories/kepler/__tests__/solve-kepler.int.spec.js.map +0 -1
  135. package/dist/categories/manoeuvres/__tests__/combine-burns-delta-v.spec.d.ts +0 -2
  136. package/dist/categories/manoeuvres/__tests__/combine-burns-delta-v.spec.d.ts.map +0 -1
  137. package/dist/categories/manoeuvres/__tests__/combine-burns-delta-v.spec.js +0 -103
  138. package/dist/categories/manoeuvres/__tests__/combine-burns-delta-v.spec.js.map +0 -1
  139. package/dist/categories/manoeuvres/__tests__/gravity-assist-turning-angle.spec.d.ts +0 -2
  140. package/dist/categories/manoeuvres/__tests__/gravity-assist-turning-angle.spec.d.ts.map +0 -1
  141. package/dist/categories/manoeuvres/__tests__/gravity-assist-turning-angle.spec.js +0 -89
  142. package/dist/categories/manoeuvres/__tests__/gravity-assist-turning-angle.spec.js.map +0 -1
  143. package/dist/categories/manoeuvres/__tests__/hohmann-transfer.spec.d.ts +0 -2
  144. package/dist/categories/manoeuvres/__tests__/hohmann-transfer.spec.d.ts.map +0 -1
  145. package/dist/categories/manoeuvres/__tests__/hohmann-transfer.spec.js +0 -101
  146. package/dist/categories/manoeuvres/__tests__/hohmann-transfer.spec.js.map +0 -1
  147. package/dist/categories/manoeuvres/__tests__/oberth-energy-gain.spec.d.ts +0 -2
  148. package/dist/categories/manoeuvres/__tests__/oberth-energy-gain.spec.d.ts.map +0 -1
  149. package/dist/categories/manoeuvres/__tests__/oberth-energy-gain.spec.js +0 -71
  150. package/dist/categories/manoeuvres/__tests__/oberth-energy-gain.spec.js.map +0 -1
  151. package/dist/categories/manoeuvres/__tests__/plane-change-delta-v.spec.d.ts +0 -2
  152. package/dist/categories/manoeuvres/__tests__/plane-change-delta-v.spec.d.ts.map +0 -1
  153. package/dist/categories/manoeuvres/__tests__/plane-change-delta-v.spec.js +0 -90
  154. package/dist/categories/manoeuvres/__tests__/plane-change-delta-v.spec.js.map +0 -1
  155. package/dist/categories/manoeuvres/__tests__/rocket-delta-v-from-isp.spec.d.ts +0 -2
  156. package/dist/categories/manoeuvres/__tests__/rocket-delta-v-from-isp.spec.d.ts.map +0 -1
  157. package/dist/categories/manoeuvres/__tests__/rocket-delta-v-from-isp.spec.js +0 -84
  158. package/dist/categories/manoeuvres/__tests__/rocket-delta-v-from-isp.spec.js.map +0 -1
  159. package/dist/categories/manoeuvres/__tests__/rocket-delta-v-from-ve.spec.d.ts +0 -2
  160. package/dist/categories/manoeuvres/__tests__/rocket-delta-v-from-ve.spec.d.ts.map +0 -1
  161. package/dist/categories/manoeuvres/__tests__/rocket-delta-v-from-ve.spec.js +0 -65
  162. package/dist/categories/manoeuvres/__tests__/rocket-delta-v-from-ve.spec.js.map +0 -1
  163. package/dist/categories/orbits/__tests__/atmospheric-drag-acceleration.spec.d.ts +0 -2
  164. package/dist/categories/orbits/__tests__/atmospheric-drag-acceleration.spec.d.ts.map +0 -1
  165. package/dist/categories/orbits/__tests__/atmospheric-drag-acceleration.spec.js +0 -71
  166. package/dist/categories/orbits/__tests__/atmospheric-drag-acceleration.spec.js.map +0 -1
  167. package/dist/categories/orbits/__tests__/characteristic-energy-c3.spec.d.ts +0 -2
  168. package/dist/categories/orbits/__tests__/characteristic-energy-c3.spec.d.ts.map +0 -1
  169. package/dist/categories/orbits/__tests__/characteristic-energy-c3.spec.js +0 -36
  170. package/dist/categories/orbits/__tests__/characteristic-energy-c3.spec.js.map +0 -1
  171. package/dist/categories/orbits/__tests__/circular-speed.spec.d.ts +0 -2
  172. package/dist/categories/orbits/__tests__/circular-speed.spec.d.ts.map +0 -1
  173. package/dist/categories/orbits/__tests__/circular-speed.spec.js +0 -69
  174. package/dist/categories/orbits/__tests__/circular-speed.spec.js.map +0 -1
  175. package/dist/categories/orbits/__tests__/cw-hill-derivatives.spec.d.ts +0 -2
  176. package/dist/categories/orbits/__tests__/cw-hill-derivatives.spec.d.ts.map +0 -1
  177. package/dist/categories/orbits/__tests__/cw-hill-derivatives.spec.js +0 -64
  178. package/dist/categories/orbits/__tests__/cw-hill-derivatives.spec.js.map +0 -1
  179. package/dist/categories/orbits/__tests__/escape-speed.spec.d.ts +0 -2
  180. package/dist/categories/orbits/__tests__/escape-speed.spec.d.ts.map +0 -1
  181. package/dist/categories/orbits/__tests__/escape-speed.spec.js +0 -79
  182. package/dist/categories/orbits/__tests__/escape-speed.spec.js.map +0 -1
  183. package/dist/categories/orbits/__tests__/flight-path-angle-from-true-anomaly.spec.d.ts +0 -2
  184. package/dist/categories/orbits/__tests__/flight-path-angle-from-true-anomaly.spec.d.ts.map +0 -1
  185. package/dist/categories/orbits/__tests__/flight-path-angle-from-true-anomaly.spec.js +0 -92
  186. package/dist/categories/orbits/__tests__/flight-path-angle-from-true-anomaly.spec.js.map +0 -1
  187. package/dist/categories/orbits/__tests__/hyperbolic-periapsis-speed.spec.d.ts +0 -2
  188. package/dist/categories/orbits/__tests__/hyperbolic-periapsis-speed.spec.d.ts.map +0 -1
  189. package/dist/categories/orbits/__tests__/hyperbolic-periapsis-speed.spec.js +0 -64
  190. package/dist/categories/orbits/__tests__/hyperbolic-periapsis-speed.spec.js.map +0 -1
  191. package/dist/categories/orbits/__tests__/j2-argument-of-perigee-rate.spec.d.ts +0 -2
  192. package/dist/categories/orbits/__tests__/j2-argument-of-perigee-rate.spec.d.ts.map +0 -1
  193. package/dist/categories/orbits/__tests__/j2-argument-of-perigee-rate.spec.js +0 -104
  194. package/dist/categories/orbits/__tests__/j2-argument-of-perigee-rate.spec.js.map +0 -1
  195. package/dist/categories/orbits/__tests__/j2-nodal-precession-rate.spec.d.ts +0 -2
  196. package/dist/categories/orbits/__tests__/j2-nodal-precession-rate.spec.d.ts.map +0 -1
  197. package/dist/categories/orbits/__tests__/j2-nodal-precession-rate.spec.js +0 -99
  198. package/dist/categories/orbits/__tests__/j2-nodal-precession-rate.spec.js.map +0 -1
  199. package/dist/categories/orbits/__tests__/kepler-period.spec.d.ts +0 -2
  200. package/dist/categories/orbits/__tests__/kepler-period.spec.d.ts.map +0 -1
  201. package/dist/categories/orbits/__tests__/kepler-period.spec.js +0 -72
  202. package/dist/categories/orbits/__tests__/kepler-period.spec.js.map +0 -1
  203. package/dist/categories/orbits/__tests__/mean-motion.spec.d.ts +0 -2
  204. package/dist/categories/orbits/__tests__/mean-motion.spec.d.ts.map +0 -1
  205. package/dist/categories/orbits/__tests__/mean-motion.spec.js +0 -48
  206. package/dist/categories/orbits/__tests__/mean-motion.spec.js.map +0 -1
  207. package/dist/categories/orbits/__tests__/peri-apoapsis-radii.spec.d.ts +0 -2
  208. package/dist/categories/orbits/__tests__/peri-apoapsis-radii.spec.d.ts.map +0 -1
  209. package/dist/categories/orbits/__tests__/peri-apoapsis-radii.spec.js +0 -74
  210. package/dist/categories/orbits/__tests__/peri-apoapsis-radii.spec.js.map +0 -1
  211. package/dist/categories/orbits/__tests__/specific-angular-momentum-from-elements.spec.d.ts +0 -2
  212. package/dist/categories/orbits/__tests__/specific-angular-momentum-from-elements.spec.d.ts.map +0 -1
  213. package/dist/categories/orbits/__tests__/specific-angular-momentum-from-elements.spec.js +0 -86
  214. package/dist/categories/orbits/__tests__/specific-angular-momentum-from-elements.spec.js.map +0 -1
  215. package/dist/categories/orbits/__tests__/specific-angular-momentum.spec.d.ts +0 -2
  216. package/dist/categories/orbits/__tests__/specific-angular-momentum.spec.d.ts.map +0 -1
  217. package/dist/categories/orbits/__tests__/specific-angular-momentum.spec.js +0 -85
  218. package/dist/categories/orbits/__tests__/specific-angular-momentum.spec.js.map +0 -1
  219. package/dist/categories/orbits/__tests__/specific-mechanical-energy.spec.d.ts +0 -2
  220. package/dist/categories/orbits/__tests__/specific-mechanical-energy.spec.d.ts.map +0 -1
  221. package/dist/categories/orbits/__tests__/specific-mechanical-energy.spec.js +0 -81
  222. package/dist/categories/orbits/__tests__/specific-mechanical-energy.spec.js.map +0 -1
  223. package/dist/categories/orbits/__tests__/sphere-of-influence-radius.spec.d.ts +0 -2
  224. package/dist/categories/orbits/__tests__/sphere-of-influence-radius.spec.d.ts.map +0 -1
  225. package/dist/categories/orbits/__tests__/sphere-of-influence-radius.spec.js +0 -62
  226. package/dist/categories/orbits/__tests__/sphere-of-influence-radius.spec.js.map +0 -1
  227. package/dist/categories/orbits/__tests__/vis-viva-speed.spec.d.ts +0 -2
  228. package/dist/categories/orbits/__tests__/vis-viva-speed.spec.d.ts.map +0 -1
  229. package/dist/categories/orbits/__tests__/vis-viva-speed.spec.js +0 -80
  230. 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,2 +0,0 @@
1
- export {};
2
- //# sourceMappingURL=wrap-angle.spec.d.ts.map
@@ -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,2 +0,0 @@
1
- export {};
2
- //# sourceMappingURL=eccentric-to-true-anomaly.spec.d.ts.map
@@ -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,2 +0,0 @@
1
- export {};
2
- //# sourceMappingURL=mean-to-eccentric-anomaly.spec.d.ts.map
@@ -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,2 +0,0 @@
1
- export {};
2
- //# sourceMappingURL=true-to-mean-anomaly.spec.d.ts.map
@@ -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"}
@@ -1,2 +0,0 @@
1
- export {};
2
- //# sourceMappingURL=body-fixed-from-inertial-dcm-iau.int.spec.d.ts.map
@@ -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
@@ -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"}
@@ -1,2 +0,0 @@
1
- export {};
2
- //# sourceMappingURL=eccentricity-squared-oblate-spheroid.spec.d.ts.map
@@ -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