1234567891011121314151617181920212223242526272829303132333435363738394041424344454647484950515253545556575859606162636465666768 |
- /*
- [auto_generated]
- boost/numeric/odeint/external/vexcl/vexcl_norm_inf.hpp
- [begin_description]
- vector_space_norm_inf specialization for vexcl
- [end_description]
- Copyright 2009-2013 Karsten Ahnert
- Copyright 2009-2013 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_VEXCL_VEXCL_NORM_INF_HPP_DEFINED
- #define BOOST_NUMERIC_ODEINT_EXTERNAL_VEXCL_VEXCL_NORM_INF_HPP_DEFINED
- #include <map>
- #include <algorithm>
- #include <vexcl/vector.hpp>
- #include <vexcl/multivector.hpp>
- #include <vexcl/reductor.hpp>
- #include <boost/numeric/odeint/algebra/vector_space_algebra.hpp>
- namespace boost {
- namespace numeric {
- namespace odeint {
- // specialization for vexcl vector
- template <typename T>
- struct vector_space_norm_inf< vex::vector<T> > {
- typedef T result_type;
- T operator()( const vex::vector<T> &x ) const {
- const auto &max = vex::get_reductor<T, vex::MAX>(x.queue_list());
- return max( fabs(x) );
- }
- };
- // specialization for vexcl multivector
- template <typename T, size_t N>
- struct vector_space_norm_inf< vex::multivector<T, N> > {
- typedef T result_type;
- T operator()( const vex::multivector<T, N> &x ) const {
- const auto &max = vex::get_reductor<T, vex::MAX>(x.queue_list());
- // Reducing a multivector results in std::array<T, N>:
- auto m = max( fabs(x) );
- // We will need to reduce it even further:
- return *std::max_element(m.begin(), m.end());
- }
- };
- } // namespace odeint
- } // namespace numeric
- } // namespace boost
- #endif // BOOST_NUMERIC_ODEINT_EXTERNAL_VEXCL_VEXCL_NORM_INF_HPP_DEFINED
|