|
sl@0
|
1 |
/*
|
|
sl@0
|
2 |
* © Portions copyright (c) 2006-2007 Nokia Corporation. All rights reserved.
|
|
sl@0
|
3 |
*
|
|
sl@0
|
4 |
* Copyright (c) 1999
|
|
sl@0
|
5 |
* Silicon Graphics Computer Systems, Inc.
|
|
sl@0
|
6 |
*
|
|
sl@0
|
7 |
* Copyright (c) 1999
|
|
sl@0
|
8 |
* Boris Fomitchev
|
|
sl@0
|
9 |
*
|
|
sl@0
|
10 |
* This material is provided "as is", with absolutely no warranty expressed
|
|
sl@0
|
11 |
* or implied. Any use is at your own risk.
|
|
sl@0
|
12 |
*
|
|
sl@0
|
13 |
* Permission to use or copy this software for any purpose is hereby granted
|
|
sl@0
|
14 |
* without fee, provided the above notices are retained on all copies.
|
|
sl@0
|
15 |
* Permission to modify the code and to distribute modified code is granted,
|
|
sl@0
|
16 |
* provided the above notices are retained, and a notice that the code was
|
|
sl@0
|
17 |
* modified is included with the above copyright notice.
|
|
sl@0
|
18 |
*
|
|
sl@0
|
19 |
*/
|
|
sl@0
|
20 |
# include "stlport_prefix.h"
|
|
sl@0
|
21 |
|
|
sl@0
|
22 |
|
|
sl@0
|
23 |
// Trigonometric and hyperbolic functions for complex<float>,
|
|
sl@0
|
24 |
// complex<double>, and complex<long double>
|
|
sl@0
|
25 |
|
|
sl@0
|
26 |
|
|
sl@0
|
27 |
#include "complex_impl.h"
|
|
sl@0
|
28 |
|
|
sl@0
|
29 |
#include <cfloat>
|
|
sl@0
|
30 |
#include <cmath>
|
|
sl@0
|
31 |
|
|
sl@0
|
32 |
#if defined(__LIBSTD_CPP_SYMBIAN32_WSD__) || defined(_STLP_LIBSTD_CPP_NO_STATIC_VAR_)
|
|
sl@0
|
33 |
#include "libstdcppwsd.h"
|
|
sl@0
|
34 |
# endif
|
|
sl@0
|
35 |
|
|
sl@0
|
36 |
_STLP_BEGIN_NAMESPACE
|
|
sl@0
|
37 |
|
|
sl@0
|
38 |
|
|
sl@0
|
39 |
//----------------------------------------------------------------------
|
|
sl@0
|
40 |
// helpers
|
|
sl@0
|
41 |
|
|
sl@0
|
42 |
#ifdef __sgi
|
|
sl@0
|
43 |
static const union { unsigned int i; float f; } float_ulimit = { 0x42b2d4fc };
|
|
sl@0
|
44 |
static const float float_limit = float_ulimit.f;
|
|
sl@0
|
45 |
static union {
|
|
sl@0
|
46 |
struct { unsigned int h; unsigned int l; } w;
|
|
sl@0
|
47 |
double d;
|
|
sl@0
|
48 |
} double_ulimit = { 0x408633ce, 0x8fb9f87d };
|
|
sl@0
|
49 |
static const double double_limit = double_ulimit.d;
|
|
sl@0
|
50 |
static union {
|
|
sl@0
|
51 |
struct { unsigned int h[2]; unsigned int l[2]; } w;
|
|
sl@0
|
52 |
long double ld;
|
|
sl@0
|
53 |
} ldouble_ulimit = {0x408633ce, 0x8fb9f87e, 0xbd23b659, 0x4e9bd8b1};
|
|
sl@0
|
54 |
# ifndef _STLP_NO_LONG_DOUBLE
|
|
sl@0
|
55 |
static const long double ldouble_limit = ldouble_ulimit.ld;
|
|
sl@0
|
56 |
# endif
|
|
sl@0
|
57 |
#else
|
|
sl@0
|
58 |
#if defined(__LIBSTD_CPP_SYMBIAN32_WSD__) || defined(_STLP_LIBSTD_CPP_NO_STATIC_VAR_)
|
|
sl@0
|
59 |
void complex_trig_limit_init()
|
|
sl@0
|
60 |
{
|
|
sl@0
|
61 |
get_complex_trig_float_limit() = _STLP_LOGF(FLT_MAX);
|
|
sl@0
|
62 |
get_complex_trig_double_limit() = _STLP_DO_LOG(double)(DBL_MAX);
|
|
sl@0
|
63 |
}
|
|
sl@0
|
64 |
# else
|
|
sl@0
|
65 |
static const float float_limit = _STLP_LOGF(FLT_MAX);
|
|
sl@0
|
66 |
static const double double_limit = _STLP_DO_LOG(double)(DBL_MAX);
|
|
sl@0
|
67 |
# endif //__LIBSTD_CPP_SYMBIAN32_WSD__
|
|
sl@0
|
68 |
# ifndef _STLP_NO_LONG_DOUBLE
|
|
sl@0
|
69 |
static const long double ldouble_limit = _STLP_LOGL(LDBL_MAX);
|
|
sl@0
|
70 |
# endif
|
|
sl@0
|
71 |
#endif
|
|
sl@0
|
72 |
|
|
sl@0
|
73 |
|
|
sl@0
|
74 |
//----------------------------------------------------------------------
|
|
sl@0
|
75 |
// sin
|
|
sl@0
|
76 |
|
|
sl@0
|
77 |
_STLP_EXP_DECLSPEC complex<float> _STLP_CALL sin(const complex<float>& z) {
|
|
sl@0
|
78 |
return complex<float>(_STLP_SINF(z._M_re) * _STLP_COSHF(z._M_im),
|
|
sl@0
|
79 |
_STLP_COSF(z._M_re) * _STLP_SINHF(z._M_im));
|
|
sl@0
|
80 |
}
|
|
sl@0
|
81 |
|
|
sl@0
|
82 |
_STLP_EXP_DECLSPEC complex<double> _STLP_CALL sin(const complex<double>& z) {
|
|
sl@0
|
83 |
return complex<double>(_STLP_SIN(z._M_re) * _STLP_COSH(z._M_im),
|
|
sl@0
|
84 |
_STLP_COS(z._M_re) * _STLP_SINH(z._M_im));
|
|
sl@0
|
85 |
}
|
|
sl@0
|
86 |
|
|
sl@0
|
87 |
#ifndef _STLP_NO_LONG_DOUBLE
|
|
sl@0
|
88 |
_STLP_EXP_DECLSPEC complex<long double> _STLP_CALL sin(const complex<long double>& z) {
|
|
sl@0
|
89 |
return complex<long double>(_STLP_SINL(z._M_re) * _STLP_COSHL(z._M_im),
|
|
sl@0
|
90 |
_STLP_COSL(z._M_re) * _STLP_SINHL(z._M_im));
|
|
sl@0
|
91 |
}
|
|
sl@0
|
92 |
#endif
|
|
sl@0
|
93 |
|
|
sl@0
|
94 |
//----------------------------------------------------------------------
|
|
sl@0
|
95 |
// cos
|
|
sl@0
|
96 |
|
|
sl@0
|
97 |
_STLP_EXP_DECLSPEC complex<float> _STLP_CALL cos(const complex<float>& z) {
|
|
sl@0
|
98 |
return complex<float>(_STLP_COSF(z._M_re) * _STLP_COSHF(z._M_im),
|
|
sl@0
|
99 |
-_STLP_SINF(z._M_re) * _STLP_SINHF(z._M_im));
|
|
sl@0
|
100 |
}
|
|
sl@0
|
101 |
|
|
sl@0
|
102 |
_STLP_EXP_DECLSPEC complex<double> _STLP_CALL cos(const complex<double>& z) {
|
|
sl@0
|
103 |
return complex<double>(_STLP_COS(z._M_re) * _STLP_COSH(z._M_im),
|
|
sl@0
|
104 |
-_STLP_SIN(z._M_re) * _STLP_SINH(z._M_im));
|
|
sl@0
|
105 |
}
|
|
sl@0
|
106 |
|
|
sl@0
|
107 |
#ifndef _STLP_NO_LONG_DOUBLE
|
|
sl@0
|
108 |
_STLP_EXP_DECLSPEC complex<long double> _STLP_CALL cos(const complex<long double>& z) {
|
|
sl@0
|
109 |
return complex<long double>(_STLP_COSL(z._M_re) * _STLP_COSHL(z._M_im),
|
|
sl@0
|
110 |
-_STLP_SINL(z._M_re) * _STLP_SINHL(z._M_im));
|
|
sl@0
|
111 |
}
|
|
sl@0
|
112 |
# endif
|
|
sl@0
|
113 |
|
|
sl@0
|
114 |
//----------------------------------------------------------------------
|
|
sl@0
|
115 |
// tan
|
|
sl@0
|
116 |
|
|
sl@0
|
117 |
_STLP_EXP_DECLSPEC complex<float> _STLP_CALL tan(const complex<float>& z) {
|
|
sl@0
|
118 |
float re2 = 2.f * z._M_re;
|
|
sl@0
|
119 |
float im2 = 2.f * z._M_im;
|
|
sl@0
|
120 |
|
|
sl@0
|
121 |
#if defined(__LIBSTD_CPP_SYMBIAN32_WSD__) || defined(_STLP_LIBSTD_CPP_NO_STATIC_VAR_)
|
|
sl@0
|
122 |
if (_STLP_ABSF(im2) > get_complex_trig_float_limit())
|
|
sl@0
|
123 |
# else
|
|
sl@0
|
124 |
if (_STLP_ABSF(im2) > float_limit)
|
|
sl@0
|
125 |
# endif
|
|
sl@0
|
126 |
return complex<float>(0.f, (im2 > 0 ? 1.f : -1.f));
|
|
sl@0
|
127 |
else {
|
|
sl@0
|
128 |
float den = _STLP_COSF(re2) + _STLP_COSHF(im2);
|
|
sl@0
|
129 |
return complex<float>(_STLP_SINF(re2) / den, _STLP_SINHF(im2) / den);
|
|
sl@0
|
130 |
}
|
|
sl@0
|
131 |
}
|
|
sl@0
|
132 |
|
|
sl@0
|
133 |
_STLP_EXP_DECLSPEC complex<double> _STLP_CALL tan(const complex<double>& z) {
|
|
sl@0
|
134 |
double re2 = 2. * z._M_re;
|
|
sl@0
|
135 |
double im2 = 2. * z._M_im;
|
|
sl@0
|
136 |
|
|
sl@0
|
137 |
#if defined(__LIBSTD_CPP_SYMBIAN32_WSD__) || defined(_STLP_LIBSTD_CPP_NO_STATIC_VAR_)
|
|
sl@0
|
138 |
if (fabs(float(im2)) > get_complex_trig_double_limit())
|
|
sl@0
|
139 |
# else
|
|
sl@0
|
140 |
if (fabs(float(im2)) > double_limit)
|
|
sl@0
|
141 |
# endif //__LIBSTD_CPP_SYMBIAN32_WSD__
|
|
sl@0
|
142 |
return complex<double>(0., (im2 > 0 ? 1. : -1.));
|
|
sl@0
|
143 |
else {
|
|
sl@0
|
144 |
double den = _STLP_COS(re2) + _STLP_COSH(im2);
|
|
sl@0
|
145 |
return complex<double>(_STLP_SIN(re2) / den, _STLP_SINH(im2) / den);
|
|
sl@0
|
146 |
}
|
|
sl@0
|
147 |
}
|
|
sl@0
|
148 |
|
|
sl@0
|
149 |
#ifndef _STLP_NO_LONG_DOUBLE
|
|
sl@0
|
150 |
_STLP_EXP_DECLSPEC complex<long double> _STLP_CALL tan(const complex<long double>& z) {
|
|
sl@0
|
151 |
long double re2 = 2.l * z._M_re;
|
|
sl@0
|
152 |
long double im2 = 2.l * z._M_im;
|
|
sl@0
|
153 |
if (_STLP_ABSL(im2) > ldouble_limit)
|
|
sl@0
|
154 |
return complex<long double>(0.l, (im2 > 0 ? 1.l : -1.l));
|
|
sl@0
|
155 |
else {
|
|
sl@0
|
156 |
long double den = _STLP_COSL(re2) + _STLP_COSHL(im2);
|
|
sl@0
|
157 |
return complex<long double>(_STLP_SINL(re2) / den, _STLP_SINHL(im2) / den);
|
|
sl@0
|
158 |
}
|
|
sl@0
|
159 |
}
|
|
sl@0
|
160 |
|
|
sl@0
|
161 |
# endif
|
|
sl@0
|
162 |
|
|
sl@0
|
163 |
//----------------------------------------------------------------------
|
|
sl@0
|
164 |
// sinh
|
|
sl@0
|
165 |
|
|
sl@0
|
166 |
_STLP_EXP_DECLSPEC complex<float> _STLP_CALL sinh(const complex<float>& z) {
|
|
sl@0
|
167 |
return complex<float>(_STLP_SINHF(z._M_re) * _STLP_COSF(z._M_im),
|
|
sl@0
|
168 |
_STLP_COSHF(z._M_re) * _STLP_SINF(z._M_im));
|
|
sl@0
|
169 |
}
|
|
sl@0
|
170 |
|
|
sl@0
|
171 |
_STLP_EXP_DECLSPEC complex<double> _STLP_CALL sinh(const complex<double>& z) {
|
|
sl@0
|
172 |
return complex<double>(_STLP_SINH(z._M_re) * _STLP_COS(z._M_im),
|
|
sl@0
|
173 |
_STLP_COSH(z._M_re) * _STLP_SIN(z._M_im));
|
|
sl@0
|
174 |
}
|
|
sl@0
|
175 |
|
|
sl@0
|
176 |
#ifndef _STLP_NO_LONG_DOUBLE
|
|
sl@0
|
177 |
_STLP_EXP_DECLSPEC complex<long double> _STLP_CALL sinh(const complex<long double>& z) {
|
|
sl@0
|
178 |
return complex<long double>(_STLP_SINHL(z._M_re) * _STLP_COSL(z._M_im),
|
|
sl@0
|
179 |
_STLP_COSHL(z._M_re) * _STLP_SINL(z._M_im));
|
|
sl@0
|
180 |
}
|
|
sl@0
|
181 |
#endif
|
|
sl@0
|
182 |
|
|
sl@0
|
183 |
//----------------------------------------------------------------------
|
|
sl@0
|
184 |
// cosh
|
|
sl@0
|
185 |
|
|
sl@0
|
186 |
_STLP_EXP_DECLSPEC complex<float> _STLP_CALL cosh(const complex<float>& z) {
|
|
sl@0
|
187 |
return complex<float>(_STLP_COSHF(z._M_re) * _STLP_COSF(z._M_im),
|
|
sl@0
|
188 |
_STLP_SINHF(z._M_re) * _STLP_SINF(z._M_im));
|
|
sl@0
|
189 |
}
|
|
sl@0
|
190 |
|
|
sl@0
|
191 |
_STLP_EXP_DECLSPEC complex<double> _STLP_CALL cosh(const complex<double>& z) {
|
|
sl@0
|
192 |
return complex<double>(_STLP_COSH(z._M_re) * _STLP_COS(z._M_im),
|
|
sl@0
|
193 |
_STLP_SINH(z._M_re) * _STLP_SIN(z._M_im));
|
|
sl@0
|
194 |
}
|
|
sl@0
|
195 |
|
|
sl@0
|
196 |
#ifndef _STLP_NO_LONG_DOUBLE
|
|
sl@0
|
197 |
_STLP_EXP_DECLSPEC complex<long double> _STLP_CALL cosh(const complex<long double>& z) {
|
|
sl@0
|
198 |
return complex<long double>(_STLP_COSHL(z._M_re) * _STLP_COSL(z._M_im),
|
|
sl@0
|
199 |
_STLP_SINHL(z._M_re) * _STLP_SINL(z._M_im));
|
|
sl@0
|
200 |
}
|
|
sl@0
|
201 |
#endif
|
|
sl@0
|
202 |
|
|
sl@0
|
203 |
//----------------------------------------------------------------------
|
|
sl@0
|
204 |
// tanh
|
|
sl@0
|
205 |
|
|
sl@0
|
206 |
_STLP_EXP_DECLSPEC complex<float> _STLP_CALL tanh(const complex<float>& z) {
|
|
sl@0
|
207 |
float re2 = 2.f * z._M_re;
|
|
sl@0
|
208 |
float im2 = 2.f * z._M_im;
|
|
sl@0
|
209 |
#if defined(__LIBSTD_CPP_SYMBIAN32_WSD__) || defined(_STLP_LIBSTD_CPP_NO_STATIC_VAR_)
|
|
sl@0
|
210 |
if (_STLP_ABSF(re2) > get_complex_trig_float_limit())
|
|
sl@0
|
211 |
# else
|
|
sl@0
|
212 |
if (_STLP_ABSF(re2) > float_limit)
|
|
sl@0
|
213 |
# endif //__LIBSTD_CPP_SYMBIAN32_WSD__
|
|
sl@0
|
214 |
return complex<float>((re2 > 0 ? 1.f : -1.f), 0.f);
|
|
sl@0
|
215 |
else {
|
|
sl@0
|
216 |
float den = _STLP_COSHF(re2) + _STLP_COSF(im2);
|
|
sl@0
|
217 |
return complex<float>(_STLP_SINHF(re2) / den, _STLP_SINF(im2) / den);
|
|
sl@0
|
218 |
}
|
|
sl@0
|
219 |
}
|
|
sl@0
|
220 |
|
|
sl@0
|
221 |
_STLP_EXP_DECLSPEC complex<double> _STLP_CALL tanh(const complex<double>& z) {
|
|
sl@0
|
222 |
double re2 = 2. * z._M_re;
|
|
sl@0
|
223 |
double im2 = 2. * z._M_im;
|
|
sl@0
|
224 |
|
|
sl@0
|
225 |
#if defined(__LIBSTD_CPP_SYMBIAN32_WSD__) || defined(_STLP_LIBSTD_CPP_NO_STATIC_VAR_)
|
|
sl@0
|
226 |
if (fabs(float(re2)) > get_complex_trig_double_limit())
|
|
sl@0
|
227 |
# else
|
|
sl@0
|
228 |
if (fabs(float(re2)) > double_limit)
|
|
sl@0
|
229 |
# endif //__LIBSTD_CPP_SYMBIAN32_WSD__
|
|
sl@0
|
230 |
return complex<double>((re2 > 0 ? 1. : -1.), 0.);
|
|
sl@0
|
231 |
else {
|
|
sl@0
|
232 |
double den = _STLP_COSH(re2) + _STLP_COS(im2);
|
|
sl@0
|
233 |
return complex<double>(_STLP_SINH(re2) / den, _STLP_SIN(im2) / den);
|
|
sl@0
|
234 |
}
|
|
sl@0
|
235 |
}
|
|
sl@0
|
236 |
|
|
sl@0
|
237 |
#ifndef _STLP_NO_LONG_DOUBLE
|
|
sl@0
|
238 |
_STLP_EXP_DECLSPEC complex<long double> _STLP_CALL tanh(const complex<long double>& z) {
|
|
sl@0
|
239 |
long double re2 = 2.l * z._M_re;
|
|
sl@0
|
240 |
long double im2 = 2.l * z._M_im;
|
|
sl@0
|
241 |
if (_STLP_ABSL(re2) > ldouble_limit)
|
|
sl@0
|
242 |
return complex<long double>((re2 > 0 ? 1.l : -1.l), 0.l);
|
|
sl@0
|
243 |
else {
|
|
sl@0
|
244 |
long double den = _STLP_COSHL(re2) + _STLP_COSL(im2);
|
|
sl@0
|
245 |
return complex<long double>(_STLP_SINHL(re2) / den, _STLP_SINL(im2) / den);
|
|
sl@0
|
246 |
}
|
|
sl@0
|
247 |
}
|
|
sl@0
|
248 |
#endif
|
|
sl@0
|
249 |
_STLP_END_NAMESPACE
|