gonzui

t2ex

bsd_source/lib/libc/src_bsd/math/s_cbrt.c - 2.6KB - 91 lines

42: 
43:         GET_HIGH_WORD(hx,x);
44:         sign=hx&0x80000000;            /* sign= sign(x) */
54:           {SET_HIGH_WORD(t,0x43500000);        /* set t= 2**54 */
55:            t*=x; GET_HIGH_WORD(high,t); SET_HIGH_WORD(t,high/3+B2);
56:           }
66:     /* chopped to 20 bits and make it larger than cbrt(x) */ 
67:         GET_HIGH_WORD(high,t);
68:         INSERT_WORDS(t,high+0x00000001,0);
78:     /* retore the sign bit */
79:         GET_HIGH_WORD(high,t);
80:         SET_HIGH_WORD(t,high|sign);