gonzui

t2ex

bsd_source/lib/libc/src_bsd/math/s_expm1.c - 8.9KB - 230 lines

200:                 GET_HIGH_WORD(high,y);
201:                 SET_HIGH_WORD(y,high+(k<<20));        /* add k to y's exponent */
202:                 return y-one;
206:                 u_int32_t high;
207:                 SET_HIGH_WORD(t,0x3ff00000 - (0x200000>>k));  /* t=1-2^-k */
208:                        y = t-(e-x);
209:                 GET_HIGH_WORD(high,y);
210:                 SET_HIGH_WORD(y,high+(k<<20));        /* add k to y's exponent */
211:            } else {
212:                 u_int32_t high;
213:                 SET_HIGH_WORD(t,((0x3ff-k)<<20));     /* 2^-k */
214:                        y = x-(e+t);
216:                 GET_HIGH_WORD(high,y);
217:                 SET_HIGH_WORD(y,high+(k<<20));        /* add k to y's exponent */
218:             }