c++ - 使用 OpenMP 在 odeint 中创建受控步进器

Question

我尝试使用 openmp_range_algebra 使用 boost::odeint 构建受控步进器

typedef vector< complex< double > > state_type;    
typedef runge_kutta_dopri5< state_type > error_stepper_type;
typedef controlled_runge_kutta< error_stepper_type > controlled_stepper_type;
controlled_stepper_type controlled_stepper(default_error_checker< double, openmp_range_algebra >;

但是，odeint 中不存在这样的构造函数，因此代码无法编译。

我的问题:如何创建一个受控的 runge_kutte-dopri5 步进器，以便我可以将它与 OpenMP 一起使用？

我真的很想并行化自适应步进器，因为这是我程序中最耗时的部分，因为状态 vector 很长(长度:2^20)。

非常感谢您的帮助

Answer 1

您需要使用范围代数对步进器进行参数化:

// Disclaimer: Not tested if it compiles
typedef runge_kutta_dopri5<
    state_type, double,
    state_type, double,
    openmp_range_algebra
> stepper_type;
typedef controlled_runge_kutta< stepper_type > controlled_stepper_type;
typedef controlled_stepper_type::error_checker_type error_checker_type;
const double eps_absolute = 1.0e-6;
const double eps_relative = 1.0e-6;
controlled_stepper_type stepper( error_checker_type( eps_absolute , eps_relative ) );

而且使用make_controlled工厂函数更简单。上面几行代码可以简化为

typedef runge_kutta_dopri5<
    state_type, double,
    state_type, double,
    openmp_range_algebra
> stepper_type;
auto stepper = make_controlled( 1.0e-6 , 1.0e-6 , stepper_type() );