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
66
67
68
69
70
71
72
73
74
|
/*
[auto_generated]
libs/numeric/odeint/test/prepare_stepper_testing.hpp
[begin_description]
This file defines some helper functions for the stepper tests.
[end_description]
Copyright 2011-2012 Karsten Ahnert
Copyright 2012-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 PREPARE_STEPPER_TESTING_HPP_
#define PREPARE_STEPPER_TESTING_HPP_
#include <boost/array.hpp>
#include <vector>
#include <complex>
#include <boost/fusion/sequence.hpp>
#include <boost/mpl/vector.hpp>
#include <boost/multiprecision/cpp_dec_float.hpp>
#include <boost/numeric/odeint/algebra/vector_space_algebra.hpp>
#include <boost/numeric/odeint/algebra/array_algebra.hpp>
#include <boost/numeric/odeint/algebra/algebra_dispatcher.hpp>
namespace mpl = boost::mpl;
namespace fusion = boost::fusion;
using namespace boost::numeric::odeint;
/* the state types that will be tested */
typedef std::vector< double > vector_type;
typedef std::vector< std::complex<double> > complex_vector_type;
typedef double vector_space_type;
typedef boost::array< double , 1 > array_type;
typedef boost::array< std::complex<double> , 1 > complex_array_type;
typedef boost::multiprecision::cpp_dec_float_50 mp_type;
typedef boost::array< mp_type , 1 > mp_array_type;
typedef mpl::vector< vector_type , complex_vector_type , vector_space_type ,
array_type , complex_array_type , mp_type , mp_array_type
>::type container_types;
namespace boost { namespace numeric { namespace odeint {
// mp_type forms a vector space
template<>
struct algebra_dispatcher< mp_type >
{
typedef vector_space_algebra algebra_type;
};
// add norm computation
template<>
struct vector_space_norm_inf< mp_type >
{
typedef mp_type result_type;
mp_type operator()( mp_type x ) const
{
using std::abs;
return abs(x);
}
};
} } }
#endif /* PREPARE_STEPPER_TESTING_HPP_ */
|
