| // (C) Copyright John Maddock 2005. |
| // Use, modification and distribution are subject to 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_MATH_COMPLEX_DETAILS_INCLUDED |
| #define BOOST_MATH_COMPLEX_DETAILS_INCLUDED |
| // |
| // This header contains all the support code that is common to the |
| // inverse trig complex functions, it also contains all the includes |
| // that we need to implement all these functions. |
| // |
| #include <boost/config.hpp> |
| #include <boost/detail/workaround.hpp> |
| #include <boost/config/no_tr1/complex.hpp> |
| #include <boost/limits.hpp> |
| #include <math.h> // isnan where available |
| #include <boost/config/no_tr1/cmath.hpp> |
| |
| #ifdef BOOST_NO_STDC_NAMESPACE |
| namespace std{ using ::sqrt; } |
| #endif |
| |
| namespace boost{ namespace math{ namespace detail{ |
| |
| template <class T> |
| inline bool test_is_nan(T t) |
| { |
| // Comparisons with Nan's always fail: |
| return std::numeric_limits<T>::has_infinity && (!(t <= std::numeric_limits<T>::infinity()) || !(t >= -std::numeric_limits<T>::infinity())); |
| } |
| #ifdef isnan |
| template<> inline bool test_is_nan<float>(float t) { return isnan(t); } |
| template<> inline bool test_is_nan<double>(double t) { return isnan(t); } |
| template<> inline bool test_is_nan<long double>(long double t) { return isnan(t); } |
| #endif |
| |
| template <class T> |
| inline T mult_minus_one(const T& t) |
| { |
| return test_is_nan(t) ? t : -t; |
| } |
| |
| template <class T> |
| inline std::complex<T> mult_i(const std::complex<T>& t) |
| { |
| return std::complex<T>(mult_minus_one(t.imag()), t.real()); |
| } |
| |
| template <class T> |
| inline std::complex<T> mult_minus_i(const std::complex<T>& t) |
| { |
| return std::complex<T>(t.imag(), mult_minus_one(t.real())); |
| } |
| |
| template <class T> |
| inline T safe_max(T t) |
| { |
| return std::sqrt((std::numeric_limits<T>::max)()) / t; |
| } |
| inline long double safe_max(long double t) |
| { |
| // long double sqrt often returns infinity due to |
| // insufficient internal precision: |
| return std::sqrt((std::numeric_limits<double>::max)()) / t; |
| } |
| #if BOOST_WORKAROUND(__BORLANDC__, BOOST_TESTED_AT(0x564)) |
| // workaround for type deduction bug: |
| inline float safe_max(float t) |
| { |
| return std::sqrt((std::numeric_limits<float>::max)()) / t; |
| } |
| inline double safe_max(double t) |
| { |
| return std::sqrt((std::numeric_limits<double>::max)()) / t; |
| } |
| #endif |
| template <class T> |
| inline T safe_min(T t) |
| { |
| return std::sqrt((std::numeric_limits<T>::min)()) * t; |
| } |
| inline long double safe_min(long double t) |
| { |
| // long double sqrt often returns zero due to |
| // insufficient internal precision: |
| return std::sqrt((std::numeric_limits<double>::min)()) * t; |
| } |
| #if BOOST_WORKAROUND(__BORLANDC__, BOOST_TESTED_AT(0x564)) |
| // type deduction workaround: |
| inline double safe_min(double t) |
| { |
| return std::sqrt((std::numeric_limits<double>::min)()) * t; |
| } |
| inline float safe_min(float t) |
| { |
| return std::sqrt((std::numeric_limits<float>::min)()) * t; |
| } |
| #endif |
| |
| } } } // namespaces |
| |
| #endif // BOOST_MATH_COMPLEX_DETAILS_INCLUDED |
| |