softfloat/bits64/softfloat.c

897 `aSig0' and `aSig1'.  The normalized exponent is stored at the location
907      bits64 aSig1,
916         shiftCount = countLeadingZeros64( aSig1 ) - 15;
918             *zSig0Ptr = aSig1>>( - shiftCount );
919             *zSig1Ptr = aSig1<<( shiftCount & 63 );
922             *zSig0Ptr = aSig1<<shiftCount;
929         shortShift128Left( aSig0, aSig1, shiftCount, zSig0Ptr, zSig1Ptr );
3939     bits64 aSig0, aSig1, bSig;
3977     aSig1 = 0;
3980         shift128Right( aSig0, 0, 1, &aSig0, &aSig1 );
3987         q = estimateDiv128To64( aSig0, aSig1, bSig );
3990         sub128( aSig0, aSig1, term0, term1, &aSig0, &aSig1 );
3991         shortShift128Left( aSig0, aSig1, 62, &aSig0, &aSig1 );
3996         q = estimateDiv128To64( aSig0, aSig1, bSig );
4000         sub128( aSig0, aSig1, term0, term1, &aSig0, &aSig1 );
4002         while ( le128( term0, term1, aSig0, aSig1 ) ) {
4004             sub128( aSig0, aSig1, term0, term1, &aSig0, &aSig1 );
4011     sub128( term0, term1, aSig0, aSig1, &alternateASig0, &alternateASig1 );
4012     if (    lt128( alternateASig0, alternateASig1, aSig0, aSig1 )
4013          || (    eq128( alternateASig0, alternateASig1, aSig0, aSig1 )
4017         aSig1 = alternateASig1;
4022             80, zSign, bExp + expDiff, aSig0, aSig1 );
4037     bits64 aSig0, aSig1, zSig0, zSig1, doubleZSig0;
4063     shift128Right( aSig0, 0, 2 + ( aExp & 1 ), &aSig0, &aSig1 );
4064     zSig0 = estimateDiv128To64( aSig0, aSig1, zSig0<<32 ) + ( zSig0<<30 );
4067     sub128( aSig0, aSig1, term0, term1, &rem0, &rem1 );
4317     bits64 aSig0, aSig1;
4319     aSig1 = extractFloat128Frac1( a );
4323     if ( ( aExp == 0x7FFF ) && ( aSig0 | aSig1 ) ) aSign = 0;
4325     aSig0 |= ( aSig1 != 0 );
4347     bits64 aSig0, aSig1, savedASig;
4350     aSig1 = extractFloat128Frac1( a );
4354     aSig0 |= ( aSig1 != 0 );
4396     bits64 aSig0, aSig1;
4398     aSig1 = extractFloat128Frac1( a );
4409                       && ( aSig1 || ( aSig0 != LIT64( 0x0001000000000000 ) ) )
4416         shortShift128Left( aSig0, aSig1, - shiftCount, &aSig0, &aSig1 );
4419         shift64ExtraRightJamming( aSig0, aSig1, shiftCount, &aSig0, &aSig1 );
4421     return roundAndPackInt64( aSign, aSig0, aSig1 );
4440     bits64 aSig0, aSig1;
4443     aSig1 = extractFloat128Frac1( a );
4453                  && ( aSig1 < LIT64( 0x0002000000000000 ) ) ) {
4454                 if ( aSig1 ) float_exception_flags |= float_flag_inexact;
4458                 if ( ! aSign || ( ( aExp == 0x7FFF ) && ( aSig0 | aSig1 ) ) ) {
4464         z = ( aSig0<<shiftCount ) | ( aSig1>>( ( - shiftCount ) & 63 ) );
4465         if ( (bits64) ( aSig1<<shiftCount ) ) {
4471             if ( aExp | aSig0 | aSig1 ) {
4477         if (    aSig1
4496     bits64 aSig0, aSig1;
4499     aSig1 = extractFloat128Frac1( a );
4509                  && ( aSig1 < LIT64( 0x0002000000000000 ) ) ) {
4510                 if ( aSig1 ) float_exception_flags |= float_flag_inexact;
4517         z = ( aSig0<<shiftCount ) | ( aSig1>>( ( - shiftCount ) & 63 ) );
4518         if ( (bits64) ( aSig1<<shiftCount ) ) {
4524             if ( aExp | aSig0 | aSig1 ) {
4530         if (aSig1 || ( shiftCount && (bits64) ( aSig0<<( shiftCount & 63 ) ) ) ) {
4552     bits64 aSig0, aSig1;
4555     aSig1 = extractFloat128Frac1( a );
4560         if ( aSig0 | aSig1 ) {
4565     aSig0 |= ( aSig1 != 0 );
4588     bits64 aSig0, aSig1;
4590     aSig1 = extractFloat128Frac1( a );
4595         if ( aSig0 | aSig1 ) {
4600     shortShift128Left( aSig0, aSig1, 14, &aSig0, &aSig1 );
4601     aSig0 |= ( aSig1 != 0 );
4624     bits64 aSig0, aSig1;
4626     aSig1 = extractFloat128Frac1( a );
4631         if ( aSig0 | aSig1 ) {
4637         if ( ( aSig0 | aSig1 ) == 0 ) return packFloatx80( aSign, 0, 0 );
4638         normalizeFloat128Subnormal( aSig0, aSig1, &aExp, &aSig0, &aSig1 );
4643     shortShift128Left( aSig0, aSig1, 15, &aSig0, &aSig1 );
4644     return roundAndPackFloatx80( 80, aSign, aExp, aSig0, aSig1 );
4768     bits64 aSig0, aSig1, bSig0, bSig1, zSig0, zSig1, zSig2;
4771     aSig1 = extractFloat128Frac1( a );
4780             if ( aSig0 | aSig1 ) return propagateFloat128NaN( a, b );
4805             aSig0, aSig1, 0, - expDiff, &aSig0, &aSig1, &zSig2 );
4810             if ( aSig0 | aSig1 | bSig0 | bSig1 ) {
4815         add128( aSig0, aSig1, bSig0, bSig1, &zSig0, &zSig1 );
4823     add128( aSig0, aSig1, bSig0, bSig1, &zSig0, &zSig1 );
4847     bits64 aSig0, aSig1, bSig0, bSig1, zSig0, zSig1;
4851     aSig1 = extractFloat128Frac1( a );
4858     shortShift128Left( aSig0, aSig1, 14, &aSig0, &aSig1 );
4863         if ( aSig0 | aSig1 | bSig0 | bSig1 ) {
4877     if ( bSig1 < aSig1 ) goto aBigger;
4878     if ( aSig1 < bSig1 ) goto bBigger;
4891     shift128RightJamming( aSig0, aSig1, - expDiff, &aSig0, &aSig1 );
4894     sub128( bSig0, bSig1, aSig0, aSig1, &zSig0, &zSig1 );
4900         if ( aSig0 | aSig1 ) return propagateFloat128NaN( a, b );
4912     sub128( aSig0, aSig1, bSig0, bSig1, &zSig0, &zSig1 );
4975     bits64 aSig0, aSig1, bSig0, bSig1, zSig0, zSig1, zSig2, zSig3;
4978     aSig1 = extractFloat128Frac1( a );
4988         if (    ( aSig0 | aSig1 )
4997         if ( ( aExp | aSig0 | aSig1 ) == 0 ) {
5007         if ( ( aSig0 | aSig1 ) == 0 ) return packFloat128( zSign, 0, 0, 0 );
5008         normalizeFloat128Subnormal( aSig0, aSig1, &aExp, &aSig0, &aSig1 );
5017     mul128To256( aSig0, aSig1, bSig0, bSig1, &zSig0, &zSig1, &zSig2, &zSig3 );
5018     add128( zSig0, zSig1, aSig0, aSig1, &zSig0, &zSig1 );
5040     bits64 aSig0, aSig1, bSig0, bSig1, zSig0, zSig1, zSig2;
5044     aSig1 = extractFloat128Frac1( a );
5054         if ( aSig0 | aSig1 ) return propagateFloat128NaN( a, b );
5067             if ( ( aExp | aSig0 | aSig1 ) == 0 ) {
5080         if ( ( aSig0 | aSig1 ) == 0 ) return packFloat128( zSign, 0, 0, 0 );
5081         normalizeFloat128Subnormal( aSig0, aSig1, &aExp, &aSig0, &aSig1 );
5085         aSig0 | LIT64( 0x0001000000000000 ), aSig1, 15, &aSig0, &aSig1 );
5088     if ( le128( bSig0, bSig1, aSig0, aSig1 ) ) {
5089         shift128Right( aSig0, aSig1, 1, &aSig0, &aSig1 );
5092     zSig0 = estimateDiv128To64( aSig0, aSig1, bSig0 );
5094     sub192( aSig0, aSig1, 0, term0, term1, term2, &rem0, &rem1, &rem2 );
5125     bits64 aSig0, aSig1, bSig0, bSig1, q, term0, term1, term2;
5130     aSig1 = extractFloat128Frac1( a );
5139         if (    ( aSig0 | aSig1 )
5160         if ( ( aSig0 | aSig1 ) == 0 ) return a;
5161         normalizeFloat128Subnormal( aSig0, aSig1, &aExp, &aSig0, &aSig1 );
5167         aSig1,
5170         &aSig1
5174     q = le128( bSig0, bSig1, aSig0, aSig1 );
5175     if ( q ) sub128( aSig0, aSig1, bSig0, bSig1, &aSig0, &aSig1 );
5178         q = estimateDiv128To64( aSig0, aSig1, bSig0 );
5182         shortShift128Left( aSig0, aSig1, 61, &aSig0, &allZero );
5183         sub128( aSig0, 0, term1, term2, &aSig0, &aSig1 );
5187         q = estimateDiv128To64( aSig0, aSig1, bSig0 );
5193             shift128Right( aSig0, aSig1, - expDiff, &aSig0, &aSig1 );
5196             shortShift128Left( aSig0, aSig1, expDiff, &aSig0, &aSig1 );
5199         sub128( aSig0, aSig1, term1, term2, &aSig0, &aSig1 );
5202         shift128Right( aSig0, aSig1, 12, &aSig0, &aSig1 );
5207         alternateASig1 = aSig1;
5209         sub128( aSig0, aSig1, bSig0, bSig1, &aSig0, &aSig1 );
5212         aSig0, aSig1, alternateASig0, alternateASig1, (bits64 *)&sigMean0, &sigMean1 );
5216         aSig1 = alternateASig1;
5219     if ( zSign ) sub128( 0, 0, aSig0, aSig1, &aSig0, &aSig1 );
5221         normalizeRoundAndPackFloat128( aSign ^ zSign, bExp - 4, aSig0, aSig1 );
5236     bits64 aSig0, aSig1, zSig0, zSig1, zSig2, doubleZSig0;
5240     aSig1 = extractFloat128Frac1( a );
5245         if ( aSig0 | aSig1 ) return propagateFloat128NaN( a, a );
5250         if ( ( aExp | aSig0 | aSig1 ) == 0 ) return a;
5258         if ( ( aSig0 | aSig1 ) == 0 ) return packFloat128( 0, 0, 0, 0 );
5259         normalizeFloat128Subnormal( aSig0, aSig1, &aExp, &aSig0, &aSig1 );
5264     shortShift128Left( aSig0, aSig1, 13 - ( aExp & 1 ), &aSig0, &aSig1 );
5265     zSig0 = estimateDiv128To64( aSig0, aSig1, zSig0<<32 ) + ( zSig0<<30 );
5268     sub128( aSig0, aSig1, term0, term1, &rem0, &rem1 );