|
@@ -278,7 +278,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]))));
|
|
@@ -373,7 +373,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])))));
|