demo.cpp

Go to the documentation of this file.

 
#include <cmath>
#include <cstdio>
#include <siderust/siderust.hpp>
 
int main() {
    using namespace siderust;
    using namespace siderust::frames;
    using namespace qtty::literals;
 
    std::printf("=== siderust-cpp demo ===\n\n");
 
    // --- Time ---
    auto jd = JulianDate::from_utc({2026, 7, 15, 22, 0, 0});
    std::printf("JD for 2026-07-15 22:00 UTC: %.6f\n", jd.value());
    std::printf("Julian centuries since J2000: %.10f\n", jd.julian_centuries());
 
    auto mjd = MJD::from_jd(jd);
    std::printf("MJD: %.6f\n\n", mjd.value());
 
    // --- Observatory ---
    auto obs = ROQUE_DE_LOS_MUCHACHOS;
    std::printf("Roque de los Muchachos: lon=%.4f  lat=%.4f  h=%.0f m\n\n",
                obs.lon.value(), obs.lat.value(), obs.height.value());
 
    // --- Sun altitude ---
    qtty::Radian sun_alt = sun::altitude_at(obs, mjd);
    std::printf("Sun altitude: %.4f rad (%.2f deg)\n\n",
                sun_alt.value(), sun_alt.to<qtty::Degree>().value());
 
    // --- Star catalog ---
    const auto& vega = VEGA;
    std::printf("Star: %s, d=%.2f ly, L=%.2f Lsun\n",
                vega.name().c_str(), vega.distance_ly(),
                vega.luminosity_solar());
 
    // --- Star altitude ---
    qtty::Radian star_alt = star_altitude::altitude_at(vega, obs, mjd);
    std::printf("Vega altitude: %.4f rad (%.2f deg)\n\n",
                star_alt.value(), star_alt.to<qtty::Degree>().value());
 
    // =================================================================
    // TYPED COORDINATE & EPHEMERIS API
    // =================================================================
    std::printf("--- Typed Coordinate API ---\n\n");
 
    // Compile-time typed ICRS direction
    spherical::direction::ICRS vega_icrs(279.23473, 38.78369);
 
    // Template-targeted transform: ICRS → EclipticMeanJ2000
    auto ecl = vega_icrs.to_frame<EclipticMeanJ2000>(jd);
    std::printf("Typed ICRS (%.4f, %.4f) -> EclMeanJ2000 (%.4f, %.4f)\n",
                vega_icrs.ra().value(), vega_icrs.dec().value(),
                ecl.lon().value(), ecl.lat().value());
 
    // Shorthand .to<>() syntax
    auto eq_j2000 = vega_icrs.to<EquatorialMeanJ2000>(jd);
    std::printf("Typed ICRS -> EquatorialJ2000 (%.4f, %.4f)\n",
                eq_j2000.ra().value(), eq_j2000.dec().value());
 
    // Horizontal transform
    auto hor = vega_icrs.to_horizontal(jd, obs);
    std::printf("Typed Horizontal: az=%.4f  alt=%.4f deg\n\n",
                hor.az().value(), hor.al().value());
 
    // Roundtrip: ICRS → Ecliptic → ICRS
    auto back = ecl.to_frame<ICRS>(jd);
    std::printf("Roundtrip: (%.6f, %.6f) -> (%.6f, %.6f) -> (%.6f, %.6f)\n",
                vega_icrs.ra().value(), vega_icrs.dec().value(),
                ecl.lon().value(), ecl.lat().value(),
                back.ra().value(), back.dec().value());
 
    // qtty unit-safe angle conversion
    qtty::Radian ra_rad = vega_icrs.ra().to<qtty::Radian>();
    std::printf("Vega RA: %.6f deg = %.6f rad\n\n", vega_icrs.ra().value(), ra_rad.value());
 
    // --- Typed Ephemeris ---
    std::printf("--- Typed Ephemeris ---\n\n");
 
    auto earth = ephemeris::earth_heliocentric(jd);
    std::printf("Earth heliocentric (typed AU): (%.8f, %.8f, %.8f)\n",
                earth.x().value(), earth.y().value(), earth.z().value());
 
    // Unit conversion: AU → km
    qtty::Kilometer x_km = earth.comp_x.to<qtty::Kilometer>();
    qtty::Kilometer y_km = earth.comp_y.to<qtty::Kilometer>();
    std::printf("Earth heliocentric (km): (%.2f, %.2f, ...)\n\n",
                x_km.value(), y_km.value());
 
    auto moon = ephemeris::moon_geocentric(jd);
    std::printf("Moon geocentric (typed km): (%.2f, %.2f, %.2f)\n",
                moon.x().value(), moon.y().value(), moon.z().value());
    auto moon_r = std::sqrt(moon.x().value() * moon.x().value() + moon.y().value() * moon.y().value() + moon.z().value() * moon.z().value());
    std::printf("Moon distance: %.2f km\n\n", moon_r);
 
    // --- Planets ---
    auto mars_data = MARS;
    std::printf("Mars: mass=%.4e kg, radius=%.2f km\n",
                mars_data.mass_kg, mars_data.radius_km);
    std::printf("  orbit: a=%.6f AU, e=%.6f\n\n",
                mars_data.orbit.semi_major_axis_au,
                mars_data.orbit.eccentricity);
 
    // --- Night periods (sun below -18°) ---
    auto night_start = mjd;
    auto night_end   = mjd + 1.0;
    auto nights      = sun::below_threshold(obs, night_start, night_end, -18.0_deg);
    std::printf("Astronomical night periods (sun < -18 deg):\n");
    for (auto& p : nights) {
        std::printf("  MJD %.6f – %.6f  (%.2f hours)\n",
                    p.start_mjd(), p.end_mjd(),
                    p.duration_days() * 24.0);
    }
 
    std::printf("\nDone.\n");
    return 0;
}