Attempt to represent the S^2->S^3 header reorganisation as a series of "hg rename" operations
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)
6 #ifndef BOOST_MATH_COMPLEX_DETAILS_INCLUDED
7 #define BOOST_MATH_COMPLEX_DETAILS_INCLUDED
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.
13 #include <boost/detail/workaround.hpp>
14 #include <boost/config.hpp>
15 #include <boost/config/no_tr1/complex.hpp>
16 #include <boost/limits.hpp>
17 #include <math.h> // isnan where available
20 #ifdef BOOST_NO_STDC_NAMESPACE
21 namespace std{ using ::sqrt; }
24 namespace boost{ namespace math{ namespace detail{
27 inline bool test_is_nan(T t)
29 // Comparisons with Nan's always fail:
30 return std::numeric_limits<T>::has_infinity && (!(t <= std::numeric_limits<T>::infinity()) || !(t >= -std::numeric_limits<T>::infinity()));
33 template<> inline bool test_is_nan<float>(float t) { return isnan(t); }
34 template<> inline bool test_is_nan<double>(double t) { return isnan(t); }
35 template<> inline bool test_is_nan<long double>(long double t) { return isnan(t); }
39 inline T mult_minus_one(const T& t)
41 return test_is_nan(t) ? t : -t;
45 inline std::complex<T> mult_i(const std::complex<T>& t)
47 return std::complex<T>(mult_minus_one(t.imag()), t.real());
51 inline std::complex<T> mult_minus_i(const std::complex<T>& t)
53 return std::complex<T>(t.imag(), mult_minus_one(t.real()));
57 inline T safe_max(T t)
59 return std::sqrt((std::numeric_limits<T>::max)()) / t;
61 inline long double safe_max(long double t)
63 // long double sqrt often returns infinity due to
64 // insufficient internal precision:
65 return std::sqrt((std::numeric_limits<double>::max)()) / t;
67 #if BOOST_WORKAROUND(__BORLANDC__, BOOST_TESTED_AT(0x564))
68 // workaround for type deduction bug:
69 inline float safe_max(float t)
71 return std::sqrt((std::numeric_limits<float>::max)()) / t;
73 inline double safe_max(double t)
75 return std::sqrt((std::numeric_limits<double>::max)()) / t;
79 inline T safe_min(T t)
81 return std::sqrt((std::numeric_limits<T>::min)()) * t;
83 inline long double safe_min(long double t)
85 // long double sqrt often returns zero due to
86 // insufficient internal precision:
87 return std::sqrt((std::numeric_limits<double>::min)()) * t;
89 #if BOOST_WORKAROUND(__BORLANDC__, BOOST_TESTED_AT(0x564))
90 // type deduction workaround:
91 inline double safe_min(double t)
93 return std::sqrt((std::numeric_limits<double>::min)()) * t;
95 inline float safe_min(float t)
97 return std::sqrt((std::numeric_limits<float>::min)()) * t;
103 #endif // BOOST_MATH_COMPLEX_DETAILS_INCLUDED