使用BOOST ODEint时integrate_const()参数类型匹配问题求助

问题描述 #

本人是C++新手，尝试使用BOOST的ODEint库求解常微分方程，但在给一维ODE应用步进器时遇到诸多问题。已确认依赖版本正确，怀疑是denseStepper的类型与integrate_const要求不匹配，恳请帮助指明方向。

原代码 #

#define ODEINT boost::numeric::odeint

// Global variables to store time and state values
std::vector<double> times;
std::vector<double> states;

// Define initial condition and time range
double x0 = 1.0; // initial condition
double t_start = 0.0;
double t_end = 5.0;
double dt = 0.1; // time step
int numSteps = (t_start - t_end) / dt;

void odeSystem(const double x, double& dxdt, const double t) {
    // Define simple ODE
    dxdt = -x;
}

float odeSolver(double X, double T) {

    auto denseStepper = ODEINT::make_dense_output<ODEINT::runge_kutta_cash_karp54<double>>(1.0e-6, 1.0e-6);

    ODEINT::integrate_const(
        denseStepper,
        odeSystem,
        x0,
        t_start,
        t_end,
        dt,
        std::function<void(const double&, const double)>([](const double& x, const double t) {
            times.push_back(t);
            states.push_back(x);
            })
    );

    int currentPt = (X - t_start) / dt;

    return states[currentPt];
}

编译错误 #

1>D:\BOOST\boost_1_84_0\boost\numeric\odeint\stepper\generation\make_dense_output.hpp(60,55): 
error C2039: 'type': is not a member of 'boost::numeric::odeint::get_dense_output<Stepper>'    
with 
[           
Stepper=boost::numeric::odeint::runge_kutta_cash_karp54<double,double,double,double,boost::numeric::odeint::algebra_dispatcher_sfinae<double,void>::algebra_type,boost::numeric::odeint::operations_dispatcher_sfinae<double,void>::operations_type,boost::numeric::odeint::initially_resizer>    
]

2>D:\BOOST\boost_1_84_0\boost\numeric\odeint\stepper\generation\make_dense_output.hpp(60,26):
see declaration of 'boost::numeric::odeint::get_dense_output<Stepper>'
with
[
Stepper=boost::numeric::odeint::runge_kutta_cash_karp54<double,double,double,double,boost::numeric::odeint::algebra_dispatcher_sfinae<double,void>::algebra_type,boost::numeric::odeint::operations_dispatcher_sfinae<double,void>::operations_type,boost::numeric::odeint::initially_resizer>
]

3>D:\BOOST\boost_1_84_0\boost\numeric\odeint\stepper\generation\make_dense_output.hpp(60,55):
the template instantiation context (the oldest one first) is
D:\Visual Studio\source\repos\TestFiles\main.cpp(48,106):
see reference to class template instantiation 'boost::numeric::odeint::result_of::make_dense_output<boost::numeric::odeint::runge_kutta_cash_karp54<double,double,State,Value,boost::numeric::odeint::algebra_dispatcher_sfinae<StateType,void>::algebra_type,boost::numeric::odeint::operations_dispatcher_sfinae<StateType,void>::operations_type,boost::numeric::odeint::initially_resizer>>' being compiled
with
[
State=double,
Value=double,
StateType=double
]

4>D:\Visual Studio\source\repos\TestFiles\main.cpp(48,23): 
error C2514: 'boost::type': class template cannot be constructed

5>D:\BOOST\boost_1_84_0\boost\type.hpp(14,3):
see declaration of 'boost::type'

6>D:\Visual Studio\source\repos\TestFiles\main.cpp(50,13): error C2664: 
'size_t boost::numeric::odeint::integrate_const<boost::type,void(__cdecl *)(double,double &,double),double,double,std::function<void (const double &,double)>>(Stepper,System,State &,Time,Time,Time,Observer)': cannot convert argument 1 from 'boost::type' to 'Stepper'
with
[
Stepper=boost::type,
System=void (__cdecl *)(double,double &,double),
State=double,
Time=double,
Observer=std::function<void (const double &,double)>
 ]
and
[
Stepper=boost::type
]

7>D:\Visual Studio\source\repos\TestFiles\main.cpp(51,9):
The target type has no constructors

8>D:\BOOST\boost_1_84_0\boost\numeric\odeint\integrate\integrate_const.hpp(104,8):
see declaration of 'boost::numeric::odeint::integrate_const'

9>D:\Visual Studio\source\repos\TestFiles\main.cpp(50,13):
while trying to match the argument list '(boost::type, overloaded-function, double, double, double, double, std::function<void (const double &,double)>)'

错误原因与解决方案 #

• 核心错误：runge_kutta_cash_karp54<double>本身不是支持稠密输出的步进器，make_dense_output需要的是基础步进器的模板，而非实例化后的类型。正确的做法是使用专门的稠密输出步进器包装，或者直接用预定义的稠密输出步进器。

• 修正步骤：

  1. 将make_dense_output的模板参数调整为正确的基础步进器类型，或直接使用runge_kutta_cash_karp54_dense_output<double>预定义步进器。
  2. 确保integrate_const传入的是合法的步进器实例，类型匹配。
  3. 全局变量times和states每次调用odeSolver前需清空，避免累积旧数据。
  4. 修正numSteps的计算逻辑，(t_end - t_start)/dt才是正确的步数计算方式。

修正后的代码 #

#define ODEINT boost::numeric::odeint

#include <vector>
#include <functional>
#include <boost/numeric/odeint.hpp>

std::vector<double> times;
std::vector<double> states;

// Define initial condition and time range
double x0 = 1.0; // initial condition
double t_start = 0.0;
double t_end = 5.0;
double dt = 0.1; // time step
int numSteps = (t_end - t_start) / dt; // 修正步数计算

void odeSystem(const double x, double& dxdt, const double t) {
    dxdt = -x;
}

float odeSolver(double X, double T) {
    // 清空历史数据
    times.clear();
    states.clear();

    // 方案1：使用make_dense_output创建稠密输出步进器
    auto denseStepper = ODEINT::make_dense_output<ODEINT::runge_kutta_cash_karp54<double>>(1.0e-6, 1.0e-6);

    // 方案2：直接使用预定义的稠密输出步进器
    // ODEINT::runge_kutta_cash_karp54_dense_output<double> denseStepper(1.0e-6, 1.0e-6);

    ODEINT::integrate_const(
        denseStepper,
        odeSystem,
        x0,
        t_start,
        t_end,
        dt,
        [](const double& x, const double t) {
            times.push_back(t);
            states.push_back(x);
        }
    );

    // 索引越界检查
    int currentPt = static_cast<int>((X - t_start) / dt);
    if (currentPt < 0 || currentPt >= static_cast<int>(states.size())) {
        return NAN; // 处理越界情况
    }

    return static_cast<float>(states[currentPt]);
}

额外说明 #

• 若仅需要固定时间步的解，无需使用稠密输出步进器，直接用普通runge_kutta_cash_karp54<double>即可。
• 尽量避免全局变量存储结果，可改为将容器作为参数传入函数，提升代码复用性与线程安全性。

内容的提问来源于stack exchange，提问作者Fredl0309