Something went wrong. Try again.
This repository has no description
Something went wrong. Try again.
on-restart-benchmarks software gecode_base gecode third-party boost numeric interval rounded_arith.hpp
5.2 kB · 120 lines
C++
at develop
123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121/* Boost interval/rounded_arith.hpp template implementation 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_ROUNDED_ARITH_HPP#define GECODE_BOOST_NUMERIC_INTERVAL_ROUNDED_ARITH_HPP
#include <gecode/third-party/boost/numeric/interval/rounding.hpp>#include <gecode/third-party/boost/numeric/interval/detail/bugs.hpp>#include <gecode/third-party/boost/config/no_tr1/cmath.hpp>
namespace gecode_boost {namespace numeric {namespace interval_lib {
/* * Three classes of rounding: exact, std, opp * See documentation for details. */
template<class T, class Rounding>struct rounded_arith_exact: Rounding { void init() { } template<class U> T conv_down(U const &v) { return v; } template<class U> T conv_up (U const &v) { return v; } T add_down (const T& x, const T& y) { return x + y; } T add_up (const T& x, const T& y) { return x + y; } T sub_down (const T& x, const T& y) { return x - y; } T sub_up (const T& x, const T& y) { return x - y; } T mul_down (const T& x, const T& y) { return x * y; } T mul_up (const T& x, const T& y) { return x * y; } T div_down (const T& x, const T& y) { return x / y; } T div_up (const T& x, const T& y) { return x / y; } T median (const T& x, const T& y) { return (x + y) / 2; } T sqrt_down(const T& x) { GECODE_BOOST_NUMERIC_INTERVAL_using_math(sqrt); return sqrt(x); } T sqrt_up (const T& x) { GECODE_BOOST_NUMERIC_INTERVAL_using_math(sqrt); return sqrt(x); } T int_down (const T& x) { GECODE_BOOST_NUMERIC_INTERVAL_using_math(floor); return floor(x); } T int_up (const T& x) { GECODE_BOOST_NUMERIC_INTERVAL_using_math(ceil); return ceil(x); }};
template<class T, class Rounding>struct rounded_arith_std: Rounding {# define GECODE_BOOST_DN(EXPR) this->downward(); return this->force_rounding(EXPR)# define GECODE_BOOST_NR(EXPR) this->to_nearest(); return this->force_rounding(EXPR)# define GECODE_BOOST_UP(EXPR) this->upward(); return this->force_rounding(EXPR) void init() { } template<class U> T conv_down(U const &v) { GECODE_BOOST_DN(v); } template<class U> T conv_up (U const &v) { GECODE_BOOST_UP(v); } T add_down(const T& x, const T& y) { GECODE_BOOST_DN(x + y); } T sub_down(const T& x, const T& y) { GECODE_BOOST_DN(x - y); } T mul_down(const T& x, const T& y) { GECODE_BOOST_DN(x * y); } T div_down(const T& x, const T& y) { GECODE_BOOST_DN(x / y); } T add_up (const T& x, const T& y) { GECODE_BOOST_UP(x + y); } T sub_up (const T& x, const T& y) { GECODE_BOOST_UP(x - y); } T mul_up (const T& x, const T& y) { GECODE_BOOST_UP(x * y); } T div_up (const T& x, const T& y) { GECODE_BOOST_UP(x / y); } T median(const T& x, const T& y) { GECODE_BOOST_NR((x + y) / 2); } T sqrt_down(const T& x) { GECODE_BOOST_NUMERIC_INTERVAL_using_math(sqrt); GECODE_BOOST_DN(sqrt(x)); } T sqrt_up (const T& x) { GECODE_BOOST_NUMERIC_INTERVAL_using_math(sqrt); GECODE_BOOST_UP(sqrt(x)); } T int_down(const T& x) { this->downward(); return this->to_int(x); } T int_up (const T& x) { this->upward(); return this->to_int(x); }# undef GECODE_BOOST_DN# undef GECODE_BOOST_NR# undef GECODE_BOOST_UP};
template<class T, class Rounding>struct rounded_arith_opp: Rounding { void init() { this->upward(); }# define GECODE_BOOST_DN(EXPR) \ this->downward(); \ T r = this->force_rounding(EXPR); \ this->upward(); \ return r# define GECODE_BOOST_NR(EXPR) \ this->to_nearest(); \ T r = this->force_rounding(EXPR); \ this->upward(); \ return r# define GECODE_BOOST_UP(EXPR) return this->force_rounding(EXPR)# define GECODE_BOOST_UP_NEG(EXPR) return -this->force_rounding(EXPR) template<class U> T conv_down(U const &v) { GECODE_BOOST_UP_NEG(-v); } template<class U> T conv_up (U const &v) { GECODE_BOOST_UP(v); } T add_down(const T& x, const T& y) { GECODE_BOOST_UP_NEG((-x) - y); } T sub_down(const T& x, const T& y) { GECODE_BOOST_UP_NEG(y - x); } T mul_down(const T& x, const T& y) { GECODE_BOOST_UP_NEG(x * (-y)); } T div_down(const T& x, const T& y) { GECODE_BOOST_UP_NEG(x / (-y)); } T add_up (const T& x, const T& y) { GECODE_BOOST_UP(x + y); } T sub_up (const T& x, const T& y) { GECODE_BOOST_UP(x - y); } T mul_up (const T& x, const T& y) { GECODE_BOOST_UP(x * y); } T div_up (const T& x, const T& y) { GECODE_BOOST_UP(x / y); } T median (const T& x, const T& y) { GECODE_BOOST_NR((x + y) / 2); } T sqrt_down(const T& x) { GECODE_BOOST_NUMERIC_INTERVAL_using_math(sqrt); GECODE_BOOST_DN(sqrt(x)); } T sqrt_up (const T& x) { GECODE_BOOST_NUMERIC_INTERVAL_using_math(sqrt); GECODE_BOOST_UP(sqrt(x)); } T int_down(const T& x) { return -this->to_int(-x); } T int_up (const T& x) { return this->to_int(x); }# undef GECODE_BOOST_DN# undef GECODE_BOOST_NR# undef GECODE_BOOST_UP# undef GECODE_BOOST_UP_NEG};
} // namespace interval_lib} // namespace numeric} // namespace gecode_boost
#endif // GECODE_BOOST_NUMERIC_INTERVAL_ROUNDED_ARITH_HPP