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

Question

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

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

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 );

Answer 2

我遇到了类似的问题。我所做的是定义全局变量，以便我的函数可以访问我通过 argvs 传递的变量。如果您需要示例，请告诉我。