JplDe.hpp

浏览该文件的文档.

 
 
#pragma once
 
#include "AstGlobal.h"
#include "AstCore/RunTime.hpp"
#include "AstUtil/OrderedMap.hpp"   // for OrderedMap
#include <stdint.h>                 // for uint32_t
#include <cstdio>                   // for FILE
#include <string>                   // for std::string
#include <mutex>                    // for std::mutex
 
AST_NAMESPACE_BEGIN
 
 
class AST_CORE_API JplDe
{
public:
    // 与GMAT的DeFile类中的天体编号一致
    enum EDataCode
    {
        eNotAvailable = -1,     
        eMercury = 0,           
        eVenus = 1,             
        eEarth = 2,             
        eMars = 3,              
        eJupiter = 4,           
        eSaturn = 5,            
        eUranus = 6,            
        eNeptune = 7,           
        ePluto = 8,             
        eMoon = 9,              
        eSun = 10,              
        eSSBarycenter = 11,     
        eEMBarycenter = 12,     
    };
 
    using ConstantList = std::vector<std::pair<std::string, double>>;
    using ConstantMap = OrderedMap<std::string, double>;
 
    JplDe();
    ~JplDe();
 
    JplDe(const JplDe&) = delete;
    JplDe& operator=(const JplDe&) = delete;
 
    errc_t open(const char* filepath);
 
 
    errc_t openDefault();
 
    bool isOpen() const;
 
 
    void close();
 
 
    errc_t getPosVelICRF(
        const TimePoint& time,
        EDataCode target,
        EDataCode referenceBody,
        Vector3d& pos,
        Vector3d& vel
    );
 
    errc_t getPosVelICRF(
        const TimePoint& time,
        EDataCode target,
        EDataCode referenceBody,
        Vector3d& pos,
        Vector3d* vel
    );
 
 
    errc_t getPosICRF(
        const TimePoint& time,
        EDataCode target,
        EDataCode referenceBody,
        Vector3d& pos
    );
    
 
    errc_t  getNutation(
        const TimePoint& time,
        double& nutLong,
        double& nutObl
    );
 
    
    errc_t      getLibration(
        const TimePoint& time,
        Euler& ang,
        Euler& angRate
    );
    
 
    errc_t       getLibration(
        const TimePoint& time,
        Euler& ang
    );
    
    uint32_t getEphemVersion() const{ return ephemVersion_; }
 
 
    uint32_t getNumConstants() const { return numConstants_; }
 
 
    const ConstantList& getConstants() const;
 
    errc_t getConstant(size_t index, double& value, std::string* name = nullptr) const;
 
 
    errc_t getConstant(const std::string& name, double& value) const;
 
    double getConstant(const std::string& name) const;
 
    double getAU() const{ return au_ * 1000; }
 
    double getEarthMoonMassRatio() const{ return emMassRatio_; }
 
    double getBodyGM(EDataCode bodyCode) const;
    errc_t getBodyGM(EDataCode bodyCode, double& value) const;
 
    errc_t getInterval(TimeInterval& interval) const;
 
private:
    errc_t         loadConstants() const;
 
    errc_t         readDataBlock(size_t idx);
    errc_t         getStateTDB(const JulianDate& jdTDB, int dataid, double pos[], double vel[]);
    errc_t         getStateTDB(const JulianDate& jdTDB, int datalist[11], double statelist[11][6]);
 
    errc_t         getState(const TimePoint& time, int dataid, double pos[], double vel[]);
    errc_t         getState(const TimePoint& time, int datalist[11], double statelist[11][6]);
protected:
 
    bool        isSameEndian_;        
    uint32_t    ephemVersion_{ 0 };   
    uint32_t    numConstants_{};      
    uint32_t    numDataBlock_{0};     
    uint32_t    numCoeff_{};          
    uint32_t    ipt_[15][3]{};        
    double      ephemStart_{0};       
    double      ephemEnd_{0};         
    double      ephemStep_{0};        
    double      au_{};                
    double      emMassRatio_{};       
    FILE*       deFile_{ NULL };      
    double**    dataBlocks_{NULL};    
    mutable ConstantMap constants_{}; 
    mutable std::mutex  mutex_{};     
};
 
 
 
AST_NAMESPACE_END