1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
|
/*@@
@header cctk_Complex.h
@date Tue Dec 14 12:28:05 1999
@author Tom Goodale
@desc
Prototypes for complex numbers.
@enddesc
@version $Header$
@@*/
#ifndef _CCTK_COMPLEX_H_
#define _CCTK_COMPLEX_H_
#ifdef __cplusplus
# define EXTERN_C_BEGIN extern "C" {
# define EXTERN_C_END }
#else
# define EXTERN_C_BEGIN
# define EXTERN_C_END
#endif
#ifndef DEFINE_CCTK_COMPLEX_EXTERN_FUNCTIONS
# define PREFIX static inline
#else
# define PREFIX
#endif
#ifdef __cplusplus
/* declare C++ overloaded operators */
# include <complex>
# include <ostream>
# define DECLARE_CMPLX_CXX_OPERATORS(CCTK_Cmplx, cctk_real, cctk_complex) \
PREFIX cctk_complex operator+ (cctk_complex const & a); \
PREFIX cctk_complex operator- (cctk_complex const & a); \
PREFIX cctk_complex conj (cctk_complex const & a); \
PREFIX cctk_real abs (cctk_complex const & a); \
PREFIX cctk_real arg (cctk_complex const & a); \
PREFIX cctk_real norm (cctk_complex const & a); \
PREFIX cctk_complex operator+ (cctk_complex const & a, cctk_complex const & b); \
PREFIX cctk_complex operator+ (cctk_real const & a, cctk_complex const & b); \
PREFIX cctk_complex operator+ (cctk_complex const & a, cctk_real const & b); \
PREFIX cctk_complex operator- (cctk_complex const & a, cctk_complex const & b); \
PREFIX cctk_complex operator- (cctk_real const & a, cctk_complex const & b); \
PREFIX cctk_complex operator- (cctk_complex const & a, cctk_real const & b); \
PREFIX cctk_complex operator* (cctk_complex const & a, cctk_complex const & b); \
PREFIX cctk_complex operator* (cctk_real const & a, cctk_complex const & b); \
PREFIX cctk_complex operator* (cctk_complex const & a, cctk_real const & b); \
PREFIX cctk_complex operator/ (cctk_complex const & a, cctk_complex const & b); \
PREFIX cctk_complex operator/ (cctk_real const & a, cctk_complex const & b); \
PREFIX cctk_complex operator/ (cctk_complex const & a, cctk_real const & b); \
PREFIX cctk_complex operator+= (cctk_complex & a, cctk_complex const & b); \
PREFIX cctk_complex operator+= (cctk_complex & a, cctk_real const & b); \
PREFIX cctk_complex operator-= (cctk_complex & a, cctk_complex const & b); \
PREFIX cctk_complex operator-= (cctk_complex & a, cctk_real const & b); \
PREFIX cctk_complex operator*= (cctk_complex & a, cctk_complex const & b); \
PREFIX cctk_complex operator*= (cctk_complex & a, cctk_real const & b); \
PREFIX cctk_complex operator/= (cctk_complex & a, cctk_complex const & b); \
PREFIX cctk_complex operator/= (cctk_complex & a, cctk_real const & b); \
PREFIX cctk_complex pow (cctk_complex const & a, int i); \
PREFIX cctk_complex pow (cctk_complex const & a, cctk_real const & r); \
PREFIX cctk_complex pow (cctk_complex const & a, cctk_complex const & b); \
PREFIX cctk_complex sin (cctk_complex const & a); \
PREFIX cctk_complex cos (cctk_complex const & a); \
PREFIX cctk_complex exp (cctk_complex const & a); \
PREFIX cctk_complex log (cctk_complex const & a); \
PREFIX cctk_complex sqrt (cctk_complex const & a); \
PREFIX std::ostream & operator << (std::ostream & os, cctk_complex const & a);
#else
/* declare no C++ overloaded operators */
# define DECLARE_CMPLX_CXX_OPERATORS(CCTK_Cmplx, cctk_real, cctk_complex)
#endif
/* macro to declare a set of complex functions for a given precision */
#define DECLARE_CMPLX_FUNCTIONS(CCTK_Cmplx, cctk_real, cctk_complex) \
EXTERN_C_BEGIN \
PREFIX cctk_complex CCTK_Cmplx (cctk_real Re, cctk_real Im); \
PREFIX cctk_real CCTK_Cmplx##Real (cctk_complex a); \
PREFIX cctk_real CCTK_Cmplx##Imag (cctk_complex a); \
PREFIX cctk_complex CCTK_Cmplx##Neg (cctk_complex a); \
PREFIX cctk_complex CCTK_Cmplx##Conjg (cctk_complex a); \
PREFIX cctk_real CCTK_Cmplx##Abs (cctk_complex a); \
PREFIX cctk_real CCTK_Cmplx##Arg (cctk_complex a); \
PREFIX cctk_real CCTK_Cmplx##Norm (cctk_complex a); \
PREFIX cctk_complex CCTK_Cmplx##Add (cctk_complex a, cctk_complex b); \
PREFIX cctk_complex CCTK_Cmplx##Sub (cctk_complex a, cctk_complex b); \
PREFIX cctk_complex CCTK_Cmplx##Mul (cctk_complex a, cctk_complex b); \
PREFIX cctk_complex CCTK_Cmplx##Div (cctk_complex a, cctk_complex b); \
PREFIX cctk_complex CCTK_Cmplx##CPow (cctk_complex a, cctk_complex b); \
PREFIX cctk_complex CCTK_Cmplx##Sin (cctk_complex a); \
PREFIX cctk_complex CCTK_Cmplx##Cos (cctk_complex a); \
PREFIX cctk_complex CCTK_Cmplx##Exp (cctk_complex a); \
PREFIX cctk_complex CCTK_Cmplx##Log (cctk_complex a); \
PREFIX cctk_complex CCTK_Cmplx##Sqrt (cctk_complex a); \
PREFIX cctk_complex CCTK_Cmplx##Pow (cctk_complex a, cctk_real b); \
PREFIX cctk_complex CCTK_Cmplx##IPow (cctk_complex a, int b); \
EXTERN_C_END \
DECLARE_CMPLX_CXX_OPERATORS(CCTK_Cmplx, cctk_real, cctk_complex)
/* declare complex functions for all available precisions */
#ifdef HAVE_CCTK_REAL4
DECLARE_CMPLX_FUNCTIONS (CCTK_Cmplx8, CCTK_REAL4, CCTK_COMPLEX8)
#endif
#ifdef HAVE_CCTK_REAL8
DECLARE_CMPLX_FUNCTIONS (CCTK_Cmplx16, CCTK_REAL8, CCTK_COMPLEX16)
#endif
#ifdef HAVE_CCTK_REAL16
DECLARE_CMPLX_FUNCTIONS (CCTK_Cmplx32, CCTK_REAL16, CCTK_COMPLEX32)
#endif
#undef PREFIX
/* declare the default precision complex functions as #define'd macros */
#ifdef CCTK_REAL_PRECISION_4
#define CCTK_Cmplx CCTK_Cmplx8
#define CCTK_CmplxReal CCTK_Cmplx8Real
#define CCTK_CmplxImag CCTK_Cmplx8Imag
#define CCTK_CmplxNeg CCTK_Cmplx8Neg
#define CCTK_CmplxConjg CCTK_Cmplx8Conjg
#define CCTK_CmplxAbs CCTK_Cmplx8Abs
#define CCTK_CmplxArg CCTK_Cmplx8Arg
#define CCTK_CmplxNorm CCTK_Cmplx8Norm
#define CCTK_CmplxAdd CCTK_Cmplx8Add
#define CCTK_CmplxSub CCTK_Cmplx8Sub
#define CCTK_CmplxMul CCTK_Cmplx8Mul
#define CCTK_CmplxDiv CCTK_Cmplx8Div
#define CCTK_CmplxCPow CCTK_Cmplx8CPow
#define CCTK_CmplxSin CCTK_Cmplx8Sin
#define CCTK_CmplxCos CCTK_Cmplx8Cos
#define CCTK_CmplxExp CCTK_Cmplx8Exp
#define CCTK_CmplxLog CCTK_Cmplx8Log
#define CCTK_CmplxSqrt CCTK_Cmplx8Sqrt
#define CCTK_CmplxPow CCTK_Cmplx8Pow
#define CCTK_CmplxIPow CCTK_Cmplx8IPow
#elif CCTK_REAL_PRECISION_8
#define CCTK_Cmplx CCTK_Cmplx16
#define CCTK_CmplxReal CCTK_Cmplx16Real
#define CCTK_CmplxImag CCTK_Cmplx16Imag
#define CCTK_CmplxNeg CCTK_Cmplx16Neg
#define CCTK_CmplxConjg CCTK_Cmplx16Conjg
#define CCTK_CmplxAbs CCTK_Cmplx16Abs
#define CCTK_CmplxArg CCTK_Cmplx16Arg
#define CCTK_CmplxNorm CCTK_Cmplx16Norm
#define CCTK_CmplxAdd CCTK_Cmplx16Add
#define CCTK_CmplxSub CCTK_Cmplx16Sub
#define CCTK_CmplxMul CCTK_Cmplx16Mul
#define CCTK_CmplxDiv CCTK_Cmplx16Div
#define CCTK_CmplxCPow CCTK_Cmplx16CPow
#define CCTK_CmplxSin CCTK_Cmplx16Sin
#define CCTK_CmplxCos CCTK_Cmplx16Cos
#define CCTK_CmplxExp CCTK_Cmplx16Exp
#define CCTK_CmplxLog CCTK_Cmplx16Log
#define CCTK_CmplxSqrt CCTK_Cmplx16Sqrt
#define CCTK_CmplxPow CCTK_Cmplx16Pow
#define CCTK_CmplxIPow CCTK_Cmplx16IPow
#elif CCTK_REAL_PRECISION_16
#define CCTK_Cmplx CCTK_Cmplx32
#define CCTK_CmplxReal CCTK_Cmplx32Real
#define CCTK_CmplxImag CCTK_Cmplx32Imag
#define CCTK_CmplxNeg CCTK_Cmplx32Neg
#define CCTK_CmplxConjg CCTK_Cmplx32Conjg
#define CCTK_CmplxAbs CCTK_Cmplx32Abs
#define CCTK_CmplxArg CCTK_Cmplx32Arg
#define CCTK_CmplxNorm CCTK_Cmplx32Norm
#define CCTK_CmplxAdd CCTK_Cmplx32Add
#define CCTK_CmplxSub CCTK_Cmplx32Sub
#define CCTK_CmplxMul CCTK_Cmplx32Mul
#define CCTK_CmplxDiv CCTK_Cmplx32Div
#define CCTK_CmplxCPow CCTK_Cmplx32CPow
#define CCTK_CmplxSin CCTK_Cmplx32Sin
#define CCTK_CmplxCos CCTK_Cmplx32Cos
#define CCTK_CmplxExp CCTK_Cmplx32Exp
#define CCTK_CmplxLog CCTK_Cmplx32Log
#define CCTK_CmplxSqrt CCTK_Cmplx32Sqrt
#define CCTK_CmplxPow CCTK_Cmplx32Pow
#define CCTK_CmplxIPow CCTK_Cmplx32IPow
#endif
#ifndef DEFINE_CCTK_COMPLEX_EXTERN_FUNCTIONS
# define DEFINE_CCTK_COMPLEX_INLINE_FUNCTIONS
# include "../main/Complex.c"
# undef DEFINE_CCTK_COMPLEX_INLINE_FUNCTIONS
#endif
#endif /* _CCTK_COMPLEX_H_ */
|
