]>
Commit | Line | Data |
---|---|---|
6599da04 JM |
1 | // The template and inlines for the -*- C++ -*- complex number classes. |
2 | // Copyright (C) 1994 Free Software Foundation | |
3 | ||
4 | // This file is part of the GNU ANSI C++ Library. This library is free | |
5 | // software; you can redistribute it and/or modify it under the terms of | |
6 | // the GNU General Public License as published by the Free Software | |
7 | // Foundation; either version 2, or (at your option) any later version. | |
8 | ||
9 | // This library is distributed in the hope that it will be useful, | |
10 | // but WITHOUT ANY WARRANTY; without even the implied warranty of | |
11 | // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the | |
12 | // GNU General Public License for more details. | |
13 | ||
14 | // You should have received a copy of the GNU General Public License | |
15 | // along with this library; see the file COPYING. If not, write to the Free | |
16 | // Software Foundation, 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA. | |
17 | ||
18 | // As a special exception, if you link this library with files compiled | |
19 | // with a GNU compiler to produce an executable, this does not cause the | |
20 | // resulting executable to be covered by the GNU General Public License. | |
21 | // This exception does not however invalidate any other reasons why the | |
22 | // executable file might be covered by the GNU General Public License. | |
23 | ||
24 | // Written by Jason Merrill based upon the specification in the 27 May 1994 | |
25 | // C++ working paper, ANSI document X3J16/94-0098. | |
26 | ||
27 | #ifndef __COMPLEXT__ | |
28 | #define __COMPLEXT__ | |
29 | ||
30 | #ifdef __GNUG__ | |
31 | #pragma interface | |
32 | #endif | |
33 | ||
34 | #include <cmath> | |
35 | ||
36 | #if ! defined (__GNUG__) && ! defined (__attribute__) | |
37 | #define __attribute__(foo) /* Ignore. */ | |
38 | #endif | |
39 | ||
40 | class istream; | |
41 | class ostream; | |
42 | ||
43 | extern "C++" { | |
905a881f JM |
44 | template <class _FLT> class complex; |
45 | template <class _FLT> complex<_FLT>& | |
46 | __doapl (complex<_FLT>* ths, const complex<_FLT>& r); | |
47 | template <class _FLT> complex<_FLT>& | |
48 | __doami (complex<_FLT>* ths, const complex<_FLT>& r); | |
49 | template <class _FLT> complex<_FLT>& | |
50 | __doaml (complex<_FLT>* ths, const complex<_FLT>& r); | |
51 | template <class _FLT> complex<_FLT>& | |
52 | __doadv (complex<_FLT>* ths, const complex<_FLT>& r); | |
53 | ||
6599da04 JM |
54 | template <class _FLT> |
55 | class complex | |
56 | { | |
57 | public: | |
58 | complex (_FLT r = 0, _FLT i = 0): re (r), im (i) { } | |
59 | complex& operator += (const complex&); | |
60 | complex& operator -= (const complex&); | |
61 | complex& operator *= (const complex&); | |
62 | complex& operator /= (const complex&); | |
63 | _FLT real () const { return re; } | |
64 | _FLT imag () const { return im; } | |
65 | private: | |
66 | _FLT re, im; | |
67 | ||
905a881f JM |
68 | friend complex& __doapl<> (complex *, const complex&); |
69 | friend complex& __doami<> (complex *, const complex&); | |
70 | friend complex& __doaml<> (complex *, const complex&); | |
71 | friend complex& __doadv<> (complex *, const complex&); | |
6599da04 JM |
72 | }; |
73 | ||
74 | // Declare specializations. | |
75 | class complex<float>; | |
76 | class complex<double>; | |
77 | class complex<long double>; | |
78 | ||
79 | template <class _FLT> | |
80 | inline complex<_FLT>& | |
81 | __doapl (complex<_FLT>* ths, const complex<_FLT>& r) | |
82 | { | |
83 | ths->re += r.re; | |
84 | ths->im += r.im; | |
85 | return *ths; | |
86 | } | |
87 | template <class _FLT> | |
88 | inline complex<_FLT>& | |
89 | complex<_FLT>::operator += (const complex<_FLT>& r) | |
90 | { | |
91 | return __doapl (this, r); | |
92 | } | |
93 | ||
94 | template <class _FLT> | |
95 | inline complex<_FLT>& | |
96 | __doami (complex<_FLT>* ths, const complex<_FLT>& r) | |
97 | { | |
98 | ths->re -= r.re; | |
99 | ths->im -= r.im; | |
100 | return *ths; | |
101 | } | |
102 | template <class _FLT> | |
103 | inline complex<_FLT>& | |
104 | complex<_FLT>::operator -= (const complex<_FLT>& r) | |
105 | { | |
106 | return __doami (this, r); | |
107 | } | |
108 | ||
109 | template <class _FLT> | |
110 | inline complex<_FLT>& | |
111 | __doaml (complex<_FLT>* ths, const complex<_FLT>& r) | |
112 | { | |
113 | _FLT f = ths->re * r.re - ths->im * r.im; | |
114 | ths->im = ths->re * r.im + ths->im * r.re; | |
115 | ths->re = f; | |
116 | return *ths; | |
117 | } | |
118 | template <class _FLT> | |
119 | inline complex<_FLT>& | |
120 | complex<_FLT>::operator *= (const complex<_FLT>& r) | |
121 | { | |
122 | return __doaml (this, r); | |
123 | } | |
124 | ||
6599da04 JM |
125 | template <class _FLT> |
126 | inline complex<_FLT>& | |
127 | complex<_FLT>::operator /= (const complex<_FLT>& r) | |
128 | { | |
129 | return __doadv (this, r); | |
130 | } | |
131 | ||
132 | template <class _FLT> inline _FLT | |
133 | imag (const complex<_FLT>& x) __attribute__ ((const)); | |
134 | ||
135 | template <class _FLT> inline _FLT | |
136 | imag (const complex<_FLT>& x) | |
137 | { | |
138 | return x.imag (); | |
139 | } | |
140 | ||
141 | template <class _FLT> inline _FLT | |
142 | real (const complex<_FLT>& x) __attribute__ ((const)); | |
143 | ||
144 | template <class _FLT> inline _FLT | |
145 | real (const complex<_FLT>& x) | |
146 | { | |
147 | return x.real (); | |
148 | } | |
149 | ||
150 | template <class _FLT> inline complex<_FLT> | |
151 | operator + (const complex<_FLT>& x, const complex<_FLT>& y) __attribute__ ((const)); | |
152 | ||
153 | template <class _FLT> inline complex<_FLT> | |
154 | operator + (const complex<_FLT>& x, const complex<_FLT>& y) | |
155 | { | |
156 | return complex<_FLT> (real (x) + real (y), imag (x) + imag (y)); | |
157 | } | |
158 | ||
159 | template <class _FLT> inline complex<_FLT> | |
160 | operator + (const complex<_FLT>& x, _FLT y) __attribute__ ((const)); | |
161 | ||
162 | template <class _FLT> inline complex<_FLT> | |
163 | operator + (const complex<_FLT>& x, _FLT y) | |
164 | { | |
165 | return complex<_FLT> (real (x) + y, imag (x)); | |
166 | } | |
167 | ||
168 | template <class _FLT> inline complex<_FLT> | |
169 | operator + (_FLT x, const complex<_FLT>& y) __attribute__ ((const)); | |
170 | ||
171 | template <class _FLT> inline complex<_FLT> | |
172 | operator + (_FLT x, const complex<_FLT>& y) | |
173 | { | |
174 | return complex<_FLT> (x + real (y), imag (y)); | |
175 | } | |
176 | ||
177 | template <class _FLT> inline complex<_FLT> | |
178 | operator - (const complex<_FLT>& x, const complex<_FLT>& y) __attribute__ ((const)); | |
179 | ||
180 | template <class _FLT> inline complex<_FLT> | |
181 | operator - (const complex<_FLT>& x, const complex<_FLT>& y) | |
182 | { | |
183 | return complex<_FLT> (real (x) - real (y), imag (x) - imag (y)); | |
184 | } | |
185 | ||
186 | template <class _FLT> inline complex<_FLT> | |
187 | operator - (const complex<_FLT>& x, _FLT y) __attribute__ ((const)); | |
188 | ||
189 | template <class _FLT> inline complex<_FLT> | |
190 | operator - (const complex<_FLT>& x, _FLT y) | |
191 | { | |
192 | return complex<_FLT> (real (x) - y, imag (x)); | |
193 | } | |
194 | ||
195 | template <class _FLT> inline complex<_FLT> | |
196 | operator - (_FLT x, const complex<_FLT>& y) __attribute__ ((const)); | |
197 | ||
198 | template <class _FLT> inline complex<_FLT> | |
199 | operator - (_FLT x, const complex<_FLT>& y) | |
200 | { | |
201 | return complex<_FLT> (x - real (y), - imag (y)); | |
202 | } | |
203 | ||
204 | template <class _FLT> inline complex<_FLT> | |
205 | operator * (const complex<_FLT>& x, const complex<_FLT>& y) __attribute__ ((const)); | |
206 | ||
207 | template <class _FLT> inline complex<_FLT> | |
208 | operator * (const complex<_FLT>& x, const complex<_FLT>& y) | |
209 | { | |
210 | return complex<_FLT> (real (x) * real (y) - imag (x) * imag (y), | |
211 | real (x) * imag (y) + imag (x) * real (y)); | |
212 | } | |
213 | ||
214 | template <class _FLT> inline complex<_FLT> | |
215 | operator * (const complex<_FLT>& x, _FLT y) __attribute__ ((const)); | |
216 | ||
217 | template <class _FLT> inline complex<_FLT> | |
218 | operator * (const complex<_FLT>& x, _FLT y) | |
219 | { | |
220 | return complex<_FLT> (real (x) * y, imag (x) * y); | |
221 | } | |
222 | ||
223 | template <class _FLT> inline complex<_FLT> | |
224 | operator * (_FLT x, const complex<_FLT>& y) __attribute__ ((const)); | |
225 | ||
226 | template <class _FLT> inline complex<_FLT> | |
227 | operator * (_FLT x, const complex<_FLT>& y) | |
228 | { | |
229 | return complex<_FLT> (x * real (y), x * imag (y)); | |
230 | } | |
231 | ||
232 | template <class _FLT> complex<_FLT> | |
233 | operator / (const complex<_FLT>& x, _FLT y) __attribute__ ((const)); | |
234 | ||
235 | template <class _FLT> complex<_FLT> | |
236 | operator / (const complex<_FLT>& x, _FLT y) | |
237 | { | |
238 | return complex<_FLT> (real (x) / y, imag (x) / y); | |
239 | } | |
240 | ||
241 | template <class _FLT> inline complex<_FLT> | |
242 | operator + (const complex<_FLT>& x) __attribute__ ((const)); | |
243 | ||
244 | template <class _FLT> inline complex<_FLT> | |
245 | operator + (const complex<_FLT>& x) | |
246 | { | |
247 | return x; | |
248 | } | |
249 | ||
250 | template <class _FLT> inline complex<_FLT> | |
251 | operator - (const complex<_FLT>& x) __attribute__ ((const)); | |
252 | ||
253 | template <class _FLT> inline complex<_FLT> | |
254 | operator - (const complex<_FLT>& x) | |
255 | { | |
256 | return complex<_FLT> (-real (x), -imag (x)); | |
257 | } | |
258 | ||
259 | template <class _FLT> inline bool | |
260 | operator == (const complex<_FLT>& x, const complex<_FLT>& y) __attribute__ ((const)); | |
261 | ||
262 | template <class _FLT> inline bool | |
263 | operator == (const complex<_FLT>& x, const complex<_FLT>& y) | |
264 | { | |
265 | return real (x) == real (y) && imag (x) == imag (y); | |
266 | } | |
267 | ||
268 | template <class _FLT> inline bool | |
269 | operator == (const complex<_FLT>& x, _FLT y) __attribute__ ((const)); | |
270 | ||
271 | template <class _FLT> inline bool | |
272 | operator == (const complex<_FLT>& x, _FLT y) | |
273 | { | |
274 | return real (x) == y && imag (x) == 0; | |
275 | } | |
276 | ||
277 | template <class _FLT> inline bool | |
278 | operator == (_FLT x, const complex<_FLT>& y) __attribute__ ((const)); | |
279 | ||
280 | template <class _FLT> inline bool | |
281 | operator == (_FLT x, const complex<_FLT>& y) | |
282 | { | |
283 | return x == real (y) && imag (y) == 0; | |
284 | } | |
285 | ||
286 | template <class _FLT> inline bool | |
287 | operator != (const complex<_FLT>& x, const complex<_FLT>& y) __attribute__ ((const)); | |
288 | ||
289 | template <class _FLT> inline bool | |
290 | operator != (const complex<_FLT>& x, const complex<_FLT>& y) | |
291 | { | |
292 | return real (x) != real (y) || imag (x) != imag (y); | |
293 | } | |
294 | ||
295 | template <class _FLT> inline bool | |
296 | operator != (const complex<_FLT>& x, _FLT y) __attribute__ ((const)); | |
297 | ||
298 | template <class _FLT> inline bool | |
299 | operator != (const complex<_FLT>& x, _FLT y) | |
300 | { | |
301 | return real (x) != y || imag (x) != 0; | |
302 | } | |
303 | ||
304 | template <class _FLT> inline bool | |
305 | operator != (_FLT x, const complex<_FLT>& y) __attribute__ ((const)); | |
306 | ||
307 | template <class _FLT> inline bool | |
308 | operator != (_FLT x, const complex<_FLT>& y) | |
309 | { | |
310 | return x != real (y) || imag (y) != 0; | |
311 | } | |
312 | ||
313 | // Some targets don't provide a prototype for hypot when -ansi. | |
314 | extern "C" double hypot (double, double) __attribute__ ((const)); | |
315 | ||
316 | template <class _FLT> inline _FLT | |
317 | abs (const complex<_FLT>& x) __attribute__ ((const)); | |
318 | ||
319 | template <class _FLT> inline _FLT | |
320 | abs (const complex<_FLT>& x) | |
321 | { | |
322 | return hypot (real (x), imag (x)); | |
323 | } | |
324 | ||
325 | template <class _FLT> inline _FLT | |
326 | arg (const complex<_FLT>& x) __attribute__ ((const)); | |
327 | ||
328 | template <class _FLT> inline _FLT | |
329 | arg (const complex<_FLT>& x) | |
330 | { | |
331 | return atan2 (imag (x), real (x)); | |
332 | } | |
333 | ||
334 | template <class _FLT> inline complex<_FLT> | |
335 | polar (_FLT r, _FLT t) __attribute__ ((const)); | |
336 | ||
337 | template <class _FLT> inline complex<_FLT> | |
338 | polar (_FLT r, _FLT t) | |
339 | { | |
340 | return complex<_FLT> (r * cos (t), r * sin (t)); | |
341 | } | |
342 | ||
343 | template <class _FLT> inline complex<_FLT> | |
344 | conj (const complex<_FLT>& x) __attribute__ ((const)); | |
345 | ||
346 | template <class _FLT> inline complex<_FLT> | |
347 | conj (const complex<_FLT>& x) | |
348 | { | |
349 | return complex<_FLT> (real (x), -imag (x)); | |
350 | } | |
351 | ||
352 | template <class _FLT> inline _FLT | |
353 | norm (const complex<_FLT>& x) __attribute__ ((const)); | |
354 | ||
355 | template <class _FLT> inline _FLT | |
356 | norm (const complex<_FLT>& x) | |
357 | { | |
358 | return real (x) * real (x) + imag (x) * imag (x); | |
359 | } | |
360 | ||
361 | // Declarations of templates in complext.ccI | |
362 | ||
363 | template <class _FLT> complex<_FLT> | |
364 | operator / (const complex<_FLT>&, const complex<_FLT>&) __attribute__ ((const)); | |
365 | template <class _FLT> complex<_FLT> | |
366 | operator / (_FLT, const complex<_FLT>&) __attribute__ ((const)); | |
367 | template <class _FLT> complex<_FLT> | |
368 | cos (const complex<_FLT>&) __attribute__ ((const)); | |
369 | template <class _FLT> complex<_FLT> | |
370 | cosh (const complex<_FLT>&) __attribute__ ((const)); | |
371 | template <class _FLT> complex<_FLT> | |
372 | exp (const complex<_FLT>&) __attribute__ ((const)); | |
373 | template <class _FLT> complex<_FLT> | |
374 | log (const complex<_FLT>&) __attribute__ ((const)); | |
375 | template <class _FLT> complex<_FLT> | |
376 | pow (const complex<_FLT>&, const complex<_FLT>&) __attribute__ ((const)); | |
377 | template <class _FLT> complex<_FLT> | |
378 | pow (const complex<_FLT>&, _FLT) __attribute__ ((const)); | |
379 | template <class _FLT> complex<_FLT> | |
380 | pow (const complex<_FLT>&, int) __attribute__ ((const)); | |
381 | template <class _FLT> complex<_FLT> | |
382 | pow (_FLT, const complex<_FLT>&) __attribute__ ((const)); | |
383 | template <class _FLT> complex<_FLT> | |
384 | sin (const complex<_FLT>&) __attribute__ ((const)); | |
385 | template <class _FLT> complex<_FLT> | |
386 | sinh (const complex<_FLT>&) __attribute__ ((const)); | |
387 | template <class _FLT> complex<_FLT> | |
388 | sqrt (const complex<_FLT>&) __attribute__ ((const)); | |
389 | ||
390 | template <class _FLT> istream& operator >> (istream&, complex<_FLT>&); | |
391 | template <class _FLT> ostream& operator << (ostream&, const complex<_FLT>&); | |
392 | } // extern "C++" | |
393 | ||
394 | // Specializations and such | |
395 | ||
396 | #include <std/fcomplex.h> | |
397 | #include <std/dcomplex.h> | |
398 | #include <std/ldcomplex.h> | |
399 | ||
400 | #endif |