ceph/src/boost/boost/math/complex/details.hpp

   1 //  (C) Copyright John Maddock 2005.
   2 //  Use, modification and distribution are subject to the
   3 //  Boost Software License, Version 1.0. (See accompanying file
   4 //  LICENSE_1_0.txt or copy at http://www.boost.org/LICENSE_1_0.txt)
   5
   6 #ifndef BOOST_MATH_COMPLEX_DETAILS_INCLUDED
   7 #define BOOST_MATH_COMPLEX_DETAILS_INCLUDED
   8 //
   9 // This header contains all the support code that is common to the
  10 // inverse trig complex functions, it also contains all the includes
  11 // that we need to implement all these functions.
  12 //
  13
  14 #include <cmath>
  15 #include <complex>
  16 #include <limits>
  17 #include <boost/math/special_functions/sign.hpp>
  18 #include <boost/math/special_functions/fpclassify.hpp>
  19 #include <boost/math/special_functions/sign.hpp>
  20 #include <boost/math/constants/constants.hpp>
  21
  22 namespace boost{ namespace math{ namespace detail{
  23
  24 template <class T>
  25 inline T mult_minus_one(const T& t)
  26 {
  27    return (boost::math::isnan)(t) ? t : (boost::math::changesign)(t);
  28 }
  29
  30 template <class T>
  31 inline std::complex<T> mult_i(const std::complex<T>& t)
  32 {
  33    return std::complex<T>(mult_minus_one(t.imag()), t.real());
  34 }
  35
  36 template <class T>
  37 inline std::complex<T> mult_minus_i(const std::complex<T>& t)
  38 {
  39    return std::complex<T>(t.imag(), mult_minus_one(t.real()));
  40 }
  41
  42 template <class T>
  43 inline T safe_max(T t)
  44 {
  45    return std::sqrt((std::numeric_limits<T>::max)()) / t;
  46 }
  47 inline long double safe_max(long double t)
  48 {
  49    // long double sqrt often returns infinity due to
  50    // insufficient internal precision:
  51    return std::sqrt((std::numeric_limits<double>::max)()) / t;
  52 }
  53
  54 template <class T>
  55 inline T safe_min(T t)
  56 {
  57    return std::sqrt((std::numeric_limits<T>::min)()) * t;
  58 }
  59 inline long double safe_min(long double t)
  60 {
  61    // long double sqrt often returns zero due to
  62    // insufficient internal precision:
  63    return std::sqrt((std::numeric_limits<double>::min)()) * t;
  64 }
  65
  66 } } } // namespaces
  67
  68 #endif // BOOST_MATH_COMPLEX_DETAILS_INCLUDED
  69