Blame view

3rdparty/boost_1_81_0/boost/math/complex/details.hpp 1.87 KB
63e88f80   Hu Chunming   提交三方库
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
  //  (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 <cmath>
  #include <complex>
  #include <limits>
  #include <boost/math/special_functions/sign.hpp>
  #include <boost/math/special_functions/fpclassify.hpp>
  #include <boost/math/special_functions/sign.hpp>
  #include <boost/math/constants/constants.hpp>
  
  namespace boost{ namespace math{ namespace detail{
  
  template <class T>
  inline T mult_minus_one(const T& t)
  {
     return (boost::math::isnan)(t) ? t : (boost::math::changesign)(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;
  }
  
  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;
  }
  
  } } } // namespaces
  
  #endif // BOOST_MATH_COMPLEX_DETAILS_INCLUDED