1234567891011121314151617181920212223242526272829303132333435363738394041424344454647484950515253545556575859606162636465 |
- /*
- [auto_generated]
- boost/numeric/odeint/external/compute/compute_algebra.hpp
- [begin_description]
- An algebra for Boost.Compute vectors.
- [end_description]
- Copyright 2009-2011 Karsten Ahnert
- Copyright 2009-2011 Mario Mulansky
- 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 BOOST_NUMERIC_ODEINT_EXTERNAL_COMPUTE_COMPUTE_ALGEBRA_HPP_DEFINED
- #define BOOST_NUMERIC_ODEINT_EXTERNAL_COMPUTE_COMPUTE_ALGEBRA_HPP_DEFINED
- #include <boost/preprocessor/repetition.hpp>
- #include <boost/compute.hpp>
- namespace boost {
- namespace numeric {
- namespace odeint {
- struct compute_algebra
- {
- #define BOOST_ODEINT_COMPUTE_STATE_PARAM(z, n, unused) \
- StateType ## n &s ## n,
- #define BOOST_ODEINT_COMPUTE_ALGEBRA(z, n, unused) \
- template< BOOST_PP_ENUM_PARAMS(n, class StateType), class Operation > \
- static void for_each ## n( \
- BOOST_PP_REPEAT(n, BOOST_ODEINT_COMPUTE_STATE_PARAM, ~) \
- Operation op \
- ) \
- { \
- op( BOOST_PP_ENUM_PARAMS(n, s) ); \
- }
- BOOST_PP_REPEAT_FROM_TO(3, 9, BOOST_ODEINT_COMPUTE_ALGEBRA, ~)
- #undef BOOST_ODEINT_COMPUTE_ALGEBRA
- #undef BOOST_ODEINT_COMPUTE_STATE_PARAM
- template < class S >
- static typename S::value_type norm_inf( const S &s ) {
- typedef typename S::value_type value_type;
- BOOST_COMPUTE_FUNCTION(value_type, max_abs, (value_type, value_type),
- {
- return max(_1, fabs(_2));
- });
- return boost::compute::accumulate(s.begin(), s.end(), value_type(), max_abs);
- }
- };
- } // odeint
- } // numeric
- } // boost
- #endif
|