Program Listing for File spice_interface.h

↰ Return to documentation for file (include/lupnt/physics/spice_interface.h)

#pragma once

#include <string.h>

#include <map>
#include <tuple>

#include "lupnt/core/constants.h"
#include "lupnt/physics/cheby.h"
#include "lupnt/physics/frame_converter.h"

namespace lupnt {

  // map from CoordSytem to string
  const std::map<Frame, std::string> frame2string = {
      {Frame::ITRF, "ITRF93"},   {Frame::GCRF, "J2000"},          {Frame::MOON_PA, "MOON_PA"},
      {Frame::MOON_CI, "J2000"}, {Frame::MARS_FIXED, "IAU_MARS"}, {Frame::VENUS_FIXED, "IAU_VENUS"},
  };

  const std::map<Time, std::string> time2string = {
      {Time::UT1, "UT1"},     {Time::UTC, "UTC"},      {Time::TAI, "TAI"}, {Time::TDB, "TDB"},
      {Time::TT, "TT"},       {Time::TCG, "TCG"},      {Time::TCB, "TCB"}, {Time::GPS, "GPS"},
      {Time::JD_TT, "JDTDT"}, {Time::JD_TDB, "JDTDB"},
  };

  namespace spice {
    static segment_t* cheby_s;
    static long cheby_n;
    void LoadSpiceKernel(void);
    void ExtractPckCoeffs(void);
    Mat6d GetFrameConversionMat(Real t_tai, Frame from_frame, Frame to_frame);

    Real String2TDB(const std::string& str);
    Real String2TAI(const std::string& str);

    std::string TAItoStringUTC(Real t_tai, int prec);
    std::string TDBtoStringUTC(Real t_tdb, int prec);

    Real ConvertTime(Real t, Time from_time, Time to_time);

    Vec6 GetBodyPosVel(const Real t_tai, NaifId center, NaifId target);
    MatX6 GetBodyPosVel(const VecX& t_tai, NaifId center, NaifId target);

    Vec3d GetBodyPosSpice(Real t_tai, NaifId obs, NaifId target, Frame refFrame = Frame::GCRF,
                          const std::string& abCorrection = "NONE");
    Vec6d GetBodyPosVelSpice(Real t_tai, NaifId obs, NaifId target, Frame refFrame = Frame::GCRF,
                             const std::string& abCorrection = "NONE");
  }  // namespace spice
}  // namespace lupnt