libc/softfloat/softfloat.c

789     bits32 aSig, bSig, zSig;
794     bSig = extractFloat32Frac( b );
798     bSig <<= 6;
808             bSig |= 0x20000000;
810         shift32RightJamming( bSig, expDiff, &bSig );
815             if ( bSig ) return propagateFloat32NaN( a, b );
829             if ( aSig | bSig ) return propagateFloat32NaN( a, b );
832         if ( aExp == 0 ) return packFloat32( zSign, 0, ( aSig + bSig )>>6 );
833         zSig = 0x40000000 + aSig + bSig;
838     zSig = ( aSig + bSig )<<1;
841         zSig = aSig + bSig;
861     bits32 aSig, bSig, zSig;
866     bSig = extractFloat32Frac( b );
870     bSig <<= 7;
874         if ( aSig | bSig ) return propagateFloat32NaN( a, b );
882     if ( bSig < aSig ) goto aBigger;
883     if ( aSig < bSig ) goto bBigger;
887         if ( bSig ) return propagateFloat32NaN( a, b );
897     bSig |= 0x40000000;
899     zSig = bSig - aSig;
912         bSig |= 0x40000000;
914     shift32RightJamming( bSig, expDiff, &bSig );
917     zSig = aSig - bSig;
980     bits32 aSig, bSig, zSig0, zSig1;
985     bSig = extractFloat32Frac( b );
990         if ( aSig || ( ( bExp == 0xFF ) && bSig ) ) {
993         if ( ( bExp | bSig ) == 0 ) {
1000         if ( bSig ) return propagateFloat32NaN( a, b );
1012         if ( bSig == 0 ) return packFloat32( zSign, 0, 0 );
1013         normalizeFloat32Subnormal( bSig, &bExp, &bSig );
1017     bSig = ( bSig | 0x00800000 )<<8;
1018     mul32To64( aSig, bSig, &zSig0, &zSig1 );
1039     bits32 aSig, bSig, zSig, rem0, rem1, term0, term1;
1044     bSig = extractFloat32Frac( b );
1051             if ( bSig ) return propagateFloat32NaN( a, b );
1058         if ( bSig ) return propagateFloat32NaN( a, b );
1062         if ( bSig == 0 ) {
1070         normalizeFloat32Subnormal( bSig, &bExp, &bSig );
1078     bSig = ( bSig | 0x00800000 )<<8;
1079     if ( bSig <= ( aSig + aSig ) ) {
1083     zSig = estimateDiv64To32( aSig, 0, bSig );
1085         mul32To64( bSig, zSig, &term0, &term1 );
1089             add64( rem0, rem1, 0, bSig, &rem0, &rem1 );
1109     bits32 aSig, bSig, q, allZero, alternateASig;
1115     bSig = extractFloat32Frac( b );
1119         if ( aSig || ( ( bExp == 0xFF ) && bSig ) ) {
1126         if ( bSig ) return propagateFloat32NaN( a, b );
1130         if ( bSig == 0 ) {
1134         normalizeFloat32Subnormal( bSig, &bExp, &bSig );
1142     bSig = ( bSig | 0x00800000 )<<8;
1147     q = ( bSig <= aSig );
1148     if ( q ) aSig -= bSig;
1151         q = estimateDiv64To32( aSig, 0, bSig );
1153         aSig = - ( ( bSig>>2 ) * q );
1158         q = estimateDiv64To32( aSig, 0, bSig );
1161         bSig >>= 2;
1162         aSig = ( ( aSig>>1 )<<( expDiff - 1 ) ) - bSig * q;
1166         bSig >>= 2;
1171         aSig -= bSig;