792     int16 expDiff;  in addFloat32Sigs()  local
798     expDiff = aExp - bExp;  in addFloat32Sigs()
801     if ( 0 < expDiff ) {  in addFloat32Sigs()
807             --expDiff;  in addFloat32Sigs()
812         shift32RightJamming( bSig, expDiff, &bSig );  in addFloat32Sigs()
815     else if ( expDiff < 0 ) {  in addFloat32Sigs()
821             ++expDiff;  in addFloat32Sigs()
826         shift32RightJamming( aSig, - expDiff, &aSig );  in addFloat32Sigs()
864     int16 expDiff;  in subFloat32Sigs()  local
870     expDiff = aExp - bExp;  in subFloat32Sigs()
873     if ( 0 < expDiff ) goto aExpBigger;  in subFloat32Sigs()
874     if ( expDiff < 0 ) goto bExpBigger;  in subFloat32Sigs()
893         ++expDiff;  in subFloat32Sigs()
898     shift32RightJamming( aSig, - expDiff, &aSig );  in subFloat32Sigs()
911         --expDiff;  in subFloat32Sigs()
916     shift32RightJamming( bSig, expDiff, &bSig );  in subFloat32Sigs()
1110     int16 aExp, bExp, expDiff;  in float32_rem()  local
1142     expDiff = aExp - bExp;  in float32_rem()
1145     if ( expDiff < 0 ) {  in float32_rem()
1146         if ( expDiff < -1 ) return a;  in float32_rem()
1151     expDiff -= 32;  in float32_rem()
1152     while ( 0 < expDiff ) {  in float32_rem()
1156         expDiff -= 30;  in float32_rem()
1158     expDiff += 32;  in float32_rem()
1159     if ( 0 < expDiff ) {  in float32_rem()
1162         q >>= 32 - expDiff;  in float32_rem()
1164         aSig = ( ( aSig>>1 )<<( expDiff - 1 ) ) - bSig * q;  in float32_rem()
1671     int16 expDiff;  in addFloat64Sigs()  local
1679     expDiff = aExp - bExp;  in addFloat64Sigs()
1680     if ( 0 < expDiff ) {  in addFloat64Sigs()
1686             --expDiff;  in addFloat64Sigs()
1692             bSig0, bSig1, 0, expDiff, &bSig0, &bSig1, &zSig2 );  in addFloat64Sigs()
1695     else if ( expDiff < 0 ) {  in addFloat64Sigs()
1701             ++expDiff;  in addFloat64Sigs()
1707             aSig0, aSig1, 0, - expDiff, &aSig0, &aSig1, &zSig2 );  in addFloat64Sigs()
1749     int16 expDiff;  in subFloat64Sigs()  local
1757     expDiff = aExp - bExp;  in subFloat64Sigs()
1760     if ( 0 < expDiff ) goto aExpBigger;  in subFloat64Sigs()
1761     if ( expDiff < 0 ) goto bExpBigger;  in subFloat64Sigs()
1784         ++expDiff;  in subFloat64Sigs()
1789     shift64RightJamming( aSig0, aSig1, - expDiff, &aSig0, &aSig1 );  in subFloat64Sigs()
1802         --expDiff;  in subFloat64Sigs()
1807     shift64RightJamming( bSig0, bSig1, expDiff, &bSig0, &bSig1 );  in subFloat64Sigs()
2015     int16 aExp, bExp, expDiff;  in float64_rem()  local
2052     expDiff = aExp - bExp;  in float64_rem()
2053     if ( expDiff < -1 ) return a;  in float64_rem()
2055         aSig0 | 0x00100000, aSig1, 11 - ( expDiff < 0 ), &aSig0, &aSig1 );  in float64_rem()
2059     expDiff -= 32;  in float64_rem()
2060     while ( 0 < expDiff ) {  in float64_rem()
2067         expDiff -= 29;  in float64_rem()
2069     if ( -32 < expDiff ) {  in float64_rem()
2072         q >>= - expDiff;  in float64_rem()
2074         expDiff += 24;  in float64_rem()
2075         if ( expDiff < 0 ) {  in float64_rem()
2076             shift64Right( aSig0, aSig1, - expDiff, &aSig0, &aSig1 );  in float64_rem()
2079             shortShift64Left( aSig0, aSig1, expDiff, &aSig0, &aSig1 );  in float64_rem()