lox-space 0.1.0-alpha.48

The Lox toolbox for space mission analysis and design
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

// SPDX-FileCopyrightText: 2026 Helge Eichhorn <git@helgeeichhorn.de>
//
// SPDX-License-Identifier: MPL-2.0

//! Analytical Sun ephemeris adapter.
//!
//! Wraps [`lox_earth::ephemeris::apparent_sun_position`] behind the
//! [`Ephemeris`](lox_ephem::Ephemeris) trait so that Earth-centred analyses can run without an
//! SPK file.

use std::error::Error;
use std::fmt::{self, Display, Formatter};

use lox_bodies::{NaifId, Origin};
use lox_core::coords::Cartesian;
use lox_core::glam::DVec3;
use lox_earth::ephemeris::apparent_sun_position;
use lox_ephem::Ephemeris;
use lox_time::{Time, time_scales::Tdb};

/// Error returned by [`AnalyticalSunEphemeris`] for unsupported body pairs.
#[derive(Debug)]
pub struct AnalyticalSunEphemerisError {
    origin_id: NaifId,
    target_id: NaifId,
}

impl Display for AnalyticalSunEphemerisError {
    fn fmt(&self, f: &mut Formatter<'_>) -> fmt::Result {
        write!(
            f,
            "analytical ephemeris only supports Earth<->Sun queries (got {}->{})",
            self.origin_id, self.target_id
        )
    }
}

impl Error for AnalyticalSunEphemerisError {}

/// Ephemeris adapter that uses [`apparent_sun_position`] for Earth↔Sun
/// queries and returns an error for anything else.
///
/// This allows power-budget and eclipse analyses to run without loading
/// an SPK kernel, at the cost of being limited to Earth-centred scenarios.
pub struct AnalyticalSunEphemeris;

impl Ephemeris for AnalyticalSunEphemeris {
    type Error = AnalyticalSunEphemerisError;

    fn state<O1: Origin, O2: Origin>(
        &self,
        time: Time<Tdb>,
        origin: O1,
        target: O2,
    ) -> Result<Cartesian, Self::Error> {
        let origin_id = origin.id();
        let target_id = target.id();
        // Earth = 399, Sun = 10
        if origin_id == NaifId(399) && target_id == NaifId(10) {
            let pos = apparent_sun_position(time);
            Ok(Cartesian::from_vecs(pos, DVec3::ZERO))
        } else if origin_id == NaifId(10) && target_id == NaifId(399) {
            let pos = apparent_sun_position(time);
            Ok(Cartesian::from_vecs(-pos, DVec3::ZERO))
        } else {
            Err(AnalyticalSunEphemerisError {
                origin_id,
                target_id,
            })
        }
    }
}