EQ2EMu/EQ2/source/depends/boost_1_72_0/libs/numeric/odeint/examples/integrate_times.cpp

/* Boost libs/numeric/odeint/examples/integrate_times.cpp

 Copyright 2009-2014 Karsten Ahnert
 Copyright 2009-2014 Mario Mulansky

 example for the use of integrate_times

 Distributed under the Boost Software License, Version 1.0.
(See accompanying file LICENSE_1_0.txt or
 copy at http://www.boost.org/LICENSE_1_0.txt)
 */


#include <iostream>
#include <boost/numeric/odeint.hpp>

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


/*
 * simple 1D ODE
 */

void rhs( const double x , double &dxdt , const double t )
{
    dxdt = 3.0/(2.0*t*t) + x/(2.0*t);
}

void write_cout( const double &x , const double t )
{
    cout << t << '\t' << x << endl;
}

// state_type = double
typedef runge_kutta_dopri5< double > stepper_type;

const double dt = 0.1;

int main()
{
    // create a vector with observation time points
    std::vector<double> times( 91 );
    for( size_t i=0 ; i<times.size() ; ++i )
        times[i] = 1.0 + dt*i;

    double x = 0.0; //initial value x(1) = 0
    // we can provide the observation time as a boost range (i.e. the vector)
    integrate_times( make_controlled( 1E-12 , 1E-12 , stepper_type() ) , rhs ,
                     x , times , dt , write_cout );
    // or as two iterators
    //integrate_times( make_controlled( 1E-12 , 1E-12 , stepper_type() ) , rhs ,
    //                   x , times.begin() , times.end() , dt , write_cout );
}