gonzui

t2ex

bsd_source/lib/libc/src_bsd/math/s_cbrtf.c - 1.9KB - 67 lines

47: 
48:         SET_FLOAT_WORD(x,hx);  /* x <- |x| */
49:     /* rough cbrt to 5 bits */
50:         if(hx<0x00800000)              /* subnormal number */
51:           {SET_FLOAT_WORD(t,0x4b800000); /* set t= 2**24 */
52:            t*=x; GET_FLOAT_WORD(high,t); SET_FLOAT_WORD(t,high/3+B2);
53:           }
54:         else
55:           SET_FLOAT_WORD(t,hx/3+B1);
56: 
64:         GET_FLOAT_WORD(high,t);
65:         SET_FLOAT_WORD(t,high|sign);
66:         return(t);