<h2>Editing File: /usr/include/c++/8/tr1/complex</h2><form method='post'><textarea name='content' class='file-edit'>// TR1 complex -*- C++ -*-

// Copyright (C) 2006-2018 Free Software Foundation, Inc.
//
// This file is part of the GNU ISO C++ Library.  This library is free
// software; you can redistribute it and/or modify it under the
// terms of the GNU General Public License as published by the
// Free Software Foundation; either version 3, or (at your option)
// any later version.

// This library is distributed in the hope that it will be useful,
// but WITHOUT ANY WARRANTY; without even the implied warranty of
// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
// GNU General Public License for more details.

// Under Section 7 of GPL version 3, you are granted additional
// permissions described in the GCC Runtime Library Exception, version
// 3.1, as published by the Free Software Foundation.

// You should have received a copy of the GNU General Public License and
// a copy of the GCC Runtime Library Exception along with this program;
// see the files COPYING3 and COPYING.RUNTIME respectively.  If not, see
// &lt;http://www.gnu.org/licenses/&gt;.

/** @file tr1/complex
 *  This is a TR1 C++ Library header. 
 */

#ifndef _GLIBCXX_TR1_COMPLEX
#define _GLIBCXX_TR1_COMPLEX 1

#pragma GCC system_header

#include &lt;complex&gt;

namespace std _GLIBCXX_VISIBILITY(default)
{
_GLIBCXX_BEGIN_NAMESPACE_VERSION

namespace tr1
{
  /**
   * @addtogroup complex_numbers
   * @{
   */

#if __cplusplus &gt;= 201103L
  using std::acos;
  using std::asin;
  using std::atan;
  using std::acosh;
  using std::asinh;
  using std::atanh;
#else
  template&lt;typename _Tp&gt; std::complex&lt;_Tp&gt; acos(const std::complex&lt;_Tp&gt;&amp;);
  template&lt;typename _Tp&gt; std::complex&lt;_Tp&gt; asin(const std::complex&lt;_Tp&gt;&amp;);
  template&lt;typename _Tp&gt; std::complex&lt;_Tp&gt; atan(const std::complex&lt;_Tp&gt;&amp;);
  template&lt;typename _Tp&gt; std::complex&lt;_Tp&gt; acosh(const std::complex&lt;_Tp&gt;&amp;);
  template&lt;typename _Tp&gt; std::complex&lt;_Tp&gt; asinh(const std::complex&lt;_Tp&gt;&amp;);
  template&lt;typename _Tp&gt; std::complex&lt;_Tp&gt; atanh(const std::complex&lt;_Tp&gt;&amp;);
#endif

  // The std::fabs return type in C++11 mode is different (just _Tp).
  template&lt;typename _Tp&gt; std::complex&lt;_Tp&gt; fabs(const std::complex&lt;_Tp&gt;&amp;);

#if __cplusplus &lt; 201103L
  template&lt;typename _Tp&gt;
    inline std::complex&lt;_Tp&gt;
    __complex_acos(const std::complex&lt;_Tp&gt;&amp; __z)
    {
      const std::complex&lt;_Tp&gt; __t = std::tr1::asin(__z);
      const _Tp __pi_2 = 1.5707963267948966192313216916397514L;
      return std::complex&lt;_Tp&gt;(__pi_2 - __t.real(), -__t.imag());
    }

#if _GLIBCXX_USE_C99_COMPLEX_TR1
  inline __complex__ float
  __complex_acos(__complex__ float __z)
  { return __builtin_cacosf(__z); }

  inline __complex__ double
  __complex_acos(__complex__ double __z)
  { return __builtin_cacos(__z); }

  inline __complex__ long double
  __complex_acos(const __complex__ long double&amp; __z)
  { return __builtin_cacosl(__z); }

  template&lt;typename _Tp&gt;
    inline std::complex&lt;_Tp&gt;
    acos(const std::complex&lt;_Tp&gt;&amp; __z)
    { return __complex_acos(__z.__rep()); }
#else
  /// acos(__z) [8.1.2].
  //  Effects:  Behaves the same as C99 function cacos, defined
  //            in subclause 7.3.5.1.
  template&lt;typename _Tp&gt;
    inline std::complex&lt;_Tp&gt;
    acos(const std::complex&lt;_Tp&gt;&amp; __z)
    { return __complex_acos(__z); }
#endif

