0

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

我有以下 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);

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

4

1 回答 1

1

使用

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

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

于 2013-01-23T16:32:08.980 回答