ephemeris.hpp

Go to the documentation of this file.

#pragma once
 
#include "coordinates.hpp"
#include "ffi_core.hpp"
#include "time.hpp"
 
namespace siderust {
 
namespace ephemeris {
 
// ============================================================================
// Typed API (new) — return cartesian::Position<C, F, U>
// ============================================================================
 
inline cartesian::position::HelioBarycentric<qtty::AstronomicalUnit> sun_barycentric(const JulianDate& jd) {
    siderust_cartesian_pos_t out;
    check_status(siderust_vsop87_sun_barycentric(jd.value(), &out),
                 "ephemeris::sun_barycentric");
    return cartesian::position::HelioBarycentric<qtty::AstronomicalUnit>::from_c(out);
}
 
inline cartesian::position::GeoBarycentric<qtty::AstronomicalUnit> earth_barycentric(const JulianDate& jd) {
    siderust_cartesian_pos_t out;
    check_status(siderust_vsop87_earth_barycentric(jd.value(), &out),
                 "ephemeris::earth_barycentric");
    return cartesian::position::GeoBarycentric<qtty::AstronomicalUnit>::from_c(out);
}
 
inline cartesian::position::EclipticMeanJ2000<qtty::AstronomicalUnit> earth_heliocentric(const JulianDate& jd) {
    siderust_cartesian_pos_t out;
    check_status(siderust_vsop87_earth_heliocentric(jd.value(), &out),
                 "ephemeris::earth_heliocentric");
    return cartesian::position::EclipticMeanJ2000<qtty::AstronomicalUnit>::from_c(out);
}
 
inline cartesian::position::MoonGeocentric<qtty::Kilometer> moon_geocentric(const JulianDate& jd) {
    siderust_cartesian_pos_t out;
    check_status(siderust_vsop87_moon_geocentric(jd.value(), &out),
                 "ephemeris::moon_geocentric");
    return cartesian::position::MoonGeocentric<qtty::Kilometer>::from_c(out);
}
 
} // namespace ephemeris
 
} // namespace siderust