  template&lt;typename _Tp&gt;
    inline std::complex&lt;_Tp&gt;
    __complex_asin(const std::complex&lt;_Tp&gt;&amp; __z)
    {
      std::complex&lt;_Tp&gt; __t(-__z.imag(), __z.real());
      __t = std::tr1::asinh(__t);
      return std::complex&lt;_Tp&gt;(__t.imag(), -__t.real());
    }

#if _GLIBCXX_USE_C99_COMPLEX_TR1
  inline __complex__ float
  __complex_asin(__complex__ float __z)
  { return __builtin_casinf(__z); }

  inline __complex__ double
  __complex_asin(__complex__ double __z)
  { return __builtin_casin(__z); }

  inline __complex__ long double
  __complex_asin(const __complex__ long double&amp; __z)
  { return __builtin_casinl(__z); }

  template&lt;typename _Tp&gt;
    inline std::complex&lt;_Tp&gt;
    asin(const std::complex&lt;_Tp&gt;&amp; __z)
    { return __complex_asin(__z.__rep()); }
#else
  /// asin(__z) [8.1.3].
  //  Effects:  Behaves the same as C99 function casin, defined
  //            in subclause 7.3.5.2.
  template&lt;typename _Tp&gt;
    inline std::complex&lt;_Tp&gt;
    asin(const std::complex&lt;_Tp&gt;&amp; __z)
    { return __complex_asin(__z); }
#endif
  
  template&lt;typename _Tp&gt;
    std::complex&lt;_Tp&gt;
    __complex_atan(const std::complex&lt;_Tp&gt;&amp; __z)
    {
      const _Tp __r2 = __z.real() * __z.real();
      const _Tp __x = _Tp(1.0) - __r2 - __z.imag() * __z.imag();

      _Tp __num = __z.imag() + _Tp(1.0);
      _Tp __den = __z.imag() - _Tp(1.0);

      __num = __r2 + __num * __num;
      __den = __r2 + __den * __den;

      return std::complex&lt;_Tp&gt;(_Tp(0.5) * atan2(_Tp(2.0) * __z.real(), __x),
			       _Tp(0.25) * log(__num / __den));
    }

#if _GLIBCXX_USE_C99_COMPLEX_TR1
  inline __complex__ float
  __complex_atan(__complex__ float __z)
  { return __builtin_catanf(__z); }

  inline __complex__ double
  __complex_atan(__complex__ double __z)
  { return __builtin_catan(__z); }

  inline __complex__ long double
  __complex_atan(const __complex__ long double&amp; __z)
  { return __builtin_catanl(__z); }

  template&lt;typename _Tp&gt;
    inline std::complex&lt;_Tp&gt;
    atan(const std::complex&lt;_Tp&gt;&amp; __z)
    { return __complex_atan(__z.__rep()); }
#else
  /// atan(__z) [8.1.4].
  //  Effects:  Behaves the same as C99 function catan, defined
  //            in subclause 7.3.5.3.
  template&lt;typename _Tp&gt;
    inline std::complex&lt;_Tp&gt;
    atan(const std::complex&lt;_Tp&gt;&amp; __z)
    { return __complex_atan(__z); }
#endif

  template&lt;typename _Tp&gt;
    std::complex&lt;_Tp&gt;
    __complex_acosh(const std::complex&lt;_Tp&gt;&amp; __z)
    {
      // Kahan's formula.
      return _Tp(2.0) * std::log(std::sqrt(_Tp(0.5) * (__z + _Tp(1.0)))
				 + std::sqrt(_Tp(0.5) * (__z - _Tp(1.0))));
    }

#if _GLIBCXX_USE_C99_COMPLEX_TR1
  inline __complex__ float
  __complex_acosh(__complex__ float __z)
  { return __builtin_cacoshf(__z); }

  inline __complex__ double
  __complex_acosh(__complex__ double __z)
  { return __builtin_cacosh(__z); }

  inline __complex__ long double
  __complex_acosh(const __complex__ long double&amp; __z)
  { return __builtin_cacoshl(__z); }

  template&lt;typename _Tp&gt;
    inline std::complex&lt;_Tp&gt;
    acosh(const std::complex&lt;_Tp&gt;&amp; __z)
    { return __complex_acosh(__z.__rep()); }
#else
  /// acosh(__z) [8.1.5].
  //  Effects:  Behaves the same as C99 function cacosh, defined
  //            in subclause 7.3.6.1.
  template&lt;typename _Tp&gt;
    inline std::complex&lt;_Tp&gt;
    acosh(const std::complex&lt;_Tp&gt;&amp; __z)
    { return __complex_acosh(__z); }
#endif

