788     bits32 aSig, bSig, zSig;  in addFloat32Sigs()  local
793     bSig = extractFloat32Frac( b );  in addFloat32Sigs()
797     bSig <<= 6;  in addFloat32Sigs()
807             bSig |= 0x20000000;  in addFloat32Sigs()
809         shift32RightJamming( bSig, expDiff, &bSig );  in addFloat32Sigs()
814             if ( bSig ) return propagateFloat32NaN( a, b );  in addFloat32Sigs()
828             if ( aSig | bSig ) return propagateFloat32NaN( a, b );  in addFloat32Sigs()
831         if ( aExp == 0 ) return packFloat32( zSign, 0, ( aSig + bSig )>>6 );  in addFloat32Sigs()
832         zSig = 0x40000000 + aSig + bSig;  in addFloat32Sigs()
837     zSig = ( aSig + bSig )<<1;  in addFloat32Sigs()
840         zSig = aSig + bSig;  in addFloat32Sigs()
860     bits32 aSig, bSig, zSig;  in subFloat32Sigs()  local
865     bSig = extractFloat32Frac( b );  in subFloat32Sigs()
869     bSig <<= 7;  in subFloat32Sigs()
873         if ( aSig | bSig ) return propagateFloat32NaN( a, b );  in subFloat32Sigs()
881     if ( bSig < aSig ) goto aBigger;  in subFloat32Sigs()
882     if ( aSig < bSig ) goto bBigger;  in subFloat32Sigs()
886         if ( bSig ) return propagateFloat32NaN( a, b );  in subFloat32Sigs()
896     bSig |= 0x40000000;  in subFloat32Sigs()
898     zSig = bSig - aSig;  in subFloat32Sigs()
911         bSig |= 0x40000000;  in subFloat32Sigs()
913     shift32RightJamming( bSig, expDiff, &bSig );  in subFloat32Sigs()
916     zSig = aSig - bSig;  in subFloat32Sigs()
979     bits32 aSig, bSig, zSig0, zSig1;  in float32_mul()  local
984     bSig = extractFloat32Frac( b );  in float32_mul()
989         if ( aSig || ( ( bExp == 0xFF ) && bSig ) ) {  in float32_mul()
992         if ( ( bExp | bSig ) == 0 ) {  in float32_mul()
999         if ( bSig ) return propagateFloat32NaN( a, b );  in float32_mul()
1011         if ( bSig == 0 ) return packFloat32( zSign, 0, 0 );  in float32_mul()
1012         normalizeFloat32Subnormal( bSig, &bExp, &bSig );  in float32_mul()
1016     bSig = ( bSig | 0x00800000 )<<8;  in float32_mul()
1017     mul32To64( aSig, bSig, &zSig0, &zSig1 );  in float32_mul()
1038     bits32 aSig, bSig, zSig, rem0, rem1, term0, term1;  in float32_div()  local
1043     bSig = extractFloat32Frac( b );  in float32_div()
1050             if ( bSig ) return propagateFloat32NaN( a, b );  in float32_div()
1057         if ( bSig ) return propagateFloat32NaN( a, b );  in float32_div()
1061         if ( bSig == 0 ) {  in float32_div()
1069         normalizeFloat32Subnormal( bSig, &bExp, &bSig );  in float32_div()
1077     bSig = ( bSig | 0x00800000 )<<8;  in float32_div()
1078     if ( bSig <= ( aSig + aSig ) ) {  in float32_div()
1082     zSig = estimateDiv64To32( aSig, 0, bSig );  in float32_div()
1084         mul32To64( bSig, zSig, &term0, &term1 );  in float32_div()
1088             add64( rem0, rem1, 0, bSig, &rem0, &rem1 );  in float32_div()
1108     bits32 aSig, bSig, q, allZero, alternateASig;  in float32_rem()  local
1114     bSig = extractFloat32Frac( b );  in float32_rem()
1118         if ( aSig || ( ( bExp == 0xFF ) && bSig ) ) {  in float32_rem()
1125         if ( bSig ) return propagateFloat32NaN( a, b );  in float32_rem()
1129         if ( bSig == 0 ) {  in float32_rem()
1133         normalizeFloat32Subnormal( bSig, &bExp, &bSig );  in float32_rem()
1141     bSig = ( bSig | 0x00800000 )<<8;  in float32_rem()
1146     q = ( bSig <= aSig );  in float32_rem()
1147     if ( q ) aSig -= bSig;  in float32_rem()
1150         q = estimateDiv64To32( aSig, 0, bSig );  in float32_rem()
1152         aSig = - ( ( bSig>>2 ) * q );  in float32_rem()
1157         q = estimateDiv64To32( aSig, 0, bSig );  in float32_rem()
1160         bSig >>= 2;  in float32_rem()
1161         aSig = ( ( aSig>>1 )<<( expDiff - 1 ) ) - bSig * q;  in float32_rem()
1165         bSig >>= 2;  in float32_rem()
1170         aSig -= bSig;  in float32_rem()