838     int16 aExp, shiftCount;  in float32_to_int32()  local
843     aExp = extractFloat32Exp( a );  in float32_to_int32()
845     if ( ( aExp == 0x7FF ) && aSig ) aSign = 0;  in float32_to_int32()
846     if ( aExp ) aSig |= 0x00800000;  in float32_to_int32()
847     shiftCount = 0xAF - aExp;  in float32_to_int32()
869     int16 aExp, shiftCount;  in float32_to_int32_round_to_zero()  local
874     aExp = extractFloat32Exp( a );  in float32_to_int32_round_to_zero()
876     shiftCount = aExp - 0x9E;  in float32_to_int32_round_to_zero()
880         if ( ! aSign || ( ( aExp == 0xFF ) && aSig ) ) return 0x7FFFFFFF;  in float32_to_int32_round_to_zero()
883     else if ( aExp <= 0x7E ) {  in float32_to_int32_round_to_zero()
884         if ( aExp | aSig ) float_raise( float_flag_inexact );  in float32_to_int32_round_to_zero()
907     int16 aExp;  in float32_to_float64()  local
911     aExp = extractFloat32Exp( a );  in float32_to_float64()
913     if ( aExp == 0xFF ) {  in float32_to_float64()
917     if ( aExp == 0 ) {  in float32_to_float64()
919         normalizeFloat32Subnormal( aSig, &aExp, &aSig );  in float32_to_float64()
920         --aExp;  in float32_to_float64()
922     return packFloat64( aSign, aExp + 0x380, ( (bits64) aSig )<<29 );  in float32_to_float64()
939     int16 aExp;  in float32_to_floatx80()  local
943     aExp = extractFloat32Exp( a );  in float32_to_floatx80()
945     if ( aExp == 0xFF ) {  in float32_to_floatx80()
949     if ( aExp == 0 ) {  in float32_to_floatx80()
951         normalizeFloat32Subnormal( aSig, &aExp, &aSig );  in float32_to_floatx80()
954     return packFloatx80( aSign, aExp + 0x3F80, ( (bits64) aSig )<<40 );  in float32_to_floatx80()
971     int16 aExp;  in float32_round_to_int()  local
976     aExp = extractFloat32Exp( a );  in float32_round_to_int()
977     if ( 0x96 <= aExp ) {  in float32_round_to_int()
978         if ( ( aExp == 0xFF ) && extractFloat32Frac( a ) ) {  in float32_round_to_int()
984     if ( aExp <= 0x7E ) {  in float32_round_to_int()
990             if ( ( aExp == 0x7E ) && extractFloat32Frac( a ) ) {  in float32_round_to_int()
1002     lastBitMask <<= 0x96 - aExp;  in float32_round_to_int()
1031     int16 aExp, bExp, zExp;  in addFloat32Sigs()  local
1036     aExp = extractFloat32Exp( a );  in addFloat32Sigs()
1039     expDiff = aExp - bExp;  in addFloat32Sigs()
1043         if ( aExp == 0xFF ) {  in addFloat32Sigs()
1054         zExp = aExp;  in addFloat32Sigs()
1061         if ( aExp == 0 ) {  in addFloat32Sigs()
1071         if ( aExp == 0xFF ) {  in addFloat32Sigs()
1075         if ( aExp == 0 ) return packFloat32( zSign, 0, ( aSig + bSig )>>6 );  in addFloat32Sigs()
1077         zExp = aExp;  in addFloat32Sigs()
1103     int16 aExp, bExp, zExp;  in subFloat32Sigs()  local
1108     aExp = extractFloat32Exp( a );  in subFloat32Sigs()
1111     expDiff = aExp - bExp;  in subFloat32Sigs()
1116     if ( aExp == 0xFF ) {  in subFloat32Sigs()
1121     if ( aExp == 0 ) {  in subFloat32Sigs()
1122         aExp = 1;  in subFloat32Sigs()
1133     if ( aExp == 0 ) {  in subFloat32Sigs()
1147     if ( aExp == 0xFF ) {  in subFloat32Sigs()
1161     zExp = aExp;  in subFloat32Sigs()
1222     int16 aExp, bExp, zExp;  in float32_mul()  local
1228     aExp = extractFloat32Exp( a );  in float32_mul()
1234     if ( aExp == 0xFF ) {  in float32_mul()
1246         if ( ( aExp | aSig ) == 0 ) {  in float32_mul()
1252     if ( aExp == 0 ) {  in float32_mul()
1254         normalizeFloat32Subnormal( aSig, &aExp, &aSig );  in float32_mul()
1260     zExp = aExp + bExp - 0x7F;  in float32_mul()
1283     int16 aExp, bExp, zExp;  in float32_div()  local
1287     aExp = extractFloat32Exp( a );  in float32_div()
1293     if ( aExp == 0xFF ) {  in float32_div()
1308             if ( ( aExp | aSig ) == 0 ) {  in float32_div()
1317     if ( aExp == 0 ) {  in float32_div()
1319         normalizeFloat32Subnormal( aSig, &aExp, &aSig );  in float32_div()
1321     zExp = aExp - bExp + 0x7D;  in float32_div()
1350     int16 aExp, bExp, expDiff;  in float32_rem()  local
1358     aExp = extractFloat32Exp( a );  in float32_rem()
1363     if ( aExp == 0xFF ) {  in float32_rem()
1381     if ( aExp == 0 ) {  in float32_rem()
1383         normalizeFloat32Subnormal( aSig, &aExp, &aSig );  in float32_rem()
1385     expDiff = aExp - bExp;  in float32_rem()
1453     int16 aExp, zExp;  in float32_sqrt()  local
1458     aExp = extractFloat32Exp( a );  in float32_sqrt()
1460     if ( aExp == 0xFF ) {  in float32_sqrt()
1467         if ( ( aExp | aSig ) == 0 ) return a;  in float32_sqrt()
1471     if ( aExp == 0 ) {  in float32_sqrt()
1473         normalizeFloat32Subnormal( aSig, &aExp, &aSig );  in float32_sqrt()
1475     zExp = ( ( aExp - 0x7F )>>1 ) + 0x7E;  in float32_sqrt()
1477     zSig = estimateSqrt32( aExp, aSig ) + 2;  in float32_sqrt()
1483             aSig >>= aExp & 1;  in float32_sqrt()
1655     int16 aExp, shiftCount;  in float64_to_int32()  local
1659     aExp = extractFloat64Exp( a );  in float64_to_int32()
1661     if ( ( aExp == 0x7FF ) && aSig ) aSign = 0;  in float64_to_int32()
1662     if ( aExp ) aSig |= LIT64( 0x0010000000000000 );  in float64_to_int32()
1663     shiftCount = 0x42C - aExp;  in float64_to_int32()
1683     int16 aExp, shiftCount;  in float64_to_int32_round_to_zero()  local
1688     aExp = extractFloat64Exp( a );  in float64_to_int32_round_to_zero()
1690     shiftCount = 0x433 - aExp;  in float64_to_int32_round_to_zero()
1692         if ( ( aExp == 0x7FF ) && aSig ) aSign = 0;  in float64_to_int32_round_to_zero()
1696         if ( aExp || aSig ) float_raise( float_flag_inexact );  in float64_to_int32_round_to_zero()
1730     int16 aExp, shiftCount;  in float64_to_uint32()  local
1734     aExp = extractFloat64Exp( a );  in float64_to_uint32()
1737     if ( aExp ) aSig |= LIT64( 0x0010000000000000 );  in float64_to_uint32()
1738     shiftCount = 0x42C - aExp;  in float64_to_uint32()
1756     int16 aExp, shiftCount;  in float64_to_uint32_round_to_zero()  local
1761     aExp = extractFloat64Exp( a );  in float64_to_uint32_round_to_zero()
1763     shiftCount = 0x433 - aExp;  in float64_to_uint32_round_to_zero()
1765         if ( ( aExp == 0x7FF ) && aSig ) aSign = 0;  in float64_to_uint32_round_to_zero()
1769         if ( aExp || aSig ) float_raise( float_flag_inexact );  in float64_to_uint32_round_to_zero()
1799     int16 aExp;  in float64_to_float32()  local
1804     aExp = extractFloat64Exp( a );  in float64_to_float32()
1806     if ( aExp == 0x7FF ) {  in float64_to_float32()
1812     if ( aExp || zSig ) {  in float64_to_float32()
1814         aExp -= 0x381;  in float64_to_float32()
1816     return roundAndPackFloat32( roundData, aSign, aExp, zSig );  in float64_to_float32()
1833     int16 aExp;  in float64_to_floatx80()  local
1837     aExp = extractFloat64Exp( a );  in float64_to_floatx80()
1839     if ( aExp == 0x7FF ) {  in float64_to_floatx80()
1843     if ( aExp == 0 ) {  in float64_to_floatx80()
1845         normalizeFloat64Subnormal( aSig, &aExp, &aSig );  in float64_to_floatx80()
1849             aSign, aExp + 0x3C00, ( aSig | LIT64( 0x0010000000000000 ) )<<11 );  in float64_to_floatx80()
1866     int16 aExp;  in float64_round_to_int()  local
1871     aExp = extractFloat64Exp( a );  in float64_round_to_int()
1872     if ( 0x433 <= aExp ) {  in float64_round_to_int()
1873         if ( ( aExp == 0x7FF ) && extractFloat64Frac( a ) ) {  in float64_round_to_int()
1878     if ( aExp <= 0x3FE ) {  in float64_round_to_int()
1884             if ( ( aExp == 0x3FE ) && extractFloat64Frac( a ) ) {  in float64_round_to_int()
1897     lastBitMask <<= 0x433 - aExp;  in float64_round_to_int()
1927     int16 aExp, bExp, zExp;  in addFloat64Sigs()  local
1932     aExp = extractFloat64Exp( a );  in addFloat64Sigs()
1935     expDiff = aExp - bExp;  in addFloat64Sigs()
1939         if ( aExp == 0x7FF ) {  in addFloat64Sigs()
1950         zExp = aExp;  in addFloat64Sigs()
1957         if ( aExp == 0 ) {  in addFloat64Sigs()
1967         if ( aExp == 0x7FF ) {  in addFloat64Sigs()
1971         if ( aExp == 0 ) return packFloat64( zSign, 0, ( aSig + bSig )>>9 );  in addFloat64Sigs()
1973         zExp = aExp;  in addFloat64Sigs()
1999     int16 aExp, bExp, zExp;  in subFloat64Sigs()  local
2004     aExp = extractFloat64Exp( a );  in subFloat64Sigs()
2007     expDiff = aExp - bExp;  in subFloat64Sigs()
2012     if ( aExp == 0x7FF ) {  in subFloat64Sigs()
2017     if ( aExp == 0 ) {  in subFloat64Sigs()
2018         aExp = 1;  in subFloat64Sigs()
2029     if ( aExp == 0 ) {  in subFloat64Sigs()
2043     if ( aExp == 0x7FF ) {  in subFloat64Sigs()
2057     zExp = aExp;  in subFloat64Sigs()
2118     int16 aExp, bExp, zExp;  in float64_mul()  local
2122     aExp = extractFloat64Exp( a );  in float64_mul()
2128     if ( aExp == 0x7FF ) {  in float64_mul()
2140         if ( ( aExp | aSig ) == 0 ) {  in float64_mul()
2146     if ( aExp == 0 ) {  in float64_mul()
2148         normalizeFloat64Subnormal( aSig, &aExp, &aSig );  in float64_mul()
2154     zExp = aExp + bExp - 0x3FF;  in float64_mul()
2177     int16 aExp, bExp, zExp;  in float64_div()  local
2183     aExp = extractFloat64Exp( a );  in float64_div()
2189     if ( aExp == 0x7FF ) {  in float64_div()
2204             if ( ( aExp | aSig ) == 0 ) {  in float64_div()
2213     if ( aExp == 0 ) {  in float64_div()
2215         normalizeFloat64Subnormal( aSig, &aExp, &aSig );  in float64_div()
2217     zExp = aExp - bExp + 0x3FD;  in float64_div()
2248     int16 aExp, bExp, expDiff;  in float64_rem()  local
2254     aExp = extractFloat64Exp( a );  in float64_rem()
2259     if ( aExp == 0x7FF ) {  in float64_rem()
2277     if ( aExp == 0 ) {  in float64_rem()
2279         normalizeFloat64Subnormal( aSig, &aExp, &aSig );  in float64_rem()
2281     expDiff = aExp - bExp;  in float64_rem()
2334     int16 aExp, zExp;  in float64_sqrt()  local
2340     aExp = extractFloat64Exp( a );  in float64_sqrt()
2342     if ( aExp == 0x7FF ) {  in float64_sqrt()
2349         if ( ( aExp | aSig ) == 0 ) return a;  in float64_sqrt()
2353     if ( aExp == 0 ) {  in float64_sqrt()
2355         normalizeFloat64Subnormal( aSig, &aExp, &aSig );  in float64_sqrt()
2357     zExp = ( ( aExp - 0x3FF )>>1 ) + 0x3FE;  in float64_sqrt()
2359     zSig = estimateSqrt32( aExp, aSig>>21 );  in float64_sqrt()
2361     aSig <<= 9 - ( aExp & 1 );  in float64_sqrt()
2544     int32 aExp, shiftCount;  in floatx80_to_int32()  local
2548     aExp = extractFloatx80Exp( a );  in floatx80_to_int32()
2550     if ( ( aExp == 0x7FFF ) && (bits64) ( aSig<<1 ) ) aSign = 0;  in floatx80_to_int32()
2551     shiftCount = 0x4037 - aExp;  in floatx80_to_int32()
2572     int32 aExp, shiftCount;  in floatx80_to_int32_round_to_zero()  local
2577     aExp = extractFloatx80Exp( a );  in floatx80_to_int32_round_to_zero()
2579     shiftCount = 0x403E - aExp;  in floatx80_to_int32_round_to_zero()
2581         if ( ( aExp == 0x7FFF ) && (bits64) ( aSig<<1 ) ) aSign = 0;  in floatx80_to_int32_round_to_zero()
2585         if ( aExp || aSig ) float_raise( float_flag_inexact );  in floatx80_to_int32_round_to_zero()
2615     int32 aExp;  in floatx80_to_float32()  local
2619     aExp = extractFloatx80Exp( a );  in floatx80_to_float32()
2621     if ( aExp == 0x7FFF ) {  in floatx80_to_float32()
2628     if ( aExp || aSig ) aExp -= 0x3F81;  in floatx80_to_float32()
2629     return roundAndPackFloat32( roundData, aSign, aExp, aSig );  in floatx80_to_float32()
2644     int32 aExp;  in floatx80_to_float64()  local
2648     aExp = extractFloatx80Exp( a );  in floatx80_to_float64()
2650     if ( aExp == 0x7FFF ) {  in floatx80_to_float64()
2657     if ( aExp || aSig ) aExp -= 0x3C01;  in floatx80_to_float64()
2658     return roundAndPackFloat64( roundData, aSign, aExp, zSig );  in floatx80_to_float64()
2673     int32 aExp;  in floatx80_round_to_int()  local
2678     aExp = extractFloatx80Exp( a );  in floatx80_round_to_int()
2679     if ( 0x403E <= aExp ) {  in floatx80_round_to_int()
2680         if ( ( aExp == 0x7FFF ) && (bits64) ( extractFloatx80Frac( a )<<1 ) ) {  in floatx80_round_to_int()
2685     if ( aExp <= 0x3FFE ) {  in floatx80_round_to_int()
2686         if (    ( aExp == 0 )  in floatx80_round_to_int()
2694             if ( ( aExp == 0x3FFE ) && (bits64) ( extractFloatx80Frac( a )<<1 )  in floatx80_round_to_int()
2713     lastBitMask <<= 0x403E - aExp;  in floatx80_round_to_int()
2747     int32 aExp, bExp, zExp;  in addFloatx80Sigs()  local
2752     aExp = extractFloatx80Exp( a );  in addFloatx80Sigs()
2755     expDiff = aExp - bExp;  in addFloatx80Sigs()
2757         if ( aExp == 0x7FFF ) {  in addFloatx80Sigs()
2763         zExp = aExp;  in addFloatx80Sigs()
2770         if ( aExp == 0 ) ++expDiff;  in addFloatx80Sigs()
2775         if ( aExp == 0x7FFF ) {  in addFloatx80Sigs()
2783         if ( aExp == 0 ) {  in addFloatx80Sigs()
2787         zExp = aExp;  in addFloatx80Sigs()
2816     int32 aExp, bExp, zExp;  in subFloatx80Sigs()  local
2822     aExp = extractFloatx80Exp( a );  in subFloatx80Sigs()
2825     expDiff = aExp - bExp;  in subFloatx80Sigs()
2828     if ( aExp == 0x7FFF ) {  in subFloatx80Sigs()
2838     if ( aExp == 0 ) {  in subFloatx80Sigs()
2839         aExp = 1;  in subFloatx80Sigs()
2851     if ( aExp == 0 ) ++expDiff;  in subFloatx80Sigs()
2859     if ( aExp == 0x7FFF ) {  in subFloatx80Sigs()
2867     zExp = aExp;  in subFloatx80Sigs()
2929     int32 aExp, bExp, zExp;  in floatx80_mul()  local
2934     aExp = extractFloatx80Exp( a );  in floatx80_mul()
2940     if ( aExp == 0x7FFF ) {  in floatx80_mul()
2950         if ( ( aExp | aSig ) == 0 ) {  in floatx80_mul()
2960     if ( aExp == 0 ) {  in floatx80_mul()
2962         normalizeFloatx80Subnormal( aSig, &aExp, &aSig );  in floatx80_mul()
2968     zExp = aExp + bExp - 0x3FFE;  in floatx80_mul()
2990     int32 aExp, bExp, zExp;  in floatx80_div()  local
2996     aExp = extractFloatx80Exp( a );  in floatx80_div()
3002     if ( aExp == 0x7FFF ) {  in floatx80_div()
3016             if ( ( aExp | aSig ) == 0 ) {  in floatx80_div()
3029     if ( aExp == 0 ) {  in floatx80_div()
3031         normalizeFloatx80Subnormal( aSig, &aExp, &aSig );  in floatx80_div()
3033     zExp = aExp - bExp + 0x3FFE;  in floatx80_div()
3072     int32 aExp, bExp, expDiff;  in floatx80_rem()  local
3078     aExp = extractFloatx80Exp( a );  in floatx80_rem()
3083     if ( aExp == 0x7FFF ) {  in floatx80_rem()
3105     if ( aExp == 0 ) {  in floatx80_rem()
3107         normalizeFloatx80Subnormal( aSig0, &aExp, &aSig0 );  in floatx80_rem()
3111     expDiff = aExp - bExp;  in floatx80_rem()
3172     int32 aExp, zExp;  in floatx80_sqrt()  local
3179     aExp = extractFloatx80Exp( a );  in floatx80_sqrt()
3181     if ( aExp == 0x7FFF ) {  in floatx80_sqrt()
3187         if ( ( aExp | aSig0 ) == 0 ) return a;  in floatx80_sqrt()
3195     if ( aExp == 0 ) {  in floatx80_sqrt()
3197         normalizeFloatx80Subnormal( aSig0, &aExp, &aSig0 );  in floatx80_sqrt()
3199     zExp = ( ( aExp - 0x3FFF )>>1 ) + 0x3FFF;  in floatx80_sqrt()
3200     zSig0 = estimateSqrt32( aExp, aSig0>>32 );  in floatx80_sqrt()
3203     shift128Right( aSig0, 0, ( aExp & 1 ) + 2, &aSig0, &aSig1 );  in floatx80_sqrt()