small cleanup patch for e_pow.c

Steve Kargl sgk at troutmask.apl.washington.edu
Sun May 10 00:38:35 UTC 2015


In reading, e_pow.c I found a small piece of code that
can be remove.  Anyone object?

Index: src/e_pow.c
===================================================================
--- src/e_pow.c	(revision 1603)
+++ src/e_pow.c	(working copy)
@@ -187,10 +187,6 @@ __ieee754_pow(double x, double y)
 
     /* |y| is huge */
 	if(iy>0x41e00000) { /* if |y| > 2**31 */
-	    if(iy>0x43f00000){	/* if |y| > 2**64, must o/uflow */
-		if(ix<=0x3fefffff) return (hy<0)? huge*huge:tiny*tiny;
-		if(ix>=0x3ff00000) return (hy>0)? huge*huge:tiny*tiny;
-	    }
 	/* over/underflow if x is not close to one */
 	    if(ix<0x3fefffff) return (hy<0)? s*huge*huge:s*tiny*tiny;
 	    if(ix>0x3ff00000) return (hy>0)? s*huge*huge:s*tiny*tiny;
-- 
Steve


More information about the freebsd-numerics mailing list