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 Time<TT, JD> &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 Time<TT, JD> &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 Time<TT, JD> &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::EclipticMeanJ2000<qtty::AstronomicalUnit>
mars_heliocentric(const Time<TT, JD> &jd) {
  siderust_cartesian_pos_t out;
  check_status(siderust_vsop87_mars_heliocentric(jd.value(), &out), "ephemeris::mars_heliocentric");
  return cartesian::position::EclipticMeanJ2000<qtty::AstronomicalUnit>::from_c(out);
}
 
inline cartesian::position::HelioBarycentric<qtty::AstronomicalUnit>
mars_barycentric(const Time<TT, JD> &jd) {
  siderust_cartesian_pos_t out;
  check_status(siderust_vsop87_mars_barycentric(jd.value(), &out), "ephemeris::mars_barycentric");
  return cartesian::position::HelioBarycentric<qtty::AstronomicalUnit>::from_c(out);
}
 
inline cartesian::position::EclipticMeanJ2000<qtty::AstronomicalUnit>
venus_heliocentric(const Time<TT, JD> &jd) {
  siderust_cartesian_pos_t out;
  check_status(siderust_vsop87_venus_heliocentric(jd.value(), &out),
               "ephemeris::venus_heliocentric");
  return cartesian::position::EclipticMeanJ2000<qtty::AstronomicalUnit>::from_c(out);
}
 
inline cartesian::position::EclipticMeanJ2000<qtty::AstronomicalUnit>
mercury_heliocentric(const Time<TT, JD> &jd) {
  siderust_cartesian_pos_t out;
  check_status(siderust_vsop87_mercury_heliocentric(jd.value(), &out),
               "ephemeris::mercury_heliocentric");
  return cartesian::position::EclipticMeanJ2000<qtty::AstronomicalUnit>::from_c(out);
}
 
inline cartesian::position::HelioBarycentric<qtty::AstronomicalUnit>
mercury_barycentric(const Time<TT, JD> &jd) {
  siderust_cartesian_pos_t out;
  check_status(siderust_vsop87_mercury_barycentric(jd.value(), &out),
               "ephemeris::mercury_barycentric");
  return cartesian::position::HelioBarycentric<qtty::AstronomicalUnit>::from_c(out);
}
 
inline cartesian::position::HelioBarycentric<qtty::AstronomicalUnit>
venus_barycentric(const Time<TT, JD> &jd) {
  siderust_cartesian_pos_t out;
  check_status(siderust_vsop87_venus_barycentric(jd.value(), &out), "ephemeris::venus_barycentric");
  return cartesian::position::HelioBarycentric<qtty::AstronomicalUnit>::from_c(out);
}
 
inline cartesian::position::EclipticMeanJ2000<qtty::AstronomicalUnit>
jupiter_heliocentric(const Time<TT, JD> &jd) {
  siderust_cartesian_pos_t out;
  check_status(siderust_vsop87_jupiter_heliocentric(jd.value(), &out),
               "ephemeris::jupiter_heliocentric");
  return cartesian::position::EclipticMeanJ2000<qtty::AstronomicalUnit>::from_c(out);
}
 
inline cartesian::position::HelioBarycentric<qtty::AstronomicalUnit>
jupiter_barycentric(const Time<TT, JD> &jd) {
  siderust_cartesian_pos_t out;
  check_status(siderust_vsop87_jupiter_barycentric(jd.value(), &out),
               "ephemeris::jupiter_barycentric");
  return cartesian::position::HelioBarycentric<qtty::AstronomicalUnit>::from_c(out);
}
 
inline cartesian::position::EclipticMeanJ2000<qtty::AstronomicalUnit>
saturn_heliocentric(const Time<TT, JD> &jd) {
  siderust_cartesian_pos_t out;
  check_status(siderust_vsop87_saturn_heliocentric(jd.value(), &out),
               "ephemeris::saturn_heliocentric");
  return cartesian::position::EclipticMeanJ2000<qtty::AstronomicalUnit>::from_c(out);
}
 
inline cartesian::position::HelioBarycentric<qtty::AstronomicalUnit>
saturn_barycentric(const Time<TT, JD> &jd) {
  siderust_cartesian_pos_t out;
  check_status(siderust_vsop87_saturn_barycentric(jd.value(), &out),
               "ephemeris::saturn_barycentric");
  return cartesian::position::HelioBarycentric<qtty::AstronomicalUnit>::from_c(out);
}
 
inline cartesian::position::EclipticMeanJ2000<qtty::AstronomicalUnit>
uranus_heliocentric(const Time<TT, JD> &jd) {
  siderust_cartesian_pos_t out;
  check_status(siderust_vsop87_uranus_heliocentric(jd.value(), &out),
               "ephemeris::uranus_heliocentric");
  return cartesian::position::EclipticMeanJ2000<qtty::AstronomicalUnit>::from_c(out);
}
 
inline cartesian::position::HelioBarycentric<qtty::AstronomicalUnit>
uranus_barycentric(const Time<TT, JD> &jd) {
  siderust_cartesian_pos_t out;
  check_status(siderust_vsop87_uranus_barycentric(jd.value(), &out),
               "ephemeris::uranus_barycentric");
  return cartesian::position::HelioBarycentric<qtty::AstronomicalUnit>::from_c(out);
}
 
inline cartesian::position::EclipticMeanJ2000<qtty::AstronomicalUnit>
neptune_heliocentric(const Time<TT, JD> &jd) {
  siderust_cartesian_pos_t out;
  check_status(siderust_vsop87_neptune_heliocentric(jd.value(), &out),
               "ephemeris::neptune_heliocentric");
  return cartesian::position::EclipticMeanJ2000<qtty::AstronomicalUnit>::from_c(out);
}
 
inline cartesian::position::HelioBarycentric<qtty::AstronomicalUnit>
neptune_barycentric(const Time<TT, JD> &jd) {
  siderust_cartesian_pos_t out;
  check_status(siderust_vsop87_neptune_barycentric(jd.value(), &out),
               "ephemeris::neptune_barycentric");
  return cartesian::position::HelioBarycentric<qtty::AstronomicalUnit>::from_c(out);
}
 
inline cartesian::position::MoonGeocentric<qtty::Kilometer>
moon_geocentric(const Time<TT, JD> &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