Lines Matching refs:x
105 __ieee754_exp(double x) /* default IEEE double exp */ in __ieee754_exp() argument
111 GET_HIGH_WORD(hx,x); in __ieee754_exp()
119 GET_LOW_WORD(lx,x); in __ieee754_exp()
121 return x+x; /* NaN */ in __ieee754_exp()
122 else return (xsb==0)? x:0.0; /* exp(+-inf)={inf,0} */ in __ieee754_exp()
124 if(x > o_threshold) return huge*huge; /* overflow */ in __ieee754_exp()
125 if(x < u_threshold) return twom1000*twom1000; /* underflow */ in __ieee754_exp()
131 hi = x-ln2HI[xsb]; lo=ln2LO[xsb]; k = 1-xsb-xsb; in __ieee754_exp()
133 k = (int)(invln2*x+halF[xsb]); in __ieee754_exp()
135 hi = x - t*ln2HI[0]; /* t*ln2HI is exact here */ in __ieee754_exp()
138 STRICT_ASSIGN(double, x, hi - lo); in __ieee754_exp()
141 if(huge+x>one) return one+x;/* trigger inexact */ in __ieee754_exp()
146 t = x*x; in __ieee754_exp()
151 c = x - t*(P1+t*(P2+t*(P3+t*(P4+t*P5)))); in __ieee754_exp()
152 if(k==0) return one-((x*c)/(c-2.0)-x); in __ieee754_exp()
153 else y = one-((lo-(x*c)/(2.0-c))-hi); in __ieee754_exp()