siderust-js 0.1.0

package/siderust-web/src/phase.rs

@@ -0,0 +1,190 @@
+// SPDX-License-Identifier: AGPL-3.0-only
+// Copyright (C) 2026 Vallés Puig, Ramon
+
+//! Lunar phase queries — geometry, events, and illumination periods.
+
+use serde::Serialize;
+use wasm_bindgen::prelude::*;
+
+use crate::body::to_js;
+use crate::events::{make_window, MjdPeriod};
+use crate::observer::Observer;
+
+use qtty::*;
+use siderust::calculus::ephemeris::Vsop87Ephemeris;
+use siderust::calculus::lunar::phase::{
+    find_phase_events, illumination_above, illumination_below, illumination_range,
+    moon_phase_geocentric, moon_phase_topocentric, MoonPhaseLabel, PhaseSearchOpts,
+};
+use siderust::time::JulianDate;
+
+// ─────────────────────────────────────────────────────────────────────────────
+// Result types
+// ─────────────────────────────────────────────────────────────────────────────
+
+#[derive(Serialize)]
+#[serde(rename_all = "camelCase")]
+pub struct MoonPhase {
+    pub phase_angle_deg: f64,
+    pub illuminated_fraction: f64,
+    pub elongation_deg: f64,
+    pub waxing: bool,
+    pub label: String,
+}
+
+#[derive(Serialize)]
+#[serde(rename_all = "camelCase")]
+pub struct PhaseEvent {
+    pub mjd: f64,
+    pub kind: String,
+}
+
+// ─────────────────────────────────────────────────────────────────────────────
+// Phase geometry
+// ─────────────────────────────────────────────────────────────────────────────
+
+/// Compute geocentric Moon phase geometry at a Julian Date.
+#[wasm_bindgen(js_name = "moonPhase")]
+pub fn moon_phase(jd: f64) -> Result<JsValue, JsError> {
+    if !jd.is_finite() {
+        return Err(JsError::new("jd must be finite"));
+    }
+    let geom = moon_phase_geocentric::<Vsop87Ephemeris>(JulianDate::new(jd));
+    let label = geom.label();
+    to_js(&MoonPhase {
+        phase_angle_deg: geom.phase_angle.to::<Degree>().value(),
+        illuminated_fraction: geom.illuminated_fraction,
+        elongation_deg: geom.elongation.to::<Degree>().value(),
+        waxing: geom.waxing,
+        label: label_to_string(label),
+    })
+}
+
+/// Compute topocentric Moon phase geometry at a Julian Date for an observer.
+#[wasm_bindgen(js_name = "moonPhaseTopocentric")]
+pub fn moon_phase_topo(jd: f64, observer: &Observer) -> Result<JsValue, JsError> {
+    if !jd.is_finite() {
+        return Err(JsError::new("jd must be finite"));
+    }
+    let geom =
+        moon_phase_topocentric::<Vsop87Ephemeris>(JulianDate::new(jd), observer.inner.clone());
+    let label = geom.label();
+    to_js(&MoonPhase {
+        phase_angle_deg: geom.phase_angle.to::<Degree>().value(),
+        illuminated_fraction: geom.illuminated_fraction,
+        elongation_deg: geom.elongation.to::<Degree>().value(),
+        waxing: geom.waxing,
+        label: label_to_string(label),
+    })
+}
+
+// ─────────────────────────────────────────────────────────────────────────────
+// Phase event finding
+// ─────────────────────────────────────────────────────────────────────────────
+
+/// Find all principal lunar phase events in a time window.
+#[wasm_bindgen(js_name = "findPhaseEvents")]
+pub fn find_phase_events_js(start_mjd: f64, end_mjd: f64) -> Result<JsValue, JsError> {
+    let window = make_window(start_mjd, end_mjd)?;
+    let events = find_phase_events::<Vsop87Ephemeris>(window, PhaseSearchOpts::default());
+    let out: Vec<PhaseEvent> = events
+        .into_iter()
+        .map(|e| PhaseEvent {
+            mjd: e.mjd.value(),
+            kind: phase_kind_to_string(e.kind),
+        })
+        .collect();
+    to_js(&out)
+}
+
+// ─────────────────────────────────────────────────────────────────────────────
+// Illumination periods
+// ─────────────────────────────────────────────────────────────────────────────
+
+/// Find periods where geocentric Moon illumination is above `kMin` ∈ [0, 1].
+#[wasm_bindgen(js_name = "moonIlluminationAbove")]
+pub fn moon_illumination_above(
+    start_mjd: f64,
+    end_mjd: f64,
+    k_min: f64,
+) -> Result<JsValue, JsError> {
+    let window = make_window(start_mjd, end_mjd)?;
+    let periods = illumination_above::<Vsop87Ephemeris>(window, k_min, PhaseSearchOpts::default());
+    let out: Vec<MjdPeriod> = periods
+        .into_iter()
+        .map(|p| MjdPeriod {
+            start_mjd: p.start.value(),
+            end_mjd: p.end.value(),
+        })
+        .collect();
+    to_js(&out)
+}
+
+/// Find periods where geocentric Moon illumination is below `kMax` ∈ [0, 1].
+#[wasm_bindgen(js_name = "moonIlluminationBelow")]
+pub fn moon_illumination_below(
+    start_mjd: f64,
+    end_mjd: f64,
+    k_max: f64,
+) -> Result<JsValue, JsError> {
+    let window = make_window(start_mjd, end_mjd)?;
+    let periods = illumination_below::<Vsop87Ephemeris>(window, k_max, PhaseSearchOpts::default());
+    let out: Vec<MjdPeriod> = periods
+        .into_iter()
+        .map(|p| MjdPeriod {
+            start_mjd: p.start.value(),
+            end_mjd: p.end.value(),
+        })
+        .collect();
+    to_js(&out)
+}
+
+/// Find periods where geocentric Moon illumination is in `[kMin, kMax]`.
+#[wasm_bindgen(js_name = "moonIlluminationRange")]
+pub fn moon_illumination_range(
+    start_mjd: f64,
+    end_mjd: f64,
+    k_min: f64,
+    k_max: f64,
+) -> Result<JsValue, JsError> {
+    let window = make_window(start_mjd, end_mjd)?;
+    let periods =
+        illumination_range::<Vsop87Ephemeris>(window, k_min, k_max, PhaseSearchOpts::default());
+    let out: Vec<MjdPeriod> = periods
+        .into_iter()
+        .map(|p| MjdPeriod {
+            start_mjd: p.start.value(),
+            end_mjd: p.end.value(),
+        })
+        .collect();
+    to_js(&out)
+}
+
+// ─────────────────────────────────────────────────────────────────────────────
+// Helpers
+// ─────────────────────────────────────────────────────────────────────────────
+
+fn label_to_string(label: MoonPhaseLabel) -> String {
+    match label {
+        MoonPhaseLabel::NewMoon => "NewMoon",
+        MoonPhaseLabel::WaxingCrescent => "WaxingCrescent",
+        MoonPhaseLabel::FirstQuarter => "FirstQuarter",
+        MoonPhaseLabel::WaxingGibbous => "WaxingGibbous",
+        MoonPhaseLabel::FullMoon => "FullMoon",
+        MoonPhaseLabel::WaningGibbous => "WaningGibbous",
+        MoonPhaseLabel::LastQuarter => "LastQuarter",
+        MoonPhaseLabel::WaningCrescent => "WaningCrescent",
+    }
+    .to_string()
+}
+
+fn phase_kind_to_string(kind: siderust::calculus::lunar::phase::PhaseKind) -> String {
+    use siderust::calculus::lunar::phase::PhaseKind;
+    match kind {
+        PhaseKind::NewMoon => "NewMoon",
+        PhaseKind::FirstQuarter => "FirstQuarter",
+        PhaseKind::FullMoon => "FullMoon",
+        PhaseKind::LastQuarter => "LastQuarter",
+    }
+    .to_string()
+}

