odeint stiff equation solverの出力を保存

Question

の出力をブーストマトリクスまたはベクトルに格納する方法があるかどうかは疑問でした。上記

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 < < ENDL文を通過します。これらのステートメントは結果を出力します。私はこれらの渡されたステートメントを変更して結果を格納する方法を見つけていません。

Answer 1

ベクトルxを格納状態の値が異なる時間に[0]

a[0] = x[0]; 

のようなので、単純な代入文を、あなたがあなたの出力を保存するベクトルで作業する場所。