Lines Matching refs:bSig
1022 bits32 aSig, bSig, zSig;
1027 bSig = extractFloat32Frac( b );
1031 bSig <<= 6;
1041 bSig |= 0x20000000;
1043 shift32RightJamming( bSig, expDiff, &bSig );
1048 if ( bSig ) return propagateFloat32NaN( a, b );
1062 if ( aSig | bSig ) return propagateFloat32NaN( a, b );
1065 if ( aExp == 0 ) return packFloat32( zSign, 0, ( aSig + bSig )>>6 );
1066 zSig = 0x40000000 + aSig + bSig;
1071 zSig = ( aSig + bSig )<<1;
1074 zSig = aSig + bSig;
1094 bits32 aSig, bSig, zSig;
1099 bSig = extractFloat32Frac( b );
1103 bSig <<= 7;
1107 if ( aSig | bSig ) return propagateFloat32NaN( a, b );
1115 if ( bSig < aSig ) goto aBigger;
1116 if ( aSig < bSig ) goto bBigger;
1120 if ( bSig ) return propagateFloat32NaN( a, b );
1130 bSig |= 0x40000000;
1132 zSig = bSig - aSig;
1145 bSig |= 0x40000000;
1147 shift32RightJamming( bSig, expDiff, &bSig );
1150 zSig = aSig - bSig;
1213 bits32 aSig, bSig;
1220 bSig = extractFloat32Frac( b );
1225 if ( aSig || ( ( bExp == 0xFF ) && bSig ) ) {
1228 if ( ( bExp | bSig ) == 0 ) {
1235 if ( bSig ) return propagateFloat32NaN( a, b );
1247 if ( bSig == 0 ) return packFloat32( zSign, 0, 0 );
1248 normalizeFloat32Subnormal( bSig, &bExp, &bSig );
1252 bSig = ( bSig | 0x00800000 )<<8;
1253 shift64RightJamming( ( (bits64) aSig ) * bSig, 32, &zSig64 );
1274 bits32 aSig, bSig, zSig;
1279 bSig = extractFloat32Frac( b );
1286 if ( bSig ) return propagateFloat32NaN( a, b );
1293 if ( bSig ) return propagateFloat32NaN( a, b );
1297 if ( bSig == 0 ) {
1305 normalizeFloat32Subnormal( bSig, &bExp, &bSig );
1313 bSig = ( bSig | 0x00800000 )<<8;
1314 if ( bSig <= ( aSig + aSig ) ) {
1318 zSig = ( ( (bits64) aSig )<<32 ) / bSig;
1320 zSig |= ( ( (bits64) bSig ) * zSig != ( (bits64) aSig )<<32 );
1337 bits32 aSig, bSig;
1346 bSig = extractFloat32Frac( b );
1350 if ( aSig || ( ( bExp == 0xFF ) && bSig ) ) {
1357 if ( bSig ) return propagateFloat32NaN( a, b );
1361 if ( bSig == 0 ) {
1365 normalizeFloat32Subnormal( bSig, &bExp, &bSig );
1373 bSig |= 0x00800000;
1376 bSig <<= 8;
1381 q = ( bSig <= aSig );
1382 if ( q ) aSig -= bSig;
1384 q = ( ( (bits64) aSig )<<32 ) / bSig;
1386 bSig >>= 2;
1387 aSig = ( ( aSig>>1 )<<( expDiff - 1 ) ) - bSig * q;
1391 bSig >>= 2;
1395 if ( bSig <= aSig ) aSig -= bSig;
1397 bSig64 = ( (bits64) bSig )<<40;
1402 aSig64 = - ( ( bSig * q64 )<<38 );
1409 bSig <<= 6;
1410 aSig = ( ( aSig64>>33 )<<( expDiff - 1 ) ) - bSig * q;
1415 aSig -= bSig;
1916 bits64 aSig, bSig, zSig;
1921 bSig = extractFloat64Frac( b );
1925 bSig <<= 9;
1935 bSig |= LIT64( 0x2000000000000000 );
1937 shift64RightJamming( bSig, expDiff, &bSig );
1942 if ( bSig ) return propagateFloat64NaN( a, b );
1956 if ( aSig | bSig ) return propagateFloat64NaN( a, b );
1959 if ( aExp == 0 ) return packFloat64( zSign, 0, ( aSig + bSig )>>9 );
1960 zSig = LIT64( 0x4000000000000000 ) + aSig + bSig;
1965 zSig = ( aSig + bSig )<<1;
1968 zSig = aSig + bSig;
1988 bits64 aSig, bSig, zSig;
1993 bSig = extractFloat64Frac( b );
1997 bSig <<= 10;
2001 if ( aSig | bSig ) return propagateFloat64NaN( a, b );
2009 if ( bSig < aSig ) goto aBigger;
2010 if ( aSig < bSig ) goto bBigger;
2014 if ( bSig ) return propagateFloat64NaN( a, b );
2024 bSig |= LIT64( 0x4000000000000000 );
2026 zSig = bSig - aSig;
2039 bSig |= LIT64( 0x4000000000000000 );
2041 shift64RightJamming( bSig, expDiff, &bSig );
2044 zSig = aSig - bSig;
2107 bits64 aSig, bSig, zSig0, zSig1;
2112 bSig = extractFloat64Frac( b );
2117 if ( aSig || ( ( bExp == 0x7FF ) && bSig ) ) {
2120 if ( ( bExp | bSig ) == 0 ) {
2127 if ( bSig ) return propagateFloat64NaN( a, b );
2139 if ( bSig == 0 ) return packFloat64( zSign, 0, 0 );
2140 normalizeFloat64Subnormal( bSig, &bExp, &bSig );
2144 bSig = ( bSig | LIT64( 0x0010000000000000 ) )<<11;
2145 mul64To128( aSig, bSig, &zSig0, &zSig1 );
2166 bits64 aSig, bSig, zSig;
2173 bSig = extractFloat64Frac( b );
2180 if ( bSig ) return propagateFloat64NaN( a, b );
2187 if ( bSig ) return propagateFloat64NaN( a, b );
2191 if ( bSig == 0 ) {
2199 normalizeFloat64Subnormal( bSig, &bExp, &bSig );
2207 bSig = ( bSig | LIT64( 0x0010000000000000 ) )<<11;
2208 if ( bSig <= ( aSig + aSig ) ) {
2212 zSig = estimateDiv128To64( aSig, 0, bSig );
2214 mul64To128( bSig, zSig, &term0, &term1 );
2218 add128( rem0, rem1, 0, bSig, &rem0, &rem1 );
2237 bits64 aSig, bSig;
2244 bSig = extractFloat64Frac( b );
2248 if ( aSig || ( ( bExp == 0x7FF ) && bSig ) ) {
2255 if ( bSig ) return propagateFloat64NaN( a, b );
2259 if ( bSig == 0 ) {
2263 normalizeFloat64Subnormal( bSig, &bExp, &bSig );
2271 bSig = ( bSig | LIT64( 0x0010000000000000 ) )<<11;
2276 q = ( bSig <= aSig );
2277 if ( q ) aSig -= bSig;
2280 q = estimateDiv128To64( aSig, 0, bSig );
2282 aSig = - ( ( bSig>>2 ) * q );
2287 q = estimateDiv128To64( aSig, 0, bSig );
2290 bSig >>= 2;
2291 aSig = ( ( aSig>>1 )<<( expDiff - 1 ) ) - bSig * q;
2295 bSig >>= 2;
2300 aSig -= bSig;
2740 bits64 aSig, bSig, zSig0, zSig1;
2745 bSig = extractFloatx80Frac( b );
2754 shift64ExtraRightJamming( bSig, 0, expDiff, &bSig, &zSig1 );
2759 if ( (bits64) ( bSig<<1 ) ) return propagateFloatx80NaN( a, b );
2768 if ( (bits64) ( ( aSig | bSig )<<1 ) ) {
2774 zSig0 = aSig + bSig;
2783 zSig0 = aSig + bSig;
2809 bits64 aSig, bSig, zSig0, zSig1;
2815 bSig = extractFloatx80Frac( b );
2821 if ( (bits64) ( ( aSig | bSig )<<1 ) ) {
2834 if ( bSig < aSig ) goto aBigger;
2835 if ( aSig < bSig ) goto bBigger;
2839 if ( (bits64) ( bSig<<1 ) ) return propagateFloatx80NaN( a, b );
2845 sub128( bSig, 0, aSig, zSig1, &zSig0, &zSig1 );
2855 shift128RightJamming( bSig, 0, expDiff, &bSig, &zSig1 );
2857 sub128( aSig, 0, bSig, zSig1, &zSig0, &zSig1 );
2921 bits64 aSig, bSig, zSig0, zSig1;
2927 bSig = extractFloatx80Frac( b );
2933 || ( ( bExp == 0x7FFF ) && (bits64) ( bSig<<1 ) ) ) {
2936 if ( ( bExp | bSig ) == 0 ) goto invalid;
2940 if ( (bits64) ( bSig<<1 ) ) return propagateFloatx80NaN( a, b );
2955 if ( bSig == 0 ) return packFloatx80( zSign, 0, 0 );
2956 normalizeFloatx80Subnormal( bSig, &bExp, &bSig );
2959 mul64To128( aSig, bSig, &zSig0, &zSig1 );
2981 bits64 aSig, bSig, zSig0, zSig1;
2988 bSig = extractFloatx80Frac( b );
2995 if ( (bits64) ( bSig<<1 ) ) return propagateFloatx80NaN( a, b );
3001 if ( (bits64) ( bSig<<1 ) ) return propagateFloatx80NaN( a, b );
3005 if ( bSig == 0 ) {
3016 normalizeFloatx80Subnormal( bSig, &bExp, &bSig );
3024 if ( bSig <= aSig ) {
3028 zSig0 = estimateDiv128To64( aSig, rem1, bSig );
3029 mul64To128( bSig, zSig0, &term0, &term1 );
3033 add128( rem0, rem1, 0, bSig, &rem0, &rem1 );
3035 zSig1 = estimateDiv128To64( rem1, 0, bSig );
3037 mul64To128( bSig, zSig1, &term1, &term2 );
3041 add128( rem1, rem2, 0, bSig, &rem1, &rem2 );
3062 bits64 aSig0, aSig1, bSig;
3069 bSig = extractFloatx80Frac( b );
3074 || ( ( bExp == 0x7FFF ) && (bits64) ( bSig<<1 ) ) ) {
3080 if ( (bits64) ( bSig<<1 ) ) return propagateFloatx80NaN( a, b );
3084 if ( bSig == 0 ) {
3091 normalizeFloatx80Subnormal( bSig, &bExp, &bSig );
3097 bSig |= LIT64( 0x8000000000000000 );
3106 q = ( bSig <= aSig0 );
3107 if ( q ) aSig0 -= bSig;
3110 q = estimateDiv128To64( aSig0, aSig1, bSig );
3112 mul64To128( bSig, q, &term0, &term1 );
3119 q = estimateDiv128To64( aSig0, aSig1, bSig );
3122 mul64To128( bSig, q<<( 64 - expDiff ), &term0, &term1 );
3124 shortShift128Left( 0, bSig, 64 - expDiff, &term0, &term1 );
3132 term0 = bSig;