[Boost-commit] svn:boost r60475 - in sandbox/odeint/libs/numeric/odeint/stuff: . test

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
-*/