libc/softfloat/softfloat.c

325 `aSig1'.  The normalized exponent is stored at the location pointed to by
335      bits32 aSig1,
344         shiftCount = countLeadingZeros32( aSig1 ) - 11;
346             *zSig0Ptr = aSig1>>( - shiftCount );
347             *zSig1Ptr = aSig1<<( shiftCount & 31 );
350             *zSig0Ptr = aSig1<<shiftCount;
357         shortShift64Left( aSig0, aSig1, shiftCount, zSig0Ptr, zSig1Ptr );
1409     bits32 aSig0, aSig1, absZ, aSigExtra;
1413     aSig1 = extractFloat64Frac1( a );
1420             if ( ( aExp == 0x7FF ) && ( aSig0 | aSig1 ) ) aSign = 0;
1424             aSig0 | 0x00100000, aSig1, shiftCount, &absZ, &aSigExtra );
1428         aSig1 = ( aSig1 != 0 );
1430             aSigExtra = aExp | aSig0 | aSig1;
1435             aSigExtra = ( aSig0<<( shiftCount & 31 ) ) | aSig1;
1483     bits32 aSig0, aSig1, absZ, aSigExtra;
1486     aSig1 = extractFloat64Frac1( a );
1493             if ( ( aExp == 0x7FF ) && ( aSig0 | aSig1 ) ) aSign = 0;
1497             aSig0 | 0x00100000, aSig1, shiftCount, &absZ, &aSigExtra );
1501             if ( aExp | aSig0 | aSig1 ) {
1507         aSigExtra = ( aSig0<<( shiftCount & 31 ) ) | aSig1;
1533     bits32 aSig0, aSig1, zSig;
1536     aSig1 = extractFloat64Frac1( a );
1541         if ( aSig0 | aSig1 ) {
1546     shift64RightJamming( aSig0, aSig1, 22, &allZero, &zSig );
1668     bits32 aSig0, aSig1, bSig0, bSig1, zSig0, zSig1, zSig2;
1671     aSig1 = extractFloat64Frac1( a );
1680             if ( aSig0 | aSig1 ) return propagateFloat64NaN( a, b );
1705             aSig0, aSig1, 0, - expDiff, &aSig0, &aSig1, &zSig2 );
1710             if ( aSig0 | aSig1 | bSig0 | bSig1 ) {
1715         add64( aSig0, aSig1, bSig0, bSig1, &zSig0, &zSig1 );
1723     add64( aSig0, aSig1, bSig0, bSig1, &zSig0, &zSig1 );
1746     bits32 aSig0, aSig1, bSig0, bSig1, zSig0, zSig1;
1749     aSig1 = extractFloat64Frac1( a );
1756     shortShift64Left( aSig0, aSig1, 10, &aSig0, &aSig1 );
1761         if ( aSig0 | aSig1 | bSig0 | bSig1 ) {
1773     if ( bSig1 < aSig1 ) goto aBigger;
1774     if ( aSig1 < bSig1 ) goto bBigger;
1787     shift64RightJamming( aSig0, aSig1, - expDiff, &aSig0, &aSig1 );
1790     sub64( bSig0, bSig1, aSig0, aSig1, &zSig0, &zSig1 );
1796         if ( aSig0 | aSig1 ) return propagateFloat64NaN( a, b );
1808     sub64( aSig0, aSig1, bSig0, bSig1, &zSig0, &zSig1 );
1871     bits32 aSig0, aSig1, bSig0, bSig1, zSig0, zSig1, zSig2, zSig3;
1873     aSig1 = extractFloat64Frac1( a );
1883         if (    ( aSig0 | aSig1 )
1892         if ( ( aExp | aSig0 | aSig1 ) == 0 ) {
1900         if ( ( aSig0 | aSig1 ) == 0 ) return packFloat64( zSign, 0, 0, 0 );
1901         normalizeFloat64Subnormal( aSig0, aSig1, &aExp, &aSig0, &aSig1 );
1910     mul64To128( aSig0, aSig1, bSig0, bSig1, &zSig0, &zSig1, &zSig2, &zSig3 );
1911     add64( zSig0, zSig1, aSig0, aSig1, &zSig0, &zSig1 );
1933     bits32 aSig0, aSig1, bSig0, bSig1, zSig0, zSig1, zSig2;
1936     aSig1 = extractFloat64Frac1( a );
1946         if ( aSig0 | aSig1 ) return propagateFloat64NaN( a, b );
1959             if ( ( aExp | aSig0 | aSig1 ) == 0 ) {
1970         if ( ( aSig0 | aSig1 ) == 0 ) return packFloat64( zSign, 0, 0, 0 );
1971         normalizeFloat64Subnormal( aSig0, aSig1, &aExp, &aSig0, &aSig1 );
1974     shortShift64Left( aSig0 | 0x00100000, aSig1, 11, &aSig0, &aSig1 );
1976     if ( le64( bSig0, bSig1, aSig0, aSig1 ) ) {
1977         shift64Right( aSig0, aSig1, 1, &aSig0, &aSig1 );
1980     zSig0 = estimateDiv64To32( aSig0, aSig1, bSig0 );
1982     sub96( aSig0, aSig1, 0, term0, term1, term2, &rem0, &rem1, &rem2 );
2014     bits32 aSig0, aSig1, bSig0, bSig1, q, term0, term1, term2;
2019     aSig1 = extractFloat64Frac1( a );
2028         if (    ( aSig0 | aSig1 )
2047         if ( ( aSig0 | aSig1 ) == 0 ) return a;
2048         normalizeFloat64Subnormal( aSig0, aSig1, &aExp, &aSig0, &aSig1 );
2053         aSig0 | 0x00100000, aSig1, 11 - ( expDiff < 0 ), &aSig0, &aSig1 );
2055     q = le64( bSig0, bSig1, aSig0, aSig1 );
2056     if ( q ) sub64( aSig0, aSig1, bSig0, bSig1, &aSig0, &aSig1 );
2059         q = estimateDiv64To32( aSig0, aSig1, bSig0 );
2063         shortShift64Left( aSig0, aSig1, 29, &aSig0, &allZero );
2064         sub64( aSig0, 0, term1, term2, &aSig0, &aSig1 );
2068         q = estimateDiv64To32( aSig0, aSig1, bSig0 );
2074             shift64Right( aSig0, aSig1, - expDiff, &aSig0, &aSig1 );
2077             shortShift64Left( aSig0, aSig1, expDiff, &aSig0, &aSig1 );
2080         sub64( aSig0, aSig1, term1, term2, &aSig0, &aSig1 );
2083         shift64Right( aSig0, aSig1, 8, &aSig0, &aSig1 );
2088         alternateASig1 = aSig1;
2090         sub64( aSig0, aSig1, bSig0, bSig1, &aSig0, &aSig1 );
2093         aSig0, aSig1, alternateASig0, alternateASig1, &sigMean0, &sigMean1 );
2097         aSig1 = alternateASig1;
2100     if ( zSign ) sub64( 0, 0, aSig0, aSig1, &aSig0, &aSig1 );
2102         normalizeRoundAndPackFloat64( aSign ^ zSign, bExp - 4, aSig0, aSig1 );
2119     bits32 aSig0, aSig1, zSig0, zSig1, zSig2, doubleZSig0;
2123     aSig1 = extractFloat64Frac1( a );
2128         if ( aSig0 | aSig1 ) return propagateFloat64NaN( a, a );
2133         if ( ( aExp | aSig0 | aSig1 ) == 0 ) return a;
2139         if ( ( aSig0 | aSig1 ) == 0 ) return packFloat64( 0, 0, 0, 0 );
2140         normalizeFloat64Subnormal( aSig0, aSig1, &aExp, &aSig0, &aSig1 );
2144     shortShift64Left( aSig0, aSig1, 11, &term0, &term1 );
2148     shortShift64Left( aSig0, aSig1, 9 - ( aExp & 1 ), &aSig0, &aSig1 );
2150     sub64( aSig0, aSig1, term0, term1, &rem0, &rem1 );