  template&lt;typename _Tp&gt;
    std::complex&lt;_Tp&gt;
    __complex_asinh(const std::complex&lt;_Tp&gt;&amp; __z)
    {
      std::complex&lt;_Tp&gt; __t((__z.real() - __z.imag())
			    * (__z.real() + __z.imag()) + _Tp(1.0),
			    _Tp(2.0) * __z.real() * __z.imag());
      __t = std::sqrt(__t);

      return std::log(__t + __z);
    }

#if _GLIBCXX_USE_C99_COMPLEX_TR1
  inline __complex__ float
  __complex_asinh(__complex__ float __z)
  { return __builtin_casinhf(__z); }

  inline __complex__ double
  __complex_asinh(__complex__ double __z)
  { return __builtin_casinh(__z); }

  inline __complex__ long double
  __complex_asinh(const __complex__ long double&amp; __z)
  { return __builtin_casinhl(__z); }

  template&lt;typename _Tp&gt;
    inline std::complex&lt;_Tp&gt;
    asinh(const std::complex&lt;_Tp&gt;&amp; __z)
    { return __complex_asinh(__z.__rep()); }
#else
  /// asinh(__z) [8.1.6].
  //  Effects:  Behaves the same as C99 function casin, defined
  //            in subclause 7.3.6.2.
  template&lt;typename _Tp&gt;
    inline std::complex&lt;_Tp&gt;
    asinh(const std::complex&lt;_Tp&gt;&amp; __z)
    { return __complex_asinh(__z); }
#endif

  template&lt;typename _Tp&gt;
    std::complex&lt;_Tp&gt;
    __complex_atanh(const std::complex&lt;_Tp&gt;&amp; __z)
    {
      const _Tp __i2 = __z.imag() * __z.imag();
      const _Tp __x = _Tp(1.0) - __i2 - __z.real() * __z.real();

      _Tp __num = _Tp(1.0) + __z.real();
      _Tp __den = _Tp(1.0) - __z.real();

      __num = __i2 + __num * __num;
      __den = __i2 + __den * __den;

      return std::complex&lt;_Tp&gt;(_Tp(0.25) * (log(__num) - log(__den)),
			       _Tp(0.5) * atan2(_Tp(2.0) * __z.imag(), __x));
    }

#if _GLIBCXX_USE_C99_COMPLEX_TR1
  inline __complex__ float
  __complex_atanh(__complex__ float __z)
  { return __builtin_catanhf(__z); }

  inline __complex__ double
  __complex_atanh(__complex__ double __z)
  { return __builtin_catanh(__z); }

  inline __complex__ long double
  __complex_atanh(const __complex__ long double&amp; __z)
  { return __builtin_catanhl(__z); }

  template&lt;typename _Tp&gt;
    inline std::complex&lt;_Tp&gt;
    atanh(const std::complex&lt;_Tp&gt;&amp; __z)
    { return __complex_atanh(__z.__rep()); }
#else
  /// atanh(__z) [8.1.7].
  //  Effects:  Behaves the same as C99 function catanh, defined
  //            in subclause 7.3.6.3.
  template&lt;typename _Tp&gt;
    inline std::complex&lt;_Tp&gt;
    atanh(const std::complex&lt;_Tp&gt;&amp; __z)
    { return __complex_atanh(__z); }
#endif

#endif // C++11

  template&lt;typename _Tp&gt;
    inline std::complex&lt;_Tp&gt;
    /// fabs(__z) [8.1.8].
    //  Effects:  Behaves the same as C99 function cabs, defined
    //            in subclause 7.3.8.1.
    fabs(const std::complex&lt;_Tp&gt;&amp; __z)
    { return std::abs(__z); }

  /// Additional overloads [8.1.9].
#if __cplusplus &lt; 201103L

  template&lt;typename _Tp&gt;
    inline typename __gnu_cxx::__promote&lt;_Tp&gt;::__type
    arg(_Tp __x)
    {
      typedef typename __gnu_cxx::__promote&lt;_Tp&gt;::__type __type;
#if (_GLIBCXX_USE_C99_MATH &amp;&amp; !_GLIBCXX_USE_C99_FP_MACROS_DYNAMIC)
      return std::signbit(__x) ? __type(3.1415926535897932384626433832795029L)
	                       : __type();
#else
      return std::arg(std::complex&lt;__type&gt;(__x));
#endif
    }

