softfloat/bits32/softfloat.c

791     bits32 aSig, bSig, zSig;
796     bSig = extractFloat32Frac( b );
800     bSig <<= 6;
810             bSig |= 0x20000000;
812         shift32RightJamming( bSig, expDiff, &bSig );
817             if ( bSig ) return propagateFloat32NaN( a, b );
831             if ( aSig | bSig ) return propagateFloat32NaN( a, b );
834         if ( aExp == 0 ) return packFloat32( zSign, 0, ( aSig + bSig )>>6 );
835         zSig = 0x40000000 + aSig + bSig;
840     zSig = ( aSig + bSig )<<1;
843         zSig = aSig + bSig;
863     bits32 aSig, bSig, zSig;
868     bSig = extractFloat32Frac( b );
872     bSig <<= 7;
876         if ( aSig | bSig ) return propagateFloat32NaN( a, b );
884     if ( bSig < aSig ) goto aBigger;
885     if ( aSig < bSig ) goto bBigger;
889         if ( bSig ) return propagateFloat32NaN( a, b );
899     bSig |= 0x40000000;
901     zSig = bSig - aSig;
914         bSig |= 0x40000000;
916     shift32RightJamming( bSig, expDiff, &bSig );
919     zSig = aSig - bSig;
982     bits32 aSig, bSig, zSig0, zSig1;
987     bSig = extractFloat32Frac( b );
992         if ( aSig || ( ( bExp == 0xFF ) && bSig ) ) {
995         if ( ( bExp | bSig ) == 0 ) {
1002         if ( bSig ) return propagateFloat32NaN( a, b );
1014         if ( bSig == 0 ) return packFloat32( zSign, 0, 0 );
1015         normalizeFloat32Subnormal( bSig, &bExp, &bSig );
1019     bSig = ( bSig | 0x00800000 )<<8;
1020     mul32To64( aSig, bSig, &zSig0, &zSig1 );
1041     bits32 aSig, bSig, zSig, rem0, rem1, term0, term1;
1046     bSig = extractFloat32Frac( b );
1053             if ( bSig ) return propagateFloat32NaN( a, b );
1060         if ( bSig ) return propagateFloat32NaN( a, b );
1064         if ( bSig == 0 ) {
1072         normalizeFloat32Subnormal( bSig, &bExp, &bSig );
1080     bSig = ( bSig | 0x00800000 )<<8;
1081     if ( bSig <= ( aSig + aSig ) ) {
1085     zSig = estimateDiv64To32( aSig, 0, bSig );
1087         mul32To64( bSig, zSig, &term0, &term1 );
1091             add64( rem0, rem1, 0, bSig, &rem0, &rem1 );
1111     bits32 aSig, bSig, q, allZero, alternateASig;
1117     bSig = extractFloat32Frac( b );
1121         if ( aSig || ( ( bExp == 0xFF ) && bSig ) ) {
1128         if ( bSig ) return propagateFloat32NaN( a, b );
1132         if ( bSig == 0 ) {
1136         normalizeFloat32Subnormal( bSig, &bExp, &bSig );
1144     bSig = ( bSig | 0x00800000 )<<8;
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 );
1163         bSig >>= 2;
1164         aSig = ( ( aSig>>1 )<<( expDiff - 1 ) ) - bSig * q;
1168         bSig >>= 2;
1173         aSig -= bSig;