Something went wrong. Try again.
This repository has no description
Something went wrong. Try again.
12 kB · 451 lines
C++
at develop
123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229230231232233234235236237238239240241242243244245246247248249250251252253254255256257258259260261262263264265266267268269270271272273274275276277278279280281282283284285286287288289290291292293294295296297298299300301302303304305306307308309310311312313314315316317318319320321322323324325326327328329330331332333334335336337338339340341342343344345346347348349350351352353354355356357358359360361362363364365366367368369370371372373374375376377378379380381382383384385386387388389390391392393394395396397398399400401402403404405406407408409410411412413414415416417418419420421422423424425426427428429430431432433434435436437438439440441442443444445446447448449450451452/* 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 GECODE_BOOST_NUMERIC_INTERVAL_INTERVAL_HPP#define GECODE_BOOST_NUMERIC_INTERVAL_INTERVAL_HPP
#include <stdexcept>#include <string>#include <gecode/third-party/boost/numeric/interval/detail/interval_prototype.hpp>
namespace gecode_boost {namespace numeric {
namespace interval_lib { class comparison_error : public std::runtime_error {public: comparison_error() : std::runtime_error("gecode_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> inlineinterval<T, Policies>::interval(): low(static_cast<T>(0)), up(static_cast<T>(0)){}
template<class T, class Policies> inlineinterval<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> inlineinterval<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> inlineinterval<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> inlineinterval<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> inlineinterval<T, Policies>::interval(interval<T, Policies> const &r): low(r.lower()), up(r.upper()){}
template<class T, class Policies> template<class Policies1> inlineinterval<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> inlineinterval<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> inlineinterval<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> inlineinterval<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> inlineinterval<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> inlineinterval<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> inlineinterval<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> inlinevoid 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> inlinevoid interval<T, Policies>::set(const T& l, const T& u){ low = l; up = u;}
template<class T, class Policies> inlinevoid interval<T, Policies>::set_empty(){ low = checking::empty_lower(); up = checking::empty_upper();}
template<class T, class Policies> inlinevoid interval<T, Policies>::set_whole(){ low = checking::neg_inf(); up = checking::pos_inf();}
template<class T, class Policies> inlineinterval<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> inlineinterval<T, Policies> interval<T, Policies>::empty(){ return interval<T, Policies>(checking::empty_lower(), checking::empty_upper(), true);}
template<class T, class Policies> inlineinterval<T, Policies> interval<T, Policies>::whole(){ return interval<T, Policies>(checking::neg_inf(), checking::pos_inf(), true);}
template<class T, class Policies> inlineconst T& interval<T, Policies>::lower() const{ return low;}
template<class T, class Policies> inlineconst T& interval<T, Policies>::upper() const{ return up;}
/* * interval/interval comparisons */
template<class T, class Policies> inlinebool 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> inlinebool 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> inlinebool 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> inlinebool 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> inlinebool interval<T, Policies>::operator== (const interval_holder& r) const{ if (!checking::is_empty(low, up)) {//Boost bug ? if (up == r.low && low == r.up) return true; if (up == r.up && low == r.low) return true; else if (up < r.low || low > r.up) return false; } throw interval_lib::comparison_error();}
template<class T, class Policies> inlinebool 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> inlinebool 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> inlinebool 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> inlinebool 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> inlinebool 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> inlinebool 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> inlinebool 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 gecode_boost
#endif // GECODE_BOOST_NUMERIC_INTERVAL_INTERVAL_HPP