- GET_FLOAT_WORD(ha,x);
- ha &= 0x7fffffff;
- GET_FLOAT_WORD(hb,y);
- hb &= 0x7fffffff;
- if (hb > ha) {
- a = y;
- b = x;
- j=ha; ha=hb; hb=j;
- } else {
- a = x;
- b = y;
- }
- a = fabsf(a);
- b = fabsf(b);
- if (ha - hb > 0xf000000) /* x/y > 2**30 */
- return a+b;
- k = 0;
- if (ha > 0x58800000) { /* a > 2**50 */
- if(ha >= 0x7f800000) { /* Inf or NaN */
- /* Use original arg order iff result is NaN; quieten sNaNs. */
- w = fabsf(x+0.0f) - fabsf(y+0.0f);
- if (ha == 0x7f800000) w = a;
- if (hb == 0x7f800000) w = b;
- return w;
- }
- /* scale a and b by 2**-68 */
- ha -= 0x22000000; hb -= 0x22000000; k += 68;
- SET_FLOAT_WORD(a, ha);
- SET_FLOAT_WORD(b, hb);
- }
- if (hb < 0x26800000) { /* b < 2**-50 */
- if (hb <= 0x007fffff) { /* subnormal b or 0 */
- if (hb == 0)
- return a;
- SET_FLOAT_WORD(t1, 0x7e800000); /* t1 = 2^126 */
- b *= t1;
- a *= t1;
- k -= 126;
- } else { /* scale a and b by 2^68 */
- ha += 0x22000000; /* a *= 2^68 */
- hb += 0x22000000; /* b *= 2^68 */
- k -= 68;
- SET_FLOAT_WORD(a, ha);
- SET_FLOAT_WORD(b, hb);
- }