Boost дает странный результат для интеграции оды - PullRequest
Новая
       7

Boost дает странный результат для интеграции оды

0 голосов
/ 24 апреля 2019

Я запустил следующий код, который решает простое дифференциальное уравнение. Результат, кажется, зависит от шагов. Результат, который я получаю для tsteps = 100, равен 9.688438503116524e-15, но для tsteps = 1000 ответ 7.124585369895499e-218 гораздо ближе к ожидаемому результату. 

#include <iostream>
#include <iomanip>
#include <boost/numeric/odeint.hpp> // odeint function definitions
using namespace std;
using namespace boost::numeric::odeint;
typedef std::vector< double> state_type;
int tsize = 1;
void my_observer( const state_type &x, const double t );
void initarrays(state_type &x)
{
    x[0]=1.0e0;
}
void my_system( const state_type &x , state_type &dxdt , const double t )
{
    dxdt[0]=-x[0];
}
void my_observer( const state_type &x, const double t )
{
    std::cout<<t<<" ";
    for(int i=0;i<tsize;i++)
    {
    std::cout<<x[i]<<" ";
    }
    std::cout<<std::endl;
}
int main()
{
    std::cout.setf ( std::ios::scientific, std::ios::floatfield );
    std::cout.precision(15);
    int size=tsize;
    state_type x0(size);
    double err_abs = 1.0e-12;
    double err_rel = 1.0e-12;
    double a_x = 1.0;
    double a_dxdt = 1.0;
    initarrays(x0);
    double t0 = 0.0e0;
    int tsteps = 1000;
    double t1 = 500.0e0;
    double dt = (t1-t0)/((double)tsteps);
    typedef runge_kutta_fehlberg78< state_type > solver;
    typedef controlled_runge_kutta< solver > controller;
    my_observer(x0,t0);
    for(int ts=0;ts<tsteps;ts++)
    {
        integrate_adaptive( make_controlled( err_abs , err_rel , solver() ), my_system, x0 , t0+ts*dt , t0+(1+ts)*dt , dt);
        my_observer(x0,t0+(1+ts)*dt);
    }
}
...