diff options
author | Lorry Tar Creator <lorry-tar-importer@baserock.org> | 2015-04-08 03:09:47 +0000 |
---|---|---|
committer | <> | 2015-05-05 14:37:32 +0000 |
commit | f2541bb90af059680aa7036f315f052175999355 (patch) | |
tree | a5b214744b256f07e1dc2bd7273035a7808c659f /libs/numeric/odeint/performance/phase_lattice_mkl.hpp | |
parent | ed232fdd34968697a68783b3195b1da4226915b5 (diff) | |
download | boost-tarball-master.tar.gz |
Imported from /home/lorry/working-area/delta_boost-tarball/boost_1_58_0.tar.bz2.HEADboost_1_58_0master
Diffstat (limited to 'libs/numeric/odeint/performance/phase_lattice_mkl.hpp')
-rw-r--r-- | libs/numeric/odeint/performance/phase_lattice_mkl.hpp | 57 |
1 files changed, 0 insertions, 57 deletions
diff --git a/libs/numeric/odeint/performance/phase_lattice_mkl.hpp b/libs/numeric/odeint/performance/phase_lattice_mkl.hpp deleted file mode 100644 index 4feea47c1..000000000 --- a/libs/numeric/odeint/performance/phase_lattice_mkl.hpp +++ /dev/null @@ -1,57 +0,0 @@ -/* - * phase_lattice_mkl.hpp - * - * 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) - */ - - - - -#ifndef PHASE_LATTICE_MKL_HPP_ -#define PHASE_LATTICE_MKL_HPP_ - -#include <cmath> - -#include <mkl_blas.h> -#include <mkl_vml_functions.h> -#include <boost/array.hpp> - -template< size_t N > -struct phase_lattice_mkl -{ - typedef double value_type; - typedef boost::array< value_type , N > state_type; - - value_type m_epsilon; - state_type m_omega; - state_type m_tmp; - - phase_lattice_mkl() : 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 int n = x.size(); - - dxdt[0] = m_omega[0] + sin( x[1] - x[0] ); - - vdSub( n-1 , &(x[1]) , &(x[0]) , &(m_tmp[0]) ); - vdSin( n-1 , &(m_tmp[0]) , &(m_tmp[0]) ); - vdAdd( n-2 , &(m_tmp[0]) , &(m_tmp[1]) , &(dxdt[1]) ); - vdAdd( n-2 , &(dxdt[1]) , &(m_omega[1]) , &(dxdt[1]) ); - - dxdt[N-1] = m_omega[N-1] + sin( x[N-1] - x[N-2] ); - } - -}; - - -#endif /* PHASE_LATTICE_MKL_HPP_ */ |