pio2
static const double pio2[] = {
fw += pio2[k] * q[i+k];
ah = pio2;
ah = -pio2;
pio2 = 1.5707963267948965580e+00,
a = pio2;
a = -pio2;
pio2 = 1.570796326794896558e+00,
return ((hy >= 0)? pio2 : -pio2);
return ((hy >= 0)? pio2 : -pio2);
return ((hy >= 0)? pio2 : -pio2);
pio2 = 1.57079632679489661923132169163975144209858469968755L,
t = pio2 - y;
*err = -(y - (pio2 - t));
t = y - pio2;
*err = y - (pio2 + t);
s = pio2 - v;
t -= (v - (pio2 - s));
s = v - pio2;
t = -(t - (v - (s + pio2)));
return ((hy >= 0)? pio2 : -pio2);
return ((hy >= 0)? pio2 : -pio2);
return ((hy >= 0)? pio2 : -pio2);
ah0 = -pio2;
ah1 = pio2;
ah1 = -pio2;
ah2 = pio2;
ah2 = -pio2;
pio2 = 1.5707963267948965580e+00,
ah0 = pio2;
ah1 = pio2;
ah1 = -pio2;
ah2 = pio2;
ah2 = -pio2;
pio2 = 1.5707963267948965580e+00,
ah0 = pio2;
ah0 = -pio2;
static const double pio2[] = {
q[0] = pio2[0] * pr[0];
q[1] = pio2[0] * pr[1] + pio2[1] * pr[0];
q[2] = pio2[0] * pr[2] + pio2[1] * pr[1] + pio2[2] * pr[0];
q[3] = pio2[0] * pr[3] + pio2[1] * pr[2] + pio2[2] * pr[1]
+ pio2[3] * pr[0];
q[i] = pio2[0] * pr[i] + pio2[1] * pr[i-1] + pio2[2] * pr[i-2]
+ pio2[3] * pr[i-3] + pio2[4] * pr[i-4];