将运行时参数传递给odeint积分器

Passing runtime parameters to odeint integrator

本文关键字：odeint 运行时 参数传递      更新时间：2023-10-16

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

我一直在关注这里列出的优秀例子

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

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

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

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

我承认，由于我缺乏C++知识，我不太理解。

我很感激了解如何接收运行时参数并将其传递给这个求解器，这样我就可以运行一堆批处理，而无需每次编译。

您可以将运行时参数传递给定义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 );

我遇到了类似的问题。我所做的是定义全局变量，这样我的函数就可以访问我通过argv传递的变量。如果你需要例子，请告诉我。

相关文章：