Question

我想在boost :: odeint中设置步进器的current_state。

我有以下MWE：

#include <queue>
#include <boost/numeric/odeint.hpp>
#include <array>

typedef std::array< double , 2 > state_type;

void eqsys(const state_type &x, state_type &dxdt, double t) {
    dxdt[0] = 0.3*x[1];
    dxdt[1] = 4*x[0];
}

int main() {
    int steps=0;
    auto stepper=make_dense_output(1.0e-6,1.0e-6,boost::numeric::odeint::runge_kutta_dopri5< state_type >() );
    state_type x={2,2};
    double stoptime=10;

    stepper.initialize(x,0,0.01);

    while(true){
        //Run stepper
        while( stepper.current_time()<stoptime && steps<10000 ) {
            stepper.do_step(eqsys);
            ++steps;
        }
        x=stepper.current_state();
        x[0]=2;
        stepper.set_current_state(x);
        stoptime+=10;
    }
}

行stepper.set_current_state(x);是不起作用的行。有没有我可以用来完成这个的命令？

我意识到我可以使用：

stepper.initialize(x, stepper.current_time(), 0.01);

但目前还不清楚这是否真的是最好的策略，或者我是否失去了步进器的某些内部状态（例如当前的步长或错误估计），最好保留它。

Answer 1

使用

stepper.initialize(x, stepper.current_time(), 0.01);

绝对可以。在这种情况下，内部状态是正确设置的。

在odeint中设置当前状态

1 个答案: