siderust-js 0.1.0

package/tempoch-js/tempoch/tempoch-core/src/julian_date_ext.rs

@@ -0,0 +1,142 @@
+// SPDX-License-Identifier: AGPL-3.0-only
+// Copyright (C) 2026 Vallés Puig, Ramon
+
+//! Julian Date (`Time<JD>`) specific extensions.
+
+use qtty::*;
+use std::ops::Add;
+
+use super::instant::Time;
+use super::scales::{JD, MJD};
+
+impl Time<JD> {
+    /// J2000.0 epoch: 2000-01-01T12:00:00 TT  (JD 2 451 545.0).
+    pub const J2000: Self = Self::new(2_451_545.0);
+
+    /// One Julian year expressed in days.
+    pub const JULIAN_YEAR: Days = Days::new(365.25);
+
+    /// One Julian century expressed in days.
+    pub const JULIAN_CENTURY: Days = Days::new(36_525.0);
+
+    /// One Julian millennium expressed in days.
+    pub const JULIAN_MILLENNIUM: Days = Days::new(365_250.0);
+
+    /// Julian millennia since J2000.0 (used by VSOP87).
+    #[inline]
+    pub fn julian_millennias(&self) -> Millennia {
+        Millennia::new(
+            ((*self - Self::J2000) / Self::JULIAN_MILLENNIUM)
+                .simplify()
+                .value(),
+        )
+    }
+
+    /// Julian centuries since J2000.0 (used by nutation, precession, sidereal time).
+    #[inline]
+    pub fn julian_centuries(&self) -> Centuries {
+        Centuries::new(
+            ((*self - Self::J2000) / Self::JULIAN_CENTURY)
+                .simplify()
+                .value(),
+        )
+    }
+
+    /// Julian years since J2000.0.
+    #[inline]
+    pub fn julian_years(&self) -> JulianYears {
+        JulianYears::new(
+            ((*self - Self::J2000) / Self::JULIAN_YEAR)
+                .simplify()
+                .value(),
+        )
+    }
+
+    /// Converts JD(TT) → JD(TDB) using the Fairhead & Bretagnon (1990)
+    /// expression for `TDB − TT`.
+    ///
+    /// The dominant term has an amplitude of ≈1.658 ms. This implementation
+    /// includes the four largest periodic terms plus a secular component,
+    /// matching the formula recommended by USNO Circular 179 (Kaplan 2005)
+    /// and consistent with IAU 2006 Resolution B3.
+    ///
+    /// Accuracy: better than 30 μs for dates within ±10 000 years of J2000.
+    ///
+    /// ## References
+    /// * Fairhead & Bretagnon (1990), A&A 229, 240
+    /// * USNO Circular 179, eq. 2.6
+    /// * SOFA `iauDtdb` (full implementation has hundreds of terms)
+    pub fn tt_to_tdb(jd_tt: Self) -> Self {
+        jd_tt + super::scales::tdb_minus_tt_days(jd_tt.quantity())
+    }
+
+    /// Convenience: MJD value corresponding to this JD.
+    ///
+    /// Kept as a convenience wrapper for `self.to::<MJD>()`.
+    #[inline]
+    pub fn to_mjd(&self) -> Time<MJD> {
+        self.to::<MJD>()
+    }
+}
+
+// ── From / Into conversions for qtty time quantities (on JD only) ────────
+
+impl Add<Years> for Time<JD> {
+    type Output = Self;
+    fn add(self, years: Years) -> Self {
+        // Treat `Years` here as Julian years for JD arithmetic stability.
+        self + Days::new(years.value() * Self::JULIAN_YEAR.value())
+    }
+}
+
+impl From<JulianYears> for Time<JD> {
+    fn from(years: JulianYears) -> Self {
+        Self::J2000 + years.to::<Day>()
+    }
+}
+
+impl From<Time<JD>> for JulianYears {
+    fn from(jd: Time<JD>) -> Self {
+        jd.julian_years()
+    }
+}
+
+impl From<Centuries> for Time<JD> {
+    fn from(centuries: Centuries) -> Self {
+        // `Centuries` are interpreted as Julian centuries relative to J2000.
+        Self::J2000 + Days::new(centuries.value() * Self::JULIAN_CENTURY.value())
+    }
+}
+
+impl From<Time<JD>> for Centuries {
+    fn from(jd: Time<JD>) -> Self {
+        jd.julian_centuries()
+    }
+}
+
+impl From<Millennia> for Time<JD> {
+    fn from(millennia: Millennia) -> Self {
+        // `Millennia` are interpreted as Julian millennia relative to J2000.
+        Self::J2000 + Days::new(millennia.value() * Self::JULIAN_MILLENNIUM.value())
+    }
+}
+
+impl From<Time<JD>> for Millennia {
+    fn from(jd: Time<JD>) -> Self {
+        jd.julian_millennias()
+    }
+}
+
+#[cfg(test)]
+mod tests {
+    use crate::{Time, JD, MJD};
+
+    #[test]
+    fn to_mjd_convenience_matches_to_generic() {
+        // to_mjd() is a convenience wrapper for .to::<MJD>().
+        let jd = Time::<JD>::J2000;
+        let via_convenience = jd.to_mjd();
+        let via_generic = jd.to::<MJD>();
+        assert_eq!(via_convenience.value(), via_generic.value());
+    }
+}

