Chris@16: /* Boost interval/rounded_transc.hpp template implementation file
Chris@16:  *
Chris@16:  * Copyright 2002-2003 Hervé Brönnimann, Guillaume Melquiond, Sylvain Pion
Chris@16:  *
Chris@16:  * Distributed under the Boost Software License, Version 1.0.
Chris@16:  * (See accompanying file LICENSE_1_0.txt or
Chris@16:  * copy at http://www.boost.org/LICENSE_1_0.txt)
Chris@16:  */
Chris@16: 
Chris@16: #ifndef BOOST_NUMERIC_INTERVAL_ROUNDED_TRANSC_HPP
Chris@16: #define BOOST_NUMERIC_INTERVAL_ROUNDED_TRANSC_HPP
Chris@16: 
Chris@16: #include <boost/numeric/interval/rounding.hpp>
Chris@16: #include <boost/numeric/interval/detail/bugs.hpp>
Chris@16: #include <boost/config/no_tr1/cmath.hpp>
Chris@16: 
Chris@16: namespace boost {
Chris@16: namespace numeric {
Chris@16: namespace interval_lib {
Chris@16: 
Chris@16: template<class T, class Rounding>
Chris@16: struct rounded_transc_exact: Rounding
Chris@16: {
Chris@16: # define BOOST_NUMERIC_INTERVAL_new_func(f) \
Chris@16:     T f##_down(const T& x) { BOOST_NUMERIC_INTERVAL_using_math(f); return f(x); } \
Chris@16:     T f##_up  (const T& x) { BOOST_NUMERIC_INTERVAL_using_math(f); return f(x); }
Chris@16:   BOOST_NUMERIC_INTERVAL_new_func(exp)
Chris@16:   BOOST_NUMERIC_INTERVAL_new_func(log)
Chris@16:   BOOST_NUMERIC_INTERVAL_new_func(sin)
Chris@16:   BOOST_NUMERIC_INTERVAL_new_func(cos)
Chris@16:   BOOST_NUMERIC_INTERVAL_new_func(tan)
Chris@16:   BOOST_NUMERIC_INTERVAL_new_func(asin)
Chris@16:   BOOST_NUMERIC_INTERVAL_new_func(acos)
Chris@16:   BOOST_NUMERIC_INTERVAL_new_func(atan)
Chris@16:   BOOST_NUMERIC_INTERVAL_new_func(sinh)
Chris@16:   BOOST_NUMERIC_INTERVAL_new_func(cosh)
Chris@16:   BOOST_NUMERIC_INTERVAL_new_func(tanh)
Chris@16: # undef BOOST_NUMERIC_INTERVAL_new_func
Chris@16: # define BOOST_NUMERIC_INTERVAL_new_func(f) \
Chris@16:     T f##_down(const T& x) { BOOST_NUMERIC_INTERVAL_using_ahyp(f); return f(x); } \
Chris@16:     T f##_up  (const T& x) { BOOST_NUMERIC_INTERVAL_using_ahyp(f); return f(x); }
Chris@16:   BOOST_NUMERIC_INTERVAL_new_func(asinh)
Chris@16:   BOOST_NUMERIC_INTERVAL_new_func(acosh)
Chris@16:   BOOST_NUMERIC_INTERVAL_new_func(atanh)
Chris@16: # undef BOOST_NUMERIC_INTERVAL_new_func
Chris@16: };
Chris@16:   
Chris@16: template<class T, class Rounding>
Chris@16: struct rounded_transc_std: Rounding
Chris@16: {
Chris@16: # define BOOST_NUMERIC_INTERVAL_new_func(f) \
Chris@16:     T f##_down(const T& x) \
Chris@16:     { BOOST_NUMERIC_INTERVAL_using_math(f); \
Chris@16:       this->downward(); return this->force_rounding(f(x)); } \
Chris@16:     T f##_up  (const T& x) \
Chris@16:     { BOOST_NUMERIC_INTERVAL_using_math(f); \
Chris@16:       this->upward(); return this->force_rounding(f(x)); }
Chris@16:   BOOST_NUMERIC_INTERVAL_new_func(exp)
Chris@16:   BOOST_NUMERIC_INTERVAL_new_func(log)
Chris@16:   BOOST_NUMERIC_INTERVAL_new_func(sin)
Chris@16:   BOOST_NUMERIC_INTERVAL_new_func(cos)
Chris@16:   BOOST_NUMERIC_INTERVAL_new_func(tan)
Chris@16:   BOOST_NUMERIC_INTERVAL_new_func(asin)
Chris@16:   BOOST_NUMERIC_INTERVAL_new_func(acos)
Chris@16:   BOOST_NUMERIC_INTERVAL_new_func(atan)
Chris@16:   BOOST_NUMERIC_INTERVAL_new_func(sinh)
Chris@16:   BOOST_NUMERIC_INTERVAL_new_func(cosh)
Chris@16:   BOOST_NUMERIC_INTERVAL_new_func(tanh)
Chris@16: # undef BOOST_NUMERIC_INTERVAL_new_func
Chris@16: # define BOOST_NUMERIC_INTERVAL_new_func(f) \
Chris@16:     T f##_down(const T& x) \
Chris@16:     { BOOST_NUMERIC_INTERVAL_using_ahyp(f); \
Chris@16:       this->downward(); return this->force_rounding(f(x)); } \
Chris@16:     T f##_up  (const T& x) \
Chris@16:     { BOOST_NUMERIC_INTERVAL_using_ahyp(f); \
Chris@16:       this->upward(); return this->force_rounding(f(x)); }
Chris@16:   BOOST_NUMERIC_INTERVAL_new_func(asinh)
Chris@16:   BOOST_NUMERIC_INTERVAL_new_func(acosh)
Chris@16:   BOOST_NUMERIC_INTERVAL_new_func(atanh)
Chris@16: # undef BOOST_NUMERIC_INTERVAL_new_func
Chris@16: };
Chris@16: 
Chris@16: template<class T, class Rounding>
Chris@16: struct rounded_transc_opp: Rounding
Chris@16: {
Chris@16: # define BOOST_NUMERIC_INTERVAL_new_func(f) \
Chris@16:     T f##_down(const T& x) \
Chris@16:     { BOOST_NUMERIC_INTERVAL_using_math(f); \
Chris@16:       this->downward(); T y = this->force_rounding(f(x)); \
Chris@16:       this->upward(); return y; } \
Chris@16:     T f##_up  (const T& x) \
Chris@16:     { BOOST_NUMERIC_INTERVAL_using_math(f); \
Chris@16:       return this->force_rounding(f(x)); }
Chris@16:   BOOST_NUMERIC_INTERVAL_new_func(exp)
Chris@16:   BOOST_NUMERIC_INTERVAL_new_func(log)
Chris@16:   BOOST_NUMERIC_INTERVAL_new_func(cos)
Chris@16:   BOOST_NUMERIC_INTERVAL_new_func(acos)
Chris@16:   BOOST_NUMERIC_INTERVAL_new_func(cosh)
Chris@16: # undef BOOST_NUMERIC_INTERVAL_new_func
Chris@16: # define BOOST_NUMERIC_INTERVAL_new_func(f) \
Chris@16:     T f##_down(const T& x) \
Chris@16:     { BOOST_NUMERIC_INTERVAL_using_math(f); \
Chris@16:       return -this->force_rounding(-f(x)); } \
Chris@16:     T f##_up  (const T& x) \
Chris@16:     { BOOST_NUMERIC_INTERVAL_using_math(f); \
Chris@16:       return this->force_rounding(f(x)); }
Chris@16:   BOOST_NUMERIC_INTERVAL_new_func(sin)
Chris@16:   BOOST_NUMERIC_INTERVAL_new_func(tan)
Chris@16:   BOOST_NUMERIC_INTERVAL_new_func(asin)
Chris@16:   BOOST_NUMERIC_INTERVAL_new_func(atan)
Chris@16:   BOOST_NUMERIC_INTERVAL_new_func(sinh)
Chris@16:   BOOST_NUMERIC_INTERVAL_new_func(tanh)
Chris@16: # undef BOOST_NUMERIC_INTERVAL_new_func
Chris@16: # define BOOST_NUMERIC_INTERVAL_new_func(f) \
Chris@16:     T f##_down(const T& x) \
Chris@16:     { BOOST_NUMERIC_INTERVAL_using_ahyp(f); \
Chris@16:       this->downward(); T y = this->force_rounding(f(x)); \
Chris@16:       this->upward(); return y; } \
Chris@16:     T f##_up  (const T& x) \
Chris@16:     { BOOST_NUMERIC_INTERVAL_using_ahyp(f); \
Chris@16:       return this->force_rounding(f(x)); }
Chris@16:   BOOST_NUMERIC_INTERVAL_new_func(asinh)
Chris@16:   BOOST_NUMERIC_INTERVAL_new_func(atanh)
Chris@16: # undef BOOST_NUMERIC_INTERVAL_new_func
Chris@16: # define BOOST_NUMERIC_INTERVAL_new_func(f) \
Chris@16:     T f##_down(const T& x) \
Chris@16:     { BOOST_NUMERIC_INTERVAL_using_ahyp(f); \
Chris@16:       return -this->force_rounding(-f(x)); } \
Chris@16:     T f##_up  (const T& x) \
Chris@16:     { BOOST_NUMERIC_INTERVAL_using_ahyp(f); \
Chris@16:       return this->force_rounding(f(x)); }
Chris@16:   BOOST_NUMERIC_INTERVAL_new_func(acosh)
Chris@16: # undef BOOST_NUMERIC_INTERVAL_new_func
Chris@16: };
Chris@16:   
Chris@16: } // namespace interval_lib
Chris@16: } // namespace numeric
Chris@16: } // namespace boost
Chris@16: 
Chris@16: #endif // BOOST_NUMERIC_INTERVAL_ROUNDED_TRANSC_HPP