scalar2.c

#include <math.h>

#include "anicca.h"
#include "scalar2.h"

SF2(bdivide, D, B, *x / (D)*y)
SF2(idivide, D, I, *x / (D)*y)
SF2(ddivide, D, D, *x / *y)
static Z zdivide(Z x, Z y) { Z z; ZR(z)=ZR(x)/ZR(y); ZI(z)=ZI(x)/ZI(y); R z; }
SF2(jdivide, Z, Z, zdivide(*x,*y))

SF2(btimes, B, B, *x * *y)
SF2(itimes, I, I, *x * *y)
SF2(dtimes, D, D, *x * *y)
static Z ztimes(Z x, Z y) { D a=ZR(x),b=ZR(y),c=ZI(x),d=ZI(y);
    Z z; ZR(z)=a*c-b*d; ZI(z)=b*c+a*d;
    R z;
}
SF2(jtimes, Z, Z, ztimes(*x,*y))

SF2(bplus, I, B, *x + *y)
SF2(iplus, I, I, *x + *y)
SF2(dplus, D, D, *x + *y)
static Z zplus(Z x, Z y) { Z z; ZR(z)=ZR(x)+ZR(y); ZI(z)=ZI(x)+ZI(y); R z; }
SF2(jplus, Z, Z, zplus(*x,*y))

SF2(bminus, I, B, *x - *y)
SF2(iminus, I, I, *x - *y)
SF2(dminus, D, D, *x - *y)
static Z zminus(Z x, Z y) { Z z; ZR(z)=ZR(x)-ZR(y); ZI(z)=ZI(x)-ZI(y); R z; }
SF2(jminus, Z, Z, zminus(*x,*y))

SF2(blthan, B, B, *x < *y)
SF2(ilthan, B, I, *x < *y)
SF2(dlthan, B, D, *x < *y)
/*SF2(jplus, Z, Z, zplus(x,y))*/

SF2(bequal, B, B, *x == *y)
SF2(iequal, B, I, *x == *y)
SF2(dequal, B, D, *x == *y)
/*SF2(jplus, Z, Z, zplus(x,y))*/

SF2(bgthan, B, B, *x > *y)
SF2(igthan, B, I, *x > *y)
SF2(dgthan, B, D, *x > *y)
/*SF2(jplus, Z, Z, zplus(x,y))*/

SF2(bpower, B, B, *x >= *y)
SF2(ipower, D, I, pow((D)*x,(D)*y))
SF2(dpower, D, D, pow(*x,*y))

SF2(bresidue, B, B, *x % *y)
SF2(iresidue, I, I, *x % *y)
SF2(dresidue, D, D, *x + *y)
/*SF2(jplus, Z, Z, zplus(x,y))*/