Chris@16
|
1 /* Boost interval/checking.hpp template implementation file
|
Chris@16
|
2 *
|
Chris@16
|
3 * Copyright 2002 Hervé Brönnimann, Guillaume Melquiond, Sylvain Pion
|
Chris@16
|
4 *
|
Chris@16
|
5 * Distributed under the Boost Software License, Version 1.0.
|
Chris@16
|
6 * (See accompanying file LICENSE_1_0.txt or
|
Chris@16
|
7 * copy at http://www.boost.org/LICENSE_1_0.txt)
|
Chris@16
|
8 */
|
Chris@16
|
9
|
Chris@16
|
10 #ifndef BOOST_NUMERIC_INTERVAL_CHECKING_HPP
|
Chris@16
|
11 #define BOOST_NUMERIC_INTERVAL_CHECKING_HPP
|
Chris@16
|
12
|
Chris@16
|
13 #include <stdexcept>
|
Chris@16
|
14 #include <string>
|
Chris@16
|
15 #include <cassert>
|
Chris@16
|
16 #include <boost/limits.hpp>
|
Chris@16
|
17
|
Chris@16
|
18 namespace boost {
|
Chris@16
|
19 namespace numeric {
|
Chris@16
|
20 namespace interval_lib {
|
Chris@16
|
21
|
Chris@16
|
22 struct exception_create_empty
|
Chris@16
|
23 {
|
Chris@16
|
24 void operator()()
|
Chris@16
|
25 {
|
Chris@16
|
26 throw std::runtime_error("boost::interval: empty interval created");
|
Chris@16
|
27 }
|
Chris@16
|
28 };
|
Chris@16
|
29
|
Chris@16
|
30 struct exception_invalid_number
|
Chris@16
|
31 {
|
Chris@16
|
32 void operator()()
|
Chris@16
|
33 {
|
Chris@16
|
34 throw std::invalid_argument("boost::interval: invalid number");
|
Chris@16
|
35 }
|
Chris@16
|
36 };
|
Chris@16
|
37
|
Chris@16
|
38 template<class T>
|
Chris@16
|
39 struct checking_base
|
Chris@16
|
40 {
|
Chris@16
|
41 static T pos_inf()
|
Chris@16
|
42 {
|
Chris@16
|
43 assert(std::numeric_limits<T>::has_infinity);
|
Chris@16
|
44 return std::numeric_limits<T>::infinity();
|
Chris@16
|
45 }
|
Chris@16
|
46 static T neg_inf()
|
Chris@16
|
47 {
|
Chris@16
|
48 assert(std::numeric_limits<T>::has_infinity);
|
Chris@16
|
49 return -std::numeric_limits<T>::infinity();
|
Chris@16
|
50 }
|
Chris@16
|
51 static T nan()
|
Chris@16
|
52 {
|
Chris@16
|
53 assert(std::numeric_limits<T>::has_quiet_NaN);
|
Chris@16
|
54 return std::numeric_limits<T>::quiet_NaN();
|
Chris@16
|
55 }
|
Chris@16
|
56 static bool is_nan(const T& x)
|
Chris@16
|
57 {
|
Chris@16
|
58 return std::numeric_limits<T>::has_quiet_NaN && (x != x);
|
Chris@16
|
59 }
|
Chris@16
|
60 static T empty_lower()
|
Chris@16
|
61 {
|
Chris@16
|
62 return (std::numeric_limits<T>::has_quiet_NaN ?
|
Chris@16
|
63 std::numeric_limits<T>::quiet_NaN() : static_cast<T>(1));
|
Chris@16
|
64 }
|
Chris@16
|
65 static T empty_upper()
|
Chris@16
|
66 {
|
Chris@16
|
67 return (std::numeric_limits<T>::has_quiet_NaN ?
|
Chris@16
|
68 std::numeric_limits<T>::quiet_NaN() : static_cast<T>(0));
|
Chris@16
|
69 }
|
Chris@16
|
70 static bool is_empty(const T& l, const T& u)
|
Chris@16
|
71 {
|
Chris@16
|
72 return !(l <= u); // safety for partial orders
|
Chris@16
|
73 }
|
Chris@16
|
74 };
|
Chris@16
|
75
|
Chris@16
|
76 template<class T, class Checking = checking_base<T>,
|
Chris@16
|
77 class Exception = exception_create_empty>
|
Chris@16
|
78 struct checking_no_empty: Checking
|
Chris@16
|
79 {
|
Chris@16
|
80 static T nan()
|
Chris@16
|
81 {
|
Chris@16
|
82 assert(false);
|
Chris@16
|
83 return Checking::nan();
|
Chris@16
|
84 }
|
Chris@16
|
85 static T empty_lower()
|
Chris@16
|
86 {
|
Chris@16
|
87 Exception()();
|
Chris@16
|
88 return Checking::empty_lower();
|
Chris@16
|
89 }
|
Chris@16
|
90 static T empty_upper()
|
Chris@16
|
91 {
|
Chris@16
|
92 Exception()();
|
Chris@16
|
93 return Checking::empty_upper();
|
Chris@16
|
94 }
|
Chris@16
|
95 static bool is_empty(const T&, const T&)
|
Chris@16
|
96 {
|
Chris@16
|
97 return false;
|
Chris@16
|
98 }
|
Chris@16
|
99 };
|
Chris@16
|
100
|
Chris@16
|
101 template<class T, class Checking = checking_base<T> >
|
Chris@16
|
102 struct checking_no_nan: Checking
|
Chris@16
|
103 {
|
Chris@16
|
104 static bool is_nan(const T&)
|
Chris@16
|
105 {
|
Chris@16
|
106 return false;
|
Chris@16
|
107 }
|
Chris@16
|
108 };
|
Chris@16
|
109
|
Chris@16
|
110 template<class T, class Checking = checking_base<T>,
|
Chris@16
|
111 class Exception = exception_invalid_number>
|
Chris@16
|
112 struct checking_catch_nan: Checking
|
Chris@16
|
113 {
|
Chris@16
|
114 static bool is_nan(const T& x)
|
Chris@16
|
115 {
|
Chris@16
|
116 if (Checking::is_nan(x)) Exception()();
|
Chris@16
|
117 return false;
|
Chris@16
|
118 }
|
Chris@16
|
119 };
|
Chris@16
|
120
|
Chris@16
|
121 template<class T>
|
Chris@16
|
122 struct checking_strict:
|
Chris@16
|
123 checking_no_nan<T, checking_no_empty<T> >
|
Chris@16
|
124 {};
|
Chris@16
|
125
|
Chris@16
|
126 } // namespace interval_lib
|
Chris@16
|
127 } // namespace numeric
|
Chris@16
|
128 } // namespace boost
|
Chris@16
|
129
|
Chris@16
|
130 #endif // BOOST_NUMERIC_INTERVAL_CHECKING_HPP
|