softfloat/bits32/softfloat.c

592     int16 aExp, shiftCount;
598     aExp = extractFloat32Exp( a );
600     shiftCount = aExp - 0x96;
602         if ( 0x9E <= aExp ) {
605                 if ( ! aSign || ( ( aExp == 0xFF ) && aSig ) ) {
615         if ( aExp < 0x7E ) {
616             aSigExtra = aExp | aSig;
663     int16 aExp, shiftCount;
668     aExp = extractFloat32Exp( a );
670     shiftCount = aExp - 0x9E;
674             if ( ! aSign || ( ( aExp == 0xFF ) && aSig ) ) return 0x7FFFFFFF;
678     else if ( aExp <= 0x7E ) {
679         if ( aExp | aSig ) set_float_exception_inexact_flag();
703     int16 aExp;
707     aExp = extractFloat32Exp( a );
709     if ( aExp == 0xFF ) {
713     if ( aExp == 0 ) {
715         normalizeFloat32Subnormal( aSig, &aExp, &aSig );
716         --aExp;
719     return packFloat64( aSign, aExp + 0x380, zSig0, zSig1 );
735     int16 aExp;
740     aExp = extractFloat32Exp( a );
741     if ( 0x96 <= aExp ) {
742         if ( ( aExp == 0xFF ) && extractFloat32Frac( a ) ) {
747     if ( aExp <= 0x7E ) {
753             if ( ( aExp == 0x7E ) && extractFloat32Frac( a ) ) {
767     lastBitMask <<= 0x96 - aExp;
798     int16 aExp, bExp, zExp;
803     aExp = extractFloat32Exp( a );
806     expDiff = aExp - bExp;
810         if ( aExp == 0xFF ) {
821         zExp = aExp;
828         if ( aExp == 0 ) {
838         if ( aExp == 0xFF ) {
842         if ( aExp == 0 ) return packFloat32( zSign, 0, ( aSig + bSig )>>6 );
844         zExp = aExp;
870     int16 aExp, bExp, zExp;
875     aExp = extractFloat32Exp( a );
878     expDiff = aExp - bExp;
883     if ( aExp == 0xFF ) {
888     if ( aExp == 0 ) {
889         aExp = 1;
900     if ( aExp == 0 ) {
914     if ( aExp == 0xFF ) {
928     zExp = aExp;
989     int16 aExp, bExp, zExp;
993     aExp = extractFloat32Exp( a );
999     if ( aExp == 0xFF ) {
1011         if ( ( aExp | aSig ) == 0 ) {
1017     if ( aExp == 0 ) {
1019         normalizeFloat32Subnormal( aSig, &aExp, &aSig );
1025     zExp = aExp + bExp - 0x7F;
1048     int16 aExp, bExp, zExp;
1052     aExp = extractFloat32Exp( a );
1058     if ( aExp == 0xFF ) {
1073             if ( ( aExp | aSig ) == 0 ) {
1082     if ( aExp == 0 ) {
1084         normalizeFloat32Subnormal( aSig, &aExp, &aSig );
1086     zExp = aExp - bExp + 0x7D;
1118     int16 aExp, bExp, expDiff;
1123     aExp = extractFloat32Exp( a );
1128     if ( aExp == 0xFF ) {
1146     if ( aExp == 0 ) {
1148         normalizeFloat32Subnormal( aSig, &aExp, &aSig );
1150     expDiff = aExp - bExp;
1205     int16 aExp, zExp;
1209     aExp = extractFloat32Exp( a );
1211     if ( aExp == 0xFF ) {
1218         if ( ( aExp | aSig ) == 0 ) return a;
1222     if ( aExp == 0 ) {
1224         normalizeFloat32Subnormal( aSig, &aExp, &aSig );
1226     zExp = ( ( aExp - 0x7F )>>1 ) + 0x7E;
1228     zSig = estimateSqrt32( aExp, aSig ) + 2;
1235             aSig >>= aExp & 1;
1358     int16 aExp, bExp;
1418     int16 aExp, shiftCount;
1425     aExp = extractFloat64Exp( a );
1427     shiftCount = aExp - 0x413;
1429         if ( 0x41E < aExp ) {
1430             if ( ( aExp == 0x7FF ) && ( aSig0 | aSig1 ) ) aSign = 0;
1439         if ( aExp < 0x3FE ) {
1440             aSigExtra = aExp | aSig0 | aSig1;
1492     int16 aExp, shiftCount;
1498     aExp = extractFloat64Exp( a );
1500     shiftCount = aExp - 0x413;
1502         if ( 0x41E < aExp ) {
1503             if ( ( aExp == 0x7FF ) && ( aSig0 | aSig1 ) ) aSign = 0;
1510         if ( aExp < 0x3FF ) {
1511             if ( aExp | aSig0 | aSig1 ) {
1542     int16 aExp;
1548     aExp = extractFloat64Exp( a );
1550     if ( aExp == 0x7FF ) {
1557     if ( aExp ) zSig |= 0x40000000;
1558     return roundAndPackFloat32( aSign, aExp - 0x381, zSig );
1574     int16 aExp;
1579     aExp = extractFloat64Exp( a );
1580     if ( 0x413 <= aExp ) {
1581         if ( 0x433 <= aExp ) {
1582             if (    ( aExp == 0x7FF )
1589         lastBitMask = ( lastBitMask<<( 0x432 - aExp ) )<<1;
1614         if ( aExp <= 0x3FE ) {
1620                 if (    ( aExp == 0x3FE )
1638         lastBitMask <<= 0x413 - aExp;
1677     int16 aExp, bExp, zExp;
1683     aExp = extractFloat64Exp( a );
1687     expDiff = aExp - bExp;
1689         if ( aExp == 0x7FF ) {
1701         zExp = aExp;
1708         if ( aExp == 0 ) {
1719         if ( aExp == 0x7FF ) {
1726         if ( aExp == 0 ) return packFloat64( zSign, 0, zSig0, zSig1 );
1729         zExp = aExp;
1755     int16 aExp, bExp, zExp;
1761     aExp = extractFloat64Exp( a );
1765     expDiff = aExp - bExp;
1770     if ( aExp == 0x7FF ) {
1777     if ( aExp == 0 ) {
1778         aExp = 1;
1791     if ( aExp == 0 ) {
1805     if ( aExp == 0x7FF ) {
1819     zExp = aExp;
1880     int16 aExp, bExp, zExp;
1885     aExp = extractFloat64Exp( a );
1892     if ( aExp == 0x7FF ) {
1902         if ( ( aExp | aSig0 | aSig1 ) == 0 ) {
1909     if ( aExp == 0 ) {
1911         normalizeFloat64Subnormal( aSig0, aSig1, &aExp, &aSig0, &aSig1 );
1917     zExp = aExp + bExp - 0x400;
1942     int16 aExp, bExp, zExp;
1948     aExp = extractFloat64Exp( a );
1955     if ( aExp == 0x7FF ) {
1969             if ( ( aExp | aSig0 | aSig1 ) == 0 ) {
1979     if ( aExp == 0 ) {
1981         normalizeFloat64Subnormal( aSig0, aSig1, &aExp, &aSig0, &aSig1 );
1983     zExp = aExp - bExp + 0x3FD;
2023     int16 aExp, bExp, expDiff;
2031     aExp = extractFloat64Exp( a );
2037     if ( aExp == 0x7FF ) {
2056     if ( aExp == 0 ) {
2058         normalizeFloat64Subnormal( aSig0, aSig1, &aExp, &aSig0, &aSig1 );
2060     expDiff = aExp - bExp;
2128     int16 aExp, zExp;
2135     aExp = extractFloat64Exp( a );
2137     if ( aExp == 0x7FF ) {
2143         if ( ( aExp | aSig0 | aSig1 ) == 0 ) return a;
2148     if ( aExp == 0 ) {
2150         normalizeFloat64Subnormal( aSig0, aSig1, &aExp, &aSig0, &aSig1 );
2152     zExp = ( ( aExp - 0x3FF )>>1 ) + 0x3FE;
2155     zSig0 = ( estimateSqrt32( aExp, term0 )>>1 ) + 1;
2158     shortShift64Left( aSig0, aSig1, 9 - ( aExp & 1 ), &aSig0, &aSig1 );