package/siderust-web/src/position.rs

@@ -0,0 +1,291 @@
+// SPDX-License-Identifier: AGPL-3.0-only
+// Copyright (C) 2026 Vallés Puig, Ramon
+
+//! Cartesian position type for WebAssembly bindings.
+
+use serde::Serialize;
+use wasm_bindgen::prelude::*;
+
+// Frame and center string constants
+pub const FRAME_ECL: &str = "EclipticMeanJ2000";
+pub const FRAME_EQ: &str = "EquatorialMeanJ2000";
+pub const FRAME_ICRS: &str = "ICRS";
+
+pub const CENTER_BARY: &str = "Barycentric";
+pub const CENTER_HELIO: &str = "Heliocentric";
+pub const CENTER_GEO: &str = "Geocentric";
+
+pub const UNIT_AU: &str = "au";
+pub const UNIT_KM: &str = "km";
+
+/// A 3D cartesian position with frame, center, and unit metadata.
+#[wasm_bindgen]
+#[derive(Clone, Serialize)]
+pub struct Position {
+    x: f64,
+    y: f64,
+    z: f64,
+    frame: String,
+    center: String,
+    unit: String,
+}
+
+#[wasm_bindgen]
+impl Position {
+    /// Create a new position with explicit coordinates, frame, center, and unit.
+    #[wasm_bindgen(constructor)]
+    pub fn new(
+        x: f64,
+        y: f64,
+        z: f64,
+        frame: String,
+        center: String,
+        unit: String,
+    ) -> Result<Position, JsError> {
+        validate_frame(&frame)?;
+        validate_center(&center)?;
+        validate_unit(&unit)?;
+        Ok(Self { x, y, z, frame, center, unit })
+    }
+
+    /// X coordinate.
+    #[wasm_bindgen(getter)]
+    pub fn x(&self) -> f64 {
+        self.x
+    }
+
+    /// Y coordinate.
+    #[wasm_bindgen(getter)]
+    pub fn y(&self) -> f64 {
+        self.y
+    }
+
+    /// Z coordinate.
+    #[wasm_bindgen(getter)]
+    pub fn z(&self) -> f64 {
+        self.z
+    }
+
+    /// Reference frame (e.g., "EclipticMeanJ2000", "ICRS").
+    #[wasm_bindgen(getter)]
+    pub fn frame(&self) -> String {
+        self.frame.clone()
+    }
+
+    /// Coordinate center (e.g., "Barycentric", "Heliocentric", "Geocentric").
+    #[wasm_bindgen(getter)]
+    pub fn center(&self) -> String {
+        self.center.clone()
+    }
+
+    /// Distance unit ("au" or "km").
+    #[wasm_bindgen(getter)]
+    pub fn unit(&self) -> String {
+        self.unit.clone()
+    }
+
+    /// Euclidean distance from origin.
+    #[wasm_bindgen]
+    pub fn magnitude(&self) -> f64 {
+        (self.x * self.x + self.y * self.y + self.z * self.z).sqrt()
+    }
+
+    /// Subtract another position to get a displacement vector.
+    #[wasm_bindgen]
+    pub fn subtract(&self, other: &Position) -> Result<Displacement, JsError> {
+        if self.frame != other.frame {
+            return Err(JsError::new(&format!(
+                "Cannot subtract positions in different frames: '{}' vs '{}'",
+                self.frame, other.frame
+            )));
+        }
+        if self.center != other.center {
+            return Err(JsError::new(&format!(
+                "Cannot subtract positions with different centers: '{}' vs '{}'",
+                self.center, other.center
+            )));
+        }
+        if self.unit != other.unit {
+            return Err(JsError::new(&format!(
+                "Cannot subtract positions with different units: '{}' vs '{}'",
+                self.unit, other.unit
+            )));
+        }
+        Ok(Displacement {
+            dx: self.x - other.x,
+            dy: self.y - other.y,
+            dz: self.z - other.z,
+            frame: self.frame.clone(),
+            unit: self.unit.clone(),
+        })
+    }
+
+    /// Add a displacement to get a new position.
+    #[wasm_bindgen(js_name = "addDisplacement")]
+    pub fn add_displacement(&self, d: &Displacement) -> Result<Position, JsError> {
+        if self.frame != d.frame {
+            return Err(JsError::new(&format!(
+                "Cannot add displacement in different frame: '{}' vs '{}'",
+                self.frame, d.frame
+            )));
+        }
+        if self.unit != d.unit {
+            return Err(JsError::new(&format!(
+                "Cannot add displacement with different unit: '{}' vs '{}'",
+                self.unit, d.unit
+            )));
+        }
+        Ok(Position {
+            x: self.x + d.dx,
+            y: self.y + d.dy,
+            z: self.z + d.dz,
+            frame: self.frame.clone(),
+            center: self.center.clone(),
+            unit: self.unit.clone(),
+        })
+    }
+}
+
+impl Position {
+    /// Internal constructor without validation (for Rust callers).
+    pub(crate) fn new_internal(
+        x: f64,
+        y: f64,
+        z: f64,
+        frame: &str,
+        center: &str,
+        unit: &str,
+    ) -> Self {
+        Self {
+            x,
+            y,
+            z,
+            frame: frame.to_string(),
+            center: center.to_string(),
+            unit: unit.to_string(),
+        }
+    }
+}
+
+/// A displacement vector (difference between two positions).
+#[wasm_bindgen]
+#[derive(Clone, Serialize)]
+pub struct Displacement {
+    dx: f64,
+    dy: f64,
+    dz: f64,
+    frame: String,
+    unit: String,
+}
+
+#[wasm_bindgen]
+impl Displacement {
+    /// Create a new displacement.
+    #[wasm_bindgen(constructor)]
+    pub fn new(
+        dx: f64,
+        dy: f64,
+        dz: f64,
+        frame: String,
+        unit: String,
+    ) -> Result<Displacement, JsError> {
+        validate_frame(&frame)?;
+        validate_unit(&unit)?;
+        Ok(Self { dx, dy, dz, frame, unit })
+    }
+
+    /// X component of displacement.
+    #[wasm_bindgen(getter)]
+    pub fn dx(&self) -> f64 {
+        self.dx
+    }
+
+    /// Y component of displacement.
+    #[wasm_bindgen(getter)]
+    pub fn dy(&self) -> f64 {
+        self.dy
+    }
+
+    /// Z component of displacement.
+    #[wasm_bindgen(getter)]
+    pub fn dz(&self) -> f64 {
+        self.dz
+    }
+
+    /// Reference frame.
+    #[wasm_bindgen(getter)]
+    pub fn frame(&self) -> String {
+        self.frame.clone()
+    }
+
+    /// Distance unit.
+    #[wasm_bindgen(getter)]
+    pub fn unit(&self) -> String {
+        self.unit.clone()
+    }
+
+    /// Euclidean magnitude.
+    #[wasm_bindgen]
+    pub fn magnitude(&self) -> f64 {
+        (self.dx * self.dx + self.dy * self.dy + self.dz * self.dz).sqrt()
+    }
+
+    /// Scale displacement by a factor.
+    #[wasm_bindgen]
+    pub fn scale(&self, factor: f64) -> Displacement {
+        Displacement {
+            dx: self.dx * factor,
+            dy: self.dy * factor,
+            dz: self.dz * factor,
+            frame: self.frame.clone(),
+            unit: self.unit.clone(),
+        }
+    }
+
+    /// Add another displacement.
+    #[wasm_bindgen]
+    pub fn add(&self, other: &Displacement) -> Result<Displacement, JsError> {
+        if self.frame != other.frame || self.unit != other.unit {
+            return Err(JsError::new(
+                "Cannot add displacements with different frame/unit",
+            ));
+        }
+        Ok(Displacement {
+            dx: self.dx + other.dx,
+            dy: self.dy + other.dy,
+            dz: self.dz + other.dz,
+            frame: self.frame.clone(),
+            unit: self.unit.clone(),
+        })
+    }
+}
+
+fn validate_frame(f: &str) -> Result<(), JsError> {
+    match f {
+        FRAME_ECL | FRAME_EQ | FRAME_ICRS | "ICRF" | "EquatorialMeanOfDate" | "EquatorialTrueOfDate" => Ok(()),
+        _ => Err(JsError::new(&format!(
+            "Unknown frame '{}'. Valid: EclipticMeanJ2000, EquatorialMeanJ2000, ICRS, ICRF",
+            f
+        ))),
+    }
+}
+
+fn validate_center(c: &str) -> Result<(), JsError> {
+    match c {
+        CENTER_BARY | CENTER_HELIO | CENTER_GEO => Ok(()),
+        _ => Err(JsError::new(&format!(
+            "Unknown center '{}'. Valid: Barycentric, Heliocentric, Geocentric",
+            c
+        ))),
+    }
+}
+
+fn validate_unit(u: &str) -> Result<(), JsError> {
+    match u {
+        UNIT_AU | UNIT_KM => Ok(()),
+        _ => Err(JsError::new(&format!(
+            "Unknown unit '{}'. Valid: au, km",
+            u
+        ))),
+    }
+}

