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
|
// Nothing like a little CPP abuse to brighten everyone's day. Would
// have been nice to do this with template functions but as of 2.5.x,
// g++ seems to fail to resolve them properly.
#define DO_VS_OP(OP) \
int l = a.length (); \
T *result = 0; \
if (l > 0) \
{ \
result = new T [l]; \
const T *x = a.data (); \
for (int i = 0; i < l; i++) \
result[i] = x[i] OP s; \
}
#define DO_SV_OP(OP) \
int l = a.length (); \
T *result = 0; \
if (l > 0) \
{ \
result = new T [l]; \
const T *x = a.data (); \
for (int i = 0; i < l; i++) \
result[i] = s OP x[i]; \
}
#define DO_VV_OP(OP) \
T *result = 0; \
if (l > 0) \
{ \
result = new T [l]; \
const T *x = a.data (); \
const T *y = b.data (); \
for (int i = 0; i < l; i++) \
result[i] = x[i] OP y[i]; \
}
#define NEG_V \
int l = a.length (); \
T *result = 0; \
if (l > 0) \
{ \
result = new T [l]; \
const T *x = a.data (); \
for (int i = 0; i < l; i++) \
result[i] = -x[i]; \
}
#define DO_VS_OP2(OP) \
int l = a.length (); \
if (l > 0) \
{ \
T *tmp = a.fortran_vec (); \
for (int i = 0; i < l; i++) \
tmp[i] OP s; \
}
#define DO_VV_OP2(OP) \
do \
{ \
T *a_tmp = a.fortran_vec (); \
const T *b_tmp = b.data (); \
for (int i = 0; i < l; i++) \
a_tmp[i] OP b_tmp[i]; \
} \
while (0)
|
