gonzui

t2ex

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

50: 
51:         SET_HIGH_WORD(x,hx);   /* x <- |x| */
52:     /* rough cbrt to 5 bits */
53:         if(hx<0x00100000)              /* subnormal number */
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:           }
57:         else
58:           SET_HIGH_WORD(t,hx/3+B1);
59: 
79:         GET_HIGH_WORD(high,t);
80:         SET_HIGH_WORD(t,high|sign);
81:         return(t);