/* 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 #include #include namespace boost { namespace numeric { namespace interval_lib { template 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 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 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