|
@@ -272,7 +272,7 @@ static const double ps2[5] = {
|
|
|
if(ix>=0x40200000) {p = pr8; q= ps8;}
|
|
|
else if(ix>=0x40122E8B){p = pr5; q= ps5;}
|
|
|
else if(ix>=0x4006DB6D){p = pr3; q= ps3;}
|
|
|
- else if(ix>=0x40000000){p = pr2; q= ps2;}
|
|
|
+ else {p = pr2; q= ps2;}
|
|
|
z = one/(x*x);
|
|
|
r = p[0]+z*(p[1]+z*(p[2]+z*(p[3]+z*(p[4]+z*p[5]))));
|
|
|
s = one+z*(q[0]+z*(q[1]+z*(q[2]+z*(q[3]+z*q[4]))));
|
|
@@ -368,7 +368,7 @@ static const double qs2[6] = {
|
|
|
if(ix>=0x40200000) {p = qr8; q= qs8;}
|
|
|
else if(ix>=0x40122E8B){p = qr5; q= qs5;}
|
|
|
else if(ix>=0x4006DB6D){p = qr3; q= qs3;}
|
|
|
- else if(ix>=0x40000000){p = qr2; q= qs2;}
|
|
|
+ else {p = qr2; q= qs2;}
|
|
|
z = one/(x*x);
|
|
|
r = p[0]+z*(p[1]+z*(p[2]+z*(p[3]+z*(p[4]+z*p[5]))));
|
|
|
s = one+z*(q[0]+z*(q[1]+z*(q[2]+z*(q[3]+z*(q[4]+z*q[5])))));
|