Imported from /home/lorry/working-area/delta_boost-tarball/boost_1_58_0.tar.bz2.HEADboost_1_58_0master

1 files changed, 0 insertions, 84 deletions

diff --git a/libs/numeric/odeint/performance/nr_rk4_lorenz.cpp b/libs/numeric/odeint/performance/nr_rk4_lorenz.cpp
deleted file mode 100644
index 36be91df8..000000000
--- a/libs/numeric/odeint/performance/nr_rk4_lorenz.cpp
+++ /dev/null
@@ -1,84 +0,0 @@
-/*
- * nr_rk4_lorenz.cpp
- *
- * Copyright 2011 Mario Mulansky
- * Copyright 2012 Karsten Ahnert
- *
- * 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 <boost/array.hpp>
-
-#include "rk_performance_test_case.hpp"
-
-#include "lorenz.hpp"
-
-const size_t dim = 3;
-
-typedef boost::array< double , dim > state_type;
-
-
-template< class System , typename T , size_t dim >
-void rk4_step( const System sys , boost::array< T , dim > &x , const double t , const double dt )
-{   // fast rk4 implementation adapted from the book 'Numerical Recipes'
-    size_t i;
-    const double hh = dt*0.5;
-    const double h6 = dt/6.0;
-    const double th = t+hh;
-    boost::array< T , dim > dydx , dym , dyt , yt;
-
-    sys( x , dydx , t );
-
-    for( i=0 ; i<dim ; i++ )
-        yt[i] = x[i] + hh*dydx[i];
-
-    sys( yt , dyt , th );
-    for( i=0 ; i<dim ; i++ )
-        yt[i] = x[i] + hh*dyt[i];
-
-    sys( yt , dym , th );
-    for( i=0 ; i<dim ; i++ ) {
-        yt[i] = x[i] + dt*dym[i];
-        dym[i] += dyt[i];
-    }
-    sys( yt , dyt , t+dt );
-    for( i=0 ; i<dim ; i++ )
-        x[i] += h6*( dydx[i] + dyt[i] + 2.0*dym[i] );
-}
-
-
-class nr_wrapper
-{
-public:
-    void reset_init_cond()
-    {
-        m_x[0] = 10.0 * rand() / RAND_MAX;
-        m_x[1] = 10.0 * rand() / RAND_MAX;
-        m_x[2] = 10.0 * rand() / RAND_MAX;
-        m_t = 0.0;
-    }
-
-    inline void do_step( const double dt )
-    {
-        rk4_step( lorenz() , m_x , m_t , dt );
-    }
-
-    double state( const size_t i ) const
-    { return m_x[i]; }
-
-private:
-    state_type m_x;
-    double m_t;
-};
-
-
-
-int main()
-{
-    nr_wrapper stepper;
-
-    run( stepper );
-}