BlockSRP.hpp

浏览该文件的文档.

 
 
#pragma once
 
#include "AstGlobal.h"
#include "AstCore/BlockDerivative.hpp"
#include "AstMath/Vector.hpp"
#include "AstCore/ForceModel.hpp"       // for ESunPosition
#include "AstCore/EclipseCalculator.hpp"
 
AST_NAMESPACE_BEGIN
 
class CelestialBody;
class Frame;
class SimTime;
 
 
constexpr double kSolarPressureAt1AU_FromLuminosity = kSunLuminosity / (4.0 * kPI * kAU * kAU * kLightSpeed);
 
constexpr double kSolarPressureAt1AU_FromSolarConstant = kSolarConstant / kLightSpeed;
 
constexpr double kSolarPressureAt1AU = kSolarPressureAt1AU_FromLuminosity;  
 
 
 
class AST_CORE_API BlockSRP : public BlockDerivative
{
public:
    BlockSRP();
 
    BlockSRP(const BlockSRP&) = delete;
    BlockSRP& operator=(const BlockSRP&) = delete;
 
    BlockSRP(EclipseCalculator* eclipseCalculator, double cr, double srpArea, Frame* propagationFrame);
 
    ~BlockSRP() override;
 
    errc_t run(const SimTime& simTime) override;
public:
    void setSunPosition(ESunPosition sunPosition){sunPosition_ = sunPosition;}
    void setCoefficient(double cr){cr_ = cr;}
    double coefficient() const {return cr_;}
protected:
    Vector3d* position_{&vectorBuffer_};            
    Vector3d* velocity_{&vectorBuffer_};            
    Vector3d* accSRP_{&vectorBuffer_};              
    Vector3d* velocityDerivative_{&vectorBuffer_};  
    double* mass_{&doubleBuffer_};                  
    double doubleBuffer_{};                         
    Vector3d vectorBuffer_{};                       
 
protected:
    EclipseCalculator* eclipseCalculator_{};        
    double cr_{1.0};                                
    double srpArea_{20.0};                          
    Frame* propagationFrame_{};                     
    ESunPosition sunPosition_{ESunPosition::eTrue}; 
};
 
AST_NAMESPACE_END