EphemerisLagrangeVar.hpp

浏览该文件的文档.

 
 
#pragma once
 
#include "AstGlobal.h"
#include "AstCore/Ephemeris.hpp"
#include "AstCore/TimePoint.hpp"
#include "AstCore/Frame.hpp"
#include "AstMath/Vector.hpp"
#include <vector>
 
AST_NAMESPACE_BEGIN
 
class AST_CORE_API EphemerisLagrangeVar: public Ephemeris
{
public:
    EphemerisLagrangeVar() = default;
    EphemerisLagrangeVar(Frame* frame) : frame_(frame){}
    ~EphemerisLagrangeVar() = default;
public:
    Frame* getFrame() const override;
    errc_t getPos(const TimePoint& tp, Vector3d& pos) const override;
    errc_t getPosVel(const TimePoint& tp, Vector3d& pos, Vector3d& vel) const override;
    errc_t getInterval(TimeInterval& interval) const override;
public:
    void setFrame(Frame* frame){frame_ = frame;}
    void setTimes(const std::vector<double>& times);
    void setTimes(std::vector<double>&& times);
    void setPositions(const std::vector<Vector3d>& positions){positions_ = positions;}
    void setPositions(std::vector<Vector3d>&& positions){positions_ = std::move(positions);}
    void setVelocities(const std::vector<Vector3d>& velocities){velocities_ = velocities;}
    void setVelocities(std::vector<Vector3d>&& velocities){velocities_ = std::move(velocities);}
    void setEpoch(const TimePoint& epoch){epoch_ = epoch;}
    void setInterpolateOrder(int order){interpolateOrder_ = order;}
public:
    size_t size() const{return times_.size();}
private:
    int findIndex(double delta) const;
protected:
    SharedPtr<Frame>      frame_;                   
    std::vector<double>   times_;                   
    std::vector<Vector3d> positions_;               
    std::vector<Vector3d> velocities_;              
    TimePoint             epoch_{};                 
    double                averageStep_{60};         
    int                   interpolateOrder_{5};     
};
 
 
AST_NAMESPACE_END