| e_pow.c (22993) | e_pow.c (23579) |
|---|---|
| 1/* @(#)e_pow.c 5.1 93/09/24 */ 2/* 3 * ==================================================== 4 * Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved. 5 * 6 * Developed at SunPro, a Sun Microsystems, Inc. business. 7 * Permission to use, copy, modify, and distribute this 8 * software is freely granted, provided that this notice 9 * is preserved. 10 * ==================================================== 11 */ 12 13#ifndef lint | 1/* @(#)e_pow.c 5.1 93/09/24 */ 2/* 3 * ==================================================== 4 * Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved. 5 * 6 * Developed at SunPro, a Sun Microsystems, Inc. business. 7 * Permission to use, copy, modify, and distribute this 8 * software is freely granted, provided that this notice 9 * is preserved. 10 * ==================================================== 11 */ 12 13#ifndef lint |
| 14static char rcsid[] = "$Id$"; | 14static char rcsid[] = "$Id: e_pow.c,v 1.4 1997/02/22 15:10:29 peter Exp $"; |
| 15#endif 16 17/* __ieee754_pow(x,y) return x**y 18 * 19 * n 20 * Method: Let x = 2 * (1+f) 21 * 1. Compute and return log2(x) in two pieces: 22 * log2(x) = w1 + w2, --- 132 unchanged lines hidden (view full) --- 155 return (hy<0)?-y: zero; 156 } 157 if(iy==0x3ff00000) { /* y is +-1 */ 158 if(hy<0) return one/x; else return x; 159 } 160 if(hy==0x40000000) return x*x; /* y is 2 */ 161 if(hy==0x3fe00000) { /* y is 0.5 */ 162 if(hx>=0) /* x >= +0 */ | 15#endif 16 17/* __ieee754_pow(x,y) return x**y 18 * 19 * n 20 * Method: Let x = 2 * (1+f) 21 * 1. Compute and return log2(x) in two pieces: 22 * log2(x) = w1 + w2, --- 132 unchanged lines hidden (view full) --- 155 return (hy<0)?-y: zero; 156 } 157 if(iy==0x3ff00000) { /* y is +-1 */ 158 if(hy<0) return one/x; else return x; 159 } 160 if(hy==0x40000000) return x*x; /* y is 2 */ 161 if(hy==0x3fe00000) { /* y is 0.5 */ 162 if(hx>=0) /* x >= +0 */ |
| 163 return sqrt(x); | 163 return __ieee754_sqrt(x); |
| 164 } 165 } 166 167 ax = fabs(x); 168 /* special value of x */ 169 if(lx==0) { 170 if(ix==0x7ff00000||ix==0||ix==0x3ff00000){ 171 z = ax; /*x is +-0,+-inf,+-1*/ --- 141 unchanged lines hidden --- | 164 } 165 } 166 167 ax = fabs(x); 168 /* special value of x */ 169 if(lx==0) { 170 if(ix==0x7ff00000||ix==0||ix==0x3ff00000){ 171 z = ax; /*x is +-0,+-inf,+-1*/ --- 141 unchanged lines hidden --- |