1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
|
/*
[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
|
