@@ -305,6 +305,7 @@ pzero (double x)
int32_t ix;
GET_HIGH_WORD (ix, x);
ix &= 0x7fffffff;
+ /* ix >= 0x40000000 for all calls to this function. */
if (ix >= 0x41b00000)
{
return one;
@@ -321,7 +322,7 @@ pzero (double x)
{
p = pR3; q = pS3;
}
- else if (ix >= 0x40000000)
+ else
{
p = pR2; q = pS2;
}
@@ -423,6 +424,7 @@ qzero (double x)
int32_t ix;
GET_HIGH_WORD (ix, x);
ix &= 0x7fffffff;
+ /* ix >= 0x40000000 for all calls to this function. */
if (ix >= 0x41b00000)
{
return -.125 / x;
@@ -439,7 +441,7 @@ qzero (double x)
{
p = qR3; q = qS3;
}
- else if (ix >= 0x40000000)
+ else
{
p = qR2; q = qS2;
}
@@ -305,6 +305,7 @@ pone (double x)
int32_t ix;
GET_HIGH_WORD (ix, x);
ix &= 0x7fffffff;
+ /* ix >= 0x40000000 for all calls to this function. */
if (ix >= 0x41b00000)
{
return one;
@@ -321,7 +322,7 @@ pone (double x)
{
p = pr3; q = ps3;
}
- else if (ix >= 0x40000000)
+ else
{
p = pr2; q = ps2;
}
@@ -424,6 +425,7 @@ qone (double x)
int32_t ix;
GET_HIGH_WORD (ix, x);
ix &= 0x7fffffff;
+ /* ix >= 0x40000000 for all calls to this function. */
if (ix >= 0x41b00000)
{
return .375 / x;
@@ -440,7 +442,7 @@ qone (double x)
{
p = qr3; q = qs3;
}
- else if (ix >= 0x40000000)
+ else
{
p = qr2; q = qs2;
}
@@ -228,10 +228,11 @@ pzerof(float x)
int32_t ix;
GET_FLOAT_WORD(ix,x);
ix &= 0x7fffffff;
+ /* ix >= 0x40000000 for all calls to this function. */
if(ix>=0x41000000) {p = pR8; q= pS8;}
else if(ix>=0x40f71c58){p = pR5; q= pS5;}
else if(ix>=0x4036db68){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]))));
@@ -324,10 +325,11 @@ qzerof(float x)
int32_t ix;
GET_FLOAT_WORD(ix,x);
ix &= 0x7fffffff;
+ /* ix >= 0x40000000 for all calls to this function. */
if(ix>=0x41000000) {p = qR8; q= qS8;}
else if(ix>=0x40f71c58){p = qR5; q= qS5;}
else if(ix>=0x4036db68){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])))));
@@ -230,10 +230,11 @@ ponef(float x)
int32_t ix;
GET_FLOAT_WORD(ix,x);
ix &= 0x7fffffff;
+ /* ix >= 0x40000000 for all calls to this function. */
if(ix>=0x41000000) {p = pr8; q= ps8;}
else if(ix>=0x40f71c58){p = pr5; q= ps5;}
else if(ix>=0x4036db68){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]))));
@@ -327,10 +328,11 @@ qonef(float x)
int32_t ix;
GET_FLOAT_WORD(ix,x);
ix &= 0x7fffffff;
+ /* ix >= 0x40000000 for all calls to this function. */
if(ix>=0x40200000) {p = qr8; q= qs8;}
else if(ix>=0x40f71c58){p = qr5; q= qs5;}
else if(ix>=0x4036db68){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])))));
@@ -356,6 +356,7 @@ pzero (long double x)
GET_LDOUBLE_WORDS (se, i0, i1, x);
ix = se & 0x7fff;
+ /* ix >= 0x4000 for all calls to this function. */
if (ix >= 0x4002)
{
p = pR8;
@@ -374,7 +375,7 @@ pzero (long double x)
p = pR3;
q = pS3;
}
- else if (ix >= 0x4000) /* x better be >= 2 */
+ else /* x >= 2 */
{
p = pR2;
q = pS2;
@@ -493,6 +494,7 @@ qzero (long double x)
GET_LDOUBLE_WORDS (se, i0, i1, x);
ix = se & 0x7fff;
+ /* ix >= 0x4000 for all calls to this function. */
if (ix >= 0x4002) /* x >= 8 */
{
p = qR8;
@@ -511,7 +513,7 @@ qzero (long double x)
p = qR3;
q = qS3;
}
- else if (ix >= 0x4000) /* x better be >= 2 */
+ else /* x >= 2 */
{
p = qR2;
q = qS2;
@@ -359,6 +359,7 @@ pone (long double x)
GET_LDOUBLE_WORDS (se, i0, i1, x);
ix = se & 0x7fff;
+ /* ix >= 0x4000 for all calls to this function. */
if (ix >= 0x4002) /* x >= 8 */
{
p = pr8;
@@ -377,7 +378,7 @@ pone (long double x)
p = pr3;
q = ps3;
}
- else if (ix >= 0x4000) /* x better be >= 2 */
+ else /* x >= 2 */
{
p = pr2;
q = ps2;
@@ -505,6 +506,7 @@ qone (long double x)
GET_LDOUBLE_WORDS (se, i0, i1, x);
ix = se & 0x7fff;
+ /* ix >= 0x4000 for all calls to this function. */
if (ix >= 0x4002) /* x >= 8 */
{
p = qr8;
@@ -523,7 +525,7 @@ qone (long double x)
p = qr3;
q = qs3;
}
- else if (ix >= 0x4000) /* x better be >= 2 */
+ else /* x >= 2 */
{
p = qr2;
q = qs2;