Lines Matching refs:bSig
1004 bits32 aSig, bSig, zSig;
1009 bSig = extractFloat32Frac( b );
1013 bSig <<= 6;
1023 bSig |= 0x20000000;
1025 shift32RightJamming( bSig, expDiff, &bSig );
1030 if ( bSig ) return propagateFloat32NaN( a, b );
1044 if ( aSig | bSig ) return propagateFloat32NaN( a, b );
1047 if ( aExp == 0 ) return packFloat32( zSign, 0, ( aSig + bSig )>>6 );
1048 zSig = 0x40000000 + aSig + bSig;
1053 zSig = ( aSig + bSig )<<1;
1056 zSig = aSig + bSig;
1076 bits32 aSig, bSig, zSig;
1081 bSig = extractFloat32Frac( b );
1085 bSig <<= 7;
1089 if ( aSig | bSig ) return propagateFloat32NaN( a, b );
1097 if ( bSig < aSig ) goto aBigger;
1098 if ( aSig < bSig ) goto bBigger;
1102 if ( bSig ) return propagateFloat32NaN( a, b );
1112 bSig |= 0x40000000;
1114 zSig = bSig - aSig;
1127 bSig |= 0x40000000;
1129 shift32RightJamming( bSig, expDiff, &bSig );
1132 zSig = aSig - bSig;
1195 bits32 aSig, bSig;
1202 bSig = extractFloat32Frac( b );
1207 if ( aSig || ( ( bExp == 0xFF ) && bSig ) ) {
1210 if ( ( bExp | bSig ) == 0 ) {
1217 if ( bSig ) return propagateFloat32NaN( a, b );
1229 if ( bSig == 0 ) return packFloat32( zSign, 0, 0 );
1230 normalizeFloat32Subnormal( bSig, &bExp, &bSig );
1234 bSig = ( bSig | 0x00800000 )<<8;
1235 shift64RightJamming( ( (bits64) aSig ) * bSig, 32, &zSig64 );
1256 bits32 aSig, bSig, zSig;
1261 bSig = extractFloat32Frac( b );
1268 if ( bSig ) return propagateFloat32NaN( a, b );
1275 if ( bSig ) return propagateFloat32NaN( a, b );
1279 if ( bSig == 0 ) {
1287 normalizeFloat32Subnormal( bSig, &bExp, &bSig );
1295 bSig = ( bSig | 0x00800000 )<<8;
1296 if ( bSig <= ( aSig + aSig ) ) {
1302 do_div( tmp, bSig );
1306 zSig |= ( ( (bits64) bSig ) * zSig != ( (bits64) aSig )<<32 );
1323 bits32 aSig, bSig;
1332 bSig = extractFloat32Frac( b );
1336 if ( aSig || ( ( bExp == 0xFF ) && bSig ) ) {
1343 if ( bSig ) return propagateFloat32NaN( a, b );
1347 if ( bSig == 0 ) {
1351 normalizeFloat32Subnormal( bSig, &bExp, &bSig );
1359 bSig |= 0x00800000;
1362 bSig <<= 8;
1367 q = ( bSig <= aSig );
1368 if ( q ) aSig -= bSig;
1371 do_div( tmp, bSig );
1374 bSig >>= 2;
1375 aSig = ( ( aSig>>1 )<<( expDiff - 1 ) ) - bSig * q;
1379 bSig >>= 2;
1383 if ( bSig <= aSig ) aSig -= bSig;
1385 bSig64 = ( (bits64) bSig )<<40;
1390 aSig64 = - ( ( bSig * q64 )<<38 );
1397 bSig <<= 6;
1398 aSig = ( ( aSig64>>33 )<<( expDiff - 1 ) ) - bSig * q;
1403 aSig -= bSig;
1900 bits64 aSig, bSig, zSig;
1905 bSig = extractFloat64Frac( b );
1909 bSig <<= 9;
1919 bSig |= LIT64( 0x2000000000000000 );
1921 shift64RightJamming( bSig, expDiff, &bSig );
1926 if ( bSig ) return propagateFloat64NaN( a, b );
1940 if ( aSig | bSig ) return propagateFloat64NaN( a, b );
1943 if ( aExp == 0 ) return packFloat64( zSign, 0, ( aSig + bSig )>>9 );
1944 zSig = LIT64( 0x4000000000000000 ) + aSig + bSig;
1949 zSig = ( aSig + bSig )<<1;
1952 zSig = aSig + bSig;
1972 bits64 aSig, bSig, zSig;
1977 bSig = extractFloat64Frac( b );
1981 bSig <<= 10;
1985 if ( aSig | bSig ) return propagateFloat64NaN( a, b );
1993 if ( bSig < aSig ) goto aBigger;
1994 if ( aSig < bSig ) goto bBigger;
1998 if ( bSig ) return propagateFloat64NaN( a, b );
2008 bSig |= LIT64( 0x4000000000000000 );
2010 zSig = bSig - aSig;
2023 bSig |= LIT64( 0x4000000000000000 );
2025 shift64RightJamming( bSig, expDiff, &bSig );
2028 zSig = aSig - bSig;
2091 bits64 aSig, bSig, zSig0, zSig1;
2096 bSig = extractFloat64Frac( b );
2101 if ( aSig || ( ( bExp == 0x7FF ) && bSig ) ) {
2104 if ( ( bExp | bSig ) == 0 ) {
2111 if ( bSig ) return propagateFloat64NaN( a, b );
2123 if ( bSig == 0 ) return packFloat64( zSign, 0, 0 );
2124 normalizeFloat64Subnormal( bSig, &bExp, &bSig );
2128 bSig = ( bSig | LIT64( 0x0010000000000000 ) )<<11;
2129 mul64To128( aSig, bSig, &zSig0, &zSig1 );
2150 bits64 aSig, bSig, zSig;
2157 bSig = extractFloat64Frac( b );
2164 if ( bSig ) return propagateFloat64NaN( a, b );
2171 if ( bSig ) return propagateFloat64NaN( a, b );
2175 if ( bSig == 0 ) {
2183 normalizeFloat64Subnormal( bSig, &bExp, &bSig );
2191 bSig = ( bSig | LIT64( 0x0010000000000000 ) )<<11;
2192 if ( bSig <= ( aSig + aSig ) ) {
2196 zSig = estimateDiv128To64( aSig, 0, bSig );
2198 mul64To128( bSig, zSig, &term0, &term1 );
2202 add128( rem0, rem1, 0, bSig, &rem0, &rem1 );
2221 bits64 aSig, bSig;
2228 bSig = extractFloat64Frac( b );
2232 if ( aSig || ( ( bExp == 0x7FF ) && bSig ) ) {
2239 if ( bSig ) return propagateFloat64NaN( a, b );
2243 if ( bSig == 0 ) {
2247 normalizeFloat64Subnormal( bSig, &bExp, &bSig );
2255 bSig = ( bSig | LIT64( 0x0010000000000000 ) )<<11;
2260 q = ( bSig <= aSig );
2261 if ( q ) aSig -= bSig;
2264 q = estimateDiv128To64( aSig, 0, bSig );
2266 aSig = - ( ( bSig>>2 ) * q );
2271 q = estimateDiv128To64( aSig, 0, bSig );
2274 bSig >>= 2;
2275 aSig = ( ( aSig>>1 )<<( expDiff - 1 ) ) - bSig * q;
2279 bSig >>= 2;
2284 aSig -= bSig;
2720 bits64 aSig, bSig, zSig0, zSig1;
2725 bSig = extractFloatx80Frac( b );
2734 shift64ExtraRightJamming( bSig, 0, expDiff, &bSig, &zSig1 );
2739 if ( (bits64) ( bSig<<1 ) ) return propagateFloatx80NaN( a, b );
2748 if ( (bits64) ( ( aSig | bSig )<<1 ) ) {
2754 zSig0 = aSig + bSig;
2763 zSig0 = aSig + bSig;
2789 bits64 aSig, bSig, zSig0, zSig1;
2795 bSig = extractFloatx80Frac( b );
2801 if ( (bits64) ( ( aSig | bSig )<<1 ) ) {
2815 if ( bSig < aSig ) goto aBigger;
2816 if ( aSig < bSig ) goto bBigger;
2820 if ( (bits64) ( bSig<<1 ) ) return propagateFloatx80NaN( a, b );
2826 sub128( bSig, 0, aSig, zSig1, &zSig0, &zSig1 );
2836 shift128RightJamming( bSig, 0, expDiff, &bSig, &zSig1 );
2838 sub128( aSig, 0, bSig, zSig1, &zSig0, &zSig1 );
2902 bits64 aSig, bSig, zSig0, zSig1;
2908 bSig = extractFloatx80Frac( b );
2914 || ( ( bExp == 0x7FFF ) && (bits64) ( bSig<<1 ) ) ) {
2917 if ( ( bExp | bSig ) == 0 ) goto invalid;
2921 if ( (bits64) ( bSig<<1 ) ) return propagateFloatx80NaN( a, b );
2937 if ( bSig == 0 ) return packFloatx80( zSign, 0, 0 );
2938 normalizeFloatx80Subnormal( bSig, &bExp, &bSig );
2941 mul64To128( aSig, bSig, &zSig0, &zSig1 );
2963 bits64 aSig, bSig, zSig0, zSig1;
2970 bSig = extractFloatx80Frac( b );
2977 if ( (bits64) ( bSig<<1 ) ) return propagateFloatx80NaN( a, b );
2983 if ( (bits64) ( bSig<<1 ) ) return propagateFloatx80NaN( a, b );
2987 if ( bSig == 0 ) {
2999 normalizeFloatx80Subnormal( bSig, &bExp, &bSig );
3007 if ( bSig <= aSig ) {
3011 zSig0 = estimateDiv128To64( aSig, rem1, bSig );
3012 mul64To128( bSig, zSig0, &term0, &term1 );
3016 add128( rem0, rem1, 0, bSig, &rem0, &rem1 );
3018 zSig1 = estimateDiv128To64( rem1, 0, bSig );
3020 mul64To128( bSig, zSig1, &term1, &term2 );
3024 add128( rem1, rem2, 0, bSig, &rem1, &rem2 );
3045 bits64 aSig0, aSig1, bSig;
3052 bSig = extractFloatx80Frac( b );
3057 || ( ( bExp == 0x7FFF ) && (bits64) ( bSig<<1 ) ) ) {
3063 if ( (bits64) ( bSig<<1 ) ) return propagateFloatx80NaN( a, b );
3067 if ( bSig == 0 ) {
3075 normalizeFloatx80Subnormal( bSig, &bExp, &bSig );
3081 bSig |= LIT64( 0x8000000000000000 );
3090 q = ( bSig <= aSig0 );
3091 if ( q ) aSig0 -= bSig;
3094 q = estimateDiv128To64( aSig0, aSig1, bSig );
3096 mul64To128( bSig, q, &term0, &term1 );
3103 q = estimateDiv128To64( aSig0, aSig1, bSig );
3106 mul64To128( bSig, q<<( 64 - expDiff ), &term0, &term1 );
3108 shortShift128Left( 0, bSig, 64 - expDiff, &term0, &term1 );
3116 term0 = bSig;