softfloat/bits32/softfloat.c

143 by the denormalized significand `aSig'.  The normalized exponent and
149  normalizeFloat32Subnormal( bits32 aSig, int16 *zExpPtr, bits32 *zSigPtr )
153     shiftCount = countLeadingZeros32( aSig ) - 8;
154     *zSigPtr = aSig<<shiftCount;
585     bits32 aSig, aSigExtra;
589     aSig = extractFloat32Frac( a );
597                 if ( ! aSign || ( ( aExp == 0xFF ) && aSig ) ) {
603         z = ( aSig | 0x00800000 )<<shiftCount;
608             aSigExtra = aExp | aSig;
612             aSig |= 0x00800000;
613             aSigExtra = aSig<<( shiftCount & 31 );
614             z = aSig>>( - shiftCount );
656     bits32 aSig;
659     aSig = extractFloat32Frac( a );
666             if ( ! aSign || ( ( aExp == 0xFF ) && aSig ) ) return 0x7FFFFFFF;
671         if ( aExp | aSig ) float_exception_flags |= float_flag_inexact;
674     aSig = ( aSig | 0x00800000 )<<8;
675     z = aSig>>( - shiftCount );
676     if ( (bits32) ( aSig<<( shiftCount & 31 ) ) ) {
696     bits32 aSig, zSig0, zSig1;
698     aSig = extractFloat32Frac( a );
702         if ( aSig ) return commonNaNToFloat64( float32ToCommonNaN( a ) );
706         if ( aSig == 0 ) return packFloat64( aSign, 0, 0, 0 );
707         normalizeFloat32Subnormal( aSig, &aExp, &aSig );
710     shift64Right( aSig, 0, 3, &zSig0, &zSig1 );
791     bits32 aSig, bSig, zSig;
794     aSig = extractFloat32Frac( a );
799     aSig <<= 6;
803             if ( aSig ) return propagateFloat32NaN( a, b );
824             aSig |= 0x20000000;
826         shift32RightJamming( aSig, - expDiff, &aSig );
831             if ( aSig | bSig ) return propagateFloat32NaN( a, b );
834         if ( aExp == 0 ) return packFloat32( zSign, 0, ( aSig + bSig )>>6 );
835         zSig = 0x40000000 + aSig + bSig;
839     aSig |= 0x20000000;
840     zSig = ( aSig + bSig )<<1;
843         zSig = aSig + bSig;
863     bits32 aSig, bSig, zSig;
866     aSig = extractFloat32Frac( a );
871     aSig <<= 7;
876         if ( aSig | bSig ) return propagateFloat32NaN( a, b );
884     if ( bSig < aSig ) goto aBigger;
885     if ( aSig < bSig ) goto bBigger;
896         aSig |= 0x40000000;
898     shift32RightJamming( aSig, - expDiff, &aSig );
901     zSig = bSig - aSig;
907         if ( aSig ) return propagateFloat32NaN( a, b );
917     aSig |= 0x40000000;
919     zSig = aSig - bSig;
982     bits32 aSig, bSig, zSig0, zSig1;
984     aSig = extractFloat32Frac( a );
992         if ( aSig || ( ( bExp == 0xFF ) && bSig ) ) {
1003         if ( ( aExp | aSig ) == 0 ) {
1010         if ( aSig == 0 ) return packFloat32( zSign, 0, 0 );
1011         normalizeFloat32Subnormal( aSig, &aExp, &aSig );
1018     aSig = ( aSig | 0x00800000 )<<7;
1020     mul32To64( aSig, bSig, &zSig0, &zSig1 );
1041     bits32 aSig, bSig, zSig, rem0, rem1, term0, term1;
1043     aSig = extractFloat32Frac( a );
1051         if ( aSig ) return propagateFloat32NaN( a, b );
1065             if ( ( aExp | aSig ) == 0 ) {
1075         if ( aSig == 0 ) return packFloat32( zSign, 0, 0 );
1076         normalizeFloat32Subnormal( aSig, &aExp, &aSig );
1079     aSig = ( aSig | 0x00800000 )<<7;
1081     if ( bSig <= ( aSig + aSig ) ) {
1082         aSig >>= 1;
1085     zSig = estimateDiv64To32( aSig, 0, bSig );
1088         sub64( aSig, 0, term0, term1, &rem0, &rem1 );
1111     bits32 aSig, bSig, q, allZero, alternateASig;
1114     aSig = extractFloat32Frac( a );
1121         if ( aSig || ( ( bExp == 0xFF ) && bSig ) ) {
1139         if ( aSig == 0 ) return a;
1140         normalizeFloat32Subnormal( aSig, &aExp, &aSig );
1143     aSig = ( aSig | 0x00800000 )<<8;
1147         aSig >>= 1;
1149     q = ( bSig <= aSig );
1150     if ( q ) aSig -= bSig;
1153         q = estimateDiv64To32( aSig, 0, bSig );
1155         aSig = - ( ( bSig>>2 ) * q );
1160         q = estimateDiv64To32( aSig, 0, bSig );
1164         aSig = ( ( aSig>>1 )<<( expDiff - 1 ) ) - bSig * q;
1167         aSig >>= 2;
1171         alternateASig = aSig;
1173         aSig -= bSig;
1174     } while ( 0 <= (sbits32) aSig );
1175     sigMean = aSig + alternateASig;
1177         aSig = alternateASig;
1179     zSign = ( (sbits32) aSig < 0 );
1180     if ( zSign ) aSig = - aSig;
1181     return normalizeRoundAndPackFloat32( aSign ^ zSign, bExp, aSig );
1198     bits32 aSig, zSig, rem0, rem1, term0, term1;
1200     aSig = extractFloat32Frac( a );
1204         if ( aSig ) return propagateFloat32NaN( a, 0 );
1210         if ( ( aExp | aSig ) == 0 ) return a;
1215         if ( aSig == 0 ) return 0;
1216         normalizeFloat32Subnormal( aSig, &aExp, &aSig );
1219     aSig = ( aSig | 0x00800000 )<<8;
1220     zSig = estimateSqrt32( aExp, aSig ) + 2;
1227             aSig >>= aExp & 1;
1229             sub64( aSig, 0, term0, term1, &rem0, &rem1 );