c++ - 将运行时参数传递给 odeint 集成器

Question

我想使用 odeint boost 积分器查看克尔时空中的测地线曲线。这需要为各种参数值运行积分器(我有初始条件和初始动量 vector ，因此系统的角动量将根据我想如何启动它而变化)

我一直在关注此处列出的优秀示例

http://headmyshoulder.github.io/odeint-v2/examples.html

特别是洛伦兹吸引子示例。

我注意到的第一件事(在我的 Kerr 系统中，而不是 Lorenz)是对于某些初始条件，我在 t=0 时间点之后得到 NaN(尽管毫无疑问这表明某处存在更深层次的错误)。为了进一步探索这一点，我希望能够更改系统的参数，而不必每次都重新编译。这相当于想要更改，比如 R 或 sigma 而无需重新编译。据我所知，没有直接的方法可以将额外的参数传递给集成函数，除了这里的第二个代码块

http://headmyshoulder.github.io/odeint-v2/doc/boost_numeric_odeint/tutorial/harmonic_oscillator.html

我承认，由于我缺乏C++知识，我不是很明白。

我很想知道如何获取运行时参数并将它们传递给这个求解器，这样我就可以运行一堆批处理而无需每次都进行编译。

Answer 1

您可以将运行时参数传递给定义 ODE 的系统函数:

struct ode
{
    double param;
    ode( double param ) : m_param( param ) {}

    void operator()( state_type const& x , state_type& dxdt , time_type t ) const 
    {
        // your ode
    }
};

integrate_const( stepper {} , ode { 0.1 } , x , t_start , t_end , dt , observer );