package/tempoch-js/tempoch/tempoch-core/src/lib.rs

@@ -0,0 +1,81 @@
+// SPDX-License-Identifier: AGPL-3.0-only
+// Copyright (C) 2026 Vallés Puig, Ramon
+
+//! Time Module
+//!
+//! This module provides time-related types and abstractions for astronomical calculations.
+//!
+//! # Core types
+//!
+//! - [`Time<S>`] — generic instant parameterised by a [`TimeScale`] marker.
+//! - [`TimeScale`] — trait that defines a time scale (epoch offset + conversions).
+//! - [`JulianDate`] — type alias for `Time<JD>`.
+//! - [`JulianEphemerisDay`] — type alias for `Time<JDE>`.
+//! - [`ModifiedJulianDate`] — type alias for `Time<MJD>`.
+//! - [`Period<S>`] — a time interval parameterised by a [`TimeScale`] marker.
+//! - [`Interval<T>`] — a generic interval over any [`TimeInstant`].
+//! - [`TimeInstant`] — trait for points in time usable with [`Interval`].
+//!
+//! # Time scales
+//!
+//! The following markers implement [`TimeScale`]:
+//!
+//! | Marker | Scale |
+//! |--------|-------|
+//! | [`JD`] | Julian Date |
+//! | [`JDE`] | Julian Ephemeris Day |
+//! | [`MJD`] | Modified Julian Date |
+//! | [`TDB`] | Barycentric Dynamical Time |
+//! | [`TT`] | Terrestrial Time |
+//! | [`TAI`] | International Atomic Time |
+//! | [`TCG`] | Geocentric Coordinate Time |
+//! | [`TCB`] | Barycentric Coordinate Time |
+//! | [`GPS`] | GPS Time |
+//! | [`UnixTime`] | Unix / POSIX time |
+//! | [`UT`] | Universal Time (Earth rotation) |
+//!
+//! # ΔT (Delta T)
+//!
+//! The difference **ΔT = TT − UT** is applied automatically by the
+//! [`UT`] time scale.  Use `Time::<UT>::new(jd_ut)` for UT-based values,
+//! or construct any scale via `from_utc()` which routes through `UT` internally.
+//! The raw ΔT value (in seconds) is available via [`Time::<UT>::delta_t()`](Time::delta_t).
+
+mod delta_t;
+pub(crate) mod instant;
+mod julian_date_ext;
+mod period;
+mod scales;
+
+// ── Re-exports ────────────────────────────────────────────────────────────
+
+pub use instant::{NonFiniteTimeError, Time, TimeInstant, TimeScale};
+pub use period::{
+    complement_within, intersect_periods, normalize_periods, validate_period_list, ConversionError,
+    Interval, InvalidIntervalError, Period, PeriodListError, UtcPeriod,
+};
+pub use scales::{tai_minus_utc, UnixTime, GPS, JD, JDE, MJD, TAI, TCB, TCG, TDB, TT, UT};
+
+// ── Backward-compatible type aliases ──────────────────────────────────────
+
+/// Julian Date — continuous count of days since the Julian Period.
+///
+/// This is a type alias for [`Time<JD>`].  All historical call-sites
+/// (`JulianDate::new(...)`, `JulianDate::J2000`, `.julian_centuries()`, …)
+/// continue to work without modification.
+pub type JulianDate = Time<JD>;
+
+/// Julian Ephemeris Day — dynamical Julian day used by many ephemeris formulas.
+///
+/// This is a type alias for [`Time<JDE>`].
+pub type JulianEphemerisDay = Time<JDE>;
+
+/// Modified Julian Date — `JD − 2 400 000.5`.
+///
+/// This is a type alias for [`Time<MJD>`].
+pub type ModifiedJulianDate = Time<MJD>;
+
+/// Universal Time — Earth-rotation civil time scale.
+///
+/// This is a type alias for [`Time<UT>`].
+pub type UniversalTime = Time<UT>;