|
sl@0
|
1 |
/*
|
|
sl@0
|
2 |
* © Portions copyright (c) 2006-2007 Nokia Corporation. All rights reserved.
|
|
sl@0
|
3 |
* Copyright (c) 1999
|
|
sl@0
|
4 |
* Silicon Graphics Computer Systems, Inc.
|
|
sl@0
|
5 |
*
|
|
sl@0
|
6 |
* Copyright (c) 1999
|
|
sl@0
|
7 |
* Boris Fomitchev
|
|
sl@0
|
8 |
*
|
|
sl@0
|
9 |
* This material is provided "as is", with absolutely no warranty expressed
|
|
sl@0
|
10 |
* or implied. Any use is at your own risk.
|
|
sl@0
|
11 |
*
|
|
sl@0
|
12 |
* Permission to use or copy this software for any purpose is hereby granted
|
|
sl@0
|
13 |
* without fee, provided the above notices are retained on all copies.
|
|
sl@0
|
14 |
* Permission to modify the code and to distribute modified code is granted,
|
|
sl@0
|
15 |
* provided the above notices are retained, and a notice that the code was
|
|
sl@0
|
16 |
* modified is included with the above copyright notice.
|
|
sl@0
|
17 |
*
|
|
sl@0
|
18 |
*/
|
|
sl@0
|
19 |
|
|
sl@0
|
20 |
# include "stlport_prefix.h"
|
|
sl@0
|
21 |
// exp, log, pow for complex<float>, complex<double>, and complex<long double>
|
|
sl@0
|
22 |
|
|
sl@0
|
23 |
#include <numeric>
|
|
sl@0
|
24 |
#include "complex_impl.h"
|
|
sl@0
|
25 |
|
|
sl@0
|
26 |
#if defined(__LIBSTD_CPP_SYMBIAN32_WSD__) || defined(_STLP_LIBSTD_CPP_NO_STATIC_VAR_)
|
|
sl@0
|
27 |
#include "libstdcppwsd.h"
|
|
sl@0
|
28 |
# endif
|
|
sl@0
|
29 |
|
|
sl@0
|
30 |
_STLP_BEGIN_NAMESPACE
|
|
sl@0
|
31 |
|
|
sl@0
|
32 |
//----------------------------------------------------------------------
|
|
sl@0
|
33 |
// exp
|
|
sl@0
|
34 |
|
|
sl@0
|
35 |
_STLP_EXP_DECLSPEC complex<float> _STLP_CALL
|
|
sl@0
|
36 |
exp(const complex<float>& z)
|
|
sl@0
|
37 |
{
|
|
sl@0
|
38 |
float expx = _STLP_EXPF(z._M_re);
|
|
sl@0
|
39 |
return complex<float>(expx * _STLP_COSF(z._M_im),
|
|
sl@0
|
40 |
expx * _STLP_SINF(z._M_im));
|
|
sl@0
|
41 |
}
|
|
sl@0
|
42 |
|
|
sl@0
|
43 |
_STLP_EXP_DECLSPEC complex<double> _STLP_CALL exp(const complex<double>& z)
|
|
sl@0
|
44 |
{
|
|
sl@0
|
45 |
double expx = _STLP_EXP(z._M_re);
|
|
sl@0
|
46 |
return complex<double>(expx * _STLP_COS(z._M_im),
|
|
sl@0
|
47 |
expx * _STLP_SIN(z._M_im));
|
|
sl@0
|
48 |
}
|
|
sl@0
|
49 |
|
|
sl@0
|
50 |
# ifndef _STLP_NO_LONG_DOUBLE
|
|
sl@0
|
51 |
_STLP_EXP_DECLSPEC complex<long double> _STLP_CALL exp(const complex<long double>& z)
|
|
sl@0
|
52 |
{
|
|
sl@0
|
53 |
long double expx = _STLP_EXPL(z._M_re);
|
|
sl@0
|
54 |
return complex<long double>(expx * _STLP_COSL(z._M_im),
|
|
sl@0
|
55 |
expx * _STLP_SINL(z._M_im));
|
|
sl@0
|
56 |
}
|
|
sl@0
|
57 |
# endif
|
|
sl@0
|
58 |
|
|
sl@0
|
59 |
//----------------------------------------------------------------------
|
|
sl@0
|
60 |
// log10
|
|
sl@0
|
61 |
|
|
sl@0
|
62 |
_STLP_EXP_DECLSPEC complex<float> _STLP_CALL log10(const complex<float>& z)
|
|
sl@0
|
63 |
{
|
|
sl@0
|
64 |
complex<float> r;
|
|
sl@0
|
65 |
|
|
sl@0
|
66 |
#if defined(__LIBSTD_CPP_SYMBIAN32_WSD__) || defined(_STLP_LIBSTD_CPP_NO_STATIC_VAR_)
|
|
sl@0
|
67 |
get_complex_exp_float_ln10_inv() = 1.f / _STLP_LOGF(10.f);
|
|
sl@0
|
68 |
r._M_im = _STLP_ATAN2F(z._M_im, z._M_re) * get_complex_exp_float_ln10_inv();
|
|
sl@0
|
69 |
# else
|
|
sl@0
|
70 |
static float ln10_inv = 1.f / _STLP_LOGF(10.f);
|
|
sl@0
|
71 |
r._M_im = _STLP_ATAN2F(z._M_im, z._M_re) * ln10_inv;
|
|
sl@0
|
72 |
# endif //__LIBSTD_CPP_SYMBIAN32_WSD__
|
|
sl@0
|
73 |
|
|
sl@0
|
74 |
r._M_re = _STLP_LOG10F(_STLP_HYPOTF(z._M_re, z._M_im));
|
|
sl@0
|
75 |
return r;
|
|
sl@0
|
76 |
}
|
|
sl@0
|
77 |
|
|
sl@0
|
78 |
_STLP_EXP_DECLSPEC complex<double> _STLP_CALL log10(const complex<double>& z)
|
|
sl@0
|
79 |
{
|
|
sl@0
|
80 |
complex<double> r;
|
|
sl@0
|
81 |
|
|
sl@0
|
82 |
#if defined(__LIBSTD_CPP_SYMBIAN32_WSD__) || defined(_STLP_LIBSTD_CPP_NO_STATIC_VAR_)
|
|
sl@0
|
83 |
get_complex_exp_double_ln10_inv() = 1. / _STLP_LOG(10.);
|
|
sl@0
|
84 |
r._M_im = _STLP_ATAN2(z._M_im, z._M_re) * get_complex_exp_double_ln10_inv();
|
|
sl@0
|
85 |
# else
|
|
sl@0
|
86 |
static double ln10_inv = 1. / _STLP_LOG(10.);
|
|
sl@0
|
87 |
r._M_im = _STLP_ATAN2(z._M_im, z._M_re) * ln10_inv;
|
|
sl@0
|
88 |
# endif //__LIBSTD_CPP_SYMBIAN32_WSD__
|
|
sl@0
|
89 |
|
|
sl@0
|
90 |
r._M_re = _STLP_LOG10(_STLP_HYPOT(z._M_re, z._M_im));
|
|
sl@0
|
91 |
return r;
|
|
sl@0
|
92 |
}
|
|
sl@0
|
93 |
|
|
sl@0
|
94 |
#ifndef _STLP_NO_LONG_DOUBLE
|
|
sl@0
|
95 |
_STLP_EXP_DECLSPEC complex<long double> _STLP_CALL log10(const complex<long double>& z)
|
|
sl@0
|
96 |
{
|
|
sl@0
|
97 |
complex<long double> result;
|
|
sl@0
|
98 |
#if defined(__LIBSTD_CPP_SYMBIAN32_WSD__) || defined(_STLP_LIBSTD_CPP_NO_STATIC_VAR_)
|
|
sl@0
|
99 |
get_complex_exp_long_double_ln10_inv() = 1.l / _STLP_LOGL(10.l);
|
|
sl@0
|
100 |
result._M_im = _STLP_ATAN2L(z._M_im, z._M_re) * get_complex_exp_long_double_ln10_inv();
|
|
sl@0
|
101 |
# else
|
|
sl@0
|
102 |
static long double ln10_inv = 1.l / _STLP_LOGL(10.l);
|
|
sl@0
|
103 |
result._M_im = _STLP_ATAN2L(z._M_im, z._M_re) * ln10_inv;
|
|
sl@0
|
104 |
# endif //__LIBSTD_CPP_SYMBIAN32_WSD__
|
|
sl@0
|
105 |
|
|
sl@0
|
106 |
result._M_re = _STLP_LOG10L(_STLP_HYPOTL(z._M_re, z._M_im));
|
|
sl@0
|
107 |
return result;
|
|
sl@0
|
108 |
}
|
|
sl@0
|
109 |
# endif
|
|
sl@0
|
110 |
|
|
sl@0
|
111 |
//----------------------------------------------------------------------
|
|
sl@0
|
112 |
// log
|
|
sl@0
|
113 |
|
|
sl@0
|
114 |
_STLP_EXP_DECLSPEC complex<float> _STLP_CALL log(const complex<float>& z)
|
|
sl@0
|
115 |
{
|
|
sl@0
|
116 |
complex<float> r;
|
|
sl@0
|
117 |
|
|
sl@0
|
118 |
r._M_im = _STLP_ATAN2F(z._M_im, z._M_re);
|
|
sl@0
|
119 |
r._M_re = _STLP_LOGF(_STLP_HYPOTF(z._M_re, z._M_im));
|
|
sl@0
|
120 |
return r;
|
|
sl@0
|
121 |
}
|
|
sl@0
|
122 |
|
|
sl@0
|
123 |
_STLP_EXP_DECLSPEC complex<double> _STLP_CALL log(const complex<double>& z)
|
|
sl@0
|
124 |
{
|
|
sl@0
|
125 |
complex<double> r;
|
|
sl@0
|
126 |
|
|
sl@0
|
127 |
r._M_im = _STLP_ATAN2(z._M_im, z._M_re);
|
|
sl@0
|
128 |
r._M_re = _STLP_LOG(_STLP_HYPOT(z._M_re, z._M_im));
|
|
sl@0
|
129 |
return r;
|
|
sl@0
|
130 |
}
|
|
sl@0
|
131 |
|
|
sl@0
|
132 |
#ifndef _STLP_NO_LONG_DOUBLE
|
|
sl@0
|
133 |
_STLP_EXP_DECLSPEC complex<long double> _STLP_CALL log(const complex<long double>& z)
|
|
sl@0
|
134 |
{
|
|
sl@0
|
135 |
complex<long double> result;
|
|
sl@0
|
136 |
|
|
sl@0
|
137 |
result._M_im = _STLP_ATAN2L(z._M_im, z._M_re);
|
|
sl@0
|
138 |
result._M_re = _STLP_LOGL(_STLP_HYPOTL(z._M_re, z._M_im));
|
|
sl@0
|
139 |
return result;
|
|
sl@0
|
140 |
}
|
|
sl@0
|
141 |
# endif
|
|
sl@0
|
142 |
|
|
sl@0
|
143 |
//----------------------------------------------------------------------
|
|
sl@0
|
144 |
// pow
|
|
sl@0
|
145 |
|
|
sl@0
|
146 |
_STLP_EXP_DECLSPEC complex<float> _STLP_CALL pow(const float& a, const complex<float>& b) {
|
|
sl@0
|
147 |
float logr = _STLP_LOGF(a);
|
|
sl@0
|
148 |
float x = _STLP_EXPF(logr*b._M_re);
|
|
sl@0
|
149 |
float y = logr*b._M_im;
|
|
sl@0
|
150 |
|
|
sl@0
|
151 |
return complex<float>(x * _STLP_COSF(y), x * _STLP_SINF(y));
|
|
sl@0
|
152 |
}
|
|
sl@0
|
153 |
|
|
sl@0
|
154 |
_STLP_EXP_DECLSPEC complex<float> _STLP_CALL pow(const complex<float>& z_in, int n) {
|
|
sl@0
|
155 |
complex<float> z = z_in;
|
|
sl@0
|
156 |
z = __power(z, (n < 0 ? -n : n), multiplies< complex<float> >());
|
|
sl@0
|
157 |
if (n < 0)
|
|
sl@0
|
158 |
return 1.f / z;
|
|
sl@0
|
159 |
else
|
|
sl@0
|
160 |
return z;
|
|
sl@0
|
161 |
}
|
|
sl@0
|
162 |
|
|
sl@0
|
163 |
_STLP_EXP_DECLSPEC complex<float> _STLP_CALL pow(const complex<float>& a, const float& b) {
|
|
sl@0
|
164 |
float logr = _STLP_LOGF(_STLP_HYPOTF(a._M_re,a._M_im));
|
|
sl@0
|
165 |
float logi = _STLP_ATAN2F(a._M_im, a._M_re);
|
|
sl@0
|
166 |
float x = _STLP_EXPF(logr * b);
|
|
sl@0
|
167 |
float y = logi * b;
|
|
sl@0
|
168 |
|
|
sl@0
|
169 |
return complex<float>(x * _STLP_COSF(y), x * _STLP_SINF(y));
|
|
sl@0
|
170 |
}
|
|
sl@0
|
171 |
|
|
sl@0
|
172 |
_STLP_EXP_DECLSPEC complex<float> _STLP_CALL pow(const complex<float>& a, const complex<float>& b) {
|
|
sl@0
|
173 |
float logr = _STLP_LOGF(_STLP_HYPOTF(a._M_re,a._M_im));
|
|
sl@0
|
174 |
float logi = _STLP_ATAN2F(a._M_im, a._M_re);
|
|
sl@0
|
175 |
float x = _STLP_EXPF(logr*b._M_re - logi*b._M_im);
|
|
sl@0
|
176 |
float y = logr*b._M_im + logi*b._M_re;
|
|
sl@0
|
177 |
|
|
sl@0
|
178 |
return complex<float>(x * _STLP_COSF(y), x * _STLP_SINF(y));
|
|
sl@0
|
179 |
}
|
|
sl@0
|
180 |
|
|
sl@0
|
181 |
|
|
sl@0
|
182 |
_STLP_EXP_DECLSPEC complex<double> _STLP_CALL pow(const double& a, const complex<double>& b) {
|
|
sl@0
|
183 |
double logr = _STLP_LOG(a);
|
|
sl@0
|
184 |
double x = _STLP_EXP(logr*b._M_re);
|
|
sl@0
|
185 |
double y = logr*b._M_im;
|
|
sl@0
|
186 |
|
|
sl@0
|
187 |
return complex<double>(x * _STLP_COS(y), x * _STLP_SIN(y));
|
|
sl@0
|
188 |
}
|
|
sl@0
|
189 |
|
|
sl@0
|
190 |
_STLP_EXP_DECLSPEC complex<double> _STLP_CALL pow(const complex<double>& z_in, int n) {
|
|
sl@0
|
191 |
complex<double> z = z_in;
|
|
sl@0
|
192 |
z = __power(z, (n < 0 ? -n : n), multiplies< complex<double> >());
|
|
sl@0
|
193 |
if (n < 0)
|
|
sl@0
|
194 |
#if !defined(__SC__) //*TY 04/15/2000 -
|
|
sl@0
|
195 |
return 1. / z;
|
|
sl@0
|
196 |
#else //*TY 04/15/2000 - added workaround for SCpp compiler
|
|
sl@0
|
197 |
return double(1.0) / z; //*TY 04/15/2000 - it incorrectly assign long double attribute to floating point literals
|
|
sl@0
|
198 |
#endif //*TY 04/15/2000 -
|
|
sl@0
|
199 |
else
|
|
sl@0
|
200 |
return z;
|
|
sl@0
|
201 |
}
|
|
sl@0
|
202 |
|
|
sl@0
|
203 |
_STLP_EXP_DECLSPEC complex<double> _STLP_CALL pow(const complex<double>& a, const double& b) {
|
|
sl@0
|
204 |
double logr = _STLP_LOG(_STLP_HYPOT(a._M_re,a._M_im));
|
|
sl@0
|
205 |
double logi = _STLP_ATAN2(a._M_im, a._M_re);
|
|
sl@0
|
206 |
double x = _STLP_EXP(logr * b);
|
|
sl@0
|
207 |
double y = logi * b;
|
|
sl@0
|
208 |
|
|
sl@0
|
209 |
return complex<double>(x * _STLP_COS(y), x * _STLP_SIN(y));
|
|
sl@0
|
210 |
}
|
|
sl@0
|
211 |
|
|
sl@0
|
212 |
_STLP_EXP_DECLSPEC complex<double> _STLP_CALL pow(const complex<double>& a, const complex<double>& b) {
|
|
sl@0
|
213 |
double logr = _STLP_LOG(_STLP_HYPOT(a._M_re,a._M_im));
|
|
sl@0
|
214 |
double logi = _STLP_ATAN2(a._M_im, a._M_re);
|
|
sl@0
|
215 |
double x = _STLP_EXP(logr*b._M_re - logi*b._M_im);
|
|
sl@0
|
216 |
double y = logr*b._M_im + logi*b._M_re;
|
|
sl@0
|
217 |
|
|
sl@0
|
218 |
return complex<double>(x * _STLP_COS(y), x * _STLP_SIN(y));
|
|
sl@0
|
219 |
}
|
|
sl@0
|
220 |
|
|
sl@0
|
221 |
|
|
sl@0
|
222 |
_STLP_EXP_DECLSPEC complex<long double> _STLP_CALL pow(const long double& a,
|
|
sl@0
|
223 |
const complex<long double>& b) {
|
|
sl@0
|
224 |
long double logr = _STLP_LOGL(a);
|
|
sl@0
|
225 |
long double x = _STLP_EXPL(logr*b._M_re);
|
|
sl@0
|
226 |
long double y = logr*b._M_im;
|
|
sl@0
|
227 |
|
|
sl@0
|
228 |
return complex<long double>(x * _STLP_COSL(y), x * _STLP_SINL(y));
|
|
sl@0
|
229 |
}
|
|
sl@0
|
230 |
|
|
sl@0
|
231 |
_STLP_EXP_DECLSPEC complex<long double> _STLP_CALL pow(const complex<long double>& z_in, int n) {
|
|
sl@0
|
232 |
complex<long double> z = z_in;
|
|
sl@0
|
233 |
z = __power(z, (n < 0 ? -n : n), multiplies< complex<long double> >());
|
|
sl@0
|
234 |
if (n < 0)
|
|
sl@0
|
235 |
return 1.l / z;
|
|
sl@0
|
236 |
else
|
|
sl@0
|
237 |
return z;
|
|
sl@0
|
238 |
}
|
|
sl@0
|
239 |
|
|
sl@0
|
240 |
_STLP_EXP_DECLSPEC complex<long double> _STLP_CALL pow(const complex<long double>& a,
|
|
sl@0
|
241 |
const long double& b) {
|
|
sl@0
|
242 |
long double logr = _STLP_LOGL(_STLP_HYPOTL(a._M_re,a._M_im));
|
|
sl@0
|
243 |
long double logi = _STLP_ATAN2L(a._M_im, a._M_re);
|
|
sl@0
|
244 |
long double x = _STLP_EXPL(logr * b);
|
|
sl@0
|
245 |
long double y = logi * b;
|
|
sl@0
|
246 |
|
|
sl@0
|
247 |
return complex<long double>(x * _STLP_COSL(y), x * _STLP_SINL(y));
|
|
sl@0
|
248 |
}
|
|
sl@0
|
249 |
|
|
sl@0
|
250 |
_STLP_EXP_DECLSPEC complex<long double> _STLP_CALL pow(const complex<long double>& a,
|
|
sl@0
|
251 |
const complex<long double>& b) {
|
|
sl@0
|
252 |
long double logr = _STLP_LOGL(_STLP_HYPOTL(a._M_re,a._M_im));
|
|
sl@0
|
253 |
long double logi = _STLP_ATAN2L(a._M_im, a._M_re);
|
|
sl@0
|
254 |
long double x = _STLP_EXPL(logr*b._M_re - logi*b._M_im);
|
|
sl@0
|
255 |
long double y = logr*b._M_im + logi*b._M_re;
|
|
sl@0
|
256 |
|
|
sl@0
|
257 |
return complex<long double>(x * _STLP_COSL(y), x * _STLP_SINL(y));
|
|
sl@0
|
258 |
}
|
|
sl@0
|
259 |
|
|
sl@0
|
260 |
_STLP_END_NAMESPACE
|
|
sl@0
|
261 |
|