PIo2
vc(PIo2, 1.5707963267948966135E0 ,0fda,40c9,68c2,a221, 1, .C90FDAA22168C2)
ic(PIo2, 1.5707963267948965580E0 , 0, 1.921FB54442D18)
#define PIo2 vccast(PIo2)
if(x==zero) return(copysign(PIo2,signy));
if(!finite(y)) return(copysign(PIo2,signy));
hi = PIo2; lo = zero;
a = PIo2-a;
a = PIo2-a; /* rtn. sign(x)*C(PI/2-|x|) */
x = copysign(PIo2-a,x);
vc(PIo2, 1.5707963267948966135E0 ,0fda,40c9,68c2,a221, 1, .C90FDAA22168C2)
ic(PIo2, 1.5707963267948965580E0 , 0, 1.921FB54442D18)
#define PIo2 vccast(PIo2)
static const double PIo2[] = {
for(fw=0.0,k=0;k<=jp&&k<=jz-i;k++) fw += PIo2[k]*q[i+k];
for(fw=0.0,k=0;k<=jp&&k<=jz-i;k++) fw += PIo2[k]*q[i+k];
static const float PIo2[] = {