存储来自odeint刚性方程求解器的输出

Question

我想知道是否有办法将以下code的输出存储到增强矩阵或向量中。

size_t num_of_steps = integrate_const( make_dense_output< rosenbrock4< double > >( 1.0e-6 , 1.0e-6 ) ,
        make_pair( stiff_system() , stiff_system_jacobi() ) ,
        x , 0.0 , 50.0 , 0.01 ,cout << phoenix::arg_names::arg2 << " " << phoenix::arg_names::arg1[0] << "\n" );

上面的第94-97行传递cout＆lt;＆lt; “”＆lt;＆lt;将endl语句作为参数放入函数中。这些语句打印结果。我没有找到通过修改这些传递的语句来存储结果的方法。

Answer 1

向量x在不同的时间存储状态值，所以简单的赋值语句，如

 a[0] = x[0];

其中a [0]是您要存储输出的矢量，将起作用。