Lines Matching +full:39 +full:x
15 /* ld80 version of __ieee754_rem_pio2l(x,y)
17 * return the remainder of x rem pi/2 in y[0]+y[1]
31 * pio2_1: first 39 bits of pi/2
33 * pio2_2: second 39 bits of pi/2
35 * pio2_3: third 39 bits of pi/2
70 __ieee754_rem_pio2l(long double x, long double *y)
78 u.e = x;
82 /* |x| ~< 2^25*(pi/2), medium size */
83 fn = rnintl(x*invpio2);
85 r = x-fn*pio2_1;
119 if(ex==0x7fff) { /* x is inf or NaN */
120 y[0]=y[1]=x-x; return 0;
122 /* set z = scalbn(|x|,ilogb(x)-23) */
123 u1.e = x;
124 e0 = ex - BIAS - 23; /* e0 = ilogb(|x|)-23; */