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_TR1_COMPLEX_HPP_INCLUDED
Chris@16: #  define BOOST_TR1_COMPLEX_HPP_INCLUDED
Chris@16: #  include <boost/tr1/detail/config.hpp>
Chris@16: #  include <complex>
Chris@16: 
Chris@16: #ifndef BOOST_HAS_TR1_COMPLEX_INVERSE_TRIG
Chris@16: 
Chris@16: #include <boost/math/complex.hpp>
Chris@16: 
Chris@16: namespace std {
Chris@16: namespace tr1 {
Chris@16: 
Chris@16: using boost::math::acos;
Chris@16: using boost::math::asin;
Chris@16: using boost::math::atan;
Chris@16: using boost::math::acosh;
Chris@16: using boost::math::asinh;
Chris@16: using boost::math::atanh;
Chris@16: using boost::math::fabs;
Chris@16: 
Chris@16: } }
Chris@16: 
Chris@16: #else
Chris@16: 
Chris@16: #  if defined(BOOST_HAS_INCLUDE_NEXT) && !defined(BOOST_TR1_DISABLE_INCLUDE_NEXT)
Chris@16: #     include_next BOOST_TR1_HEADER(complex)
Chris@16: #  else
Chris@16: #     include <boost/tr1/detail/config_all.hpp>
Chris@16: #     include BOOST_TR1_STD_HEADER(BOOST_TR1_PATH(complex))
Chris@16: #  endif
Chris@16: 
Chris@16: #endif
Chris@16: 
Chris@16: #ifndef BOOST_HAS_TR1_COMPLEX_OVERLOADS
Chris@16: 
Chris@16: #include <boost/tr1/detail/math_overloads.hpp>
Chris@16: #include <boost/assert.hpp>
Chris@16: #include <boost/detail/workaround.hpp>
Chris@16: #include <boost/config/no_tr1/cmath.hpp>
Chris@16: 
Chris@16: namespace std{ 
Chris@16: 
Chris@16: #ifdef BOOST_NO_STDC_NAMESPACE
Chris@16:    using :: atan2;
Chris@16: #endif
Chris@16: 
Chris@16: #ifndef BOOST_NO_FUNCTION_TEMPLATE_ORDERING
Chris@16: template <class T>
Chris@16: inline BOOST_TR1_MATH_RETURN(double) arg(const T& t)
Chris@16: {
Chris@16:    return ::std::atan2(0.0, static_cast<double>(t));
Chris@16: }
Chris@16: #else
Chris@16: inline double arg(const double& t)
Chris@16: {
Chris@16:    return ::std::atan2(0.0, t);
Chris@16: }
Chris@16: #endif
Chris@16: inline long double arg(const long double& t)
Chris@16: {
Chris@16:    return ::std::atan2(0.0L, static_cast<long double>(t));
Chris@16: }
Chris@16: inline float arg(const float& t)
Chris@16: {
Chris@16:    return ::std::atan2(0.0F, static_cast<float>(t));
Chris@16: }
Chris@16: 
Chris@16: #ifndef BOOST_NO_FUNCTION_TEMPLATE_ORDERING
Chris@16: template <class T>
Chris@16: inline BOOST_TR1_MATH_RETURN(double) norm(const T& t)
Chris@16: {
Chris@16:    double r = static_cast<double>(t);
Chris@16:    return r*r;
Chris@16: }
Chris@16: #else
Chris@16: inline double norm(const double& t)
Chris@16: {
Chris@16:    return t*t;
Chris@16: }
Chris@16: #endif
Chris@16: inline long double norm(const long double& t)
Chris@16: {
Chris@16:    long double l = t;
Chris@16:    return l*l;
Chris@16: }
Chris@16: inline float norm(const float& t)
Chris@16: {
Chris@16:    float f = t;
Chris@16:    return f*f;
Chris@16: }
Chris@16: 
Chris@16: #ifndef BOOST_NO_FUNCTION_TEMPLATE_ORDERING
Chris@16: template <class T>
Chris@16: inline BOOST_TR1_MATH_RETURN(std::complex<double>) conj(const T& t)
Chris@16: {
Chris@16:    return ::std::conj(std::complex<double>(static_cast<double>(t)));
Chris@16: }
Chris@16: #else
Chris@16: inline std::complex<double> conj(const double& t)
Chris@16: {
Chris@16:    return ::std::conj(std::complex<double>(t));
Chris@16: }
Chris@16: #endif
Chris@16: inline std::complex<long double> conj(const long double& t)
Chris@16: {
Chris@16:    return ::std::conj(std::complex<long double>(t));
Chris@16: }
Chris@16: inline std::complex<float> conj(const float& t)
Chris@16: {
Chris@16:    std::complex<float> ct(t);
Chris@16:    ct = ::std::conj(ct);
Chris@16:    return ct;
Chris@16: }
Chris@16: 
Chris@101: #if !BOOST_WORKAROUND(__BORLANDC__, <=0x570)
Chris@16: inline complex<double> polar(const char& rho, const char& theta = 0)
Chris@16: { return ::std::polar(static_cast<double>(rho), static_cast<double>(theta)); }
Chris@16: inline complex<double> polar(const unsigned char& rho, const unsigned char& theta = 0)
Chris@16: { return ::std::polar(static_cast<double>(rho), static_cast<double>(theta)); }
Chris@16: inline complex<double> polar(const signed char& rho, const signed char& theta = 0)
Chris@16: { return ::std::polar(static_cast<double>(rho), static_cast<double>(theta)); }
Chris@16: inline complex<double> polar(const short& rho, const short& theta = 0)
Chris@16: { return ::std::polar(static_cast<double>(rho), static_cast<double>(theta)); }
Chris@16: inline complex<double> polar(const unsigned short& rho, const unsigned short& theta = 0)
Chris@16: { return ::std::polar(static_cast<double>(rho), static_cast<double>(theta)); }
Chris@16: inline complex<double> polar(const int& rho, const int& theta = 0)
Chris@16: { return ::std::polar(static_cast<double>(rho), static_cast<double>(theta)); }
Chris@16: inline complex<double> polar(const unsigned int& rho, const unsigned int& theta = 0)
Chris@16: { return ::std::polar(static_cast<double>(rho), static_cast<double>(theta)); }
Chris@16: inline complex<double> polar(const long& rho, const long& theta = 0)
Chris@16: { return ::std::polar(static_cast<double>(rho), static_cast<double>(theta)); }
Chris@16: inline complex<double> polar(const unsigned long& rho, const unsigned long& theta = 0)
Chris@16: { return ::std::polar(static_cast<double>(rho), static_cast<double>(theta)); }
Chris@16: #ifdef BOOST_HAS_LONG_LONG
Chris@16: inline complex<double> polar(const long long& rho, const long long& theta = 0)
Chris@16: { return ::std::polar(static_cast<double>(rho), static_cast<double>(theta)); }
Chris@16: inline complex<double> polar(const unsigned long long& rho, const unsigned long long& theta = 0)
Chris@16: { return ::std::polar(static_cast<double>(rho), static_cast<double>(theta)); }
Chris@16: #elif defined(BOOST_HAS_MS_INT64)
Chris@16: inline complex<double> polar(const __int64& rho, const __int64& theta = 0)
Chris@16: { return ::std::polar(static_cast<double>(rho), static_cast<double>(theta)); }
Chris@16: inline complex<double> polar(const unsigned __int64& rho, const unsigned __int64& theta = 0)
Chris@16: { return ::std::polar(static_cast<double>(rho), static_cast<double>(theta)); }
Chris@16: #endif
Chris@16: 
Chris@16: template<class T, class U> 
Chris@16: inline complex<typename boost::tr1_detail::promote_to_real<T, U>::type> 
Chris@16:    polar(const T& rho, const U& theta)
Chris@16: {
Chris@16:    typedef typename boost::tr1_detail::promote_to_real<T, U>::type real_type;
Chris@16:    return std::polar(static_cast<real_type>(rho), static_cast<real_type>(theta));
Chris@16: }
Chris@16: #endif
Chris@16: 
Chris@16: #ifndef BOOST_NO_FUNCTION_TEMPLATE_ORDERING
Chris@16: template <class T>
Chris@16: inline BOOST_TR1_MATH_RETURN(double) imag(const T& )
Chris@16: {
Chris@16:    return 0;
Chris@16: }
Chris@16: #else
Chris@16: inline double imag(const double& )
Chris@16: {
Chris@16:    return 0;
Chris@16: }
Chris@16: #endif
Chris@16: inline long double imag(const long double& )
Chris@16: {
Chris@16:    return 0;
Chris@16: }
Chris@16: inline float imag(const float& )
Chris@16: {
Chris@16:    return 0;
Chris@16: }
Chris@16: 
Chris@16: #ifndef BOOST_NO_FUNCTION_TEMPLATE_ORDERING
Chris@16: template <class T>
Chris@16: inline BOOST_TR1_MATH_RETURN(double) real(const T& t)
Chris@16: {
Chris@16:    return static_cast<double>(t);
Chris@16: }
Chris@16: #else
Chris@16: inline double real(const double& t)
Chris@16: {
Chris@16:    return t;
Chris@16: }
Chris@16: #endif
Chris@16: inline long double real(const long double& t)
Chris@16: {
Chris@16:    return t;
Chris@16: }
Chris@16: inline float real(const float& t)
Chris@16: {
Chris@16:    return t;
Chris@16: }
Chris@16: 
Chris@16: template<class T, class U>
Chris@16: inline complex<typename boost::tr1_detail::largest_real<T, U>::type>
Chris@16:    pow(const complex<T>& x, const complex<U>& y)
Chris@16: {
Chris@16:    typedef complex<typename boost::tr1_detail::largest_real<T, U>::type> result_type;
Chris@16:    typedef typename boost::mpl::if_<boost::is_same<result_type, complex<T> >, result_type const&, result_type>::type cast1_type;
Chris@16:    typedef typename boost::mpl::if_<boost::is_same<result_type, complex<U> >, result_type const&, result_type>::type cast2_type;
Chris@16:    cast1_type x1(x);
Chris@16:    cast2_type y1(y);
Chris@16:    return std::pow(x1, y1);
Chris@16: }
Chris@16: template<class T, class U> 
Chris@16: inline complex<typename boost::tr1_detail::promote_to_real<T, U>::type>
Chris@16:    pow (const complex<T>& x, const U& y)
Chris@16: {
Chris@16:    typedef typename boost::tr1_detail::promote_to_real<T, U>::type real_type;
Chris@16:    typedef complex<typename boost::tr1_detail::promote_to_real<T, U>::type> result_type;
Chris@16:    typedef typename boost::mpl::if_<boost::is_same<result_type, complex<T> >, result_type const&, result_type>::type cast1_type;
Chris@16:    real_type r = y;
Chris@16:    cast1_type x1(x);
Chris@16:    std::complex<real_type> y1(r);
Chris@16:    return std::pow(x1, y1);
Chris@16: }
Chris@16: 
Chris@16: template<class T, class U> 
Chris@16: inline complex<typename boost::tr1_detail::promote_to_real<T, U>::type>
Chris@16:    pow (const T& x, const complex<U>& y)
Chris@16: {
Chris@16:    typedef typename boost::tr1_detail::promote_to_real<T, U>::type real_type;
Chris@16:    typedef complex<typename boost::tr1_detail::promote_to_real<T, U>::type> result_type;
Chris@16:    typedef typename boost::mpl::if_<boost::is_same<result_type, complex<U> >, result_type const&, result_type>::type cast_type;
Chris@16:    real_type r = x;
Chris@16:    std::complex<real_type> x1(r);
Chris@16:    cast_type y1(y);
Chris@16:    return std::pow(x1, y1);
Chris@16: }
Chris@16: 
Chris@16: }
Chris@16: 
Chris@16: #endif
Chris@16: 
Chris@16: #endif
Chris@16: