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
|
/*
** -*- C -*-
**
**
** Made by Raffaele.Nutricato@tiscali.it
**
** Copyright Raffaele Nutricato
*/
/* Modified by Arnaud Torset */
#ifndef __OPLOGEQ_H__
#define __OPLOGEQ_H__
#include "floatComplex.h"
#include "doubleComplex.h"
#include "types.h"
#define s0s0OpLogEqs0(in1,in2) (float) (in1 == in2)
#define d0d0OpLogEqd0(in1,in2) (double) (in1 == in2)
#define c0c0OpLogEqs0(in1,in2) (float) ((creals(in1) == creals(in2)) && (cimags(in1) == cimags(in2)))
#define z0z0OpLogEqd0(in1,in2) (double) ((zreals(in1) == zreals(in2)) && (zimags(in1) == zimags(in2)))
#define u80u80OpLogEqu80(in1,in2) (uint8) (in1 == in2)
#define i80i80OpLogEqi80(in1,in2) (int8) (in1 == in2)
#define u160u160OpLogEqu160(in1,in2) (uint16) (in1 == in2)
#define i160i160OpLogEqi160(in1,in2) (int16) (in1 == in2)
#define s0c0OpLogEqs0(in1,in2) (float) ((in1==creals(in2)) && (0==cimags(in2)))
#define d0z0OpLogEqd0(in1,in2) (double) ((in1==zreals(in2)) && (0==zimags(in2)))
#define c0s0OpLogEqs0(in1,in2) s0c0OpLogEqs0(in2,in1)
#define z0d0OpLogEqd0(in1,in2) d0z0OpLogEqd0(in2,in1)
#define u80d0OpLogEqu80(in1,in2) (uint8) (in1 == (uint8)in2)
#define i80d0OpLogEqi80(in1,in2) (int8) (in1 == (int8)in2)
#define u160d0OpLogEqu160(in1,in2) (uint16) (in1 == (uint16)in2)
#define i160d0OpLogEqi160(in1,in2) (int16) (in1 == (int16)in2)
#define d0u80OpLogEqu80(in1,in2) (uint8) ((uint8)in1 == in2)
#define d0i80OpLogEqi80(in1,in2) (int8) ((int8)in1 == in2)
#define d0u160OpLogEqu160(in1,in2) (uint16) ((uint16)in1 == in2)
#define d0i160OpLogEqi160(in1,in2) (int16) ((int16)in1 == in2)
#define s2s0OpLogEqs2(in1,size,in2,out) {int i;\
for (i=0;i<size[0]*size[1];i++) out[i]=(float)(in1[i]==in2);\
}
#define d2d0OpLogEqd2(in1,size,in2,out) {int i;\
for (i=0;i<size[0]*size[1];i++) out[i]=(double)(in1[i]==in2);\
}
#define c2c0OpLogEqs2(in1,size,in2,out) {int i;\
for (i=0;i<size[0]*size[1];i++) out[i]=(float)((creals(in1[i])==creals(in2))&&(cimags(in1[i])==cimags(in2)));\
}
#define z2z0OpLogEqd2(in1,size,in2,out) {int i;\
for (i=0;i<size[0]*size[1];i++) out[i]=(double)((zreals(in1[i])==zreals(in2))&&(zimags(in1[i])==zimags(in2)));\
}
#define c2s0OpLogEqs2(in1,size,in2,out) c2c0OpLogEqs2(in1,size,FloatComplex(in2,0),out)
#define z2d0OpLogEqd2(in1,size,in2,out) z2z0OpLogEqd2(in1,size,DoubleComplex(in2,0),out)
#define s2c0OpLogEqs2(in1,size,in2,out) {int i;\
for (i=0;i<size[0]*size[1];i++) out[i]=(float)((in1[i]==creals(in2))&&(cimags(in2)==0));\
}
#define d2z0OpLogEqd2(in1,size,in2,out) {int i;\
for (i=0;i<size[0]*size[1];i++) out[i]=(double)((in1[i]==zreals(in2))&&(zimags(in2)==0));\
}
#define u82u80OpLogEqu82(in1,size,in2,out) {int i;\
for (i=0;i<size[0]*size[1];i++) out[i]=(uint8)(in1[i]==in2);\
}
#define i82i80OpLogEqi82(in1,size,in2,out) {int i;\
for (i=0;i<size[0]*size[1];i++) out[i]=(int8)(in1[i]==in2);\
}
#define u162u160OpLogEqu162(in1,size,in2,out) {int i;\
for (i=0;i<size[0]*size[1];i++) out[i]=(uint16)(in1[i]==in2);\
}
#define i162i160OpLogEqi162(in1,size,in2,out) {int i;\
for (i=0;i<size[0]*size[1];i++) out[i]=(int16)(in1[i]==in2);\
}
#define s0s2OpLogEqs2(in1,in2,inSize,out) s2s0OpLogEqs2(in2,inSize,in1,out)
#define c0s2OpLogEqs2(in1,in2,inSize,out) s2c0OpLogEqs2(in2,inSize,in1,out)
#define d0d2OpLogEqd2(in1,in2,inSize,out) d2d0OpLogEqd2(in2,inSize,in1,out)
#define z0d2OpLogEqd2(in1,in2,inSize,out) d2z0OpLogEqd2(in2,inSize,in1,out)
#define s0c2OpLogEqs2(in1,in2,inSize,out) c2s0OpLogEqs2(in2,inSize,in1,out)
#define c0c2OpLogEqs2(in1,in2,inSize,out) c2c0OpLogEqs2(in2,inSize,in1,out)
#define d0z2OpLogEqd2(in1,in2,inSize,out) z2d0OpLogEqd2(in2,inSize,in1,out)
#define z0z2OpLogEqd2(in1,in2,inSize,out) z2z0OpLogEqd2(in2,inSize,in1,out)
#define u80u82OpLogEqu82(in1,in2,inSize,out) u82u80OpLogEqu82(in2,inSize,in1,out)
#define i80i82OpLogEqi82(in1,in2,inSize,out) i82i80OpLogEqi82(in2,inSize,in1,out)
#define u160u162OpLogEqu162(in1,in2,inSize,out) u162u160OpLogEq162(in2,inSize,in1,out)
#define i160i162OpLogEqi162(in1,in2,inSize,out) i162i160OpLogEqi162(in2,inSize,in1,out)
/* we must have size1=size2 */
#define s2s2OpLogEqs2(in1,size1,in2,size2,out) {int i;\
for (i=0;i<size1[0]*size2[1];i++) out[i]=(float)(in1[i]==in2[i]);\
}
#define d2d2OpLogEqd2(in1,size1,in2,size2,out) {int i;\
for (i=0;i<size1[0]*size2[1];i++) out[i]=(double)(in1[i]==in2[i]);\
}
#define c2c2OpLogEqs2(in1,size1,in2,size2,out) {int i;\
for (i=0;i<size1[0]*size2[1];i++) \
out[i]=(float)((creals(in1[i])==creals(in2[i]))&&(cimags(in1[i])==cimags(in2[i])));\
}
#define z2z2OpLogEqd2(in1,size1,in2,size2,out) {int i;\
for (i=0;i<size1[0]*size2[1];i++) \
out[i]=(double)((zreals(in1[i])==zreals(in2[i]))&&(zimags(in1[i])==zimags(in2[i])));\
}
#define s2c2OpLogEqs2(in1,size1,in2,size2,out) {int i;\
for (i=0;i<size1[0]*size2[1];i++) \
out[i]=(float)((in1[i]==creals(in2[i]))&&(0==cimags(in2[i])));\
}
#define d2z2OpLogEqd2(in1,size1,in2,size2,out) {int i;\
for (i=0;i<size1[0]*size2[1];i++) \
out[i]=(double)((in1[i]==zreals(in2[i]))&&(0==zimags(in2[i])));\
}
#define c2s2OpLogEqs2(in1,size1,in2,size2,out) s2c2OpLogEqs2(in2,size2,in1,size1,out)
#define z2d2OpLogEqd2(in1,size1,in2,size2,out) d2z2OpLogEqd2(in2,size2,in1,size1,out)
#define u82u82OpLogEqu82(in1,size1,in2,size2,out) {int i;\
for (i=0;i<size1[0]*size2[1];i++) out[i]=(uint8)(in1[i]==in2[i]);\
}
#define i82i82OpLogEqi82(in1,size1,in2,size2,out) {int i;\
for (i=0;i<size1[0]*size2[1];i++) out[i]=(int8)(in1[i]==in2[i]);\
}
#define u162u162OpLogEqu162(in1,size1,in2,size2,out) {int i;\
for (i=0;i<size1[0]*size2[1];i++) out[i]=(uint16)(in1[i]==in2[i]);\
}
#define i162i162OpLogEqi162(in1,size1,in2,size2,out) {int i;\
for (i=0;i<size1[0]*size2[1];i++) out[i]=(int16)(in1[i]==in2[i]);\
}
#endif /* !__OPLOGEQ_H__ */
|
