boost/numeric/interval/rounded_transc.hpp

   1 /* Boost interval/rounded_transc.hpp template implementation file
   2  *
   3  * Copyright 2002-2003 Hervé Brönnimann, Guillaume Melquiond, Sylvain Pion
   4  *
   5  * Distributed under the Boost Software License, Version 1.0.
   6  * (See accompanying file LICENSE_1_0.txt or
   7  * copy at http://www.boost.org/LICENSE_1_0.txt)
   8  */
   9
  10 #ifndef BOOST_NUMERIC_INTERVAL_ROUNDED_TRANSC_HPP
  11 #define BOOST_NUMERIC_INTERVAL_ROUNDED_TRANSC_HPP
  12
  13 #include <boost/numeric/interval/rounding.hpp>
  14 #include <boost/numeric/interval/detail/bugs.hpp>
  15 #include <boost/config/no_tr1/cmath.hpp>
  16
  17 namespace boost {
  18 namespace numeric {
  19 namespace interval_lib {
  20
  21 template<class T, class Rounding>
  22 struct rounded_transc_exact: Rounding
  23 {
  24 # define BOOST_NUMERIC_INTERVAL_new_func(f) \
  25     T f##_down(const T& x) { BOOST_NUMERIC_INTERVAL_using_math(f); return f(x); } \
  26     T f##_up  (const T& x) { BOOST_NUMERIC_INTERVAL_using_math(f); return f(x); }
  27   BOOST_NUMERIC_INTERVAL_new_func(exp)
  28   BOOST_NUMERIC_INTERVAL_new_func(log)
  29   BOOST_NUMERIC_INTERVAL_new_func(sin)
  30   BOOST_NUMERIC_INTERVAL_new_func(cos)
  31   BOOST_NUMERIC_INTERVAL_new_func(tan)
  32   BOOST_NUMERIC_INTERVAL_new_func(asin)
  33   BOOST_NUMERIC_INTERVAL_new_func(acos)
  34   BOOST_NUMERIC_INTERVAL_new_func(atan)
  35   BOOST_NUMERIC_INTERVAL_new_func(sinh)
  36   BOOST_NUMERIC_INTERVAL_new_func(cosh)
  37   BOOST_NUMERIC_INTERVAL_new_func(tanh)
  38 # undef BOOST_NUMERIC_INTERVAL_new_func
  39 # define BOOST_NUMERIC_INTERVAL_new_func(f) \
  40     T f##_down(const T& x) { BOOST_NUMERIC_INTERVAL_using_ahyp(f); return f(x); } \
  41     T f##_up  (const T& x) { BOOST_NUMERIC_INTERVAL_using_ahyp(f); return f(x); }
  42   BOOST_NUMERIC_INTERVAL_new_func(asinh)
  43   BOOST_NUMERIC_INTERVAL_new_func(acosh)
  44   BOOST_NUMERIC_INTERVAL_new_func(atanh)
  45 # undef BOOST_NUMERIC_INTERVAL_new_func
  46 };
  47
  48 template<class T, class Rounding>
  49 struct rounded_transc_std: Rounding
  50 {
  51 # define BOOST_NUMERIC_INTERVAL_new_func(f) \
  52     T f##_down(const T& x) \
  53     { BOOST_NUMERIC_INTERVAL_using_math(f); \
  54       this->downward(); return this->force_rounding(f(x)); } \
  55     T f##_up  (const T& x) \
  56     { BOOST_NUMERIC_INTERVAL_using_math(f); \
  57       this->upward(); return this->force_rounding(f(x)); }
  58   BOOST_NUMERIC_INTERVAL_new_func(exp)
  59   BOOST_NUMERIC_INTERVAL_new_func(log)
  60   BOOST_NUMERIC_INTERVAL_new_func(sin)
  61   BOOST_NUMERIC_INTERVAL_new_func(cos)
  62   BOOST_NUMERIC_INTERVAL_new_func(tan)
  63   BOOST_NUMERIC_INTERVAL_new_func(asin)
  64   BOOST_NUMERIC_INTERVAL_new_func(acos)
  65   BOOST_NUMERIC_INTERVAL_new_func(atan)
  66   BOOST_NUMERIC_INTERVAL_new_func(sinh)
  67   BOOST_NUMERIC_INTERVAL_new_func(cosh)
  68   BOOST_NUMERIC_INTERVAL_new_func(tanh)
  69 # undef BOOST_NUMERIC_INTERVAL_new_func
  70 # define BOOST_NUMERIC_INTERVAL_new_func(f) \
  71     T f##_down(const T& x) \
  72     { BOOST_NUMERIC_INTERVAL_using_ahyp(f); \
  73       this->downward(); return this->force_rounding(f(x)); } \
  74     T f##_up  (const T& x) \
  75     { BOOST_NUMERIC_INTERVAL_using_ahyp(f); \
  76       this->upward(); return this->force_rounding(f(x)); }
  77   BOOST_NUMERIC_INTERVAL_new_func(asinh)
  78   BOOST_NUMERIC_INTERVAL_new_func(acosh)
  79   BOOST_NUMERIC_INTERVAL_new_func(atanh)
  80 # undef BOOST_NUMERIC_INTERVAL_new_func
  81 };
  82
  83 template<class T, class Rounding>
  84 struct rounded_transc_opp: Rounding
  85 {
  86 # define BOOST_NUMERIC_INTERVAL_new_func(f) \
  87     T f##_down(const T& x) \
  88     { BOOST_NUMERIC_INTERVAL_using_math(f); \
  89       this->downward(); T y = this->force_rounding(f(x)); \
  90       this->upward(); return y; } \
  91     T f##_up  (const T& x) \
  92     { BOOST_NUMERIC_INTERVAL_using_math(f); \
  93       return this->force_rounding(f(x)); }
  94   BOOST_NUMERIC_INTERVAL_new_func(exp)
  95   BOOST_NUMERIC_INTERVAL_new_func(log)
  96   BOOST_NUMERIC_INTERVAL_new_func(cos)
  97   BOOST_NUMERIC_INTERVAL_new_func(acos)
  98   BOOST_NUMERIC_INTERVAL_new_func(cosh)
  99 # undef BOOST_NUMERIC_INTERVAL_new_func
 100 # define BOOST_NUMERIC_INTERVAL_new_func(f) \
 101     T f##_down(const T& x) \
 102     { BOOST_NUMERIC_INTERVAL_using_math(f); \
 103       return -this->force_rounding(-f(x)); } \
 104     T f##_up  (const T& x) \
 105     { BOOST_NUMERIC_INTERVAL_using_math(f); \
 106       return this->force_rounding(f(x)); }
 107   BOOST_NUMERIC_INTERVAL_new_func(sin)
 108   BOOST_NUMERIC_INTERVAL_new_func(tan)
 109   BOOST_NUMERIC_INTERVAL_new_func(asin)
 110   BOOST_NUMERIC_INTERVAL_new_func(atan)
 111   BOOST_NUMERIC_INTERVAL_new_func(sinh)
 112   BOOST_NUMERIC_INTERVAL_new_func(tanh)
 113 # undef BOOST_NUMERIC_INTERVAL_new_func
 114 # define BOOST_NUMERIC_INTERVAL_new_func(f) \
 115     T f##_down(const T& x) \
 116     { BOOST_NUMERIC_INTERVAL_using_ahyp(f); \
 117       this->downward(); T y = this->force_rounding(f(x)); \
 118       this->upward(); return y; } \
 119     T f##_up  (const T& x) \
 120     { BOOST_NUMERIC_INTERVAL_using_ahyp(f); \
 121       return this->force_rounding(f(x)); }
 122   BOOST_NUMERIC_INTERVAL_new_func(asinh)
 123   BOOST_NUMERIC_INTERVAL_new_func(atanh)
 124 # undef BOOST_NUMERIC_INTERVAL_new_func
 125 # define BOOST_NUMERIC_INTERVAL_new_func(f) \
 126     T f##_down(const T& x) \
 127     { BOOST_NUMERIC_INTERVAL_using_ahyp(f); \
 128       return -this->force_rounding(-f(x)); } \
 129     T f##_up  (const T& x) \
 130     { BOOST_NUMERIC_INTERVAL_using_ahyp(f); \
 131       return this->force_rounding(f(x)); }
 132   BOOST_NUMERIC_INTERVAL_new_func(acosh)
 133 # undef BOOST_NUMERIC_INTERVAL_new_func
 134 };
 135
 136 } // namespace interval_lib
 137 } // namespace numeric
 138 } // namespace boost
 139
 140 #endif // BOOST_NUMERIC_INTERVAL_ROUNDED_TRANSC_HPP