include/std/std_complex.h

200     complex<_Tp>::complex(const complex<_Up>& __z)
201     : _M_real(__z.real()), _M_imag(__z.imag()) { }
253     complex<_Tp>::operator=(const complex<_Up>& __z)
255       _M_real = __z.real();
256       _M_imag = __z.imag();
264     complex<_Tp>::operator+=(const complex<_Up>& __z)
266       _M_real += __z.real();
267       _M_imag += __z.imag();
275     complex<_Tp>::operator-=(const complex<_Up>& __z)
277       _M_real -= __z.real();
278       _M_imag -= __z.imag();
287     complex<_Tp>::operator*=(const complex<_Up>& __z)
289       const _Tp __r = _M_real * __z.real() - _M_imag * __z.imag();
290       _M_imag = _M_real * __z.imag() + _M_imag * __z.real();
300     complex<_Tp>::operator/=(const complex<_Up>& __z)
302       const _Tp __r =  _M_real * __z.real() + _M_imag * __z.imag();
303       const _Tp __n = std::norm(__z);
304       _M_imag = (_M_imag * __z.real() - _M_real * __z.imag()) / __n;
531     real(complex<_Tp>& __z)
532     { return __z.real(); }
536     real(const complex<_Tp>& __z)
537     { return __z.real(); }
541     imag(complex<_Tp>& __z)
542     { return __z.imag(); }
546     imag(const complex<_Tp>& __z)
547     { return __z.imag(); }
549   // 26.2.7/3 abs(__z):  Returns the magnitude of __z.
552     __complex_abs(const complex<_Tp>& __z)
554       _Tp __x = __z.real();
555       _Tp __y = __z.imag();
566   __complex_abs(__complex__ float __z) { return __builtin_cabsf(__z); }
569   __complex_abs(__complex__ double __z) { return __builtin_cabs(__z); }
572   __complex_abs(const __complex__ long double& __z)
573   { return __builtin_cabsl(__z); }
577     abs(const complex<_Tp>& __z) { return __complex_abs(__z.__rep()); }
581     abs(const complex<_Tp>& __z) { return __complex_abs(__z); }
585   // 26.2.7/4: arg(__z): Returns the phase angle of __z.
588     __complex_arg(const complex<_Tp>& __z)
589     { return  atan2(__z.imag(), __z.real()); }
593   __complex_arg(__complex__ float __z) { return __builtin_cargf(__z); }
596   __complex_arg(__complex__ double __z) { return __builtin_carg(__z); }
599   __complex_arg(const __complex__ long double& __z)
600   { return __builtin_cargl(__z); }
604     arg(const complex<_Tp>& __z) { return __complex_arg(__z.__rep()); }
608     arg(const complex<_Tp>& __z) { return __complex_arg(__z); }
611   // 26.2.7/5: norm(__z) returns the squared magintude of __z.
620         static inline _Tp _S_do_it(const complex<_Tp>& __z)
622           const _Tp __x = __z.real();
623           const _Tp __y = __z.imag();
632         static inline _Tp _S_do_it(const complex<_Tp>& __z)
634           _Tp __res = std::abs(__z);
641     norm(const complex<_Tp>& __z)
644 	&& !_GLIBCXX_FAST_MATH>::_S_do_it(__z);
654     conj(const complex<_Tp>& __z)
655     { return complex<_Tp>(__z.real(), -__z.imag()); }
659   // 26.2.8/1 cos(__z):  Returns the cosine of __z.
662     __complex_cos(const complex<_Tp>& __z)
664       const _Tp __x = __z.real();
665       const _Tp __y = __z.imag();
671   __complex_cos(__complex__ float __z) { return __builtin_ccosf(__z); }
674   __complex_cos(__complex__ double __z) { return __builtin_ccos(__z); }
677   __complex_cos(const __complex__ long double& __z)
678   { return __builtin_ccosl(__z); }
682     cos(const complex<_Tp>& __z) { return __complex_cos(__z.__rep()); }
686     cos(const complex<_Tp>& __z) { return __complex_cos(__z); }
689   // 26.2.8/2 cosh(__z): Returns the hyperbolic cosine of __z.
692     __complex_cosh(const complex<_Tp>& __z)
694       const _Tp __x = __z.real();
695       const _Tp __y = __z.imag();
701   __complex_cosh(__complex__ float __z) { return __builtin_ccoshf(__z); }
704   __complex_cosh(__complex__ double __z) { return __builtin_ccosh(__z); }
707   __complex_cosh(const __complex__ long double& __z)
708   { return __builtin_ccoshl(__z); }
712     cosh(const complex<_Tp>& __z) { return __complex_cosh(__z.__rep()); }
716     cosh(const complex<_Tp>& __z) { return __complex_cosh(__z); }
719   // 26.2.8/3 exp(__z): Returns the complex base e exponential of x
722     __complex_exp(const complex<_Tp>& __z)
723     { return std::polar(exp(__z.real()), __z.imag()); }
727   __complex_exp(__complex__ float __z) { return __builtin_cexpf(__z); }
730   __complex_exp(__complex__ double __z) { return __builtin_cexp(__z); }
733   __complex_exp(const __complex__ long double& __z)
734   { return __builtin_cexpl(__z); }
738     exp(const complex<_Tp>& __z) { return __complex_exp(__z.__rep()); }
742     exp(const complex<_Tp>& __z) { return __complex_exp(__z); }
745   // 26.2.8/5 log(__z): Reurns the natural complex logaritm of __z.
749     __complex_log(const complex<_Tp>& __z)
750     { return complex<_Tp>(log(std::abs(__z)), std::arg(__z)); }
754   __complex_log(__complex__ float __z) { return __builtin_clogf(__z); }
757   __complex_log(__complex__ double __z) { return __builtin_clog(__z); }
760   __complex_log(const __complex__ long double& __z)
761   { return __builtin_clogl(__z); }
765     log(const complex<_Tp>& __z) { return __complex_log(__z.__rep()); }
769     log(const complex<_Tp>& __z) { return __complex_log(__z); }
774     log10(const complex<_Tp>& __z)
775     { return std::log(__z) / log(_Tp(10.0)); }
777   // 26.2.8/10 sin(__z): Returns the sine of __z.
780     __complex_sin(const complex<_Tp>& __z)
782       const _Tp __x = __z.real();
783       const _Tp __y = __z.imag();
789   __complex_sin(__complex__ float __z) { return __builtin_csinf(__z); }
792   __complex_sin(__complex__ double __z) { return __builtin_csin(__z); }
795   __complex_sin(const __complex__ long double& __z)
796   { return __builtin_csinl(__z); }
800     sin(const complex<_Tp>& __z) { return __complex_sin(__z.__rep()); }
804     sin(const complex<_Tp>& __z) { return __complex_sin(__z); }
807   // 26.2.8/11 sinh(__z): Returns the hyperbolic sine of __z.
810     __complex_sinh(const complex<_Tp>& __z)
812       const _Tp __x = __z.real();
813       const _Tp  __y = __z.imag();
819   __complex_sinh(__complex__ float __z) { return __builtin_csinhf(__z); }
822   __complex_sinh(__complex__ double __z) { return __builtin_csinh(__z); }
825   __complex_sinh(const __complex__ long double& __z)
826   { return __builtin_csinhl(__z); }
830     sinh(const complex<_Tp>& __z) { return __complex_sinh(__z.__rep()); }
834     sinh(const complex<_Tp>& __z) { return __complex_sinh(__z); }
837   // 26.2.8/13 sqrt(__z): Returns the complex square root of __z.
841     __complex_sqrt(const complex<_Tp>& __z)
843       _Tp __x = __z.real();
844       _Tp __y = __z.imag();
853           _Tp __t = sqrt(2 * (std::abs(__z) + abs(__x)));
863   __complex_sqrt(__complex__ float __z) { return __builtin_csqrtf(__z); }
866   __complex_sqrt(__complex__ double __z) { return __builtin_csqrt(__z); }
869   __complex_sqrt(const __complex__ long double& __z)
870   { return __builtin_csqrtl(__z); }
874     sqrt(const complex<_Tp>& __z) { return __complex_sqrt(__z.__rep()); }
878     sqrt(const complex<_Tp>& __z) { return __complex_sqrt(__z); }
881   // 26.2.8/14 tan(__z):  Return the complex tangent of __z.
885     __complex_tan(const complex<_Tp>& __z)
886     { return std::sin(__z) / std::cos(__z); }
890   __complex_tan(__complex__ float __z) { return __builtin_ctanf(__z); }
893   __complex_tan(__complex__ double __z) { return __builtin_ctan(__z); }
896   __complex_tan(const __complex__ long double& __z)
897   { return __builtin_ctanl(__z); }
901     tan(const complex<_Tp>& __z) { return __complex_tan(__z.__rep()); }
905     tan(const complex<_Tp>& __z) { return __complex_tan(__z); }
909   // 26.2.8/15 tanh(__z):  Returns the hyperbolic tangent of __z.
913     __complex_tanh(const complex<_Tp>& __z)
914     { return std::sinh(__z) / std::cosh(__z); }
918   __complex_tanh(__complex__ float __z) { return __builtin_ctanhf(__z); }
921   __complex_tanh(__complex__ double __z) { return __builtin_ctanh(__z); }
924   __complex_tanh(const __complex__ long double& __z)
925   { return __builtin_ctanhl(__z); }
929     tanh(const complex<_Tp>& __z) { return __complex_tanh(__z.__rep()); }
933     tanh(const complex<_Tp>& __z) { return __complex_tanh(__z); }
942     pow(const complex<_Tp>& __z, int __n)
943     { return std::__pow_helper(__z, __n); }
1007       complex(_ComplexT __z) : _M_value(__z) { }
1106   complex<float>::operator=(const complex<_Tp>& __z)
1108     __real__ _M_value = __z.real();
1109     __imag__ _M_value = __z.imag();
1115   complex<float>::operator+=(const complex<_Tp>& __z)
1117     __real__ _M_value += __z.real();
1118     __imag__ _M_value += __z.imag();
1124     complex<float>::operator-=(const complex<_Tp>& __z)
1126      __real__ _M_value -= __z.real();
1127      __imag__ _M_value -= __z.imag();
1133     complex<float>::operator*=(const complex<_Tp>& __z)
1136       __real__ __t = __z.real();
1137       __imag__ __t = __z.imag();
1144     complex<float>::operator/=(const complex<_Tp>& __z)
1147       __real__ __t = __z.real();
1148       __imag__ __t = __z.imag();
1161       complex(_ComplexT __z) : _M_value(__z) { }
1259     complex<double>::operator=(const complex<_Tp>& __z)
1261       __real__ _M_value = __z.real();
1262       __imag__ _M_value = __z.imag();
1268     complex<double>::operator+=(const complex<_Tp>& __z)
1270       __real__ _M_value += __z.real();
1271       __imag__ _M_value += __z.imag();
1277     complex<double>::operator-=(const complex<_Tp>& __z)
1279       __real__ _M_value -= __z.real();
1280       __imag__ _M_value -= __z.imag();
1286     complex<double>::operator*=(const complex<_Tp>& __z)
1289       __real__ __t = __z.real();
1290       __imag__ __t = __z.imag();
1297     complex<double>::operator/=(const complex<_Tp>& __z)
1300       __real__ __t = __z.real();
1301       __imag__ __t = __z.imag();
1314       complex(_ComplexT __z) : _M_value(__z) { }
1412     complex<long double>::operator=(const complex<_Tp>& __z)
1414       __real__ _M_value = __z.real();
1415       __imag__ _M_value = __z.imag();
1421     complex<long double>::operator+=(const complex<_Tp>& __z)
1423       __real__ _M_value += __z.real();
1424       __imag__ _M_value += __z.imag();
1430     complex<long double>::operator-=(const complex<_Tp>& __z)
1432       __real__ _M_value -= __z.real();
1433       __imag__ _M_value -= __z.imag();
1439     complex<long double>::operator*=(const complex<_Tp>& __z)
1442       __real__ __t = __z.real();
1443       __imag__ __t = __z.imag();
1450     complex<long double>::operator/=(const complex<_Tp>& __z)
1453       __real__ __t = __z.real();
1454       __imag__ __t = __z.imag();
1464   complex<float>::complex(const complex<double>& __z)
1465   : _M_value(__z.__rep()) { }
1468   complex<float>::complex(const complex<long double>& __z)
1469   : _M_value(__z.__rep()) { }
1472   complex<double>::complex(const complex<float>& __z)
1473   : _M_value(__z.__rep()) { }
1476   complex<double>::complex(const complex<long double>& __z)
1477     : _M_value(__z.__rep()) { }
1480   complex<long double>::complex(const complex<float>& __z)
1481   : _M_value(__z.__rep()) { }
1484   complex<long double>::complex(const complex<double>& __z)
1485   : _M_value(__z.__rep()) { }