Chris@16: /* Boost interval/transc.hpp template implementation file
Chris@16:  *
Chris@16:  * Copyright 2000 Jens Maurer
Chris@16:  * Copyright 2002 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_TRANSC_HPP
Chris@16: #define BOOST_NUMERIC_INTERVAL_TRANSC_HPP
Chris@16: 
Chris@16: #include <boost/config.hpp>
Chris@16: #include <boost/numeric/interval/detail/interval_prototype.hpp>
Chris@16: #include <boost/numeric/interval/detail/bugs.hpp>
Chris@16: #include <boost/numeric/interval/detail/test_input.hpp>
Chris@16: #include <boost/numeric/interval/rounding.hpp>
Chris@16: #include <boost/numeric/interval/constants.hpp>
Chris@16: #include <boost/numeric/interval/arith.hpp>
Chris@16: #include <boost/numeric/interval/arith2.hpp>
Chris@16: #include <algorithm>
Chris@16: 
Chris@16: namespace boost {
Chris@16: namespace numeric {
Chris@16: 
Chris@16: template<class T, class Policies> inline
Chris@16: interval<T, Policies> exp(const interval<T, Policies>& x)
Chris@16: {
Chris@16:   typedef interval<T, Policies> I;
Chris@16:   if (interval_lib::detail::test_input(x))
Chris@16:     return I::empty();
Chris@16:   typename Policies::rounding rnd;
Chris@16:   return I(rnd.exp_down(x.lower()), rnd.exp_up(x.upper()), true);
Chris@16: }
Chris@16: 
Chris@16: template<class T, class Policies> inline
Chris@16: interval<T, Policies> log(const interval<T, Policies>& x)
Chris@16: {
Chris@16:   typedef interval<T, Policies> I;
Chris@16:   if (interval_lib::detail::test_input(x) ||
Chris@16:       !interval_lib::user::is_pos(x.upper()))
Chris@16:     return I::empty();
Chris@16:   typename Policies::rounding rnd;
Chris@16:   typedef typename Policies::checking checking;
Chris@16:   T l = !interval_lib::user::is_pos(x.lower())
Chris@16:              ? checking::neg_inf() : rnd.log_down(x.lower());
Chris@16:   return I(l, rnd.log_up(x.upper()), true);
Chris@16: }
Chris@16: 
Chris@16: template<class T, class Policies> inline
Chris@16: interval<T, Policies> cos(const interval<T, Policies>& x)
Chris@16: {
Chris@16:   if (interval_lib::detail::test_input(x))
Chris@16:     return interval<T, Policies>::empty();
Chris@16:   typename Policies::rounding rnd;
Chris@16:   typedef interval<T, Policies> I;
Chris@16:   typedef typename interval_lib::unprotect<I>::type R;
Chris@16: 
Chris@16:   // get lower bound within [0, pi]
Chris@16:   const R pi2 = interval_lib::pi_twice<R>();
Chris@16:   R tmp = fmod((const R&)x, pi2);
Chris@16:   if (width(tmp) >= pi2.lower())
Chris@16:     return I(static_cast<T>(-1), static_cast<T>(1), true); // we are covering a full period
Chris@16:   if (tmp.lower() >= interval_lib::constants::pi_upper<T>())
Chris@16:     return -cos(tmp - interval_lib::pi<R>());
Chris@16:   T l = tmp.lower();
Chris@16:   T u = tmp.upper();
Chris@16: 
Chris@16:   BOOST_USING_STD_MIN();
Chris@16:   // separate into monotone subintervals
Chris@16:   if (u <= interval_lib::constants::pi_lower<T>())
Chris@16:     return I(rnd.cos_down(u), rnd.cos_up(l), true);
Chris@16:   else if (u <= pi2.lower())
Chris@16:     return I(static_cast<T>(-1), rnd.cos_up(min BOOST_PREVENT_MACRO_SUBSTITUTION(rnd.sub_down(pi2.lower(), u), l)), true);
Chris@16:   else
Chris@16:     return I(static_cast<T>(-1), static_cast<T>(1), true);
Chris@16: }
Chris@16: 
Chris@16: template<class T, class Policies> inline
Chris@16: interval<T, Policies> sin(const interval<T, Policies>& x)
Chris@16: {
Chris@16:   typedef interval<T, Policies> I;
Chris@16:   if (interval_lib::detail::test_input(x))
Chris@16:     return I::empty();
Chris@16:   typename Policies::rounding rnd;
Chris@16:   typedef typename interval_lib::unprotect<I>::type R;
Chris@16:   I r = cos((const R&)x - interval_lib::pi_half<R>());
Chris@16:   (void)&rnd;
Chris@16:   return r;
Chris@16: }
Chris@16: 
Chris@16: template<class T, class Policies> inline
Chris@16: interval<T, Policies> tan(const interval<T, Policies>& x)
Chris@16: {
Chris@16:   typedef interval<T, Policies> I;
Chris@16:   if (interval_lib::detail::test_input(x))
Chris@16:     return I::empty();
Chris@16:   typename Policies::rounding rnd;
Chris@16:   typedef typename interval_lib::unprotect<I>::type R;
Chris@16: 
Chris@16:   // get lower bound within [-pi/2, pi/2]
Chris@16:   const R pi = interval_lib::pi<R>();
Chris@16:   R tmp = fmod((const R&)x, pi);
Chris@16:   const T pi_half_d = interval_lib::constants::pi_half_lower<T>();
Chris@16:   if (tmp.lower() >= pi_half_d)
Chris@16:     tmp -= pi;
Chris@16:   if (tmp.lower() <= -pi_half_d || tmp.upper() >= pi_half_d)
Chris@16:     return I::whole();
Chris@16:   return I(rnd.tan_down(tmp.lower()), rnd.tan_up(tmp.upper()), true);
Chris@16: }
Chris@16: 
Chris@16: template<class T, class Policies> inline
Chris@16: interval<T, Policies> asin(const interval<T, Policies>& x)
Chris@16: {
Chris@16:   typedef interval<T, Policies> I;
Chris@16:   if (interval_lib::detail::test_input(x)
Chris@16:      || x.upper() < static_cast<T>(-1) || x.lower() > static_cast<T>(1))
Chris@16:     return I::empty();
Chris@16:   typename Policies::rounding rnd;
Chris@16:   T l = (x.lower() <= static_cast<T>(-1))
Chris@16:              ? -interval_lib::constants::pi_half_upper<T>()
Chris@16:              : rnd.asin_down(x.lower());
Chris@16:   T u = (x.upper() >= static_cast<T>(1) )
Chris@16:              ?  interval_lib::constants::pi_half_upper<T>()
Chris@16:              : rnd.asin_up  (x.upper());
Chris@16:   return I(l, u, true);
Chris@16: }
Chris@16: 
Chris@16: template<class T, class Policies> inline
Chris@16: interval<T, Policies> acos(const interval<T, Policies>& x)
Chris@16: {
Chris@16:   typedef interval<T, Policies> I;
Chris@16:   if (interval_lib::detail::test_input(x)
Chris@16:      || x.upper() < static_cast<T>(-1) || x.lower() > static_cast<T>(1))
Chris@16:     return I::empty();
Chris@16:   typename Policies::rounding rnd;
Chris@16:   T l = (x.upper() >= static_cast<T>(1) )
Chris@16:           ? static_cast<T>(0)
Chris@16:           : rnd.acos_down(x.upper());
Chris@16:   T u = (x.lower() <= static_cast<T>(-1))
Chris@16:           ? interval_lib::constants::pi_upper<T>()
Chris@16:           : rnd.acos_up  (x.lower());
Chris@16:   return I(l, u, true);
Chris@16: }
Chris@16: 
Chris@16: template<class T, class Policies> inline
Chris@16: interval<T, Policies> atan(const interval<T, Policies>& x)
Chris@16: {
Chris@16:   typedef interval<T, Policies> I;
Chris@16:   if (interval_lib::detail::test_input(x))
Chris@16:     return I::empty();
Chris@16:   typename Policies::rounding rnd;
Chris@16:   return I(rnd.atan_down(x.lower()), rnd.atan_up(x.upper()), true);
Chris@16: }
Chris@16: 
Chris@16: template<class T, class Policies> inline
Chris@16: interval<T, Policies> sinh(const interval<T, Policies>& x)
Chris@16: {
Chris@16:   typedef interval<T, Policies> I;
Chris@16:   if (interval_lib::detail::test_input(x))
Chris@16:     return I::empty();
Chris@16:   typename Policies::rounding rnd;
Chris@16:   return I(rnd.sinh_down(x.lower()), rnd.sinh_up(x.upper()), true);
Chris@16: }
Chris@16: 
Chris@16: template<class T, class Policies> inline
Chris@16: interval<T, Policies> cosh(const interval<T, Policies>& x)
Chris@16: {
Chris@16:   typedef interval<T, Policies> I;
Chris@16:   if (interval_lib::detail::test_input(x))
Chris@16:     return I::empty();
Chris@16:   typename Policies::rounding rnd;
Chris@16:   if (interval_lib::user::is_neg(x.upper()))
Chris@16:     return I(rnd.cosh_down(x.upper()), rnd.cosh_up(x.lower()), true);
Chris@16:   else if (!interval_lib::user::is_neg(x.lower()))
Chris@16:     return I(rnd.cosh_down(x.lower()), rnd.cosh_up(x.upper()), true);
Chris@16:   else
Chris@16:     return I(static_cast<T>(1), rnd.cosh_up(-x.lower() > x.upper() ? x.lower() : x.upper()), true);
Chris@16: }
Chris@16: 
Chris@16: template<class T, class Policies> inline
Chris@16: interval<T, Policies> tanh(const interval<T, Policies>& x)
Chris@16: {
Chris@16:   typedef interval<T, Policies> I;
Chris@16:   if (interval_lib::detail::test_input(x))
Chris@16:     return I::empty();
Chris@16:   typename Policies::rounding rnd;
Chris@16:   return I(rnd.tanh_down(x.lower()), rnd.tanh_up(x.upper()), true);
Chris@16: }
Chris@16: 
Chris@16: template<class T, class Policies> inline
Chris@16: interval<T, Policies> asinh(const interval<T, Policies>& x)
Chris@16: {
Chris@16:   typedef interval<T, Policies> I;
Chris@16:   if (interval_lib::detail::test_input(x))
Chris@16:     return I::empty();
Chris@16:   typename Policies::rounding rnd;
Chris@16:   return I(rnd.asinh_down(x.lower()), rnd.asinh_up(x.upper()), true);
Chris@16: }
Chris@16: 
Chris@16: template<class T, class Policies> inline
Chris@16: interval<T, Policies> acosh(const interval<T, Policies>& x)
Chris@16: {
Chris@16:   typedef interval<T, Policies> I;
Chris@16:   if (interval_lib::detail::test_input(x) || x.upper() < static_cast<T>(1))
Chris@16:     return I::empty();
Chris@16:   typename Policies::rounding rnd;
Chris@16:   T l = x.lower() <= static_cast<T>(1) ? static_cast<T>(0) : rnd.acosh_down(x.lower());
Chris@16:   return I(l, rnd.acosh_up(x.upper()), true);
Chris@16: }
Chris@16: 
Chris@16: template<class T, class Policies> inline
Chris@16: interval<T, Policies> atanh(const interval<T, Policies>& x)
Chris@16: {
Chris@16:   typedef interval<T, Policies> I;
Chris@16:   if (interval_lib::detail::test_input(x)
Chris@16:       || x.upper() < static_cast<T>(-1) || x.lower() > static_cast<T>(1))
Chris@16:     return I::empty();
Chris@16:   typename Policies::rounding rnd;
Chris@16:   typedef typename Policies::checking checking;
Chris@16:   T l = (x.lower() <= static_cast<T>(-1))
Chris@16:              ? checking::neg_inf() : rnd.atanh_down(x.lower());
Chris@16:   T u = (x.upper() >= static_cast<T>(1) )
Chris@16:              ? checking::pos_inf() : rnd.atanh_up  (x.upper());
Chris@16:   return I(l, u, true);
Chris@16: }
Chris@16: 
Chris@16: } // namespace numeric
Chris@16: } // namespace boost
Chris@16: 
Chris@16: #endif // BOOST_NUMERIC_INTERVAL_TRANSC_HPP