source: NonGTP/Boost/boost/numeric/interval/interval.hpp @ 857

Revision 857, 12.3 KB checked in by igarcia, 19 years ago (diff)
Line 
1/* Boost interval/interval.hpp header file
2 *
3 * Copyright 2002-2003 Hervé Brönnimann, Guillaume Melquiond, Sylvain Pion
4 *
5 * Distributed under the Boost Software License, Version 1.0.
6 * (See accompanying file LICENSE_1_0.txt or
7 * copy at http://www.boost.org/LICENSE_1_0.txt)
8 */
9
10#ifndef BOOST_NUMERIC_INTERVAL_INTERVAL_HPP
11#define BOOST_NUMERIC_INTERVAL_INTERVAL_HPP
12
13#include <stdexcept>
14#include <string>
15#include <boost/numeric/interval/detail/interval_prototype.hpp>
16
17namespace boost {
18namespace numeric {
19
20namespace interval_lib {
21   
22class comparison_error
23  : public std::runtime_error
24{
25public:
26  comparison_error()
27    : std::runtime_error("boost::interval: uncertain comparison")
28  { }
29};
30
31} // namespace interval_lib
32
33/*
34 * interval class
35 */
36
37template<class T, class Policies>
38class interval
39{
40private:
41  struct interval_holder;
42  struct number_holder;
43public:
44  typedef T base_type;
45  typedef Policies traits_type;
46
47  T const &lower() const;
48  T const &upper() const;
49
50  interval();
51  interval(T const &v);
52  template<class T1> interval(T1 const &v);
53  interval(T const &l, T const &u);
54  template<class T1, class T2> interval(T1 const &l, T2 const &u);
55  interval(interval<T, Policies> const &r);
56  template<class Policies1> interval(interval<T, Policies1> const &r);
57  template<class T1, class Policies1> interval(interval<T1, Policies1> const &r);
58
59  interval &operator=(T const &v);
60  template<class T1> interval &operator=(T1 const &v);
61  interval &operator=(interval<T, Policies> const &r);
62  template<class Policies1> interval &operator=(interval<T, Policies1> const &r);
63  template<class T1, class Policies1> interval &operator=(interval<T1, Policies1> const &r);
64 
65  void assign(const T& l, const T& u);
66
67  static interval empty();
68  static interval whole();
69  static interval hull(const T& x, const T& y);
70
71  interval& operator+= (const T& r);
72  interval& operator+= (const interval& r);
73  interval& operator-= (const T& r);
74  interval& operator-= (const interval& r);
75  interval& operator*= (const T& r);
76  interval& operator*= (const interval& r);
77  interval& operator/= (const T& r);
78  interval& operator/= (const interval& r);
79
80  bool operator< (const interval_holder& r) const;
81  bool operator> (const interval_holder& r) const;
82  bool operator<= (const interval_holder& r) const;
83  bool operator>= (const interval_holder& r) const;
84  bool operator== (const interval_holder& r) const;
85  bool operator!= (const interval_holder& r) const;
86
87  bool operator< (const number_holder& r) const;
88  bool operator> (const number_holder& r) const;
89  bool operator<= (const number_holder& r) const;
90  bool operator>= (const number_holder& r) const;
91  bool operator== (const number_holder& r) const;
92  bool operator!= (const number_holder& r) const;
93
94  // the following is for internal use only, it is not a published interface
95  // nevertheless, it's public because friends don't always work correctly.
96  interval(const T& l, const T& u, bool): low(l), up(u) {}
97  void set_empty();
98  void set_whole();
99  void set(const T& l, const T& u);
100
101private:
102  struct interval_holder {
103    template<class Policies2>
104    interval_holder(const interval<T, Policies2>& r)
105      : low(r.lower()), up(r.upper())
106    {
107      typedef typename Policies2::checking checking2;
108      if (checking2::is_empty(low, up))
109        throw interval_lib::comparison_error();
110    }
111
112    const T& low;
113    const T& up;
114  };
115
116  struct number_holder {
117    number_holder(const T& r) : val(r)
118    {
119      typedef typename Policies::checking checking;
120      if (checking::is_nan(r))
121        throw interval_lib::comparison_error();
122    }
123   
124    const T& val;
125  };
126
127  typedef typename Policies::checking checking;
128  typedef typename Policies::rounding rounding;
129
130  T low;
131  T up;
132};
133
134template<class T, class Policies> inline
135interval<T, Policies>::interval():
136  low(static_cast<T>(0)), up(static_cast<T>(0))
137{}
138
139template<class T, class Policies> inline
140interval<T, Policies>::interval(T const &v): low(v), up(v)
141{
142  if (checking::is_nan(v)) set_empty();
143}
144
145template<class T, class Policies> template<class T1> inline
146interval<T, Policies>::interval(T1 const &v)
147{
148  if (checking::is_nan(v)) set_empty();
149  else {
150    rounding rnd;
151    low = rnd.conv_down(v);
152    up  = rnd.conv_up  (v);
153  }
154}
155
156template<class T, class Policies> template<class T1, class T2> inline
157interval<T, Policies>::interval(T1 const &l, T2 const &u)
158{
159  if (checking::is_nan(l) || checking::is_nan(u) || !(l <= u)) set_empty();
160  else {
161    rounding rnd;
162    low = rnd.conv_down(l);
163    up  = rnd.conv_up  (u);
164  }
165}
166
167template<class T, class Policies> inline
168interval<T, Policies>::interval(T const &l, T const &u): low(l), up(u)
169{
170  if (checking::is_nan(l) || checking::is_nan(u) || !(l <= u))
171    set_empty();
172}
173
174
175template<class T, class Policies> inline
176interval<T, Policies>::interval(interval<T, Policies> const &r): low(r.lower()), up(r.upper())
177{}
178
179template<class T, class Policies> template<class Policies1> inline
180interval<T, Policies>::interval(interval<T, Policies1> const &r): low(r.lower()), up(r.upper())
181{
182  typedef typename Policies1::checking checking1;
183  if (checking1::is_empty(r.lower(), r.upper())) set_empty();
184}
185
186template<class T, class Policies> template<class T1, class Policies1> inline
187interval<T, Policies>::interval(interval<T1, Policies1> const &r)
188{
189  typedef typename Policies1::checking checking1;
190  if (checking1::is_empty(r.lower(), r.upper())) set_empty();
191  else {
192    rounding rnd;
193    low = rnd.conv_down(r.lower());
194    up  = rnd.conv_up  (r.upper());
195  }
196}
197
198template<class T, class Policies> inline
199interval<T, Policies> &interval<T, Policies>::operator=(T const &v)
200{
201  if (checking::is_nan(v)) set_empty();
202  else low = up = v;
203  return *this;
204}
205
206template<class T, class Policies> template<class T1> inline
207interval<T, Policies> &interval<T, Policies>::operator=(T1 const &v)
208{
209  if (checking::is_nan(v)) set_empty();
210  else {
211    rounding rnd;
212    low = rnd.conv_down(v);
213    up  = rnd.conv_up  (v);
214  }
215  return *this;
216}
217
218template<class T, class Policies> inline
219interval<T, Policies> &interval<T, Policies>::operator=(interval<T, Policies> const &r)
220{
221  low = r.lower();
222  up  = r.upper();
223  return *this;
224}
225
226template<class T, class Policies> template<class Policies1> inline
227interval<T, Policies> &interval<T, Policies>::operator=(interval<T, Policies1> const &r)
228{
229  typedef typename Policies1::checking checking1;
230  if (checking1::is_empty(r.lower(), r.upper())) set_empty();
231  else {
232    low = r.lower();
233    up  = r.upper();
234  }
235  return *this;
236}
237
238template<class T, class Policies> template<class T1, class Policies1> inline
239interval<T, Policies> &interval<T, Policies>::operator=(interval<T1, Policies1> const &r)
240{
241  typedef typename Policies1::checking checking1;
242  if (checking1::is_empty(r.lower(), r.upper())) set_empty();
243  else {
244    rounding rnd;
245    low = rnd.conv_down(r.lower());
246    up  = rnd.conv_up  (r.upper());
247  }
248  return *this;
249}
250
251template<class T, class Policies> inline
252void interval<T, Policies>::assign(const T& l, const T& u)
253{
254  if (checking::is_nan(l) || checking::is_nan(u) || !(l <= u))
255    set_empty();
256  else set(l, u);
257}
258
259template<class T, class Policies> inline
260void interval<T, Policies>::set(const T& l, const T& u)
261{
262  low = l;
263  up  = u;
264}
265
266template<class T, class Policies> inline
267void interval<T, Policies>::set_empty()
268{
269  low = checking::empty_lower();
270  up  = checking::empty_upper();
271}
272
273template<class T, class Policies> inline
274void interval<T, Policies>::set_whole()
275{
276  low = checking::neg_inf();
277  up  = checking::pos_inf();
278}
279
280template<class T, class Policies> inline
281interval<T, Policies> interval<T, Policies>::hull(const T& x, const T& y)
282{
283  bool bad_x = checking::is_nan(x);
284  bool bad_y = checking::is_nan(y);
285  if (bad_x)
286    if (bad_y) return interval::empty();
287    else       return interval(y, y, true);
288  else
289    if (bad_y) return interval(x, x, true);
290  if (x <= y) return interval(x, y, true);
291  else        return interval(y, x, true);
292}
293
294template<class T, class Policies> inline
295interval<T, Policies> interval<T, Policies>::empty()
296{
297  return interval<T, Policies>(checking::empty_lower(),
298                               checking::empty_upper(), true);
299}
300
301template<class T, class Policies> inline
302interval<T, Policies> interval<T, Policies>::whole()
303{
304  return interval<T, Policies>(checking::neg_inf(), checking::pos_inf(), true);
305}
306
307template<class T, class Policies> inline
308const T& interval<T, Policies>::lower() const
309{
310  return low;
311}
312
313template<class T, class Policies> inline
314const T& interval<T, Policies>::upper() const
315{
316  return up;
317}
318
319/*
320 * interval/interval comparisons
321 */
322
323template<class T, class Policies> inline
324bool interval<T, Policies>::operator< (const interval_holder& r) const
325{
326  if (!checking::is_empty(low, up)) {
327    if (up < r.low) return true;
328    else if (low >= r.up) return false;
329  }
330  throw interval_lib::comparison_error();
331}
332
333template<class T, class Policies> inline
334bool interval<T, Policies>::operator> (const interval_holder& r) const
335{
336  if (!checking::is_empty(low, up)) {
337    if (low > r.up) return true;
338    else if (up <= r.low) return false;
339  }
340  throw interval_lib::comparison_error();
341}
342
343template<class T, class Policies> inline
344bool interval<T, Policies>::operator<= (const interval_holder& r) const
345{
346  if (!checking::is_empty(low, up)) {
347    if (up <= r.low) return true;
348    else if (low > r.up) return false;
349  }
350  throw interval_lib::comparison_error();
351}
352
353template<class T, class Policies> inline
354bool interval<T, Policies>::operator>= (const interval_holder& r) const
355{
356  if (!checking::is_empty(low, up)) {
357    if (low >= r.up) return true;
358    else if (up < r.low) return false;
359  }
360  throw interval_lib::comparison_error();
361}
362
363template<class T, class Policies> inline
364bool interval<T, Policies>::operator== (const interval_holder& r) const
365{
366  if (!checking::is_empty(low, up)) {
367    if (up == r.low && low == r.up) return true;
368    else if (up < r.low || low > r.up) return false;
369  }
370  throw interval_lib::comparison_error();
371}
372
373template<class T, class Policies> inline
374bool interval<T, Policies>::operator!= (const interval_holder& r) const
375{
376  if (!checking::is_empty(low, up)) {
377    if (up < r.low || low > r.up) return true;
378    else if (up == r.low && low == r.up) return false;
379  }
380  throw interval_lib::comparison_error();
381}
382
383/*
384 * interval/number comparisons
385 */
386
387template<class T, class Policies> inline
388bool interval<T, Policies>::operator< (const number_holder& r) const
389{
390  if (!checking::is_empty(low, up)) {
391    if (up < r.val) return true;
392    else if (low >= r.val) return false;
393  }
394  throw interval_lib::comparison_error();
395}
396
397template<class T, class Policies> inline
398bool interval<T, Policies>::operator> (const number_holder& r) const
399{
400  if (!checking::is_empty(low, up)) {
401    if (low > r.val) return true;
402    else if (up <= r.val) return false;
403  }
404  throw interval_lib::comparison_error();
405}
406
407template<class T, class Policies> inline
408bool interval<T, Policies>::operator<= (const number_holder& r) const
409{
410  if (!checking::is_empty(low, up)) {
411    if (up <= r.val) return true;
412    else if (low > r.val) return false;
413  }
414  throw interval_lib::comparison_error();
415}
416
417template<class T, class Policies> inline
418bool interval<T, Policies>::operator>= (const number_holder& r) const
419{
420  if (!checking::is_empty(low, up)) {
421    if (low >= r.val) return true;
422    else if (up < r.val) return false;
423  }
424  throw interval_lib::comparison_error();
425}
426
427template<class T, class Policies> inline
428bool interval<T, Policies>::operator== (const number_holder& r) const
429{
430  if (!checking::is_empty(low, up)) {
431    if (up == r.val && low == r.val) return true;
432    else if (up < r.val || low > r.val) return false;
433  }
434  throw interval_lib::comparison_error();
435}
436
437template<class T, class Policies> inline
438bool interval<T, Policies>::operator!= (const number_holder& r) const
439{
440  if (!checking::is_empty(low, up)) {
441    if (up < r.val || low > r.val) return true;
442    else if (up == r.val && low == r.val) return false;
443  }
444  throw interval_lib::comparison_error();
445}
446
447} // namespace numeric
448} // namespace boost
449
450#endif // BOOST_NUMERIC_INTERVAL_INTERVAL_HPP
Note: See TracBrowser for help on using the repository browser.