friend class complex<long double>;
// friend float abs<>(const complex<float>&);
- friend float arg<>(const complex<float>&);
friend complex<float> conj<>(const complex<float>&);
friend class complex<long double>;
// friend double abs<>(const complex<double>&);
- friend double arg<>(const complex<double>&);
friend complex<double> conj<>(const complex<double>&);
friend complex<double> cos<>(const complex<double>&);
friend class complex<double>;
// friend long double abs<>(const complex<long double>&);
- friend long double arg<>(const complex<long double>&);
friend complex<long double> conj<>(const complex<long double>&);
friend complex<long double> cos<>(const complex<long double>&);
template<typename _Tp>
inline bool
operator==(const complex<_Tp>& __x, const _Tp& __y)
- { return __x.real() == __y && __x.imag() == 0; }
+ { return __x.real() == __y && __x.imag() == _Tp(); }
template<typename _Tp>
inline bool
operator==(const _Tp& __x, const complex<_Tp>& __y)
- { return __x == __y.real() && 0 == __y.imag(); }
+ { return __x == __y.real() && _Tp() == __y.imag(); }
template<typename _Tp>
inline bool
template<typename _Tp>
inline bool
operator!=(const complex<_Tp>& __x, const _Tp& __y)
- { return __x.real() != __y || __x.imag() != 0; }
+ { return __x.real() != __y || __x.imag() != _Tp(); }
template<typename _Tp>
inline bool
operator!=(const _Tp& __x, const complex<_Tp>& __y)
- { return __x != __y.real() || 0 != __y.imag(); }
+ { return __x != __y.real() || _Tp() != __y.imag(); }
template<typename _Tp, typename _CharT, class _Traits>
basic_istream<_CharT, _Traits>&
__x /= __s; __y /= __s;
return __s * sqrt(__x * __x + __y * __y);
}
-
+
+ template<typename _Tp>
+ inline _Tp
+ arg(const complex<_Tp>& __z)
+ { return atan2(__z.imag(), __z.real()); }
+
+
+ template<typename _Tp>
+ inline complex<_Tp>
+ polar(const _Tp& __rho, const _Tp& __theta)
+ { return complex<_Tp>(__rho * cos(__theta), __rho * sin(__theta)); }
+
+ template<typename _Tp>
+ inline complex<_Tp>
+ conj(const complex<_Tp>& __z)
+ { return complex<_Tp>(__z.real(), -__z.imag()); }
// We use here a few more specializations.
template<>
// return hypot (__real__ __x._M_value, __imag__ __x._M_value);
// }
- template<>
- FLT
- arg(const complex<FLT>& __x)
- { return carg(__x._M_value); }
+// template<>
+// FLT
+// arg(const complex<FLT>& __x)
+// { return carg(__x._M_value); }
- template<>
- complex<FLT>
- polar(const FLT& __rho, const FLT& __theta)
- {
-#if 0
-// XXX
-// defined(_GLIBCPP_HAVE_SINCOS) && !defined(__osf__)
- // Although sincos does exist on OSF3.2 and OSF4.0 we cannot use it
- // since the necessary types are not defined in the headers.
- FLT __sinx, __cosx;
- sincos(__theta, &__sinx, &__cosx);
- return complex<FLT>(__rho * __cosx, __rho * __sinx);
-#else
- return complex<FLT>(__rho * cos(__theta), __rho * sin(__theta));
-#endif
- }
+// template<>
+// complex<FLT>
+// polar(const FLT& __rho, const FLT& __theta)
+// {
+// #if 0
+// // XXX
+// // defined(_GLIBCPP_HAVE_SINCOS) && !defined(__osf__)
+// // Although sincos does exist on OSF3.2 and OSF4.0 we cannot use it
+// // since the necessary types are not defined in the headers.
+// FLT __sinx, __cosx;
+// sincos(__theta, &__sinx, &__cosx);
+// return complex<FLT>(__rho * __cosx, __rho * __sinx);
+// #else
+// return complex<FLT>(__rho * cos(__theta), __rho * sin(__theta));
+// #endif
+// }
template<>
complex<FLT>
{ return complex<FLT>(ctanh(__x._M_value)); }
} // namespace std
-
-
-
-
-