summaryrefslogtreecommitdiff
path: root/libs/numeric/odeint/performance/phase_lattice.hpp
diff options
context:
space:
mode:
Diffstat (limited to 'libs/numeric/odeint/performance/phase_lattice.hpp')
-rw-r--r--libs/numeric/odeint/performance/phase_lattice.hpp43
1 files changed, 0 insertions, 43 deletions
diff --git a/libs/numeric/odeint/performance/phase_lattice.hpp b/libs/numeric/odeint/performance/phase_lattice.hpp
deleted file mode 100644
index 6657edeea..000000000
--- a/libs/numeric/odeint/performance/phase_lattice.hpp
+++ /dev/null
@@ -1,43 +0,0 @@
-/*
- * 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 <cmath>
-
-#include <boost/array.hpp>
-
-template< size_t N >
-struct phase_lattice
-{
- typedef double value_type;
- typedef boost::array< value_type , N > state_type;
-
- value_type m_epsilon;
- state_type m_omega;
-
- phase_lattice() : m_epsilon( 6.0/(N*N) ) // should be < 8/N^2 to see phase locking
- {
- for( size_t i=1 ; i<N-1 ; ++i )
- m_omega[i] = m_epsilon*(N-i);
- }
-
- void inline operator()( const state_type &x , state_type &dxdt , const double t ) const
- {
- double c = 0.0;
-
- for( size_t i=0 ; i<N-1 ; ++i )
- {
- dxdt[i] = m_omega[i] + c;
- c = ( x[i+1] - x[i] );
- dxdt[i] += c;
- }
-
- //dxdt[N-1] = m_omega[N-1] + sin( x[N-1] - x[N-2] );
- }
-
-};