@rybosome/tspice 0.0.1
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/LICENSE +21 -0
- package/NOTICE +65 -0
- package/README.md +151 -0
- package/backend-contract/LICENSE +21 -0
- package/backend-contract/dist/.tsbuildinfo +1 -0
- package/backend-contract/dist/index.d.ts +237 -0
- package/backend-contract/dist/index.js +2 -0
- package/backend-fake/dist/.tsbuildinfo +1 -0
- package/backend-fake/dist/index.d.ts +4 -0
- package/backend-fake/dist/index.js +627 -0
- package/backend-node/LICENSE +21 -0
- package/backend-node/NOTICE +32 -0
- package/backend-node/dist/.tsbuildinfo +1 -0
- package/backend-node/dist/index.d.ts +4 -0
- package/backend-node/dist/index.js +389 -0
- package/backend-node/dist/native.d.ts +112 -0
- package/backend-node/dist/native.js +102 -0
- package/backend-wasm/LICENSE +21 -0
- package/backend-wasm/NOTICE +32 -0
- package/backend-wasm/dist/.tsbuildinfo +1 -0
- package/backend-wasm/dist/index.d.ts +8 -0
- package/backend-wasm/dist/index.js +1505 -0
- package/backend-wasm/dist/tspice_backend_wasm.js +19 -0
- package/backend-wasm/dist/tspice_backend_wasm.wasm +0 -0
- package/core/LICENSE +21 -0
- package/core/dist/.tsbuildinfo +1 -0
- package/core/dist/index.d.ts +6 -0
- package/core/dist/index.js +15 -0
- package/dist/.tsbuildinfo +1 -0
- package/dist/backend.d.ts +18 -0
- package/dist/backend.js +37 -0
- package/dist/errors.d.ts +8 -0
- package/dist/errors.js +16 -0
- package/dist/index.d.ts +9 -0
- package/dist/index.js +4 -0
- package/dist/spice-types.d.ts +36 -0
- package/dist/spice-types.js +2 -0
- package/dist/spice.d.ts +13 -0
- package/dist/spice.js +84 -0
- package/dist/types.d.ts +54 -0
- package/dist/types.js +2 -0
- package/package.json +40 -0
|
@@ -0,0 +1,627 @@
|
|
|
1
|
+
/**
|
|
2
|
+
* A deterministic, pure-TS "toy" backend.
|
|
3
|
+
*
|
|
4
|
+
* This backend exists for:
|
|
5
|
+
* - tests
|
|
6
|
+
* - demos (e.g. viewer smoke tests)
|
|
7
|
+
* - environments where native + WASM backends are unavailable
|
|
8
|
+
*
|
|
9
|
+
* It is intentionally **not** a physically-accurate ephemeris.
|
|
10
|
+
*
|
|
11
|
+
* Time conversion notes:
|
|
12
|
+
* - `str2et` supports ISO-8601 / RFC3339-style UTC timestamps only.
|
|
13
|
+
* - Leap seconds are ignored.
|
|
14
|
+
* - The J2000 epoch is treated as `2000-01-01T12:00:00Z`.
|
|
15
|
+
*/
|
|
16
|
+
const TWO_PI = Math.PI * 2;
|
|
17
|
+
export const FAKE_SPICE_VERSION = "tspice-fake-backend@0.0.0";
|
|
18
|
+
const J2000_UTC_MS = Date.parse("2000-01-01T12:00:00.000Z");
|
|
19
|
+
const BODY_IDS = {
|
|
20
|
+
SUN: 10,
|
|
21
|
+
EARTH: 399,
|
|
22
|
+
MOON: 301,
|
|
23
|
+
};
|
|
24
|
+
const BODY_META = [
|
|
25
|
+
{ id: BODY_IDS.SUN, name: "SUN", meanRadiusKm: 695_700 },
|
|
26
|
+
{ id: BODY_IDS.EARTH, name: "EARTH", meanRadiusKm: 6_371 },
|
|
27
|
+
{ id: BODY_IDS.MOON, name: "MOON", meanRadiusKm: 1_737.4 },
|
|
28
|
+
];
|
|
29
|
+
const ID_TO_BODY = new Map(BODY_META.map((b) => [b.id, b]));
|
|
30
|
+
const NAME_TO_ID = new Map(BODY_META.flatMap((b) => [
|
|
31
|
+
[b.name, b.id],
|
|
32
|
+
[b.name.toLowerCase(), b.id],
|
|
33
|
+
]));
|
|
34
|
+
const FRAME_CODES = {
|
|
35
|
+
J2000: 1,
|
|
36
|
+
IAU_EARTH: 10013,
|
|
37
|
+
IAU_MOON: 10020,
|
|
38
|
+
};
|
|
39
|
+
const FRAME_NAME_TO_CODE = new Map([
|
|
40
|
+
["J2000", FRAME_CODES.J2000],
|
|
41
|
+
["j2000", FRAME_CODES.J2000],
|
|
42
|
+
["IAU_EARTH", FRAME_CODES.IAU_EARTH],
|
|
43
|
+
["iau_earth", FRAME_CODES.IAU_EARTH],
|
|
44
|
+
["IAU_MOON", FRAME_CODES.IAU_MOON],
|
|
45
|
+
["iau_moon", FRAME_CODES.IAU_MOON],
|
|
46
|
+
]);
|
|
47
|
+
const FRAME_CODE_TO_NAME = new Map([
|
|
48
|
+
[FRAME_CODES.J2000, "J2000"],
|
|
49
|
+
[FRAME_CODES.IAU_EARTH, "IAU_EARTH"],
|
|
50
|
+
[FRAME_CODES.IAU_MOON, "IAU_MOON"],
|
|
51
|
+
]);
|
|
52
|
+
/**
|
|
53
|
+
* Frame spin rates (rad/s) relative to J2000.
|
|
54
|
+
*
|
|
55
|
+
* Used for deterministic `pxform`/`sxform`. These are not intended to be
|
|
56
|
+
* authoritative values.
|
|
57
|
+
*/
|
|
58
|
+
const FRAME_SPIN_RATE_RAD_PER_SEC = {
|
|
59
|
+
J2000: 0,
|
|
60
|
+
// Approx sidereal rotation (deterministic constant).
|
|
61
|
+
IAU_EARTH: TWO_PI / 86164.0905,
|
|
62
|
+
// Approx synchronous rotation with orbital period (deterministic constant).
|
|
63
|
+
IAU_MOON: TWO_PI / (27.321661 * 24 * 60 * 60),
|
|
64
|
+
};
|
|
65
|
+
function normalizeName(name) {
|
|
66
|
+
return name.trim();
|
|
67
|
+
}
|
|
68
|
+
function parseBodyRef(ref) {
|
|
69
|
+
const trimmed = normalizeName(ref);
|
|
70
|
+
// Accept numeric IDs as strings.
|
|
71
|
+
if (/^-?\d+$/.test(trimmed)) {
|
|
72
|
+
return Number(trimmed);
|
|
73
|
+
}
|
|
74
|
+
const id = NAME_TO_ID.get(trimmed) ?? NAME_TO_ID.get(trimmed.toLowerCase());
|
|
75
|
+
if (id === undefined) {
|
|
76
|
+
throw new Error(`Fake backend: unsupported body: ${JSON.stringify(ref)}`);
|
|
77
|
+
}
|
|
78
|
+
return id;
|
|
79
|
+
}
|
|
80
|
+
function parseFrameName(ref) {
|
|
81
|
+
const trimmed = normalizeName(ref);
|
|
82
|
+
const code = FRAME_NAME_TO_CODE.get(trimmed) ?? FRAME_NAME_TO_CODE.get(trimmed.toLowerCase());
|
|
83
|
+
const name = code === undefined ? undefined : FRAME_CODE_TO_NAME.get(code);
|
|
84
|
+
if (!name) {
|
|
85
|
+
throw new Error(`Fake backend: unsupported frame: ${JSON.stringify(ref)}`);
|
|
86
|
+
}
|
|
87
|
+
return name;
|
|
88
|
+
}
|
|
89
|
+
function getBodyRadiusKm(bodyId) {
|
|
90
|
+
return ID_TO_BODY.get(bodyId)?.meanRadiusKm ?? 1;
|
|
91
|
+
}
|
|
92
|
+
function vadd(a, b) {
|
|
93
|
+
return [a[0] + b[0], a[1] + b[1], a[2] + b[2]];
|
|
94
|
+
}
|
|
95
|
+
function vsub(a, b) {
|
|
96
|
+
return [a[0] - b[0], a[1] - b[1], a[2] - b[2]];
|
|
97
|
+
}
|
|
98
|
+
function vscale(s, v) {
|
|
99
|
+
return [s * v[0], s * v[1], s * v[2]];
|
|
100
|
+
}
|
|
101
|
+
function vdot(a, b) {
|
|
102
|
+
return a[0] * b[0] + a[1] * b[1] + a[2] * b[2];
|
|
103
|
+
}
|
|
104
|
+
function vcrss(a, b) {
|
|
105
|
+
return [
|
|
106
|
+
a[1] * b[2] - a[2] * b[1],
|
|
107
|
+
a[2] * b[0] - a[0] * b[2],
|
|
108
|
+
a[0] * b[1] - a[1] * b[0],
|
|
109
|
+
];
|
|
110
|
+
}
|
|
111
|
+
function vnorm(v) {
|
|
112
|
+
return Math.sqrt(vdot(v, v));
|
|
113
|
+
}
|
|
114
|
+
function vhat(v) {
|
|
115
|
+
const n = vnorm(v);
|
|
116
|
+
if (n === 0)
|
|
117
|
+
return [0, 0, 0];
|
|
118
|
+
return vscale(1 / n, v);
|
|
119
|
+
}
|
|
120
|
+
function clamp(x, lo, hi) {
|
|
121
|
+
return Math.max(lo, Math.min(hi, x));
|
|
122
|
+
}
|
|
123
|
+
function angleBetween(a, b) {
|
|
124
|
+
const na = vnorm(a);
|
|
125
|
+
const nb = vnorm(b);
|
|
126
|
+
if (na === 0 || nb === 0)
|
|
127
|
+
return 0;
|
|
128
|
+
const c = clamp(vdot(a, b) / (na * nb), -1, 1);
|
|
129
|
+
return Math.acos(c);
|
|
130
|
+
}
|
|
131
|
+
function canonicalizeZero(n) {
|
|
132
|
+
return Object.is(n, -0) ? 0 : n;
|
|
133
|
+
}
|
|
134
|
+
function rotZRowMajor(theta) {
|
|
135
|
+
const c = Math.cos(theta);
|
|
136
|
+
const s = Math.sin(theta);
|
|
137
|
+
return [
|
|
138
|
+
canonicalizeZero(c),
|
|
139
|
+
canonicalizeZero(-s),
|
|
140
|
+
0,
|
|
141
|
+
canonicalizeZero(s),
|
|
142
|
+
canonicalizeZero(c),
|
|
143
|
+
0,
|
|
144
|
+
0,
|
|
145
|
+
0,
|
|
146
|
+
1,
|
|
147
|
+
];
|
|
148
|
+
}
|
|
149
|
+
function drotZRowMajor(theta, w) {
|
|
150
|
+
// d/dt rotZ(theta) = w * d/dtheta rotZ(theta)
|
|
151
|
+
const c = Math.cos(theta);
|
|
152
|
+
const s = Math.sin(theta);
|
|
153
|
+
return [
|
|
154
|
+
canonicalizeZero(-w * s),
|
|
155
|
+
canonicalizeZero(-w * c),
|
|
156
|
+
0,
|
|
157
|
+
canonicalizeZero(w * c),
|
|
158
|
+
canonicalizeZero(-w * s),
|
|
159
|
+
0,
|
|
160
|
+
0,
|
|
161
|
+
0,
|
|
162
|
+
0,
|
|
163
|
+
];
|
|
164
|
+
}
|
|
165
|
+
function mmul3(a, b) {
|
|
166
|
+
// Row-major 3x3 multiply: out = a*b
|
|
167
|
+
const out = new Array(9).fill(0);
|
|
168
|
+
for (let r = 0; r < 3; r++) {
|
|
169
|
+
for (let c = 0; c < 3; c++) {
|
|
170
|
+
let sum = 0;
|
|
171
|
+
for (let k = 0; k < 3; k++) {
|
|
172
|
+
sum += a[r * 3 + k] * b[k * 3 + c];
|
|
173
|
+
}
|
|
174
|
+
out[r * 3 + c] = sum;
|
|
175
|
+
}
|
|
176
|
+
}
|
|
177
|
+
return out;
|
|
178
|
+
}
|
|
179
|
+
function mtx3(m) {
|
|
180
|
+
return [
|
|
181
|
+
m[0],
|
|
182
|
+
m[3],
|
|
183
|
+
m[6],
|
|
184
|
+
m[1],
|
|
185
|
+
m[4],
|
|
186
|
+
m[7],
|
|
187
|
+
m[2],
|
|
188
|
+
m[5],
|
|
189
|
+
m[8],
|
|
190
|
+
];
|
|
191
|
+
}
|
|
192
|
+
function mxv(m, v) {
|
|
193
|
+
return [
|
|
194
|
+
m[0] * v[0] + m[1] * v[1] + m[2] * v[2],
|
|
195
|
+
m[3] * v[0] + m[4] * v[1] + m[5] * v[2],
|
|
196
|
+
m[6] * v[0] + m[7] * v[1] + m[8] * v[2],
|
|
197
|
+
];
|
|
198
|
+
}
|
|
199
|
+
function mtxv(m, v) {
|
|
200
|
+
// (m^T) * v
|
|
201
|
+
return mxv(mtx3(m), v);
|
|
202
|
+
}
|
|
203
|
+
function sxformRowMajor(from, to, et) {
|
|
204
|
+
const wFrom = FRAME_SPIN_RATE_RAD_PER_SEC[from];
|
|
205
|
+
const wTo = FRAME_SPIN_RATE_RAD_PER_SEC[to];
|
|
206
|
+
const wDelta = wFrom - wTo;
|
|
207
|
+
const theta = wDelta * et;
|
|
208
|
+
const r = rotZRowMajor(theta);
|
|
209
|
+
const dr = drotZRowMajor(theta, wDelta);
|
|
210
|
+
// [ R 0 ; dR R ]
|
|
211
|
+
return [
|
|
212
|
+
// row 0
|
|
213
|
+
r[0], r[1], r[2], 0, 0, 0,
|
|
214
|
+
// row 1
|
|
215
|
+
r[3], r[4], r[5], 0, 0, 0,
|
|
216
|
+
// row 2
|
|
217
|
+
r[6], r[7], r[8], 0, 0, 0,
|
|
218
|
+
// row 3
|
|
219
|
+
dr[0], dr[1], dr[2], r[0], r[1], r[2],
|
|
220
|
+
// row 4
|
|
221
|
+
dr[3], dr[4], dr[5], r[3], r[4], r[5],
|
|
222
|
+
// row 5
|
|
223
|
+
dr[6], dr[7], dr[8], r[6], r[7], r[8],
|
|
224
|
+
];
|
|
225
|
+
}
|
|
226
|
+
function mxv6(m, v) {
|
|
227
|
+
const out = new Array(6).fill(0);
|
|
228
|
+
for (let r = 0; r < 6; r++) {
|
|
229
|
+
let sum = 0;
|
|
230
|
+
for (let c = 0; c < 6; c++) {
|
|
231
|
+
sum += m[r * 6 + c] * v[c];
|
|
232
|
+
}
|
|
233
|
+
out[r] = sum;
|
|
234
|
+
}
|
|
235
|
+
return out;
|
|
236
|
+
}
|
|
237
|
+
/**
|
|
238
|
+
* Deterministic toy ephemerides in J2000.
|
|
239
|
+
*/
|
|
240
|
+
function getAbsoluteStateInJ2000(bodyId, et) {
|
|
241
|
+
switch (bodyId) {
|
|
242
|
+
case BODY_IDS.SUN:
|
|
243
|
+
return { posKm: [0, 0, 0], velKmPerSec: [0, 0, 0] };
|
|
244
|
+
case BODY_IDS.EARTH: {
|
|
245
|
+
const rKm = 149_597_870.7; // 1 AU
|
|
246
|
+
const periodSec = 365.25 * 24 * 60 * 60;
|
|
247
|
+
const w = TWO_PI / periodSec;
|
|
248
|
+
const t = w * et;
|
|
249
|
+
const c = Math.cos(t);
|
|
250
|
+
const s = Math.sin(t);
|
|
251
|
+
const x = rKm * c;
|
|
252
|
+
const y = rKm * s;
|
|
253
|
+
const vx = -rKm * w * s;
|
|
254
|
+
const vy = rKm * w * c;
|
|
255
|
+
return { posKm: [x, y, 0], velKmPerSec: [vx, vy, 0] };
|
|
256
|
+
}
|
|
257
|
+
case BODY_IDS.MOON: {
|
|
258
|
+
const earth = getAbsoluteStateInJ2000(BODY_IDS.EARTH, et);
|
|
259
|
+
const rKm = 384_400;
|
|
260
|
+
const periodSec = 27.321661 * 24 * 60 * 60;
|
|
261
|
+
const w = TWO_PI / periodSec;
|
|
262
|
+
const t = w * et;
|
|
263
|
+
const c = Math.cos(t);
|
|
264
|
+
const s = Math.sin(t);
|
|
265
|
+
const xRel = rKm * c;
|
|
266
|
+
const yRel = rKm * s;
|
|
267
|
+
const vxRel = -rKm * w * s;
|
|
268
|
+
const vyRel = rKm * w * c;
|
|
269
|
+
return {
|
|
270
|
+
posKm: [earth.posKm[0] + xRel, earth.posKm[1] + yRel, 0],
|
|
271
|
+
velKmPerSec: [earth.velKmPerSec[0] + vxRel, earth.velKmPerSec[1] + vyRel, 0],
|
|
272
|
+
};
|
|
273
|
+
}
|
|
274
|
+
default:
|
|
275
|
+
throw new Error(`Fake backend: unsupported body id: ${bodyId}`);
|
|
276
|
+
}
|
|
277
|
+
}
|
|
278
|
+
function getRelativeStateInJ2000(target, observer, et) {
|
|
279
|
+
const targetId = parseBodyRef(target);
|
|
280
|
+
const observerId = parseBodyRef(observer);
|
|
281
|
+
const t = getAbsoluteStateInJ2000(targetId, et);
|
|
282
|
+
const o = getAbsoluteStateInJ2000(observerId, et);
|
|
283
|
+
const relPos = vsub(t.posKm, o.posKm);
|
|
284
|
+
const relVel = vsub(t.velKmPerSec, o.velKmPerSec);
|
|
285
|
+
return [relPos[0], relPos[1], relPos[2], relVel[0], relVel[1], relVel[2]];
|
|
286
|
+
}
|
|
287
|
+
function applyStateTransform(from, to, et, state) {
|
|
288
|
+
const xform = sxformRowMajor(from, to, et);
|
|
289
|
+
return mxv6(xform, state);
|
|
290
|
+
}
|
|
291
|
+
function isIso8601OrRfc3339Utcish(s) {
|
|
292
|
+
// Intentionally conservative. We only guarantee ISO/RFC3339-style parsing.
|
|
293
|
+
// Examples:
|
|
294
|
+
// - 2000-01-01T12:00:00Z
|
|
295
|
+
// - 2000-01-01T12:00:00.123Z
|
|
296
|
+
// - 2000-01-01T12:00:00+00:00
|
|
297
|
+
return /^\d{4}-\d{2}-\d{2}T\d{2}:\d{2}:\d{2}(\.\d+)?(Z|[+-]\d{2}:\d{2})$/.test(s);
|
|
298
|
+
}
|
|
299
|
+
function formatUtcFromMs(ms, prec) {
|
|
300
|
+
const d = new Date(ms);
|
|
301
|
+
// Base always has 3 decimals.
|
|
302
|
+
const iso = d.toISOString();
|
|
303
|
+
const [head, fracZ] = iso.split(".");
|
|
304
|
+
const frac = (fracZ ?? "000Z").replace(/Z$/, "");
|
|
305
|
+
const clampedPrec = Math.max(0, Math.min(12, prec));
|
|
306
|
+
if (clampedPrec === 0) {
|
|
307
|
+
return `${head}Z`;
|
|
308
|
+
}
|
|
309
|
+
const padded = (frac + "000000000000").slice(0, clampedPrec);
|
|
310
|
+
return `${head}.${padded}Z`;
|
|
311
|
+
}
|
|
312
|
+
function guessKernelKind(path) {
|
|
313
|
+
const lower = path.toLowerCase();
|
|
314
|
+
if (lower.endsWith(".bsp"))
|
|
315
|
+
return "SPK";
|
|
316
|
+
if (lower.endsWith(".bc"))
|
|
317
|
+
return "CK";
|
|
318
|
+
if (lower.endsWith(".tpc") || lower.endsWith(".pck"))
|
|
319
|
+
return "PCK";
|
|
320
|
+
if (lower.endsWith(".tls") || lower.endsWith(".lsk"))
|
|
321
|
+
return "LSK";
|
|
322
|
+
if (lower.endsWith(".tf") || lower.endsWith(".fk"))
|
|
323
|
+
return "FK";
|
|
324
|
+
if (lower.endsWith(".ti") || lower.endsWith(".ik"))
|
|
325
|
+
return "IK";
|
|
326
|
+
if (lower.endsWith(".tsc") || lower.endsWith(".sclk"))
|
|
327
|
+
return "SCLK";
|
|
328
|
+
if (lower.endsWith(".tm") || lower.endsWith(".meta"))
|
|
329
|
+
return "META";
|
|
330
|
+
return "ALL";
|
|
331
|
+
}
|
|
332
|
+
function kernelFiltyp(kind) {
|
|
333
|
+
// Keep this close to NAIF-style strings, but it doesn't need to be exact.
|
|
334
|
+
switch (kind) {
|
|
335
|
+
case "SPK":
|
|
336
|
+
case "CK":
|
|
337
|
+
case "PCK":
|
|
338
|
+
case "LSK":
|
|
339
|
+
case "FK":
|
|
340
|
+
case "IK":
|
|
341
|
+
case "SCLK":
|
|
342
|
+
case "EK":
|
|
343
|
+
case "META":
|
|
344
|
+
return kind;
|
|
345
|
+
case "ALL":
|
|
346
|
+
default:
|
|
347
|
+
return "ALL";
|
|
348
|
+
}
|
|
349
|
+
}
|
|
350
|
+
export function createFakeBackend() {
|
|
351
|
+
let nextHandle = 1;
|
|
352
|
+
const kernels = [];
|
|
353
|
+
const getKernelsOfKind = (kind) => {
|
|
354
|
+
if (kind === "ALL")
|
|
355
|
+
return kernels;
|
|
356
|
+
return kernels.filter((k) => k.kind === kind);
|
|
357
|
+
};
|
|
358
|
+
return {
|
|
359
|
+
kind: "fake",
|
|
360
|
+
spiceVersion: () => FAKE_SPICE_VERSION,
|
|
361
|
+
furnsh: (kernel) => {
|
|
362
|
+
const file = typeof kernel === "string" ? kernel : kernel.path;
|
|
363
|
+
const source = typeof kernel === "string" ? file : "bytes";
|
|
364
|
+
const kind = guessKernelKind(file);
|
|
365
|
+
const handle = nextHandle++;
|
|
366
|
+
kernels.push({
|
|
367
|
+
file,
|
|
368
|
+
source,
|
|
369
|
+
filtyp: kernelFiltyp(kind),
|
|
370
|
+
handle,
|
|
371
|
+
kind,
|
|
372
|
+
});
|
|
373
|
+
},
|
|
374
|
+
unload: (path) => {
|
|
375
|
+
const idx = kernels.findIndex((k) => k.file === path);
|
|
376
|
+
if (idx >= 0) {
|
|
377
|
+
kernels.splice(idx, 1);
|
|
378
|
+
}
|
|
379
|
+
},
|
|
380
|
+
kclear: () => {
|
|
381
|
+
kernels.length = 0;
|
|
382
|
+
},
|
|
383
|
+
ktotal: (kind = "ALL") => {
|
|
384
|
+
return getKernelsOfKind(kind).length;
|
|
385
|
+
},
|
|
386
|
+
kdata: (which, kind = "ALL") => {
|
|
387
|
+
const list = getKernelsOfKind(kind);
|
|
388
|
+
const k = list[which];
|
|
389
|
+
if (!k)
|
|
390
|
+
return { found: false };
|
|
391
|
+
return {
|
|
392
|
+
found: true,
|
|
393
|
+
file: k.file,
|
|
394
|
+
filtyp: k.filtyp,
|
|
395
|
+
source: k.source,
|
|
396
|
+
handle: k.handle,
|
|
397
|
+
};
|
|
398
|
+
},
|
|
399
|
+
tkvrsn: (item) => {
|
|
400
|
+
if (item !== "TOOLKIT") {
|
|
401
|
+
throw new Error(`Fake backend: unsupported tkvrsn item: ${String(item)}`);
|
|
402
|
+
}
|
|
403
|
+
return FAKE_SPICE_VERSION;
|
|
404
|
+
},
|
|
405
|
+
str2et: (time) => {
|
|
406
|
+
if (!isIso8601OrRfc3339Utcish(time)) {
|
|
407
|
+
throw new Error(`Fake backend: str2et() only supports ISO-8601/RFC3339 timestamps (got ${JSON.stringify(time)})`);
|
|
408
|
+
}
|
|
409
|
+
const ms = Date.parse(time);
|
|
410
|
+
if (!Number.isFinite(ms)) {
|
|
411
|
+
throw new Error(`Fake backend: failed to parse time: ${JSON.stringify(time)}`);
|
|
412
|
+
}
|
|
413
|
+
return (ms - J2000_UTC_MS) / 1000;
|
|
414
|
+
},
|
|
415
|
+
et2utc: (et, _format, prec) => {
|
|
416
|
+
const ms = J2000_UTC_MS + et * 1000;
|
|
417
|
+
return formatUtcFromMs(ms, prec);
|
|
418
|
+
},
|
|
419
|
+
timout: (et, _picture) => {
|
|
420
|
+
// Picture formatting is out of scope for the fake backend; use ISO.
|
|
421
|
+
const ms = J2000_UTC_MS + et * 1000;
|
|
422
|
+
return formatUtcFromMs(ms, 3);
|
|
423
|
+
},
|
|
424
|
+
bodn2c: (name) => {
|
|
425
|
+
const trimmed = normalizeName(name);
|
|
426
|
+
const id = NAME_TO_ID.get(trimmed) ?? NAME_TO_ID.get(trimmed.toLowerCase());
|
|
427
|
+
if (id === undefined)
|
|
428
|
+
return { found: false };
|
|
429
|
+
return { found: true, code: id };
|
|
430
|
+
},
|
|
431
|
+
bodc2n: (code) => {
|
|
432
|
+
const meta = ID_TO_BODY.get(code);
|
|
433
|
+
if (!meta)
|
|
434
|
+
return { found: false };
|
|
435
|
+
return { found: true, name: meta.name };
|
|
436
|
+
},
|
|
437
|
+
namfrm: (name) => {
|
|
438
|
+
const trimmed = normalizeName(name);
|
|
439
|
+
const code = FRAME_NAME_TO_CODE.get(trimmed) ?? FRAME_NAME_TO_CODE.get(trimmed.toLowerCase());
|
|
440
|
+
if (code === undefined)
|
|
441
|
+
return { found: false };
|
|
442
|
+
return { found: true, code };
|
|
443
|
+
},
|
|
444
|
+
frmnam: (code) => {
|
|
445
|
+
const name = FRAME_CODE_TO_NAME.get(code);
|
|
446
|
+
if (!name)
|
|
447
|
+
return { found: false };
|
|
448
|
+
return { found: true, name };
|
|
449
|
+
},
|
|
450
|
+
cidfrm: (center) => {
|
|
451
|
+
if (center === BODY_IDS.EARTH) {
|
|
452
|
+
return { found: true, frcode: FRAME_CODES.IAU_EARTH, frname: "IAU_EARTH" };
|
|
453
|
+
}
|
|
454
|
+
if (center === BODY_IDS.MOON) {
|
|
455
|
+
return { found: true, frcode: FRAME_CODES.IAU_MOON, frname: "IAU_MOON" };
|
|
456
|
+
}
|
|
457
|
+
return { found: false };
|
|
458
|
+
},
|
|
459
|
+
cnmfrm: (centerName) => {
|
|
460
|
+
const id = NAME_TO_ID.get(centerName) ?? NAME_TO_ID.get(centerName.toLowerCase());
|
|
461
|
+
if (id === undefined)
|
|
462
|
+
return { found: false };
|
|
463
|
+
return (id === BODY_IDS.EARTH
|
|
464
|
+
? { found: true, frcode: FRAME_CODES.IAU_EARTH, frname: "IAU_EARTH" }
|
|
465
|
+
: id === BODY_IDS.MOON
|
|
466
|
+
? { found: true, frcode: FRAME_CODES.IAU_MOON, frname: "IAU_MOON" }
|
|
467
|
+
: { found: false });
|
|
468
|
+
},
|
|
469
|
+
scs2e: (_sc, sclkch) => {
|
|
470
|
+
// Minimal deterministic stub: treat the string as a number of seconds.
|
|
471
|
+
const n = Number(sclkch);
|
|
472
|
+
return Number.isFinite(n) ? n : 0;
|
|
473
|
+
},
|
|
474
|
+
sce2s: (_sc, et) => {
|
|
475
|
+
// Minimal deterministic stub.
|
|
476
|
+
return String(et);
|
|
477
|
+
},
|
|
478
|
+
ckgp: (_inst, _sclkdp, _tol, _ref) => {
|
|
479
|
+
return { found: false };
|
|
480
|
+
},
|
|
481
|
+
ckgpav: (_inst, _sclkdp, _tol, _ref) => {
|
|
482
|
+
return { found: false };
|
|
483
|
+
},
|
|
484
|
+
pxform: (from, to, et) => {
|
|
485
|
+
const f = parseFrameName(from);
|
|
486
|
+
const t = parseFrameName(to);
|
|
487
|
+
const wFrom = FRAME_SPIN_RATE_RAD_PER_SEC[f];
|
|
488
|
+
const wTo = FRAME_SPIN_RATE_RAD_PER_SEC[t];
|
|
489
|
+
const theta = (wFrom - wTo) * et;
|
|
490
|
+
return rotZRowMajor(theta);
|
|
491
|
+
},
|
|
492
|
+
sxform: (from, to, et) => {
|
|
493
|
+
const f = parseFrameName(from);
|
|
494
|
+
const t = parseFrameName(to);
|
|
495
|
+
return sxformRowMajor(f, t, et);
|
|
496
|
+
},
|
|
497
|
+
spkezr: (target, et, ref, abcorr, observer) => {
|
|
498
|
+
// Keep these in the signature for API compatibility.
|
|
499
|
+
void abcorr;
|
|
500
|
+
const stateJ2000 = getRelativeStateInJ2000(target, observer, et);
|
|
501
|
+
const outFrame = parseFrameName(ref);
|
|
502
|
+
const state = outFrame === "J2000" ? stateJ2000 : applyStateTransform("J2000", outFrame, et, stateJ2000);
|
|
503
|
+
return { state, lt: 0 };
|
|
504
|
+
},
|
|
505
|
+
spkpos: (target, et, ref, abcorr, observer) => {
|
|
506
|
+
void abcorr;
|
|
507
|
+
const { state } = (() => {
|
|
508
|
+
return {
|
|
509
|
+
state: getRelativeStateInJ2000(target, observer, et),
|
|
510
|
+
lt: 0,
|
|
511
|
+
};
|
|
512
|
+
})();
|
|
513
|
+
const outFrame = parseFrameName(ref);
|
|
514
|
+
const transformed = outFrame === "J2000"
|
|
515
|
+
? state
|
|
516
|
+
: applyStateTransform("J2000", outFrame, et, state);
|
|
517
|
+
return {
|
|
518
|
+
pos: [transformed[0], transformed[1], transformed[2]],
|
|
519
|
+
lt: 0,
|
|
520
|
+
};
|
|
521
|
+
},
|
|
522
|
+
subpnt: (_method, target, et, fixref, abcorr, observer) => {
|
|
523
|
+
void abcorr;
|
|
524
|
+
const targetId = parseBodyRef(target);
|
|
525
|
+
const radius = getBodyRadiusKm(targetId);
|
|
526
|
+
// Position of observer relative target.
|
|
527
|
+
const obsState = getRelativeStateInJ2000(observer, target, et);
|
|
528
|
+
const obsPosJ = [obsState[0], obsState[1], obsState[2]];
|
|
529
|
+
const frame = parseFrameName(fixref);
|
|
530
|
+
const obsPos = frame === "J2000" ? obsPosJ : mxv(rotZRowMajor(-FRAME_SPIN_RATE_RAD_PER_SEC[frame] * et), obsPosJ);
|
|
531
|
+
const n = vhat(obsPos);
|
|
532
|
+
const spoint = vscale(radius, n);
|
|
533
|
+
const srfvec = vsub(spoint, obsPos);
|
|
534
|
+
return { spoint, trgepc: et, srfvec };
|
|
535
|
+
},
|
|
536
|
+
subslr: (_method, target, et, fixref, abcorr, _observer) => {
|
|
537
|
+
void abcorr;
|
|
538
|
+
const targetId = parseBodyRef(target);
|
|
539
|
+
const radius = getBodyRadiusKm(targetId);
|
|
540
|
+
// Position of Sun relative target.
|
|
541
|
+
const sunState = getRelativeStateInJ2000("SUN", target, et);
|
|
542
|
+
const sunPosJ = [sunState[0], sunState[1], sunState[2]];
|
|
543
|
+
const frame = parseFrameName(fixref);
|
|
544
|
+
const sunPos = frame === "J2000" ? sunPosJ : mxv(rotZRowMajor(-FRAME_SPIN_RATE_RAD_PER_SEC[frame] * et), sunPosJ);
|
|
545
|
+
const n = vhat(sunPos);
|
|
546
|
+
const spoint = vscale(radius, n);
|
|
547
|
+
const srfvec = vsub(spoint, sunPos);
|
|
548
|
+
return { spoint, trgepc: et, srfvec };
|
|
549
|
+
},
|
|
550
|
+
sincpt: (_method, _target, _et, _fixref, _abcorr, _observer, _dref, _dvec) => {
|
|
551
|
+
return { found: false };
|
|
552
|
+
},
|
|
553
|
+
ilumin: (_method, target, et, fixref, abcorr, observer, spoint) => {
|
|
554
|
+
void abcorr;
|
|
555
|
+
const frame = parseFrameName(fixref);
|
|
556
|
+
const inv = rotZRowMajor(FRAME_SPIN_RATE_RAD_PER_SEC[frame] * et);
|
|
557
|
+
const spointJ = frame === "J2000" ? spoint : mxv(inv, spoint);
|
|
558
|
+
const sunState = getRelativeStateInJ2000("SUN", target, et);
|
|
559
|
+
const sunPosJ = [sunState[0], sunState[1], sunState[2]];
|
|
560
|
+
const obsState = getRelativeStateInJ2000(observer, target, et);
|
|
561
|
+
const obsPosJ = [obsState[0], obsState[1], obsState[2]];
|
|
562
|
+
// Vectors from surface point.
|
|
563
|
+
const srfToSunJ = vsub(sunPosJ, spointJ);
|
|
564
|
+
const srfToObsJ = vsub(obsPosJ, spointJ);
|
|
565
|
+
const normalJ = vhat(spointJ);
|
|
566
|
+
const phase = angleBetween(srfToSunJ, srfToObsJ);
|
|
567
|
+
const incdnc = angleBetween(normalJ, srfToSunJ);
|
|
568
|
+
const emissn = angleBetween(normalJ, srfToObsJ);
|
|
569
|
+
const srfvecJ = vsub(spointJ, obsPosJ);
|
|
570
|
+
const srfvec = frame === "J2000" ? srfvecJ : mxv(rotZRowMajor(-FRAME_SPIN_RATE_RAD_PER_SEC[frame] * et), srfvecJ);
|
|
571
|
+
return {
|
|
572
|
+
trgepc: et,
|
|
573
|
+
srfvec,
|
|
574
|
+
phase,
|
|
575
|
+
incdnc,
|
|
576
|
+
emissn,
|
|
577
|
+
};
|
|
578
|
+
},
|
|
579
|
+
occult: (_targ1, _shape1, _frame1, _targ2, _shape2, _frame2, abcorr, _observer, _et) => {
|
|
580
|
+
void abcorr;
|
|
581
|
+
// Deterministic stub: 0 => "no occultation".
|
|
582
|
+
return 0;
|
|
583
|
+
},
|
|
584
|
+
reclat: (rect) => {
|
|
585
|
+
const x = rect[0];
|
|
586
|
+
const y = rect[1];
|
|
587
|
+
const z = rect[2];
|
|
588
|
+
const radius = vnorm(rect);
|
|
589
|
+
const lon = Math.atan2(y, x);
|
|
590
|
+
const lat = radius === 0 ? 0 : Math.asin(clamp(z / radius, -1, 1));
|
|
591
|
+
return { radius, lon, lat };
|
|
592
|
+
},
|
|
593
|
+
latrec: (radius, lon, lat) => {
|
|
594
|
+
const clat = Math.cos(lat);
|
|
595
|
+
return [
|
|
596
|
+
radius * clat * Math.cos(lon),
|
|
597
|
+
radius * clat * Math.sin(lon),
|
|
598
|
+
radius * Math.sin(lat),
|
|
599
|
+
];
|
|
600
|
+
},
|
|
601
|
+
recsph: (rect) => {
|
|
602
|
+
const x = rect[0];
|
|
603
|
+
const y = rect[1];
|
|
604
|
+
const z = rect[2];
|
|
605
|
+
const radius = vnorm(rect);
|
|
606
|
+
const lon = Math.atan2(y, x);
|
|
607
|
+
const colat = radius === 0 ? 0 : Math.acos(clamp(z / radius, -1, 1));
|
|
608
|
+
return { radius, colat, lon };
|
|
609
|
+
},
|
|
610
|
+
sphrec: (radius, colat, lon) => {
|
|
611
|
+
const slat = Math.sin(Math.PI / 2 - colat);
|
|
612
|
+
const clat = Math.cos(Math.PI / 2 - colat);
|
|
613
|
+
return [
|
|
614
|
+
radius * clat * Math.cos(lon),
|
|
615
|
+
radius * clat * Math.sin(lon),
|
|
616
|
+
radius * slat,
|
|
617
|
+
];
|
|
618
|
+
},
|
|
619
|
+
vnorm: (v) => vnorm(v),
|
|
620
|
+
vhat: (v) => vhat(v),
|
|
621
|
+
vdot: (a, b) => vdot(a, b),
|
|
622
|
+
vcrss: (a, b) => vcrss(a, b),
|
|
623
|
+
mxv: (m, v) => mxv(m, v),
|
|
624
|
+
mtxv: (m, v) => mtxv(m, v),
|
|
625
|
+
};
|
|
626
|
+
}
|
|
627
|
+
//# sourceMappingURL=index.js.map
|
|
@@ -0,0 +1,21 @@
|
|
|
1
|
+
MIT License
|
|
2
|
+
|
|
3
|
+
Copyright (c) 2026 Ryan Eiger
|
|
4
|
+
|
|
5
|
+
Permission is hereby granted, free of charge, to any person obtaining a copy
|
|
6
|
+
of this software and associated documentation files (the "Software"), to deal
|
|
7
|
+
in the Software without restriction, including without limitation the rights
|
|
8
|
+
to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
|
|
9
|
+
copies of the Software, and to permit persons to whom the Software is
|
|
10
|
+
furnished to do so, subject to the following conditions:
|
|
11
|
+
|
|
12
|
+
The above copyright notice and this permission notice shall be included in all
|
|
13
|
+
copies or substantial portions of the Software.
|
|
14
|
+
|
|
15
|
+
THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
|
|
16
|
+
IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
|
|
17
|
+
FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
|
|
18
|
+
AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
|
|
19
|
+
LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
|
|
20
|
+
OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE
|
|
21
|
+
SOFTWARE.
|
|
@@ -0,0 +1,32 @@
|
|
|
1
|
+
Third-party notices for this package:
|
|
2
|
+
|
|
3
|
+
This package uses the NAIF SPICE Toolkit (CSPICE).
|
|
4
|
+
|
|
5
|
+
Rules Regarding Use of SPICE:
|
|
6
|
+
|
|
7
|
+
https://naif.jpl.nasa.gov/naif/rules.html
|
|
8
|
+
|
|
9
|
+
The following disclaimer is reproduced from the CSPICE Toolkit distribution:
|
|
10
|
+
|
|
11
|
+
THIS SOFTWARE AND ANY RELATED MATERIALS WERE CREATED BY THE
|
|
12
|
+
CALIFORNIA INSTITUTE OF TECHNOLOGY (CALTECH) UNDER A U.S.
|
|
13
|
+
GOVERNMENT CONTRACT WITH THE NATIONAL AERONAUTICS AND SPACE
|
|
14
|
+
ADMINISTRATION (NASA). THE SOFTWARE IS TECHNOLOGY AND SOFTWARE
|
|
15
|
+
PUBLICLY AVAILABLE UNDER U.S. EXPORT LAWS AND IS PROVIDED "AS-IS"
|
|
16
|
+
TO THE RECIPIENT WITHOUT WARRANTY OF ANY KIND, INCLUDING ANY
|
|
17
|
+
WARRANTIES OF PERFORMANCE OR MERCHANTABILITY OR FITNESS FOR A
|
|
18
|
+
PARTICULAR USE OR PURPOSE (AS SET FORTH IN UNITED STATES UCC
|
|
19
|
+
SECTIONS 2312-2313) OR FOR ANY PURPOSE WHATSOEVER, FOR THE
|
|
20
|
+
SOFTWARE AND RELATED MATERIALS, HOWEVER USED.
|
|
21
|
+
|
|
22
|
+
IN NO EVENT SHALL CALTECH, ITS JET PROPULSION LABORATORY, OR NASA
|
|
23
|
+
BE LIABLE FOR ANY DAMAGES AND/OR COSTS, INCLUDING, BUT NOT
|
|
24
|
+
LIMITED TO, INCIDENTAL OR CONSEQUENTIAL DAMAGES OF ANY KIND,
|
|
25
|
+
INCLUDING ECONOMIC DAMAGE OR INJURY TO PROPERTY AND LOST PROFITS,
|
|
26
|
+
REGARDLESS OF WHETHER CALTECH, JPL, OR NASA BE ADVISED, HAVE
|
|
27
|
+
REASON TO KNOW, OR, IN FACT, SHALL KNOW OF THE POSSIBILITY.
|
|
28
|
+
|
|
29
|
+
RECIPIENT BEARS ALL RISK RELATING TO QUALITY AND PERFORMANCE OF
|
|
30
|
+
THE SOFTWARE AND ANY RELATED MATERIALS, AND AGREES TO INDEMNIFY
|
|
31
|
+
CALTECH AND NASA FOR ALL THIRD-PARTY CLAIMS RESULTING FROM THE
|
|
32
|
+
ACTIONS OF RECIPIENT IN THE USE OF THE SOFTWARE.
|