empyrean 0.8.2

Uncertainty-first orbit propagation, ephemeris, orbit determination, and event detection for asteroids and comets, powered by automatic differentiation
Documentation 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89

//! Ephemeris I/O — write across parquet, JSON, and CSV.
//!
//! Read paths are not exposed: the propagator's
//! [`Context::generate_ephemeris`](crate::Context::generate_ephemeris)
//! is the canonical producer.

use std::ffi::CString;
use std::path::Path;

use crate::ephemeris::EphemerisEntry;
use crate::error::{Error, Result};

use super::path_to_cstring;

fn ephemeris_to_ffi_array(
    entries: &[EphemerisEntry],
) -> Result<(Vec<empyrean_sys::EmpyreanEphemerisEntry>, Vec<CString>)> {
    let mut keep = Vec::with_capacity(entries.len());
    let ffi: Vec<empyrean_sys::EmpyreanEphemerisEntry> = entries
        .iter()
        .map(|e| {
            let id_c = CString::new(e.orbit_id.as_bytes()).unwrap_or_default();
            let id_ptr = id_c.as_ptr() as *mut std::ffi::c_char;
            keep.push(id_c);
            let mut obs_code = [0u8; 4];
            let bytes = e.obs_code.as_bytes();
            let n = bytes.len().min(3);
            obs_code[..n].copy_from_slice(&bytes[..n]);
            Ok(empyrean_sys::EmpyreanEphemerisEntry {
                orbit_id: id_ptr,
                epoch_mjd_tdb: e.epoch.mjd_tdb()?,
                ra_deg: e.ra_deg,
                dec_deg: e.dec_deg,
                rho_au: e.rho_au,
                vrho_au_day: e.vrho_au_day,
                vra_deg_day: e.vra_deg_day,
                vdec_deg_day: e.vdec_deg_day,
                light_time_days: e.light_time_days,
                phase_angle_deg: e.phase_angle_deg,
                elongation_deg: e.elongation_deg,
                heliocentric_distance_au: e.heliocentric_distance_au,
                mag: e.mag,
                mag_sigma: e.mag_sigma,
                zenith_angle_deg: e.zenith_angle_deg,
                azimuth_deg: e.azimuth_deg,
                hour_angle_deg: e.hour_angle_deg,
                lunar_elongation_deg: e.lunar_elongation_deg,
                position_angle_deg: e.position_angle_deg,
                sky_rate_deg_day: e.sky_rate_deg_day,
                obs_code,
            })
        })
        .collect::<Result<Vec<_>>>()?;
    Ok((ffi, keep))
}

fn write_ephemeris_via<F>(path: &Path, entries: &[EphemerisEntry], c_call: F) -> Result<()>
where
    F: FnOnce(*const std::ffi::c_char, *const empyrean_sys::EmpyreanEphemerisEntry, usize) -> i32,
{
    let path_c = path_to_cstring(path)?;
    let (ffi, _keep) = ephemeris_to_ffi_array(entries)?;
    let code = c_call(path_c.as_ptr(), ffi.as_ptr(), ffi.len());
    if code != 0 {
        return Err(Error::capture(code));
    }
    Ok(())
}

/// Write ephemeris entries to a parquet file.
pub fn write_ephemeris_parquet(path: impl AsRef<Path>, entries: &[EphemerisEntry]) -> Result<()> {
    write_ephemeris_via(path.as_ref(), entries, |p, ptr, n| unsafe {
        empyrean_sys::empyrean_ephemeris_write_parquet(p, ptr, n)
    })
}

/// Write ephemeris entries to JSON.
pub fn write_ephemeris_json(path: impl AsRef<Path>, entries: &[EphemerisEntry]) -> Result<()> {
    write_ephemeris_via(path.as_ref(), entries, |p, ptr, n| unsafe {
        empyrean_sys::empyrean_ephemeris_write_json(p, ptr, n)
    })
}

/// Write ephemeris entries to CSV.
pub fn write_ephemeris_csv(path: impl AsRef<Path>, entries: &[EphemerisEntry]) -> Result<()> {
    write_ephemeris_via(path.as_ref(), entries, |p, ptr, n| unsafe {
        empyrean_sys::empyrean_ephemeris_write_csv(p, ptr, n)
    })
}