/* Boost interval/interval.hpp header file | |
* | |
* Copyright 2002-2003 Hervé Brönnimann, Guillaume Melquiond, Sylvain Pion | |
* | |
* 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_INTERVAL_INTERVAL_HPP | |
#define BOOST_NUMERIC_INTERVAL_INTERVAL_HPP | |
#include <stdexcept> | |
#include <string> | |
#include <boost/numeric/interval/detail/interval_prototype.hpp> | |
namespace boost { | |
namespace numeric { | |
namespace interval_lib { | |
class comparison_error | |
: public std::runtime_error | |
{ | |
public: | |
comparison_error() | |
: std::runtime_error("boost::interval: uncertain comparison") | |
{ } | |
}; | |
} // namespace interval_lib | |
/* | |
* interval class | |
*/ | |
template<class T, class Policies> | |
class interval | |
{ | |
private: | |
struct interval_holder; | |
struct number_holder; | |
public: | |
typedef T base_type; | |
typedef Policies traits_type; | |
T const &lower() const; | |
T const &upper() const; | |
interval(); | |
interval(T const &v); | |
template<class T1> interval(T1 const &v); | |
interval(T const &l, T const &u); | |
template<class T1, class T2> interval(T1 const &l, T2 const &u); | |
interval(interval<T, Policies> const &r); | |
template<class Policies1> interval(interval<T, Policies1> const &r); | |
template<class T1, class Policies1> interval(interval<T1, Policies1> const &r); | |
interval &operator=(T const &v); | |
template<class T1> interval &operator=(T1 const &v); | |
interval &operator=(interval<T, Policies> const &r); | |
template<class Policies1> interval &operator=(interval<T, Policies1> const &r); | |
template<class T1, class Policies1> interval &operator=(interval<T1, Policies1> const &r); | |
void assign(const T& l, const T& u); | |
static interval empty(); | |
static interval whole(); | |
static interval hull(const T& x, const T& y); | |
interval& operator+= (const T& r); | |
interval& operator+= (const interval& r); | |
interval& operator-= (const T& r); | |
interval& operator-= (const interval& r); | |
interval& operator*= (const T& r); | |
interval& operator*= (const interval& r); | |
interval& operator/= (const T& r); | |
interval& operator/= (const interval& r); | |
bool operator< (const interval_holder& r) const; | |
bool operator> (const interval_holder& r) const; | |
bool operator<= (const interval_holder& r) const; | |
bool operator>= (const interval_holder& r) const; | |
bool operator== (const interval_holder& r) const; | |
bool operator!= (const interval_holder& r) const; | |
bool operator< (const number_holder& r) const; | |
bool operator> (const number_holder& r) const; | |
bool operator<= (const number_holder& r) const; | |
bool operator>= (const number_holder& r) const; | |
bool operator== (const number_holder& r) const; | |
bool operator!= (const number_holder& r) const; | |
// the following is for internal use only, it is not a published interface | |
// nevertheless, it's public because friends don't always work correctly. | |
interval(const T& l, const T& u, bool): low(l), up(u) {} | |
void set_empty(); | |
void set_whole(); | |
void set(const T& l, const T& u); | |
private: | |
struct interval_holder { | |
template<class Policies2> | |
interval_holder(const interval<T, Policies2>& r) | |
: low(r.lower()), up(r.upper()) | |
{ | |
typedef typename Policies2::checking checking2; | |
if (checking2::is_empty(low, up)) | |
throw interval_lib::comparison_error(); | |
} | |
const T& low; | |
const T& up; | |
}; | |
struct number_holder { | |
number_holder(const T& r) : val(r) | |
{ | |
typedef typename Policies::checking checking; | |
if (checking::is_nan(r)) | |
throw interval_lib::comparison_error(); | |
} | |
const T& val; | |
}; | |
typedef typename Policies::checking checking; | |
typedef typename Policies::rounding rounding; | |
T low; | |
T up; | |
}; | |
template<class T, class Policies> inline | |
interval<T, Policies>::interval(): | |
low(static_cast<T>(0)), up(static_cast<T>(0)) | |
{} | |
template<class T, class Policies> inline | |
interval<T, Policies>::interval(T const &v): low(v), up(v) | |
{ | |
if (checking::is_nan(v)) set_empty(); | |
} | |
template<class T, class Policies> template<class T1> inline | |
interval<T, Policies>::interval(T1 const &v) | |
{ | |
if (checking::is_nan(v)) set_empty(); | |
else { | |
rounding rnd; | |
low = rnd.conv_down(v); | |
up = rnd.conv_up (v); | |
} | |
} | |
template<class T, class Policies> template<class T1, class T2> inline | |
interval<T, Policies>::interval(T1 const &l, T2 const &u) | |
{ | |
if (checking::is_nan(l) || checking::is_nan(u) || !(l <= u)) set_empty(); | |
else { | |
rounding rnd; | |
low = rnd.conv_down(l); | |
up = rnd.conv_up (u); | |
} | |
} | |
template<class T, class Policies> inline | |
interval<T, Policies>::interval(T const &l, T const &u): low(l), up(u) | |
{ | |
if (checking::is_nan(l) || checking::is_nan(u) || !(l <= u)) | |
set_empty(); | |
} | |
template<class T, class Policies> inline | |
interval<T, Policies>::interval(interval<T, Policies> const &r): low(r.lower()), up(r.upper()) | |
{} | |
template<class T, class Policies> template<class Policies1> inline | |
interval<T, Policies>::interval(interval<T, Policies1> const &r): low(r.lower()), up(r.upper()) | |
{ | |
typedef typename Policies1::checking checking1; | |
if (checking1::is_empty(r.lower(), r.upper())) set_empty(); | |
} | |
template<class T, class Policies> template<class T1, class Policies1> inline | |
interval<T, Policies>::interval(interval<T1, Policies1> const &r) | |
{ | |
typedef typename Policies1::checking checking1; | |
if (checking1::is_empty(r.lower(), r.upper())) set_empty(); | |
else { | |
rounding rnd; | |
low = rnd.conv_down(r.lower()); | |
up = rnd.conv_up (r.upper()); | |
} | |
} | |
template<class T, class Policies> inline | |
interval<T, Policies> &interval<T, Policies>::operator=(T const &v) | |
{ | |
if (checking::is_nan(v)) set_empty(); | |
else low = up = v; | |
return *this; | |
} | |
template<class T, class Policies> template<class T1> inline | |
interval<T, Policies> &interval<T, Policies>::operator=(T1 const &v) | |
{ | |
if (checking::is_nan(v)) set_empty(); | |
else { | |
rounding rnd; | |
low = rnd.conv_down(v); | |
up = rnd.conv_up (v); | |
} | |
return *this; | |
} | |
template<class T, class Policies> inline | |
interval<T, Policies> &interval<T, Policies>::operator=(interval<T, Policies> const &r) | |
{ | |
low = r.lower(); | |
up = r.upper(); | |
return *this; | |
} | |
template<class T, class Policies> template<class Policies1> inline | |
interval<T, Policies> &interval<T, Policies>::operator=(interval<T, Policies1> const &r) | |
{ | |
typedef typename Policies1::checking checking1; | |
if (checking1::is_empty(r.lower(), r.upper())) set_empty(); | |
else { | |
low = r.lower(); | |
up = r.upper(); | |
} | |
return *this; | |
} | |
template<class T, class Policies> template<class T1, class Policies1> inline | |
interval<T, Policies> &interval<T, Policies>::operator=(interval<T1, Policies1> const &r) | |
{ | |
typedef typename Policies1::checking checking1; | |
if (checking1::is_empty(r.lower(), r.upper())) set_empty(); | |
else { | |
rounding rnd; | |
low = rnd.conv_down(r.lower()); | |
up = rnd.conv_up (r.upper()); | |
} | |
return *this; | |
} | |
template<class T, class Policies> inline | |
void interval<T, Policies>::assign(const T& l, const T& u) | |
{ | |
if (checking::is_nan(l) || checking::is_nan(u) || !(l <= u)) | |
set_empty(); | |
else set(l, u); | |
} | |
template<class T, class Policies> inline | |
void interval<T, Policies>::set(const T& l, const T& u) | |
{ | |
low = l; | |
up = u; | |
} | |
template<class T, class Policies> inline | |
void interval<T, Policies>::set_empty() | |
{ | |
low = checking::empty_lower(); | |
up = checking::empty_upper(); | |
} | |
template<class T, class Policies> inline | |
void interval<T, Policies>::set_whole() | |
{ | |
low = checking::neg_inf(); | |
up = checking::pos_inf(); | |
} | |
template<class T, class Policies> inline | |
interval<T, Policies> interval<T, Policies>::hull(const T& x, const T& y) | |
{ | |
bool bad_x = checking::is_nan(x); | |
bool bad_y = checking::is_nan(y); | |
if (bad_x) | |
if (bad_y) return interval::empty(); | |
else return interval(y, y, true); | |
else | |
if (bad_y) return interval(x, x, true); | |
if (x <= y) return interval(x, y, true); | |
else return interval(y, x, true); | |
} | |
template<class T, class Policies> inline | |
interval<T, Policies> interval<T, Policies>::empty() | |
{ | |
return interval<T, Policies>(checking::empty_lower(), | |
checking::empty_upper(), true); | |
} | |
template<class T, class Policies> inline | |
interval<T, Policies> interval<T, Policies>::whole() | |
{ | |
return interval<T, Policies>(checking::neg_inf(), checking::pos_inf(), true); | |
} | |
template<class T, class Policies> inline | |
const T& interval<T, Policies>::lower() const | |
{ | |
return low; | |
} | |
template<class T, class Policies> inline | |
const T& interval<T, Policies>::upper() const | |
{ | |
return up; | |
} | |
/* | |
* interval/interval comparisons | |
*/ | |
template<class T, class Policies> inline | |
bool interval<T, Policies>::operator< (const interval_holder& r) const | |
{ | |
if (!checking::is_empty(low, up)) { | |
if (up < r.low) return true; | |
else if (low >= r.up) return false; | |
} | |
throw interval_lib::comparison_error(); | |
} | |
template<class T, class Policies> inline | |
bool interval<T, Policies>::operator> (const interval_holder& r) const | |
{ | |
if (!checking::is_empty(low, up)) { | |
if (low > r.up) return true; | |
else if (up <= r.low) return false; | |
} | |
throw interval_lib::comparison_error(); | |
} | |
template<class T, class Policies> inline | |
bool interval<T, Policies>::operator<= (const interval_holder& r) const | |
{ | |
if (!checking::is_empty(low, up)) { | |
if (up <= r.low) return true; | |
else if (low > r.up) return false; | |
} | |
throw interval_lib::comparison_error(); | |
} | |
template<class T, class Policies> inline | |
bool interval<T, Policies>::operator>= (const interval_holder& r) const | |
{ | |
if (!checking::is_empty(low, up)) { | |
if (low >= r.up) return true; | |
else if (up < r.low) return false; | |
} | |
throw interval_lib::comparison_error(); | |
} | |
template<class T, class Policies> inline | |
bool interval<T, Policies>::operator== (const interval_holder& r) const | |
{ | |
if (!checking::is_empty(low, up)) { | |
if (up == r.low && low == r.up) return true; | |
else if (up < r.low || low > r.up) return false; | |
} | |
throw interval_lib::comparison_error(); | |
} | |
template<class T, class Policies> inline | |
bool interval<T, Policies>::operator!= (const interval_holder& r) const | |
{ | |
if (!checking::is_empty(low, up)) { | |
if (up < r.low || low > r.up) return true; | |
else if (up == r.low && low == r.up) return false; | |
} | |
throw interval_lib::comparison_error(); | |
} | |
/* | |
* interval/number comparisons | |
*/ | |
template<class T, class Policies> inline | |
bool interval<T, Policies>::operator< (const number_holder& r) const | |
{ | |
if (!checking::is_empty(low, up)) { | |
if (up < r.val) return true; | |
else if (low >= r.val) return false; | |
} | |
throw interval_lib::comparison_error(); | |
} | |
template<class T, class Policies> inline | |
bool interval<T, Policies>::operator> (const number_holder& r) const | |
{ | |
if (!checking::is_empty(low, up)) { | |
if (low > r.val) return true; | |
else if (up <= r.val) return false; | |
} | |
throw interval_lib::comparison_error(); | |
} | |
template<class T, class Policies> inline | |
bool interval<T, Policies>::operator<= (const number_holder& r) const | |
{ | |
if (!checking::is_empty(low, up)) { | |
if (up <= r.val) return true; | |
else if (low > r.val) return false; | |
} | |
throw interval_lib::comparison_error(); | |
} | |
template<class T, class Policies> inline | |
bool interval<T, Policies>::operator>= (const number_holder& r) const | |
{ | |
if (!checking::is_empty(low, up)) { | |
if (low >= r.val) return true; | |
else if (up < r.val) return false; | |
} | |
throw interval_lib::comparison_error(); | |
} | |
template<class T, class Policies> inline | |
bool interval<T, Policies>::operator== (const number_holder& r) const | |
{ | |
if (!checking::is_empty(low, up)) { | |
if (up == r.val && low == r.val) return true; | |
else if (up < r.val || low > r.val) return false; | |
} | |
throw interval_lib::comparison_error(); | |
} | |
template<class T, class Policies> inline | |
bool interval<T, Policies>::operator!= (const number_holder& r) const | |
{ | |
if (!checking::is_empty(low, up)) { | |
if (up < r.val || low > r.val) return true; | |
else if (up == r.val && low == r.val) return false; | |
} | |
throw interval_lib::comparison_error(); | |
} | |
} // namespace numeric | |
} // namespace boost | |
#endif // BOOST_NUMERIC_INTERVAL_INTERVAL_HPP |