package/siderust-web/src/star.rs

@@ -0,0 +1,178 @@
+// SPDX-License-Identifier: AGPL-3.0-only
+// Copyright (C) 2026 Vallés Puig, Ramon
+
+//! Star — catalog and custom stars for altitude/azimuth/coordinate queries.
+
+use wasm_bindgen::prelude::*;
+
+use qtty::length::nominal::SolarRadiuses;
+use qtty::*;
+use siderust::bodies::{self, Star as SiderustStar};
+use siderust::coordinates::centers::Geocentric;
+use siderust::coordinates::frames::EquatorialMeanJ2000;
+use siderust::coordinates::spherical;
+use siderust::targets::CoordinateWithPM;
+use siderust::time::JulianDate;
+use std::borrow::Cow;
+
+// ─────────────────────────────────────────────────────────────────────────────
+// Star class
+// ─────────────────────────────────────────────────────────────────────────────
+
+/// A star with physical parameters and sky coordinates.
+#[wasm_bindgen]
+pub struct Star {
+    pub(crate) inner: SiderustStar<'static>,
+}
+
+#[wasm_bindgen]
+impl Star {
+    /// Create a custom star.
+    ///
+    /// @param name        — Display name.
+    /// @param distanceLy  — Distance in light-years.
+    /// @param massSolar   — Mass in solar masses (M☉).
+    /// @param radiusSolar — Radius in solar radii (R☉).
+    /// @param luminositySolar — Luminosity in solar luminosities (L☉).
+    /// @param raDeg       — Right ascension at J2000.0 in degrees.
+    /// @param decDeg      — Declination at J2000.0 in degrees.
+    #[wasm_bindgen(constructor)]
+    pub fn new(
+        name: &str,
+        distance_ly: f64,
+        mass_solar: f64,
+        radius_solar: f64,
+        luminosity_solar: f64,
+        ra_deg: f64,
+        dec_deg: f64,
+    ) -> Result<Star, JsError> {
+        for (label, v) in [
+            ("distanceLy", distance_ly),
+            ("massSolar", mass_solar),
+            ("radiusSolar", radius_solar),
+            ("luminositySolar", luminosity_solar),
+            ("raDeg", ra_deg),
+            ("decDeg", dec_deg),
+        ] {
+            if !v.is_finite() {
+                return Err(JsError::new(&format!(
+                    "{label} must be finite (not NaN or ±infinity)"
+                )));
+            }
+        }
+
+        let pos = spherical::Position::<Geocentric, EquatorialMeanJ2000, LightYear>::new(
+            Degrees::new(ra_deg),
+            Degrees::new(dec_deg),
+            LightYears::new(distance_ly),
+        );
+        let target =
+            CoordinateWithPM::<spherical::Position<Geocentric, EquatorialMeanJ2000, LightYear>>::new_static(
+                pos,
+                JulianDate::J2000,
+            );
+
+        let star = SiderustStar::new(
+            Cow::<'static, str>::Owned(name.to_string()),
+            LightYears::new(distance_ly),
+            SolarMasses::new(mass_solar),
+            SolarRadiuses::new(radius_solar),
+            SolarLuminosities::new(luminosity_solar),
+            target,
+        );
+
+        Ok(Self { inner: star })
+    }
+
+    /// Look up a star from the built-in catalog by name.
+    ///
+    /// Supported names (case-insensitive): Vega, Sirius, Polaris, Canopus,
+    /// Arcturus, Rigel, Betelgeuse, Procyon, Aldebaran, Altair.
+    pub fn catalog(name: &str) -> Result<Star, JsError> {
+        let star: &SiderustStar<'static> = match name.to_uppercase().as_str() {
+            "VEGA" => &bodies::VEGA,
+            "SIRIUS" => &bodies::SIRIUS,
+            "POLARIS" => &bodies::POLARIS,
+            "CANOPUS" => &bodies::CANOPUS,
+            "ARCTURUS" => &bodies::ARCTURUS,
+            "RIGEL" => &bodies::RIGEL,
+            "BETELGEUSE" => &bodies::BETELGEUSE,
+            "PROCYON" => &bodies::PROCYON,
+            "ALDEBARAN" => &bodies::ALDEBARAN,
+            "ALTAIR" => &bodies::ALTAIR,
+            _ => {
+                return Err(JsError::new(&format!(
+                    "Unknown catalog star: \"{name}\". Available: Vega, Sirius, Polaris, Canopus, Arcturus, Rigel, Betelgeuse, Procyon, Aldebaran, Altair."
+                )));
+            }
+        };
+        Ok(Self {
+            inner: star.clone(),
+        })
+    }
+
+    // ── accessors ────────────────────────────────────────────────────
+
+    /// Star name.
+    #[wasm_bindgen(getter)]
+    pub fn name(&self) -> String {
+        self.inner.name.to_string()
+    }
+
+    /// Distance in light-years.
+    #[wasm_bindgen(getter, js_name = "distanceLy")]
+    pub fn distance_ly(&self) -> f64 {
+        self.inner.distance.value()
+    }
+
+    /// Mass in solar masses (M☉).
+    #[wasm_bindgen(getter, js_name = "massSolar")]
+    pub fn mass_solar(&self) -> f64 {
+        self.inner.mass.value()
+    }
+
+    /// Radius in solar radii (R☉).
+    #[wasm_bindgen(getter, js_name = "radiusSolar")]
+    pub fn radius_solar(&self) -> f64 {
+        self.inner.radius.value()
+    }
+
+    /// Luminosity in solar luminosities (L☉).
+    #[wasm_bindgen(getter, js_name = "luminositySolar")]
+    pub fn luminosity_solar(&self) -> f64 {
+        self.inner.luminosity.value()
+    }
+
+    /// Right ascension at J2000.0 in degrees.
+    #[wasm_bindgen(getter, js_name = "raDeg")]
+    pub fn ra_deg(&self) -> f64 {
+        self.inner.coordinate.get_position().azimuth.value()
+    }
+
+    /// Declination at J2000.0 in degrees.
+    #[wasm_bindgen(getter, js_name = "decDeg")]
+    pub fn dec_deg(&self) -> f64 {
+        self.inner.coordinate.get_position().polar.value()
+    }
+
+    /// Human-readable representation.
+    pub fn format(&self) -> String {
+        format!(
+            "Star({}, d={:.1} ly, RA={:.4}°, Dec={:.4}°)",
+            self.inner.name,
+            self.inner.distance.value(),
+            self.inner.coordinate.get_position().azimuth.value(),
+            self.inner.coordinate.get_position().polar.value(),
+        )
+    }
+}
+
+/// List all catalog star names.
+#[wasm_bindgen(js_name = "listCatalogStars")]
+pub fn list_catalog_stars() -> JsValue {
+    let names: Vec<&str> = vec![
+        "Vega", "Sirius", "Polaris", "Canopus", "Arcturus", "Rigel", "Betelgeuse", "Procyon",
+        "Aldebaran", "Altair",
+    ];
+    serde_wasm_bindgen::to_value(&names).unwrap()
+}