ODEVarStepIntegrator.hpp

浏览该文件的文档.

 
 
#pragma once
 
#include "AstGlobal.h"
#include "ODEFixedStepIntegrator.hpp"
 
 
AST_NAMESPACE_BEGIN
 
class AST_MATH_API ODEVarStepIntegrator : public ODEFixedStepIntegrator
{
public:
    AST_OBJECT(ODEVarStepIntegrator)
    AST_PROPERT(useMinStep)
    AST_PROPERT(useMaxStep)
    AST_PROPERT(warnOnMinStep)
    AST_PROPERT(maxStepAttempts)
    AST_PROPERT(minStepSize)
    AST_PROPERT(maxStepSize)
    AST_PROPERT(maxAbsErr)
    AST_PROPERT(maxRelErr)
    AST_PROPERT(minStepScaleFactor)
    AST_PROPERT(maxStepScaleFactor)
    AST_PROPERT(safetyCoeffLow)
    AST_PROPERT(safetyCoeffHigh)
    AST_PROPERT(errCtrPowthLow)
    AST_PROPERT(errCtrPowthHigh)
;
    
    ODEVarStepIntegrator();
    ~ODEVarStepIntegrator();
    using ODEFixedStepIntegrator::integrate;
 
    errc_t integrate(ODE& ode, double* y, double& t, double tf) final;
 
    errc_t integrateStep(ODE& ode, double* y, double& t, double tf) final;
 
    void setMaxAbsErr(double maxAbsErr){maxAbsErr_ = maxAbsErr;};
 
    void setMaxRelErr(double maxRelErr){maxRelErr_ = maxRelErr;};
 
    void setMaxIterations(int maxIterations){maxStepAttempts_ = maxIterations;};
    int getMaxIterations() const{return maxStepAttempts_;};
 
    void setMaxStepSize(double maxStepSize){maxStepSize_ = maxStepSize;};
    double getMaxStepSize() const{return maxStepSize_;};
 
    void setMinStepSize(double minStepSize){minStepSize_ = minStepSize;};
    double getMinStepSize() const{return minStepSize_;};
 
    void setSafetyCoeffLow(double safetyCoeffLow){safetyCoeffLow_ = safetyCoeffLow;};
    double getSafetyCoeffLow() const{return safetyCoeffLow_;};
 
    void setSafetyCoeffHigh(double safetyCoeffHigh){safetyCoeffHigh_ = safetyCoeffHigh;};
    double getSafetyCoeffHigh() const{return safetyCoeffHigh_;};
 
 
    void setUseMinStep(bool useMinStep){useMinStep_ = useMinStep;};
    bool getUseMinStep() const{return useMinStep_;};
 
    void setUseMaxStep(bool useMaxStep){useMaxStep_ = useMaxStep;};
    bool getUseMaxStep() const{return useMaxStep_;};
 
    // void setUseFixedStep(bool useFixedStep){useFixedStep_ = useFixedStep;};
    // bool getUseFixedStep() const{return useFixedStep_;};
 
    void setInitialStepSize(double initialStepSize){setStepSize(initialStepSize);};
 
    double getLargestStepSize() const;
 
    double getSmallestStepSize() const;
 
    using ODEFixedStepIntegrator::getNumSteps;
protected:
    bool isErrorMeet(double& absh, const double* y, const double* ynew);
PROPERTIES:
    bool useMinStep_;           
    bool useMaxStep_;           
    // bool useFixedStepSize_;     ///< 是否使用固定步长
    bool warnOnMinStep_;        
    int  maxStepAttempts_;      
    double minStepSize_;        
    double maxStepSize_;        
    double maxAbsErr_;          
    double maxRelErr_;          
    double minStepScaleFactor_; 
    double maxStepScaleFactor_; 
    double safetyCoeffLow_;     
    double safetyCoeffHigh_;    
    double errCtrPowthLow_;     
    double errCtrPowthHigh_;    
};
 
AST_NAMESPACE_END