Временной шаг управления с использованием odeint (Boost)

Question

Я пытаюсь изучить решение проблемы ODE с помощью библиотеки ODEINT в следующем примере. Однако, когда я выводя результат, шаг времени просто прыгает через 0; 1; 5,5; 25 ... Есть ли способ контролировать этот временной шаг, чтобы увеличить его на «1». Благодаря!

#include <iostream> 
#include <boost/array.hpp> 

#include <boost/numeric/odeint.hpp> 

using namespace std; 
using namespace boost::numeric::odeint; 

const double sigma = 0.0018; 


typedef std::vector<double> state_type; 

void my_ode(const state_type &x , state_type &dxdt , int t) 
{ 
    dxdt[0] = sigma * x[0]*(1 - x[0]); 

} 

void write_ode(const state_type &x , const double t) 
{ 
    cout << t << '\t' << x[0] << endl; 
} 

int main(int argc, char **argv) 
{ 
    state_type x(1); // Initial condition, vector of 1 element (scalar problem) 
    x[0] = 0.001; 
    integrate(my_ode , x , 0.0 , 6000.0 , 1.0 , write_ode); 
} 

Answer 1

Да, есть простой способ сделать это. В учебнике полно примеров.

Вам необходимо определить конкретный решатель. Хорошим выбором может быть плотный выходной решатель, который имеет контроль размера шага и позволяет вычислить решение на произвольных временных шагах. Он может использоваться как

// not tested 
state_type x(1); 
x[0] = 0.001; 
auto rkd = runge_kutta_dopri5<state_type>{}; 
auto stepper = make_dense_output(1.0e-6, 1.0e-6, rkd); 
integrate_const(stepper, x, 0.0, 6000.0, 1.0, write_ode);