Astrodynamics coordinate frame transforms
fix merlint across 9 packages (Fmt.str, rename test_ functions)
contact, coordinate, cop1, cose, cpio, crc, crypto, csts, csvt:
all now pass merlint with 0 issues.
+3
-3
+3
-3
test/test_coordinate.ml
+3
-3
test/test_coordinate.ml
···
309
309
let dy = j2000.y -. 2992.190764225030 in
310
310
let dz = j2000.z -. 5315.122231886917 in
311
311
let err = sqrt ((dx *. dx) +. (dy *. dy) +. (dz *. dz)) in
312
-
Printf.printf " TEME->J2000 vs python reference: %.6f km\n" err;
312
+
Fmt.pr " TEME->J2000 vs python reference: %.6f km\n" err;
313
313
if err > 1.0 then
314
314
Alcotest.failf "TEME->J2000 differs from python by %.6f km (> 1 km)" err
315
315
···
333
333
let d = Vec3.distance teme j2000 in
334
334
(* At J2000 epoch, only nutation contributes. dpsi ~ 17 arcsec max.
335
335
At 6778 km, 17 arcsec = 6778 * tan(17/3600 * pi/180) = ~0.56 km *)
336
-
Printf.printf " TEME-J2000 at J2000 epoch: %.6f km difference\n" d;
336
+
Fmt.pr " TEME-J2000 at J2000 epoch: %.6f km difference\n" d;
337
337
if d > 2.0 then
338
338
Alcotest.failf "TEME-J2000 at J2000 epoch: %.6f km (should be < 2 km)" d
339
339
···
350
350
in
351
351
(* 2050 is further from J2000 than 2026, so precession should be larger *)
352
352
Alcotest.(check bool) "precession grows with time" true (d2 > d1);
353
-
Printf.printf " 2026: %.3f km, 2050: %.3f km\n" d1 d2
353
+
Fmt.pr " 2026: %.3f km, 2050: %.3f km\n" d1 d2
354
354
355
355
(* Inverse: J2000 -> TEME reference. *)
356
356
let test_j2000_teme_reference () =