123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140 |
- /* 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
|