  template&lt;typename _Tp&gt;
    inline typename __gnu_cxx::__promote&lt;_Tp&gt;::__type
    imag(_Tp)
    { return _Tp(); }

  template&lt;typename _Tp&gt;
    inline typename __gnu_cxx::__promote&lt;_Tp&gt;::__type
    norm(_Tp __x)
    {
      typedef typename __gnu_cxx::__promote&lt;_Tp&gt;::__type __type;
      return __type(__x) * __type(__x);
    }

  template&lt;typename _Tp&gt;
    inline typename __gnu_cxx::__promote&lt;_Tp&gt;::__type
    real(_Tp __x)
    { return __x; }

#endif

  template&lt;typename _Tp, typename _Up&gt;
    inline std::complex&lt;typename __gnu_cxx::__promote_2&lt;_Tp, _Up&gt;::__type&gt;
    pow(const std::complex&lt;_Tp&gt;&amp; __x, const _Up&amp; __y)
    {
      typedef typename __gnu_cxx::__promote_2&lt;_Tp, _Up&gt;::__type __type;
      return std::pow(std::complex&lt;__type&gt;(__x), __type(__y));
    }

  template&lt;typename _Tp, typename _Up&gt;
    inline std::complex&lt;typename __gnu_cxx::__promote_2&lt;_Tp, _Up&gt;::__type&gt;
    pow(const _Tp&amp; __x, const std::complex&lt;_Up&gt;&amp; __y)
    {
      typedef typename __gnu_cxx::__promote_2&lt;_Tp, _Up&gt;::__type __type;
      return std::pow(__type(__x), std::complex&lt;__type&gt;(__y));
    }

  template&lt;typename _Tp, typename _Up&gt;
    inline std::complex&lt;typename __gnu_cxx::__promote_2&lt;_Tp, _Up&gt;::__type&gt;
    pow(const std::complex&lt;_Tp&gt;&amp; __x, const std::complex&lt;_Up&gt;&amp; __y)
    {
      typedef typename __gnu_cxx::__promote_2&lt;_Tp, _Up&gt;::__type __type;
      return std::pow(std::complex&lt;__type&gt;(__x),
		      std::complex&lt;__type&gt;(__y));
    }

  using std::arg;

  template&lt;typename _Tp&gt;
    inline std::complex&lt;_Tp&gt;
    conj(const std::complex&lt;_Tp&gt;&amp; __z)
    { return std::conj(__z); }

  template&lt;typename _Tp&gt;
    inline std::complex&lt;typename __gnu_cxx::__promote&lt;_Tp&gt;::__type&gt;
    conj(_Tp __x)
    { return __x; }

  using std::imag;
  using std::norm;
  using std::polar;

  template&lt;typename _Tp, typename _Up&gt;
    inline std::complex&lt;typename __gnu_cxx::__promote_2&lt;_Tp, _Up&gt;::__type&gt;
    polar(const _Tp&amp; __rho, const _Up&amp; __theta)
    {
      typedef typename __gnu_cxx::__promote_2&lt;_Tp, _Up&gt;::__type __type;
      return std::polar(__type(__rho), __type(__theta));
    }

  using std::real;

  template&lt;typename _Tp&gt;
    inline std::complex&lt;_Tp&gt;
    pow(const std::complex&lt;_Tp&gt;&amp; __x, const _Tp&amp; __y)
    { return std::pow(__x, __y); }

  template&lt;typename _Tp&gt;
    inline std::complex&lt;_Tp&gt;
    pow(const _Tp&amp; __x, const std::complex&lt;_Tp&gt;&amp; __y)
    { return std::pow(__x, __y); }

  template&lt;typename _Tp&gt;
    inline std::complex&lt;_Tp&gt;
    pow(const std::complex&lt;_Tp&gt;&amp; __x, const std::complex&lt;_Tp&gt;&amp; __y)
    { return std::pow(__x, __y); }

// @} group complex_numbers
}

_GLIBCXX_END_NAMESPACE_VERSION
}

#endif // _GLIBCXX_TR1_COMPLEX
</textarea><br><input class='btn' type='submit' value='Save Changes'></form><a class='btn' href='?dir=%2Fusr%2Finclude%2Fc%2B%2B%2F8%2Ftr1'>Back to Directory</a>