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
75
76
77
78
|
/*
[auto_generated]
boost/numeric/odeint/util/detail/less_with_sign.hpp
[begin_description]
Helper function to compare times taking into account the sign of dt
[end_description]
Copyright 2012-2013 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 BOOST_NUMERIC_ODEINT_INTEGRATE_DETAIL_LESS_WITH_SIGN_HPP_INCLUDED
#define BOOST_NUMERIC_ODEINT_INTEGRATE_DETAIL_LESS_WITH_SIGN_HPP_INCLUDED
#include <limits>
#include <boost/numeric/odeint/util/unit_helper.hpp>
namespace boost {
namespace numeric {
namespace odeint {
namespace detail {
/**
* return t1 < t2 if dt > 0 and t1 > t2 if dt < 0 with epsilon accuracy
*/
template< typename T >
bool less_with_sign( T t1 , T t2 , T dt )
{
if( get_unit_value(dt) > 0 )
//return t1 < t2;
return t2-t1 > std::numeric_limits<T>::epsilon();
else
//return t1 > t2;
return t1-t2 > std::numeric_limits<T>::epsilon();
}
/**
* return t1 <= t2 if dt > 0 and t1 => t2 if dt < 0 with epsilon accuracy
*/
template< typename T >
bool less_eq_with_sign( T t1 , T t2 , T dt )
{
if( get_unit_value(dt) > 0 )
return t1-t2 < std::numeric_limits<T>::epsilon();
else
return t2-t1 < std::numeric_limits<T>::epsilon();
}
template< typename T >
T min_abs( T t1 , T t2 )
{
BOOST_USING_STD_MIN();
BOOST_USING_STD_MAX();
if( t1>0 )
return min BOOST_PREVENT_MACRO_SUBSTITUTION ( t1 , t2 );
else
return max BOOST_PREVENT_MACRO_SUBSTITUTION ( t1 , t2 );
}
template< typename T >
T max_abs( T t1 , T t2 )
{
BOOST_USING_STD_MIN();
BOOST_USING_STD_MAX();
if( t1>0 )
return max BOOST_PREVENT_MACRO_SUBSTITUTION ( t1 , t2 );
else
return min BOOST_PREVENT_MACRO_SUBSTITUTION ( t1 , t2 );
}
} } } }
#endif
|
