softfloat/bits64/softfloat.c

627 and extended significand formed by the concatenation of `zSig0' and `zSig1',
651      int8 roundingPrecision, flag zSign, int32 zExp, bits64 zSig0, bits64 zSig1
672     zSig0 |= ( zSig1 != 0 );
728     increment = ( (sbits64) zSig1 < 0 );
735                 increment = ( roundingMode == float_round_down ) && zSig1;
738                 increment = ( roundingMode == float_round_up ) && zSig1;
766             shift64ExtraRightJamming( zSig0, zSig1, 1 - zExp, &zSig0, &zSig1 );
768             if ( isTiny && zSig1 ) float_raise( float_flag_underflow );
769             if ( zSig1 ) float_exception_flags |= float_flag_inexact;
771                 increment = ( (sbits64) zSig1 < 0 );
775                     increment = ( roundingMode == float_round_down ) && zSig1;
778                     increment = ( roundingMode == float_round_up ) && zSig1;
784                     ~ ( ( (bits64) ( zSig1<<1 ) == 0 ) & roundNearestEven );
790     if ( zSig1 ) float_exception_flags |= float_flag_inexact;
798             zSig0 &= ~ ( ( (bits64) ( zSig1<<1 ) == 0 ) & roundNearestEven );
811 `zExp', and significand formed by the concatenation of `zSig0' and `zSig1',
820      int8 roundingPrecision, flag zSign, int32 zExp, bits64 zSig0, bits64 zSig1
826         zSig0 = zSig1;
827         zSig1 = 0;
831     shortShift128Left( zSig0, zSig1, shiftCount, &zSig0, &zSig1 );
834         roundAndPackFloatx80( roundingPrecision, zSign, zExp, zSig0, zSig1 );
938 by the concatenation of `zSig0' and `zSig1' into a quadruple-precision
945 whenever `zSig0' and `zSig1' concatenated form a complete, normalized
950  packFloat128( flag zSign, int32 zExp, bits64 zSig0, bits64 zSig1 )
954     z.low = zSig1;
963 and extended significand formed by the concatenation of `zSig0', `zSig1',
984      flag zSign, int32 zExp, bits64 zSig0, bits64 zSig1, bits64 zSig2 )
1012                          zSig1
1039                        zSig1,
1044                 zSig0, zSig1, zSig2, - zExp, &zSig0, &zSig1, &zSig2 );
1062         add128( zSig0, zSig1, 0, 1, &zSig0, &zSig1 );
1063         zSig1 &= ~ ( ( zSig2 + zSig2 == 0 ) & roundNearestEven );
1066         if ( ( zSig0 | zSig1 ) == 0 ) zExp = 0;
1068     return packFloat128( zSign, zExp, zSig0, zSig1 );
1075 and significand formed by the concatenation of `zSig0' and `zSig1', and
1085      flag zSign, int32 zExp, bits64 zSig0, bits64 zSig1 )
1091         zSig0 = zSig1;
1092         zSig1 = 0;
1098         shortShift128Left( zSig0, zSig1, shiftCount, &zSig0, &zSig1 );
1102             zSig0, zSig1, 0, - shiftCount, &zSig0, &zSig1, &zSig2 );
1105     return roundAndPackFloat128( zSign, zExp, zSig0, zSig1, zSig2 );
1350     bits64 zSig0, zSig1;
1358         zSig1 = 0;
1363         zSig1 = absA;
1366     shortShift128Left( zSig0, zSig1, shiftCount, &zSig0, &zSig1 );
1367     return packFloat128( zSign, zExp, zSig0, zSig1 );
2577     bits64 aSig, zSig0, zSig1;
2591     shift128Right( aSig, 0, 4, &zSig0, &zSig1 );
2592     return packFloat128( aSign, aExp + 0x3C00, zSig0, zSig1 );
2866     bits64 aSig, bSig, zSig0, zSig1;
2904     mul64To128( aSig, bSig, &zSig0, &zSig1 );
2905     zSig0 |= ( zSig1 != 0 );
3516     bits64 aSig, zSig0, zSig1;
3524     shift128Right( aSig<<1, 0, 16, &zSig0, &zSig1 );
3525     return packFloat128( aSign, aExp, zSig0, zSig1 );
3619     bits64 aSig, bSig, zSig0, zSig1;
3633         shift64ExtraRightJamming( bSig, 0, expDiff, &bSig, &zSig1 );
3642         shift64ExtraRightJamming( aSig, 0, - expDiff, &aSig, &zSig1 );
3652         zSig1 = 0;
3664     shift64ExtraRightJamming( zSig0, zSig1, 1, &zSig0, &zSig1 );
3670             floatx80_rounding_precision, zSign, zExp, zSig0, zSig1 );
3686     bits64 aSig, bSig, zSig0, zSig1;
3710     zSig1 = 0;
3720     shift128RightJamming( aSig, 0, - expDiff, &aSig, &zSig1 );
3722     sub128( bSig, 0, aSig, zSig1, &zSig0, &zSig1 );
3732     shift128RightJamming( bSig, 0, expDiff, &bSig, &zSig1 );
3734     sub128( aSig, 0, bSig, zSig1, &zSig0, &zSig1 );
3739             floatx80_rounding_precision, zSign, zExp, zSig0, zSig1 );
3798     bits64 aSig, bSig, zSig0, zSig1;
3836     mul64To128( aSig, bSig, &zSig0, &zSig1 );
3838         shortShift128Left( zSig0, zSig1, 1, &zSig0, &zSig1 );
3843             floatx80_rounding_precision, zSign, zExp, zSig0, zSig1 );
3858     bits64 aSig, bSig, zSig0, zSig1;
3912     zSig1 = estimateDiv128To64( rem1, 0, bSig );
3913     if ( (bits64) ( zSig1<<1 ) <= 8 ) {
3914         mul64To128( bSig, zSig1, &term1, &term2 );
3917             --zSig1;
3920         zSig1 |= ( ( rem1 | rem2 ) != 0 );
3924             floatx80_rounding_precision, zSign, zExp, zSig0, zSig1 );
4037     bits64 aSig0, aSig1, zSig0, zSig1, doubleZSig0;
4073     zSig1 = estimateDiv128To64( rem1, 0, doubleZSig0 );
4074     if ( ( zSig1 & LIT64( 0x3FFFFFFFFFFFFFFF ) ) <= 5 ) {
4075         if ( zSig1 == 0 ) zSig1 = 1;
4076         mul64To128( doubleZSig0, zSig1, &term1, &term2 );
4078         mul64To128( zSig1, zSig1, &term2, &term3 );
4081             --zSig1;
4082             shortShift128Left( 0, zSig1, 1, &term2, &term3 );
4087         zSig1 |= ( ( rem1 | rem2 | rem3 ) != 0 );
4089     shortShift128Left( 0, zSig1, 1, &zSig0, &zSig1 );
4093             floatx80_rounding_precision, 0, zExp, zSig0, zSig1 );
4768     bits64 aSig0, aSig1, bSig0, bSig1, zSig0, zSig1, zSig2;
4815         add128( aSig0, aSig1, bSig0, bSig1, &zSig0, &zSig1 );
4816         if ( aExp == 0 ) return packFloat128( zSign, 0, zSig0, zSig1 );
4823     add128( aSig0, aSig1, bSig0, bSig1, &zSig0, &zSig1 );
4829         zSig0, zSig1, zSig2, 1, &zSig0, &zSig1, &zSig2 );
4831     return roundAndPackFloat128( zSign, zExp, zSig0, zSig1, zSig2 );
4847     bits64 aSig0, aSig1, bSig0, bSig1, zSig0, zSig1;
4894     sub128( bSig0, bSig1, aSig0, aSig1, &zSig0, &zSig1 );
4912     sub128( aSig0, aSig1, bSig0, bSig1, &zSig0, &zSig1 );
4916     return normalizeRoundAndPackFloat128( zSign, zExp - 14, zSig0, zSig1 );
4975     bits64 aSig0, aSig1, bSig0, bSig1, zSig0, zSig1, zSig2, zSig3;
5017     mul128To256( aSig0, aSig1, bSig0, bSig1, &zSig0, &zSig1, &zSig2, &zSig3 );
5018     add128( zSig0, zSig1, aSig0, aSig1, &zSig0, &zSig1 );
5022             zSig0, zSig1, zSig2, 1, &zSig0, &zSig1, &zSig2 );
5025     return roundAndPackFloat128( zSign, zExp, zSig0, zSig1, zSig2 );
5040     bits64 aSig0, aSig1, bSig0, bSig1, zSig0, zSig1, zSig2;
5099     zSig1 = estimateDiv128To64( rem1, rem2, bSig0 );
5100     if ( ( zSig1 & 0x3FFF ) <= 4 ) {
5101         mul128By64To192( bSig0, bSig1, zSig1, &term1, &term2, &term3 );
5104             --zSig1;
5107         zSig1 |= ( ( rem1 | rem2 | rem3 ) != 0 );
5109     shift128ExtraRightJamming( zSig0, zSig1, 0, 15, &zSig0, &zSig1, &zSig2 );
5110     return roundAndPackFloat128( zSign, zExp, zSig0, zSig1, zSig2 );
5236     bits64 aSig0, aSig1, zSig0, zSig1, zSig2, doubleZSig0;
5274     zSig1 = estimateDiv128To64( rem1, 0, doubleZSig0 );
5275     if ( ( zSig1 & 0x1FFF ) <= 5 ) {
5276         if ( zSig1 == 0 ) zSig1 = 1;
5277         mul64To128( doubleZSig0, zSig1, &term1, &term2 );
5279         mul64To128( zSig1, zSig1, &term2, &term3 );
5282             --zSig1;
5283             shortShift128Left( 0, zSig1, 1, &term2, &term3 );
5288         zSig1 |= ( ( rem1 | rem2 | rem3 ) != 0 );
5290     shift128ExtraRightJamming( zSig0, zSig1, 0, 14, &zSig0, &zSig1, &zSig2 );
5291     return roundAndPackFloat128( 0, zExp, zSig0, zSig1, zSig2 );