float FTEMP1;
int ITEMP1;
#define ELLIPSE(val,x,ps,c,xnorm,nsects)\
for(ITEMP1=0;ITEMP1<nsects;ITEMP1++) {\
FTEMP1 = x + c[ITEMP1][0] * ps[ITEMP1][0] + c[ITEMP1][2] * ps[ITEMP1][1] - c[ITEMP1][1] * ps[ITEMP1][2] - c[ITEMP1][3] * ps[ITEMP1][3];\
ps[ITEMP1][1] = ps[ITEMP1][0];\
ps[ITEMP1][0] = x;\
ps[ITEMP1][3] = ps[ITEMP1][2];\
ps[ITEMP1][2] = FTEMP1;\
x = FTEMP1; }\
val = x*xnorm;
