38 #ifndef _GLIBCXX_COMPLEX
39 #define _GLIBCXX_COMPLEX 1
41 #pragma GCC system_header
49 _GLIBCXX_BEGIN_NAMESPACE(std)
61 template<> class complex<
float>;
62 template<> class complex<
double>;
63 template<> class complex<
long double>;
66 template<typename _Tp> _Tp
abs(const complex<_Tp>&);
68 template<typename _Tp> _Tp
arg(const complex<_Tp>&);
70 template<typename _Tp> _Tp
norm(const complex<_Tp>&);
73 template<typename _Tp> complex<_Tp>
conj(const complex<_Tp>&);
75 template<typename _Tp> complex<_Tp>
polar(const _Tp&, const _Tp& = 0);
79 template<typename _Tp> complex<_Tp>
cos(const complex<_Tp>&);
81 template<typename _Tp> complex<_Tp>
cosh(const complex<_Tp>&);
83 template<typename _Tp> complex<_Tp>
exp(const complex<_Tp>&);
85 template<typename _Tp> complex<_Tp>
log(const complex<_Tp>&);
87 template<typename _Tp> complex<_Tp>
log10(const complex<_Tp>&);
88 #ifndef __GXX_EXPERIMENTAL_CXX0X__
121 template<
typename _Tp>
129 complex(
const _Tp& __r = _Tp(),
const _Tp& __i = _Tp())
130 : _M_real(__r), _M_imag(__i) { }
135 template<
typename _Up>
137 : _M_real(__z.real()), _M_imag(__z.imag()) { }
139 #ifdef __GXX_EXPERIMENTAL_CXX0X__
153 const _Tp& real()
const
161 const _Tp& imag()
const
174 complex<_Tp>& operator=(
const _Tp&);
179 operator+=(
const _Tp& __t)
188 operator-=(
const _Tp& __t)
203 template<
typename _Up>
206 template<
typename _Up>
209 template<
typename _Up>
212 template<
typename _Up>
215 template<
typename _Up>
218 const complex& __rep()
const
226 template<
typename _Tp>
236 template<
typename _Tp>
246 template<
typename _Tp>
255 template<
typename _Tp>
256 template<
typename _Up>
260 _M_real = __z.real();
261 _M_imag = __z.imag();
266 template<
typename _Tp>
267 template<
typename _Up>
271 _M_real += __z.real();
272 _M_imag += __z.imag();
277 template<
typename _Tp>
278 template<
typename _Up>
282 _M_real -= __z.real();
283 _M_imag -= __z.imag();
289 template<
typename _Tp>
290 template<
typename _Up>
294 const _Tp __r = _M_real * __z.real() - _M_imag * __z.imag();
295 _M_imag = _M_real * __z.imag() + _M_imag * __z.real();
302 template<
typename _Tp>
303 template<
typename _Up>
307 const _Tp __r = _M_real * __z.real() + _M_imag * __z.imag();
309 _M_imag = (_M_imag * __z.real() - _M_real * __z.imag()) / __n;
317 template<
typename _Tp>
326 template<
typename _Tp>
335 template<
typename _Tp>
347 template<
typename _Tp>
356 template<
typename _Tp>
365 template<
typename _Tp>
377 template<
typename _Tp>
386 template<
typename _Tp>
395 template<
typename _Tp>
407 template<
typename _Tp>
416 template<
typename _Tp>
425 template<
typename _Tp>
436 template<
typename _Tp>
442 template<
typename _Tp>
449 template<
typename _Tp>
452 {
return __x.real() == __y.real() && __x.imag() == __y.imag(); }
454 template<
typename _Tp>
457 {
return __x.real() == __y && __x.imag() == _Tp(); }
459 template<
typename _Tp>
462 {
return __x == __y.real() && _Tp() == __y.imag(); }
467 template<
typename _Tp>
470 {
return __x.real() != __y.real() || __x.imag() != __y.imag(); }
472 template<
typename _Tp>
475 {
return __x.real() != __y || __x.imag() != _Tp(); }
477 template<
typename _Tp>
480 {
return __x != __y.real() || _Tp() != __y.imag(); }
484 template<
typename _Tp,
typename _CharT,
class _Traits>
485 basic_istream<_CharT, _Traits>&
493 __is >> __re_x >> __ch;
496 __is >> __im_x >> __ch;
502 else if (__ch ==
')')
517 template<
typename _Tp,
typename _CharT,
class _Traits>
518 basic_ostream<_CharT, _Traits>&
519 operator<<(basic_ostream<_CharT, _Traits>& __os,
const complex<_Tp>& __x)
522 __s.
flags(__os.flags());
523 __s.
imbue(__os.getloc());
525 __s <<
'(' << __x.real() <<
',' << __x.imag() <<
')';
526 return __os << __s.
str();
530 #ifdef __GXX_EXPERIMENTAL_CXX0X__
531 template<
typename _Tp>
533 real(
const complex<_Tp>& __z)
534 {
return __z.real(); }
536 template<
typename _Tp>
538 imag(
const complex<_Tp>& __z)
539 {
return __z.imag(); }
541 template<
typename _Tp>
543 real(complex<_Tp>& __z)
544 {
return __z.real(); }
546 template<
typename _Tp>
548 real(
const complex<_Tp>& __z)
549 {
return __z.real(); }
551 template<
typename _Tp>
553 imag(complex<_Tp>& __z)
554 {
return __z.imag(); }
556 template<
typename _Tp>
558 imag(
const complex<_Tp>& __z)
559 {
return __z.imag(); }
563 template<
typename _Tp>
565 __complex_abs(
const complex<_Tp>& __z)
567 _Tp __x = __z.real();
568 _Tp __y = __z.imag();
574 return __s *
sqrt(__x * __x + __y * __y);
577 #if _GLIBCXX_USE_C99_COMPLEX
579 __complex_abs(__complex__
float __z) {
return __builtin_cabsf(__z); }
582 __complex_abs(__complex__
double __z) {
return __builtin_cabs(__z); }
585 __complex_abs(
const __complex__
long double& __z)
586 {
return __builtin_cabsl(__z); }
588 template<
typename _Tp>
590 abs(
const complex<_Tp>& __z) {
return __complex_abs(__z.__rep()); }
592 template<
typename _Tp>
599 template<
typename _Tp>
601 __complex_arg(
const complex<_Tp>& __z)
602 {
return atan2(__z.imag(), __z.real()); }
604 #if _GLIBCXX_USE_C99_COMPLEX
606 __complex_arg(__complex__
float __z) {
return __builtin_cargf(__z); }
609 __complex_arg(__complex__
double __z) {
return __builtin_carg(__z); }
612 __complex_arg(
const __complex__
long double& __z)
613 {
return __builtin_cargl(__z); }
615 template<
typename _Tp>
617 arg(
const complex<_Tp>& __z) {
return __complex_arg(__z.__rep()); }
619 template<
typename _Tp>
632 template<
typename _Tp>
633 static inline _Tp _S_do_it(
const complex<_Tp>& __z)
635 const _Tp __x = __z.real();
636 const _Tp __y = __z.imag();
637 return __x * __x + __y * __y;
642 struct _Norm_helper<true>
644 template<
typename _Tp>
645 static inline _Tp _S_do_it(
const complex<_Tp>& __z)
648 return __res * __res;
652 template<
typename _Tp>
656 return _Norm_helper<__is_floating<_Tp>::__value
657 && !_GLIBCXX_FAST_MATH>::_S_do_it(__z);
660 template<
typename _Tp>
662 polar(
const _Tp& __rho,
const _Tp& __theta)
665 template<
typename _Tp>
673 template<
typename _Tp>
675 __complex_cos(
const complex<_Tp>& __z)
677 const _Tp __x = __z.real();
678 const _Tp __y = __z.imag();
682 #if _GLIBCXX_USE_C99_COMPLEX
683 inline __complex__
float
684 __complex_cos(__complex__
float __z) {
return __builtin_ccosf(__z); }
686 inline __complex__
double
687 __complex_cos(__complex__
double __z) {
return __builtin_ccos(__z); }
689 inline __complex__
long double
690 __complex_cos(
const __complex__
long double& __z)
691 {
return __builtin_ccosl(__z); }
693 template<
typename _Tp>
695 cos(
const complex<_Tp>& __z) {
return __complex_cos(__z.__rep()); }
697 template<
typename _Tp>
703 template<
typename _Tp>
705 __complex_cosh(
const complex<_Tp>& __z)
707 const _Tp __x = __z.real();
708 const _Tp __y = __z.imag();
712 #if _GLIBCXX_USE_C99_COMPLEX
713 inline __complex__
float
714 __complex_cosh(__complex__
float __z) {
return __builtin_ccoshf(__z); }
716 inline __complex__
double
717 __complex_cosh(__complex__
double __z) {
return __builtin_ccosh(__z); }
719 inline __complex__
long double
720 __complex_cosh(
const __complex__
long double& __z)
721 {
return __builtin_ccoshl(__z); }
723 template<
typename _Tp>
725 cosh(
const complex<_Tp>& __z) {
return __complex_cosh(__z.__rep()); }
727 template<
typename _Tp>
733 template<
typename _Tp>
735 __complex_exp(
const complex<_Tp>& __z)
738 #if _GLIBCXX_USE_C99_COMPLEX
739 inline __complex__
float
740 __complex_exp(__complex__
float __z) {
return __builtin_cexpf(__z); }
742 inline __complex__
double
743 __complex_exp(__complex__
double __z) {
return __builtin_cexp(__z); }
745 inline __complex__
long double
746 __complex_exp(
const __complex__
long double& __z)
747 {
return __builtin_cexpl(__z); }
749 template<
typename _Tp>
751 exp(
const complex<_Tp>& __z) {
return __complex_exp(__z.__rep()); }
753 template<
typename _Tp>
760 template<
typename _Tp>
762 __complex_log(
const complex<_Tp>& __z)
765 #if _GLIBCXX_USE_C99_COMPLEX
766 inline __complex__
float
767 __complex_log(__complex__
float __z) {
return __builtin_clogf(__z); }
769 inline __complex__
double
770 __complex_log(__complex__
double __z) {
return __builtin_clog(__z); }
772 inline __complex__
long double
773 __complex_log(
const __complex__
long double& __z)
774 {
return __builtin_clogl(__z); }
776 template<
typename _Tp>
778 log(
const complex<_Tp>& __z) {
return __complex_log(__z.__rep()); }
780 template<
typename _Tp>
785 template<
typename _Tp>
791 template<
typename _Tp>
793 __complex_sin(
const complex<_Tp>& __z)
795 const _Tp __x = __z.real();
796 const _Tp __y = __z.imag();
800 #if _GLIBCXX_USE_C99_COMPLEX
801 inline __complex__
float
802 __complex_sin(__complex__
float __z) {
return __builtin_csinf(__z); }
804 inline __complex__
double
805 __complex_sin(__complex__
double __z) {
return __builtin_csin(__z); }
807 inline __complex__
long double
808 __complex_sin(
const __complex__
long double& __z)
809 {
return __builtin_csinl(__z); }
811 template<
typename _Tp>
813 sin(
const complex<_Tp>& __z) {
return __complex_sin(__z.__rep()); }
815 template<
typename _Tp>
821 template<
typename _Tp>
823 __complex_sinh(
const complex<_Tp>& __z)
825 const _Tp __x = __z.real();
826 const _Tp __y = __z.imag();
830 #if _GLIBCXX_USE_C99_COMPLEX
831 inline __complex__
float
832 __complex_sinh(__complex__
float __z) {
return __builtin_csinhf(__z); }
834 inline __complex__
double
835 __complex_sinh(__complex__
double __z) {
return __builtin_csinh(__z); }
837 inline __complex__
long double
838 __complex_sinh(
const __complex__
long double& __z)
839 {
return __builtin_csinhl(__z); }
841 template<
typename _Tp>
843 sinh(
const complex<_Tp>& __z) {
return __complex_sinh(__z.__rep()); }
845 template<
typename _Tp>
852 template<
typename _Tp>
854 __complex_sqrt(
const complex<_Tp>& __z)
856 _Tp __x = __z.real();
857 _Tp __y = __z.imag();
862 return complex<_Tp>(__t, __y < _Tp() ? -__t : __t);
869 ? complex<_Tp>(__u, __y / __t)
870 : complex<_Tp>(
abs(__y) / __t, __y < _Tp() ? -__u : __u);
874 #if _GLIBCXX_USE_C99_COMPLEX
875 inline __complex__
float
876 __complex_sqrt(__complex__
float __z) {
return __builtin_csqrtf(__z); }
878 inline __complex__
double
879 __complex_sqrt(__complex__
double __z) {
return __builtin_csqrt(__z); }
881 inline __complex__
long double
882 __complex_sqrt(
const __complex__
long double& __z)
883 {
return __builtin_csqrtl(__z); }
885 template<
typename _Tp>
887 sqrt(
const complex<_Tp>& __z) {
return __complex_sqrt(__z.__rep()); }
889 template<
typename _Tp>
896 template<
typename _Tp>
898 __complex_tan(
const complex<_Tp>& __z)
901 #if _GLIBCXX_USE_C99_COMPLEX
902 inline __complex__
float
903 __complex_tan(__complex__
float __z) {
return __builtin_ctanf(__z); }
905 inline __complex__
double
906 __complex_tan(__complex__
double __z) {
return __builtin_ctan(__z); }
908 inline __complex__
long double
909 __complex_tan(
const __complex__
long double& __z)
910 {
return __builtin_ctanl(__z); }
912 template<
typename _Tp>
914 tan(
const complex<_Tp>& __z) {
return __complex_tan(__z.__rep()); }
916 template<
typename _Tp>
924 template<
typename _Tp>
926 __complex_tanh(
const complex<_Tp>& __z)
929 #if _GLIBCXX_USE_C99_COMPLEX
930 inline __complex__
float
931 __complex_tanh(__complex__
float __z) {
return __builtin_ctanhf(__z); }
933 inline __complex__
double
934 __complex_tanh(__complex__
double __z) {
return __builtin_ctanh(__z); }
936 inline __complex__
long double
937 __complex_tanh(
const __complex__
long double& __z)
938 {
return __builtin_ctanhl(__z); }
940 template<
typename _Tp>
942 tanh(
const complex<_Tp>& __z) {
return __complex_tanh(__z.__rep()); }
944 template<
typename _Tp>
953 #ifndef __GXX_EXPERIMENTAL_CXX0X__
956 template<
typename _Tp>
958 pow(
const complex<_Tp>& __z,
int __n)
959 {
return std::__pow_helper(__z, __n); }
962 template<
typename _Tp>
966 #ifndef _GLIBCXX_USE_C99_COMPLEX
970 if (__x.imag() == _Tp() && __x.real() > _Tp())
971 return pow(__x.real(), __y);
977 template<
typename _Tp>
979 __complex_pow(
const complex<_Tp>& __x,
const complex<_Tp>& __y)
980 {
return __x == _Tp() ? _Tp() : std::
exp(__y * std::
log(__x)); }
982 #if _GLIBCXX_USE_C99_COMPLEX
983 inline __complex__
float
984 __complex_pow(__complex__
float __x, __complex__
float __y)
985 {
return __builtin_cpowf(__x, __y); }
987 inline __complex__
double
988 __complex_pow(__complex__
double __x, __complex__
double __y)
989 {
return __builtin_cpow(__x, __y); }
991 inline __complex__
long double
992 __complex_pow(
const __complex__
long double& __x,
993 const __complex__
long double& __y)
994 {
return __builtin_cpowl(__x, __y); }
996 template<
typename _Tp>
998 pow(
const complex<_Tp>& __x,
const complex<_Tp>& __y)
999 {
return __complex_pow(__x.__rep(), __y.__rep()); }
1001 template<
typename _Tp>
1004 {
return __complex_pow(__x, __y); }
1007 template<
typename _Tp>
1012 __y.imag() *
log(__x))
1019 struct complex<float>
1021 typedef float value_type;
1022 typedef __complex__
float _ComplexT;
1024 complex(_ComplexT __z) : _M_value(__z) { }
1026 complex(
float __r = 0.0f,
float __i = 0.0f)
1028 __real__ _M_value = __r;
1029 __imag__ _M_value = __i;
1032 explicit complex(
const complex<double>&);
1033 explicit complex(
const complex<long double>&);
1035 #ifdef __GXX_EXPERIMENTAL_CXX0X__
1039 {
return __real__ _M_value; }
1042 {
return __imag__ _M_value; }
1045 {
return __real__ _M_value; }
1047 const float& real()
const
1048 {
return __real__ _M_value; }
1051 {
return __imag__ _M_value; }
1053 const float& imag()
const
1054 {
return __imag__ _M_value; }
1059 void real(
float __val)
1060 { __real__ _M_value = __val; }
1062 void imag(
float __val)
1063 { __imag__ _M_value = __val; }
1066 operator=(
float __f)
1068 __real__ _M_value = __f;
1069 __imag__ _M_value = 0.0f;
1074 operator+=(
float __f)
1076 __real__ _M_value += __f;
1081 operator-=(
float __f)
1083 __real__ _M_value -= __f;
1088 operator*=(
float __f)
1095 operator/=(
float __f)
1105 template<
typename _Tp>
1107 operator=(
const complex<_Tp>& __z)
1109 __real__ _M_value = __z.real();
1110 __imag__ _M_value = __z.imag();
1114 template<
typename _Tp>
1116 operator+=(
const complex<_Tp>& __z)
1118 __real__ _M_value += __z.real();
1119 __imag__ _M_value += __z.imag();
1125 operator-=(
const complex<_Tp>& __z)
1127 __real__ _M_value -= __z.real();
1128 __imag__ _M_value -= __z.imag();
1134 operator*=(
const complex<_Tp>& __z)
1137 __real__ __t = __z.real();
1138 __imag__ __t = __z.imag();
1145 operator/=(
const complex<_Tp>& __z)
1148 __real__ __t = __z.real();
1149 __imag__ __t = __z.imag();
1154 const _ComplexT& __rep()
const {
return _M_value; }
1163 struct complex<double>
1165 typedef double value_type;
1166 typedef __complex__
double _ComplexT;
1168 complex(_ComplexT __z) : _M_value(__z) { }
1170 complex(
double __r = 0.0,
double __i = 0.0)
1172 __real__ _M_value = __r;
1173 __imag__ _M_value = __i;
1176 complex(
const complex<float>& __z)
1177 : _M_value(__z.__rep()) { }
1179 explicit complex(
const complex<long double>&);
1181 #ifdef __GXX_EXPERIMENTAL_CXX0X__
1185 {
return __real__ _M_value; }
1188 {
return __imag__ _M_value; }
1191 {
return __real__ _M_value; }
1193 const double& real()
const
1194 {
return __real__ _M_value; }
1197 {
return __imag__ _M_value; }
1199 const double& imag()
const
1200 {
return __imag__ _M_value; }
1205 void real(
double __val)
1206 { __real__ _M_value = __val; }
1208 void imag(
double __val)
1209 { __imag__ _M_value = __val; }
1212 operator=(
double __d)
1214 __real__ _M_value = __d;
1215 __imag__ _M_value = 0.0;
1220 operator+=(
double __d)
1222 __real__ _M_value += __d;
1227 operator-=(
double __d)
1229 __real__ _M_value -= __d;
1234 operator*=(
double __d)
1241 operator/=(
double __d)
1250 template<
typename _Tp>
1252 operator=(
const complex<_Tp>& __z)
1254 __real__ _M_value = __z.real();
1255 __imag__ _M_value = __z.imag();
1259 template<
typename _Tp>
1261 operator+=(
const complex<_Tp>& __z)
1263 __real__ _M_value += __z.real();
1264 __imag__ _M_value += __z.imag();
1268 template<
typename _Tp>
1270 operator-=(
const complex<_Tp>& __z)
1272 __real__ _M_value -= __z.real();
1273 __imag__ _M_value -= __z.imag();
1277 template<
typename _Tp>
1279 operator*=(
const complex<_Tp>& __z)
1282 __real__ __t = __z.real();
1283 __imag__ __t = __z.imag();
1288 template<
typename _Tp>
1290 operator/=(
const complex<_Tp>& __z)
1293 __real__ __t = __z.real();
1294 __imag__ __t = __z.imag();
1299 const _ComplexT& __rep()
const {
return _M_value; }
1308 struct complex<long double>
1310 typedef long double value_type;
1311 typedef __complex__
long double _ComplexT;
1313 complex(_ComplexT __z) : _M_value(__z) { }
1315 complex(
long double __r = 0.0L,
long double __i = 0.0L)
1317 __real__ _M_value = __r;
1318 __imag__ _M_value = __i;
1321 complex(
const complex<float>& __z)
1322 : _M_value(__z.__rep()) { }
1324 complex(
const complex<double>& __z)
1325 : _M_value(__z.__rep()) { }
1327 #ifdef __GXX_EXPERIMENTAL_CXX0X__
1330 long double real()
const
1331 {
return __real__ _M_value; }
1333 long double imag()
const
1334 {
return __imag__ _M_value; }
1337 {
return __real__ _M_value; }
1339 const long double& real()
const
1340 {
return __real__ _M_value; }
1343 {
return __imag__ _M_value; }
1345 const long double& imag()
const
1346 {
return __imag__ _M_value; }
1351 void real(
long double __val)
1352 { __real__ _M_value = __val; }
1354 void imag(
long double __val)
1355 { __imag__ _M_value = __val; }
1357 complex<long double>&
1358 operator=(
long double __r)
1360 __real__ _M_value = __r;
1361 __imag__ _M_value = 0.0L;
1365 complex<long double>&
1366 operator+=(
long double __r)
1368 __real__ _M_value += __r;
1372 complex<long double>&
1373 operator-=(
long double __r)
1375 __real__ _M_value -= __r;
1379 complex<long double>&
1380 operator*=(
long double __r)
1386 complex<long double>&
1387 operator/=(
long double __r)
1396 template<
typename _Tp>
1397 complex<long double>&
1398 operator=(
const complex<_Tp>& __z)
1400 __real__ _M_value = __z.real();
1401 __imag__ _M_value = __z.imag();
1405 template<
typename _Tp>
1406 complex<long double>&
1407 operator+=(
const complex<_Tp>& __z)
1409 __real__ _M_value += __z.real();
1410 __imag__ _M_value += __z.imag();
1414 template<
typename _Tp>
1415 complex<long double>&
1416 operator-=(
const complex<_Tp>& __z)
1418 __real__ _M_value -= __z.real();
1419 __imag__ _M_value -= __z.imag();
1423 template<
typename _Tp>
1424 complex<long double>&
1425 operator*=(
const complex<_Tp>& __z)
1428 __real__ __t = __z.real();
1429 __imag__ __t = __z.imag();
1434 template<
typename _Tp>
1435 complex<long double>&
1436 operator/=(
const complex<_Tp>& __z)
1439 __real__ __t = __z.real();
1440 __imag__ __t = __z.imag();
1445 const _ComplexT& __rep()
const {
return _M_value; }
1454 complex<float>::complex(
const complex<double>& __z)
1455 : _M_value(__z.__rep()) { }
1458 complex<float>::complex(
const complex<long double>& __z)
1459 : _M_value(__z.__rep()) { }
1462 complex<double>::complex(
const complex<long double>& __z)
1463 : _M_value(__z.__rep()) { }
1468 #if _GLIBCXX_EXTERN_TEMPLATE
1470 extern template ostream& operator<<(ostream&, const complex<float>&);
1472 extern template ostream& operator<<(ostream&, const complex<double>&);
1474 extern template ostream& operator<<(ostream&, const complex<long double>&);
1476 #ifdef _GLIBCXX_USE_WCHAR_T
1478 extern template wostream& operator<<(wostream&, const complex<float>&);
1480 extern template wostream& operator<<(wostream&, const complex<double>&);
1482 extern template wostream& operator<<(wostream&, const complex<long double>&);
1488 _GLIBCXX_END_NAMESPACE
1490 _GLIBCXX_BEGIN_NAMESPACE(__gnu_cxx)
1493 template<typename _Tp, typename _Up>
1494 struct __promote_2<std::complex<_Tp>, _Up>
1500 template<
typename _Tp,
typename _Up>
1501 struct __promote_2<_Tp, std::complex<_Up> >
1507 template<
typename _Tp,
typename _Up>
1508 struct __promote_2<std::complex<_Tp>,
std::complex<_Up> >
1514 _GLIBCXX_END_NAMESPACE
1516 #ifdef __GXX_EXPERIMENTAL_CXX0X__
1517 # if defined(_GLIBCXX_INCLUDE_AS_TR1)
1518 # error C++0x header cannot be included from TR1 header
1520 # if defined(_GLIBCXX_INCLUDE_AS_CXX0X)
1523 # define _GLIBCXX_INCLUDE_AS_CXX0X
1524 # define _GLIBCXX_BEGIN_NAMESPACE_TR1
1525 # define _GLIBCXX_END_NAMESPACE_TR1
1526 # define _GLIBCXX_TR1
1528 # undef _GLIBCXX_TR1
1529 # undef _GLIBCXX_END_NAMESPACE_TR1
1530 # undef _GLIBCXX_BEGIN_NAMESPACE_TR1
1531 # undef _GLIBCXX_INCLUDE_AS_CXX0X
1534 _GLIBCXX_BEGIN_NAMESPACE(std)
1538 template<typename _Tp> std::complex<_Tp> proj(const std::complex<_Tp>&);
1540 template<typename _Tp>
1542 __complex_proj(const std::complex<_Tp>& __z)
1544 const _Tp __den = (__z.real() * __z.real()
1545 + __z.imag() * __z.imag() + _Tp(1.0));
1548 (_Tp(2.0) * __z.imag()) / __den);
1551 #if _GLIBCXX_USE_C99_COMPLEX
1552 inline __complex__
float
1553 __complex_proj(__complex__
float __z)
1554 {
return __builtin_cprojf(__z); }
1556 inline __complex__
double
1557 __complex_proj(__complex__
double __z)
1558 {
return __builtin_cproj(__z); }
1560 inline __complex__
long double
1561 __complex_proj(
const __complex__
long double& __z)
1562 {
return __builtin_cprojl(__z); }
1564 template<
typename _Tp>
1567 {
return __complex_proj(__z.__rep()); }
1569 template<
typename _Tp>
1572 {
return __complex_proj(__z); }
1575 template<
typename _Tp>
1579 typedef typename __gnu_cxx::__promote<_Tp>::__type __type;
1583 _GLIBCXX_END_NAMESPACE