123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145 |
- /*
- * Copyright (c) 1999
- * Silicon Graphics Computer Systems, Inc.
- *
- * Copyright (c) 1999
- * Boris Fomitchev
- *
- * This material is provided "as is", with absolutely no warranty expressed
- * or implied. Any use is at your own risk.
- *
- * Permission to use or copy this software for any purpose is hereby granted
- * without fee, provided the above notices are retained on all copies.
- * Permission to modify the code and to distribute modified code is granted,
- * provided the above notices are retained, and a notice that the code was
- * modified is included with the above copyright notice.
- *
- */
- #ifndef _STLP_COMPLEX_C
- #define _STLP_COMPLEX_C
- #ifndef _STLP_INTERNAL_COMPLEX
- # include <stl/_complex.h>
- #endif
- #if !defined (_STLP_USE_NO_IOSTREAMS)
- # ifndef _STLP_INTERNAL_ISTREAM
- # include <stl/_istream.h>
- # endif
- # ifndef _STLP_INTERNAL_SSTREAM
- # include <stl/_sstream.h>
- # endif
- # ifndef _STLP_STRING_IO_H
- # include <stl/_string_io.h>
- # endif
- #endif
- _STLP_BEGIN_NAMESPACE
- // Non-inline member functions.
- template <class _Tp>
- void complex<_Tp>::_div(const _Tp& __z1_r, const _Tp& __z1_i,
- const _Tp& __z2_r, const _Tp& __z2_i,
- _Tp& __res_r, _Tp& __res_i) {
- _Tp __ar = __z2_r >= 0 ? __z2_r : -__z2_r;
- _Tp __ai = __z2_i >= 0 ? __z2_i : -__z2_i;
- if (__ar <= __ai) {
- _Tp __ratio = __z2_r / __z2_i;
- _Tp __denom = __z2_i * (1 + __ratio * __ratio);
- __res_r = (__z1_r * __ratio + __z1_i) / __denom;
- __res_i = (__z1_i * __ratio - __z1_r) / __denom;
- }
- else {
- _Tp __ratio = __z2_i / __z2_r;
- _Tp __denom = __z2_r * (1 + __ratio * __ratio);
- __res_r = (__z1_r + __z1_i * __ratio) / __denom;
- __res_i = (__z1_i - __z1_r * __ratio) / __denom;
- }
- }
- template <class _Tp>
- void complex<_Tp>::_div(const _Tp& __z1_r,
- const _Tp& __z2_r, const _Tp& __z2_i,
- _Tp& __res_r, _Tp& __res_i) {
- _Tp __ar = __z2_r >= 0 ? __z2_r : -__z2_r;
- _Tp __ai = __z2_i >= 0 ? __z2_i : -__z2_i;
- if (__ar <= __ai) {
- _Tp __ratio = __z2_r / __z2_i;
- _Tp __denom = __z2_i * (1 + __ratio * __ratio);
- __res_r = (__z1_r * __ratio) / __denom;
- __res_i = - __z1_r / __denom;
- }
- else {
- _Tp __ratio = __z2_i / __z2_r;
- _Tp __denom = __z2_r * (1 + __ratio * __ratio);
- __res_r = __z1_r / __denom;
- __res_i = - (__z1_r * __ratio) / __denom;
- }
- }
- // I/O.
- #if !defined (_STLP_USE_NO_IOSTREAMS)
- // Complex output, in the form (re,im). We use a two-step process
- // involving stringstream so that we get the padding right.
- template <class _Tp, class _CharT, class _Traits>
- basic_ostream<_CharT, _Traits>& _STLP_CALL
- operator<<(basic_ostream<_CharT, _Traits>& __os, const complex<_Tp>& __z) {
- basic_ostringstream<_CharT, _Traits, allocator<_CharT> > __tmp;
- __tmp.flags(__os.flags());
- __tmp.imbue(__os.getloc());
- __tmp.precision(__os.precision());
- __tmp << '(' << __z.real() << ',' << __z.imag() << ')';
- return __os << __tmp.str();
- }
- // Complex input from arbitrary streams. Note that results in some
- // locales may be confusing, since the decimal character varies with
- // locale and the separator between real and imaginary parts does not.
- template <class _Tp, class _CharT, class _Traits>
- basic_istream<_CharT, _Traits>& _STLP_CALL
- operator>>(basic_istream<_CharT, _Traits>& __is, complex<_Tp>& __z) {
- _Tp __re = 0;
- _Tp __im = 0;
- const ctype<_CharT>& __c_type = *__is._M_ctype_facet();
- const char __punct[4] = "(,)";
- _CharT __wpunct[3];
- __c_type.widen(__punct, __punct + 3, __wpunct);
- _CharT __c;
- __is >> __c;
- if (_Traits::eq(__c, __wpunct[0])) { // Left paren
- __is >> __re >> __c;
- if (_Traits::eq(__c, __wpunct[1])) // Comma
- __is >> __im >> __c;
- if (!_Traits::eq(__c, __wpunct[2])) // Right paren
- __is.setstate(ios_base::failbit);
- }
- else {
- __is.putback(__c);
- __is >> __re;
- }
- if (__is)
- __z = complex<_Tp>(__re, __im);
- return __is;
- }
- #endif /* _STLP_USE_NO_IOSTREAMS */
- _STLP_END_NAMESPACE
- #endif /* _STLP_COMPLEX_C */
- // Local Variables:
- // mode:C++
- // End:
|