aboutsummaryrefslogtreecommitdiff
path: root/target/m68k/softfloat.c
diff options
context:
space:
mode:
Diffstat (limited to 'target/m68k/softfloat.c')
-rw-r--r--target/m68k/softfloat.c83
1 files changed, 0 insertions, 83 deletions
diff --git a/target/m68k/softfloat.c b/target/m68k/softfloat.c
index 9f120cf15e..b6d0ed7acf 100644
--- a/target/m68k/softfloat.c
+++ b/target/m68k/softfloat.c
@@ -43,89 +43,6 @@ static floatx80 propagateFloatx80NaNOneArg(floatx80 a, float_status *status)
}
/*
- * Returns the modulo remainder of the extended double-precision floating-point
- * value `a' with respect to the corresponding value `b'.
- */
-
-floatx80 floatx80_mod(floatx80 a, floatx80 b, float_status *status)
-{
- bool aSign, zSign;
- int32_t aExp, bExp, expDiff;
- uint64_t aSig0, aSig1, bSig;
- uint64_t qTemp, term0, term1;
-
- aSig0 = extractFloatx80Frac(a);
- aExp = extractFloatx80Exp(a);
- aSign = extractFloatx80Sign(a);
- bSig = extractFloatx80Frac(b);
- bExp = extractFloatx80Exp(b);
-
- if (aExp == 0x7FFF) {
- if ((uint64_t) (aSig0 << 1)
- || ((bExp == 0x7FFF) && (uint64_t) (bSig << 1))) {
- return propagateFloatx80NaN(a, b, status);
- }
- goto invalid;
- }
- if (bExp == 0x7FFF) {
- if ((uint64_t) (bSig << 1)) {
- return propagateFloatx80NaN(a, b, status);
- }
- return a;
- }
- if (bExp == 0) {
- if (bSig == 0) {
- invalid:
- float_raise(float_flag_invalid, status);
- return floatx80_default_nan(status);
- }
- normalizeFloatx80Subnormal(bSig, &bExp, &bSig);
- }
- if (aExp == 0) {
- if ((uint64_t) (aSig0 << 1) == 0) {
- return a;
- }
- normalizeFloatx80Subnormal(aSig0, &aExp, &aSig0);
- }
- bSig |= UINT64_C(0x8000000000000000);
- zSign = aSign;
- expDiff = aExp - bExp;
- aSig1 = 0;
- if (expDiff < 0) {
- return a;
- }
- qTemp = (bSig <= aSig0);
- if (qTemp) {
- aSig0 -= bSig;
- }
- expDiff -= 64;
- while (0 < expDiff) {
- qTemp = estimateDiv128To64(aSig0, aSig1, bSig);
- qTemp = (2 < qTemp) ? qTemp - 2 : 0;
- mul64To128(bSig, qTemp, &term0, &term1);
- sub128(aSig0, aSig1, term0, term1, &aSig0, &aSig1);
- shortShift128Left(aSig0, aSig1, 62, &aSig0, &aSig1);
- expDiff -= 62;
- }
- expDiff += 64;
- if (0 < expDiff) {
- qTemp = estimateDiv128To64(aSig0, aSig1, bSig);
- qTemp = (2 < qTemp) ? qTemp - 2 : 0;
- qTemp >>= 64 - expDiff;
- mul64To128(bSig, qTemp << (64 - expDiff), &term0, &term1);
- sub128(aSig0, aSig1, term0, term1, &aSig0, &aSig1);
- shortShift128Left(0, bSig, 64 - expDiff, &term0, &term1);
- while (le128(term0, term1, aSig0, aSig1)) {
- ++qTemp;
- sub128(aSig0, aSig1, term0, term1, &aSig0, &aSig1);
- }
- }
- return
- normalizeRoundAndPackFloatx80(
- 80, zSign, bExp + expDiff, aSig0, aSig1, status);
-}
-
-/*
* Returns the mantissa of the extended double-precision floating-point
* value `a'.
*/