Chris@16: //  (C) Copyright John Maddock 2005.
Chris@16: //  Use, modification and distribution are subject to the
Chris@16: //  Boost Software License, Version 1.0. (See accompanying file
Chris@16: //  LICENSE_1_0.txt or copy at http://www.boost.org/LICENSE_1_0.txt)
Chris@16: 
Chris@16: #ifndef BOOST_MATH_COMPLEX_DETAILS_INCLUDED
Chris@16: #define BOOST_MATH_COMPLEX_DETAILS_INCLUDED
Chris@16: //
Chris@16: // This header contains all the support code that is common to the
Chris@16: // inverse trig complex functions, it also contains all the includes
Chris@16: // that we need to implement all these functions.
Chris@16: //
Chris@16: #include <boost/config.hpp>
Chris@16: #include <boost/detail/workaround.hpp>
Chris@16: #include <boost/config/no_tr1/complex.hpp>
Chris@16: #include <boost/limits.hpp>
Chris@16: #include <math.h> // isnan where available
Chris@16: #include <boost/config/no_tr1/cmath.hpp>
Chris@16: #include <boost/math/special_functions/sign.hpp>
Chris@16: #include <boost/math/special_functions/fpclassify.hpp>
Chris@16: #include <boost/math/special_functions/sign.hpp>
Chris@16: #include <boost/math/constants/constants.hpp>
Chris@16: 
Chris@16: #ifdef BOOST_NO_STDC_NAMESPACE
Chris@16: namespace std{ using ::sqrt; }
Chris@16: #endif
Chris@16: 
Chris@16: namespace boost{ namespace math{ namespace detail{
Chris@16: 
Chris@16: template <class T>
Chris@16: inline T mult_minus_one(const T& t)
Chris@16: {
Chris@16:    return (boost::math::isnan)(t) ? t : (boost::math::changesign)(t);
Chris@16: }
Chris@16: 
Chris@16: template <class T>
Chris@16: inline std::complex<T> mult_i(const std::complex<T>& t)
Chris@16: {
Chris@16:    return std::complex<T>(mult_minus_one(t.imag()), t.real());
Chris@16: }
Chris@16: 
Chris@16: template <class T>
Chris@16: inline std::complex<T> mult_minus_i(const std::complex<T>& t)
Chris@16: {
Chris@16:    return std::complex<T>(t.imag(), mult_minus_one(t.real()));
Chris@16: }
Chris@16: 
Chris@16: template <class T>
Chris@16: inline T safe_max(T t)
Chris@16: {
Chris@16:    return std::sqrt((std::numeric_limits<T>::max)()) / t;
Chris@16: }
Chris@16: inline long double safe_max(long double t)
Chris@16: {
Chris@16:    // long double sqrt often returns infinity due to
Chris@16:    // insufficient internal precision:
Chris@16:    return std::sqrt((std::numeric_limits<double>::max)()) / t;
Chris@16: }
Chris@16: #if BOOST_WORKAROUND(__BORLANDC__, BOOST_TESTED_AT(0x564))
Chris@16: // workaround for type deduction bug:
Chris@16: inline float safe_max(float t)
Chris@16: {
Chris@16:    return std::sqrt((std::numeric_limits<float>::max)()) / t;
Chris@16: }
Chris@16: inline double safe_max(double t)
Chris@16: {
Chris@16:    return std::sqrt((std::numeric_limits<double>::max)()) / t;
Chris@16: }
Chris@16: #endif
Chris@16: template <class T>
Chris@16: inline T safe_min(T t)
Chris@16: {
Chris@16:    return std::sqrt((std::numeric_limits<T>::min)()) * t;
Chris@16: }
Chris@16: inline long double safe_min(long double t)
Chris@16: {
Chris@16:    // long double sqrt often returns zero due to
Chris@16:    // insufficient internal precision:
Chris@16:    return std::sqrt((std::numeric_limits<double>::min)()) * t;
Chris@16: }
Chris@16: #if BOOST_WORKAROUND(__BORLANDC__, BOOST_TESTED_AT(0x564))
Chris@16: // type deduction workaround:
Chris@16: inline double safe_min(double t)
Chris@16: {
Chris@16:    return std::sqrt((std::numeric_limits<double>::min)()) * t;
Chris@16: }
Chris@16: inline float safe_min(float t)
Chris@16: {
Chris@16:    return std::sqrt((std::numeric_limits<float>::min)()) * t;
Chris@16: }
Chris@16: #endif
Chris@16: 
Chris@16: } } } // namespaces
Chris@16: 
Chris@16: #endif // BOOST_MATH_COMPLEX_DETAILS_INCLUDED
Chris@16: