diff options
Diffstat (limited to 'src/c/operations/interfaces/int_OpLogEq.h')
| -rw-r--r-- | src/c/operations/interfaces/int_OpLogEq.h | 94 |
1 files changed, 94 insertions, 0 deletions
diff --git a/src/c/operations/interfaces/int_OpLogEq.h b/src/c/operations/interfaces/int_OpLogEq.h new file mode 100644 index 00000000..b15a8b6a --- /dev/null +++ b/src/c/operations/interfaces/int_OpLogEq.h @@ -0,0 +1,94 @@ +/* +** -*- 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" + +#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 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 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 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) + +/* 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) +#endif /* !__OPLOGGT_H__ */ |