summaryrefslogtreecommitdiff
path: root/libs/numeric/odeint/performance/gsl_rk4_lorenz.cpp
diff options
context:
space:
mode:
Diffstat (limited to 'libs/numeric/odeint/performance/gsl_rk4_lorenz.cpp')
-rw-r--r--libs/numeric/odeint/performance/gsl_rk4_lorenz.cpp71
1 files changed, 0 insertions, 71 deletions
diff --git a/libs/numeric/odeint/performance/gsl_rk4_lorenz.cpp b/libs/numeric/odeint/performance/gsl_rk4_lorenz.cpp
deleted file mode 100644
index f5f13882a..000000000
--- a/libs/numeric/odeint/performance/gsl_rk4_lorenz.cpp
+++ /dev/null
@@ -1,71 +0,0 @@
-/*
- * gsl_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 <gsl/gsl_odeiv.h>
-
-#include "rk_performance_test_case.hpp"
-
-#include "lorenz_gsl.hpp"
-
-const size_t dim = 3;
-
-class gsl_wrapper
-{
-public:
-
- gsl_wrapper()
- {
- m_s = gsl_odeiv_step_alloc( gsl_odeiv_step_rk4 , dim);
- m_sys.function = lorenz_gsl;
- m_sys.jacobian = 0;
- m_sys.dimension = dim;
- m_sys.params = 0;
- }
-
- 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 )
- {
- gsl_odeiv_step_apply ( m_s , m_t , dt , m_x , m_x_err , 0 , 0 , &m_sys );
- //m_t += dt;
- }
-
- double state( const size_t i ) const
- { return m_x[i]; }
-
- ~gsl_wrapper()
- {
- gsl_odeiv_step_free( m_s );
- }
-
-private:
- double m_x[dim];
- double m_x_err[dim];
- double m_t;
- gsl_odeiv_step *m_s;
- gsl_odeiv_system m_sys;
-};
-
-
-
-int main()
-{
- gsl_wrapper stepper;
-
- run( stepper , 20000000 / 3 , 1E-10 * 3);
-}