Subject: [Boost-commit] svn:boost r60475 - in sandbox/odeint/libs/numeric/odeint/stuff: . test
From: karsten.ahnert_at_[hidden]
Date: 2010-03-11 03:47:46
Author: karsten
Date: 2010-03-11 03:47:45 EST (Thu, 11 Mar 2010)
New Revision: 60475
URL: http://svn.boost.org/trac/boost/changeset/60475
Log:
more changes in the dir structure
Added:
sandbox/odeint/libs/numeric/odeint/stuff/test_stepper_bs.cpp
- copied unchanged from r60427, /sandbox/odeint/libs/numeric/odeint/stuff/test/stepper_bs.cpp
sandbox/odeint/libs/numeric/odeint/stuff/test_stepper_midpoint.cpp
- copied unchanged from r60427, /sandbox/odeint/libs/numeric/odeint/stuff/test/stepper_midpoint.cpp
Removed:
sandbox/odeint/libs/numeric/odeint/stuff/test/stepper_bs.cpp
sandbox/odeint/libs/numeric/odeint/stuff/test/stepper_midpoint.cpp
Deleted: sandbox/odeint/libs/numeric/odeint/stuff/test/stepper_bs.cpp
==============================================================================
--- sandbox/odeint/libs/numeric/odeint/stuff/test/stepper_bs.cpp 2010-03-11 03:47:45 EST (Thu, 11 Mar 2010)
+++ (empty file)
@@ -1,51 +0,0 @@
-#include <iostream>
-#include <fstream>
-#include <iomanip>
-#include <vector>
-#include <cmath>
-
-#include <boost/numeric/odeint.hpp>
-
-#define tab "\t"
-
-using namespace std;
-using namespace boost::numeric::odeint;
-
-
-typedef std::tr1::array< double , 2 > state_type;
-
-
-const double gam = .15;
-
-void harmonic_oscillator(const state_type &x, state_type &dxdt, const double t)
-{
- dxdt[0] = x[1];
- dxdt[1] = -x[0] - gam*x[1];
-}
-
-
-int main( int argc , char **argv )
-{
- state_type x1 = {{ 1.0, 0.0 }};
-
- const size_t step_count = 100;
- controlled_stepper_bs< state_type > stepper_bs(0.0, 1E-11, 1.0, 1.0);
- double t = 0.0;
- double dt = 0.01;
- cout << " Everything initialized - starting time evolution " << endl;
- cout.precision(16);
- clog.precision(16);
- for( size_t step=0; step < step_count; step++ )
- {
- controlled_step_result res = stepper_bs.try_step(harmonic_oscillator, x1, t, dt);
- //clog << " ####################################################### " << endl;
- cout << t << tab << dt << tab << x1[0] << tab << x1[1] << tab << res << endl;
- //clog << " ####################################################### " << endl;
- }
-}
-
-/*
- Compile with
- g++ -Wall -O3 -I$BOOST_ROOT -I../../../../../ stepper_bs.cpp
-*/
-
Deleted: sandbox/odeint/libs/numeric/odeint/stuff/test/stepper_midpoint.cpp
==============================================================================
--- sandbox/odeint/libs/numeric/odeint/stuff/test/stepper_midpoint.cpp 2010-03-11 03:47:45 EST (Thu, 11 Mar 2010)
+++ (empty file)
@@ -1,63 +0,0 @@
-#include <iostream>
-#include <fstream>
-#include <iomanip>
-#include <vector>
-#include <cmath>
-
-#include <boost/numeric/odeint.hpp>
-
-#define tab "\t"
-
-using namespace std;
-using namespace boost::numeric::odeint;
-
-
-typedef std::tr1::array< double , 2 > state_type;
-
-
-const double gam = .15;
-
-void harmonic_oscillator(const state_type &x, state_type &dxdt, const double t)
-{
- dxdt[0] = x[1];
- dxdt[1] = -x[0] - gam*x[1];
-}
-
-
-int main( int argc , char **argv )
-{
- const double dt = 0.01;
- const size_t olen = 10000;
- state_type x1 = {{ 1.0, 0.0 }};
- state_type x2 = {{ 1.0, 0.0 }};
- state_type x3 = {{ 1.0, 0.0 }};
- state_type x4 = {{ 1.0, 0.0 }};
- state_type x5 = {{ 1.0, 0.0 }};
-
- stepper_euler< state_type > stepper_euler;
- stepper_midpoint< state_type > stepper_mp;
- stepper_rk4< state_type > stepper_rk4;
-
- double t = 0.0;
- for( size_t oi=0 ; oi<olen ; ++oi,t+=dt )
- {
- stepper_euler.do_step( harmonic_oscillator , x1 , t , dt );
- stepper_mp.set_stepcount(2);
- stepper_mp.do_step( harmonic_oscillator , x2 , t , dt );
- stepper_mp.set_stepcount(4);
- stepper_mp.do_step( harmonic_oscillator , x3 , t , dt );
- stepper_mp.set_stepcount(8);
- stepper_mp.do_step( harmonic_oscillator , x4 , t , dt );
- stepper_rk4.do_step( harmonic_oscillator , x5 , t , dt );
- cout<< t << tab << x1[0]*x1[0] + x1[1]*x1[1];
- cout<< tab << x2[0]*x2[0] + x2[1]*x2[1];
- cout<< tab << x3[0]*x3[0] + x3[1]*x3[1];
- cout<< tab << x4[0]*x4[0] + x4[1]*x4[1];
- cout<< tab << x5[0]*x5[0] + x5[1]*x5[1] << endl;
- }
-}
-
-/*
- Compile with
- g++ -Wall -O3 -I$BOOST_ROOT -I../../../../../ stepper_midpoint.cpp
-*/