/* Boost interval/rounded_transc.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 BOOST_NUMERIC_INTERVAL_ROUNDED_TRANSC_HPP | |
#define BOOST_NUMERIC_INTERVAL_ROUNDED_TRANSC_HPP | |
#include <boost/numeric/interval/rounding.hpp> | |
#include <boost/numeric/interval/detail/bugs.hpp> | |
#include <boost/config/no_tr1/cmath.hpp> | |
namespace boost { | |
namespace numeric { | |
namespace interval_lib { | |
template<class T, class Rounding> | |
struct rounded_transc_exact: Rounding | |
{ | |
# define BOOST_NUMERIC_INTERVAL_new_func(f) \ | |
T f##_down(const T& x) { BOOST_NUMERIC_INTERVAL_using_math(f); return f(x); } \ | |
T f##_up (const T& x) { BOOST_NUMERIC_INTERVAL_using_math(f); return f(x); } | |
BOOST_NUMERIC_INTERVAL_new_func(exp) | |
BOOST_NUMERIC_INTERVAL_new_func(log) | |
BOOST_NUMERIC_INTERVAL_new_func(sin) | |
BOOST_NUMERIC_INTERVAL_new_func(cos) | |
BOOST_NUMERIC_INTERVAL_new_func(tan) | |
BOOST_NUMERIC_INTERVAL_new_func(asin) | |
BOOST_NUMERIC_INTERVAL_new_func(acos) | |
BOOST_NUMERIC_INTERVAL_new_func(atan) | |
BOOST_NUMERIC_INTERVAL_new_func(sinh) | |
BOOST_NUMERIC_INTERVAL_new_func(cosh) | |
BOOST_NUMERIC_INTERVAL_new_func(tanh) | |
# undef BOOST_NUMERIC_INTERVAL_new_func | |
# define BOOST_NUMERIC_INTERVAL_new_func(f) \ | |
T f##_down(const T& x) { BOOST_NUMERIC_INTERVAL_using_ahyp(f); return f(x); } \ | |
T f##_up (const T& x) { BOOST_NUMERIC_INTERVAL_using_ahyp(f); return f(x); } | |
BOOST_NUMERIC_INTERVAL_new_func(asinh) | |
BOOST_NUMERIC_INTERVAL_new_func(acosh) | |
BOOST_NUMERIC_INTERVAL_new_func(atanh) | |
# undef BOOST_NUMERIC_INTERVAL_new_func | |
}; | |
template<class T, class Rounding> | |
struct rounded_transc_std: Rounding | |
{ | |
# define BOOST_NUMERIC_INTERVAL_new_func(f) \ | |
T f##_down(const T& x) \ | |
{ BOOST_NUMERIC_INTERVAL_using_math(f); \ | |
this->downward(); return this->force_rounding(f(x)); } \ | |
T f##_up (const T& x) \ | |
{ BOOST_NUMERIC_INTERVAL_using_math(f); \ | |
this->upward(); return this->force_rounding(f(x)); } | |
BOOST_NUMERIC_INTERVAL_new_func(exp) | |
BOOST_NUMERIC_INTERVAL_new_func(log) | |
BOOST_NUMERIC_INTERVAL_new_func(sin) | |
BOOST_NUMERIC_INTERVAL_new_func(cos) | |
BOOST_NUMERIC_INTERVAL_new_func(tan) | |
BOOST_NUMERIC_INTERVAL_new_func(asin) | |
BOOST_NUMERIC_INTERVAL_new_func(acos) | |
BOOST_NUMERIC_INTERVAL_new_func(atan) | |
BOOST_NUMERIC_INTERVAL_new_func(sinh) | |
BOOST_NUMERIC_INTERVAL_new_func(cosh) | |
BOOST_NUMERIC_INTERVAL_new_func(tanh) | |
# undef BOOST_NUMERIC_INTERVAL_new_func | |
# define BOOST_NUMERIC_INTERVAL_new_func(f) \ | |
T f##_down(const T& x) \ | |
{ BOOST_NUMERIC_INTERVAL_using_ahyp(f); \ | |
this->downward(); return this->force_rounding(f(x)); } \ | |
T f##_up (const T& x) \ | |
{ BOOST_NUMERIC_INTERVAL_using_ahyp(f); \ | |
this->upward(); return this->force_rounding(f(x)); } | |
BOOST_NUMERIC_INTERVAL_new_func(asinh) | |
BOOST_NUMERIC_INTERVAL_new_func(acosh) | |
BOOST_NUMERIC_INTERVAL_new_func(atanh) | |
# undef BOOST_NUMERIC_INTERVAL_new_func | |
}; | |
template<class T, class Rounding> | |
struct rounded_transc_opp: Rounding | |
{ | |
# define BOOST_NUMERIC_INTERVAL_new_func(f) \ | |
T f##_down(const T& x) \ | |
{ BOOST_NUMERIC_INTERVAL_using_math(f); \ | |
this->downward(); T y = this->force_rounding(f(x)); \ | |
this->upward(); return y; } \ | |
T f##_up (const T& x) \ | |
{ BOOST_NUMERIC_INTERVAL_using_math(f); \ | |
return this->force_rounding(f(x)); } | |
BOOST_NUMERIC_INTERVAL_new_func(exp) | |
BOOST_NUMERIC_INTERVAL_new_func(log) | |
BOOST_NUMERIC_INTERVAL_new_func(cos) | |
BOOST_NUMERIC_INTERVAL_new_func(acos) | |
BOOST_NUMERIC_INTERVAL_new_func(cosh) | |
# undef BOOST_NUMERIC_INTERVAL_new_func | |
# define BOOST_NUMERIC_INTERVAL_new_func(f) \ | |
T f##_down(const T& x) \ | |
{ BOOST_NUMERIC_INTERVAL_using_math(f); \ | |
return -this->force_rounding(-f(x)); } \ | |
T f##_up (const T& x) \ | |
{ BOOST_NUMERIC_INTERVAL_using_math(f); \ | |
return this->force_rounding(f(x)); } | |
BOOST_NUMERIC_INTERVAL_new_func(sin) | |
BOOST_NUMERIC_INTERVAL_new_func(tan) | |
BOOST_NUMERIC_INTERVAL_new_func(asin) | |
BOOST_NUMERIC_INTERVAL_new_func(atan) | |
BOOST_NUMERIC_INTERVAL_new_func(sinh) | |
BOOST_NUMERIC_INTERVAL_new_func(tanh) | |
# undef BOOST_NUMERIC_INTERVAL_new_func | |
# define BOOST_NUMERIC_INTERVAL_new_func(f) \ | |
T f##_down(const T& x) \ | |
{ BOOST_NUMERIC_INTERVAL_using_ahyp(f); \ | |
this->downward(); T y = this->force_rounding(f(x)); \ | |
this->upward(); return y; } \ | |
T f##_up (const T& x) \ | |
{ BOOST_NUMERIC_INTERVAL_using_ahyp(f); \ | |
return this->force_rounding(f(x)); } | |
BOOST_NUMERIC_INTERVAL_new_func(asinh) | |
BOOST_NUMERIC_INTERVAL_new_func(atanh) | |
# undef BOOST_NUMERIC_INTERVAL_new_func | |
# define BOOST_NUMERIC_INTERVAL_new_func(f) \ | |
T f##_down(const T& x) \ | |
{ BOOST_NUMERIC_INTERVAL_using_ahyp(f); \ | |
return -this->force_rounding(-f(x)); } \ | |
T f##_up (const T& x) \ | |
{ BOOST_NUMERIC_INTERVAL_using_ahyp(f); \ | |
return this->force_rounding(f(x)); } | |
BOOST_NUMERIC_INTERVAL_new_func(acosh) | |
# undef BOOST_NUMERIC_INTERVAL_new_func | |
}; | |
} // namespace interval_lib | |
} // namespace numeric | |
} // namespace boost | |
#endif // BOOST_NUMERIC_INTERVAL_ROUNDED_TRANSC_HPP |