.. _program_listing_file_data_kernels.h:

Program Listing for File kernels.h
==================================

|exhale_lsh| :ref:`Return to documentation for file <file_data_kernels.h>` (``data/kernels.h``)

.. |exhale_lsh| unicode:: U+021B0 .. UPWARDS ARROW WITH TIP LEFTWARDS

.. code-block:: cpp

   /*
   Greg Miller (gmiller@gregmiller.net) 2022
   Released as public domain
   http://www.celestialprogramming.com/
   
   Class to read binary versions of JPL's Development Ephemeris.  Files in
   the propper format can be obtained from:
   ftp://ssd.jpl.nasa.gov/pub/eph/planets/Linux
   
   #    Properties       Units          Center Description
   0    x,y,z            km             SSB    Mercury
   1    x,y,z            km             SSB    Venus
   2    x,y,z            km             SSB    Earth-Moon barycenter
   3    x,y,z            km             SSB    Mars
   4    x,y,z            km             SSB    Jupiter
   5    x,y,z            km             SSB    Saturn
   6    x,y,z            km             SSB    Uranus
   7    x,y,z            km             SSB    Neptune
   8    x,y,z            km             SSB    Pluto
   9    x,y,z            km             Earth  Moon (geocentric)
   10   x,y,z            km             SSB    Sun
   11   dPsi,dEps        radians               Earth Nutations in lon and obliquity
   12   phi,theta,psi    radians               Lunar mantle libration
   13   Ox,Oy,Oz         radians/day           Lunar mantle angular velocity
   14   t                seconds               TT-TDB (at geocenter)
   
   Example: (prints x coordinate of venus using first JD available)
   
   JPLDE.DE de = new JPLDE.DE(@"E:\Astronomy\_Ephemeris\JPLDEBinaries\jpleph.405");
   Console.WriteLine(de.getPlanet(1, de.getHeader().jdStart)[0]);
   
   24857048.3412405
   */
   #pragma once
   
   #include <fstream>
   #include <iostream>
   #include <sstream>
   #include <string>
   #include <unordered_map>
   #include <vector>
   
   #include "lupnt/core/constants.h"
   #include "lupnt/interfaces/spice.h"
   
   #define MAXCOEFF 1020
   
   namespace lupnt {
   
     double GetTtTdbDifference(double t_tdb);
   
     Vec6 GetLunarMantleData(Real t_tdb, bool compute_vel = true);
   
     MatX6 GetLunarMantleData(VecX t_tdb, bool compute_vel = true);
   
     Vec6 GetBodyPosVel(Real t_tdb, BodyId center, BodyId target, Frame frame);
   
     Vec6 GetBodyPosVel(Real t_tdb, BodyId target, Frame frame);
   
     Vec6 GetBodyPosVel(Real t_tdb, BodyId center, BodyId target, Frame frame, const UnitSystem& units,
                        CoordinateScale scale = CoordinateScale::TDB);
   
     Vec6 GetBodyPosVel(Real t_tdb, BodyId target, Frame frame, const UnitSystem& units,
                        CoordinateScale scale = CoordinateScale::TDB);
   
     MatX6 GetBodyPosVel(const VecX& t_tdb, BodyId target, Frame frame);
   
     MatX6 GetBodyPosVel(const VecX& t_tdb, BodyId center, BodyId target, Frame frame);
   
     MatX6 GetBodyPosVel(const VecX& t_tdb, BodyId target, Frame frame, const UnitSystem& units,
                         CoordinateScale scale = CoordinateScale::TDB);
   
     MatX6 GetBodyPosVel(const VecX& t_tdb, BodyId center, BodyId target, Frame frame,
                         const UnitSystem& units, CoordinateScale scale = CoordinateScale::TDB);
   
     Vec3 GetBodyPos(Real t_tdb, BodyId center, BodyId target, Frame frame);
   
     Vec3 GetBodyPos(Real t_tdb, BodyId target, Frame frame);
   
     Vec3 GetBodyPos(Real t_tdb, BodyId center, BodyId target, Frame frame, const UnitSystem& units,
                     CoordinateScale scale = CoordinateScale::TDB);
   
     Vec3 GetBodyPos(Real t_tdb, BodyId target, Frame frame, const UnitSystem& units,
                     CoordinateScale scale = CoordinateScale::TDB);
   
     MatX3 GetBodyPos(const VecX& t_tdb, BodyId target, Frame frame);
   
     MatX3 GetBodyPos(const VecX& t_tdb, BodyId center, BodyId target, Frame frame);
   
     MatX3 GetBodyPos(const VecX& t_tdb, BodyId target, Frame frame, const UnitSystem& units,
                      CoordinateScale scale = CoordinateScale::TDB);
   
     MatX3 GetBodyPos(const VecX& t_tdb, BodyId center, BodyId target, Frame frame,
                      const UnitSystem& units, CoordinateScale scale = CoordinateScale::TDB);
   
   }  // namespace lupnt