diff options
author | Srikant Patnaik | 2015-01-11 12:28:04 +0530 |
---|---|---|
committer | Srikant Patnaik | 2015-01-11 12:28:04 +0530 |
commit | 871480933a1c28f8a9fed4c4d34d06c439a7a422 (patch) | |
tree | 8718f573808810c2a1e8cb8fb6ac469093ca2784 /ANDROID_3.4.5/arch/m68k/math-emu | |
parent | 9d40ac5867b9aefe0722bc1f110b965ff294d30d (diff) | |
download | FOSSEE-netbook-kernel-source-871480933a1c28f8a9fed4c4d34d06c439a7a422.tar.gz FOSSEE-netbook-kernel-source-871480933a1c28f8a9fed4c4d34d06c439a7a422.tar.bz2 FOSSEE-netbook-kernel-source-871480933a1c28f8a9fed4c4d34d06c439a7a422.zip |
Moved, renamed, and deleted files
The original directory structure was scattered and unorganized.
Changes are basically to make it look like kernel structure.
Diffstat (limited to 'ANDROID_3.4.5/arch/m68k/math-emu')
-rw-r--r-- | ANDROID_3.4.5/arch/m68k/math-emu/Makefile | 9 | ||||
-rw-r--r-- | ANDROID_3.4.5/arch/m68k/math-emu/fp_arith.c | 701 | ||||
-rw-r--r-- | ANDROID_3.4.5/arch/m68k/math-emu/fp_arith.h | 52 | ||||
-rw-r--r-- | ANDROID_3.4.5/arch/m68k/math-emu/fp_cond.S | 334 | ||||
-rw-r--r-- | ANDROID_3.4.5/arch/m68k/math-emu/fp_decode.h | 417 | ||||
-rw-r--r-- | ANDROID_3.4.5/arch/m68k/math-emu/fp_emu.h | 146 | ||||
-rw-r--r-- | ANDROID_3.4.5/arch/m68k/math-emu/fp_entry.S | 324 | ||||
-rw-r--r-- | ANDROID_3.4.5/arch/m68k/math-emu/fp_log.c | 219 | ||||
-rw-r--r-- | ANDROID_3.4.5/arch/m68k/math-emu/fp_move.S | 244 | ||||
-rw-r--r-- | ANDROID_3.4.5/arch/m68k/math-emu/fp_movem.S | 368 | ||||
-rw-r--r-- | ANDROID_3.4.5/arch/m68k/math-emu/fp_scan.S | 478 | ||||
-rw-r--r-- | ANDROID_3.4.5/arch/m68k/math-emu/fp_trig.c | 183 | ||||
-rw-r--r-- | ANDROID_3.4.5/arch/m68k/math-emu/fp_trig.h | 32 | ||||
-rw-r--r-- | ANDROID_3.4.5/arch/m68k/math-emu/fp_util.S | 1454 | ||||
-rw-r--r-- | ANDROID_3.4.5/arch/m68k/math-emu/multi_arith.h | 289 |
15 files changed, 0 insertions, 5250 deletions
diff --git a/ANDROID_3.4.5/arch/m68k/math-emu/Makefile b/ANDROID_3.4.5/arch/m68k/math-emu/Makefile deleted file mode 100644 index 547c23c6..00000000 --- a/ANDROID_3.4.5/arch/m68k/math-emu/Makefile +++ /dev/null @@ -1,9 +0,0 @@ -# -# Makefile for the linux kernel. -# - -#asflags-y := -DFPU_EMU_DEBUG -#ccflags-y := -DFPU_EMU_DEBUG - -obj-y := fp_entry.o fp_scan.o fp_util.o fp_move.o fp_movem.o \ - fp_cond.o fp_arith.o fp_log.o fp_trig.o diff --git a/ANDROID_3.4.5/arch/m68k/math-emu/fp_arith.c b/ANDROID_3.4.5/arch/m68k/math-emu/fp_arith.c deleted file mode 100644 index 08f286db..00000000 --- a/ANDROID_3.4.5/arch/m68k/math-emu/fp_arith.c +++ /dev/null @@ -1,701 +0,0 @@ -/* - - fp_arith.c: floating-point math routines for the Linux-m68k - floating point emulator. - - Copyright (c) 1998-1999 David Huggins-Daines. - - Somewhat based on the AlphaLinux floating point emulator, by David - Mosberger-Tang. - - You may copy, modify, and redistribute this file under the terms of - the GNU General Public License, version 2, or any later version, at - your convenience. - */ - -#include "fp_emu.h" -#include "multi_arith.h" -#include "fp_arith.h" - -const struct fp_ext fp_QNaN = -{ - .exp = 0x7fff, - .mant = { .m64 = ~0 } -}; - -const struct fp_ext fp_Inf = -{ - .exp = 0x7fff, -}; - -/* let's start with the easy ones */ - -struct fp_ext * -fp_fabs(struct fp_ext *dest, struct fp_ext *src) -{ - dprint(PINSTR, "fabs\n"); - - fp_monadic_check(dest, src); - - dest->sign = 0; - - return dest; -} - -struct fp_ext * -fp_fneg(struct fp_ext *dest, struct fp_ext *src) -{ - dprint(PINSTR, "fneg\n"); - - fp_monadic_check(dest, src); - - dest->sign = !dest->sign; - - return dest; -} - -/* Now, the slightly harder ones */ - -/* fp_fadd: Implements the kernel of the FADD, FSADD, FDADD, FSUB, - FDSUB, and FCMP instructions. */ - -struct fp_ext * -fp_fadd(struct fp_ext *dest, struct fp_ext *src) -{ - int diff; - - dprint(PINSTR, "fadd\n"); - - fp_dyadic_check(dest, src); - - if (IS_INF(dest)) { - /* infinity - infinity == NaN */ - if (IS_INF(src) && (src->sign != dest->sign)) - fp_set_nan(dest); - return dest; - } - if (IS_INF(src)) { - fp_copy_ext(dest, src); - return dest; - } - - if (IS_ZERO(dest)) { - if (IS_ZERO(src)) { - if (src->sign != dest->sign) { - if (FPDATA->rnd == FPCR_ROUND_RM) - dest->sign = 1; - else - dest->sign = 0; - } - } else - fp_copy_ext(dest, src); - return dest; - } - - dest->lowmant = src->lowmant = 0; - - if ((diff = dest->exp - src->exp) > 0) - fp_denormalize(src, diff); - else if ((diff = -diff) > 0) - fp_denormalize(dest, diff); - - if (dest->sign == src->sign) { - if (fp_addmant(dest, src)) - if (!fp_addcarry(dest)) - return dest; - } else { - if (dest->mant.m64 < src->mant.m64) { - fp_submant(dest, src, dest); - dest->sign = !dest->sign; - } else - fp_submant(dest, dest, src); - } - - return dest; -} - -/* fp_fsub: Implements the kernel of the FSUB, FSSUB, and FDSUB - instructions. - - Remember that the arguments are in assembler-syntax order! */ - -struct fp_ext * -fp_fsub(struct fp_ext *dest, struct fp_ext *src) -{ - dprint(PINSTR, "fsub "); - - src->sign = !src->sign; - return fp_fadd(dest, src); -} - - -struct fp_ext * -fp_fcmp(struct fp_ext *dest, struct fp_ext *src) -{ - dprint(PINSTR, "fcmp "); - - FPDATA->temp[1] = *dest; - src->sign = !src->sign; - return fp_fadd(&FPDATA->temp[1], src); -} - -struct fp_ext * -fp_ftst(struct fp_ext *dest, struct fp_ext *src) -{ - dprint(PINSTR, "ftst\n"); - - (void)dest; - - return src; -} - -struct fp_ext * -fp_fmul(struct fp_ext *dest, struct fp_ext *src) -{ - union fp_mant128 temp; - int exp; - - dprint(PINSTR, "fmul\n"); - - fp_dyadic_check(dest, src); - - /* calculate the correct sign now, as it's necessary for infinities */ - dest->sign = src->sign ^ dest->sign; - - /* Handle infinities */ - if (IS_INF(dest)) { - if (IS_ZERO(src)) - fp_set_nan(dest); - return dest; - } - if (IS_INF(src)) { - if (IS_ZERO(dest)) - fp_set_nan(dest); - else - fp_copy_ext(dest, src); - return dest; - } - - /* Of course, as we all know, zero * anything = zero. You may - not have known that it might be a positive or negative - zero... */ - if (IS_ZERO(dest) || IS_ZERO(src)) { - dest->exp = 0; - dest->mant.m64 = 0; - dest->lowmant = 0; - - return dest; - } - - exp = dest->exp + src->exp - 0x3ffe; - - /* shift up the mantissa for denormalized numbers, - so that the highest bit is set, this makes the - shift of the result below easier */ - if ((long)dest->mant.m32[0] >= 0) - exp -= fp_overnormalize(dest); - if ((long)src->mant.m32[0] >= 0) - exp -= fp_overnormalize(src); - - /* now, do a 64-bit multiply with expansion */ - fp_multiplymant(&temp, dest, src); - - /* normalize it back to 64 bits and stuff it back into the - destination struct */ - if ((long)temp.m32[0] > 0) { - exp--; - fp_putmant128(dest, &temp, 1); - } else - fp_putmant128(dest, &temp, 0); - - if (exp >= 0x7fff) { - fp_set_ovrflw(dest); - return dest; - } - dest->exp = exp; - if (exp < 0) { - fp_set_sr(FPSR_EXC_UNFL); - fp_denormalize(dest, -exp); - } - - return dest; -} - -/* fp_fdiv: Implements the "kernel" of the FDIV, FSDIV, FDDIV and - FSGLDIV instructions. - - Note that the order of the operands is counter-intuitive: instead - of src / dest, the result is actually dest / src. */ - -struct fp_ext * -fp_fdiv(struct fp_ext *dest, struct fp_ext *src) -{ - union fp_mant128 temp; - int exp; - - dprint(PINSTR, "fdiv\n"); - - fp_dyadic_check(dest, src); - - /* calculate the correct sign now, as it's necessary for infinities */ - dest->sign = src->sign ^ dest->sign; - - /* Handle infinities */ - if (IS_INF(dest)) { - /* infinity / infinity = NaN (quiet, as always) */ - if (IS_INF(src)) - fp_set_nan(dest); - /* infinity / anything else = infinity (with approprate sign) */ - return dest; - } - if (IS_INF(src)) { - /* anything / infinity = zero (with appropriate sign) */ - dest->exp = 0; - dest->mant.m64 = 0; - dest->lowmant = 0; - - return dest; - } - - /* zeroes */ - if (IS_ZERO(dest)) { - /* zero / zero = NaN */ - if (IS_ZERO(src)) - fp_set_nan(dest); - /* zero / anything else = zero */ - return dest; - } - if (IS_ZERO(src)) { - /* anything / zero = infinity (with appropriate sign) */ - fp_set_sr(FPSR_EXC_DZ); - dest->exp = 0x7fff; - dest->mant.m64 = 0; - - return dest; - } - - exp = dest->exp - src->exp + 0x3fff; - - /* shift up the mantissa for denormalized numbers, - so that the highest bit is set, this makes lots - of things below easier */ - if ((long)dest->mant.m32[0] >= 0) - exp -= fp_overnormalize(dest); - if ((long)src->mant.m32[0] >= 0) - exp -= fp_overnormalize(src); - - /* now, do the 64-bit divide */ - fp_dividemant(&temp, dest, src); - - /* normalize it back to 64 bits and stuff it back into the - destination struct */ - if (!temp.m32[0]) { - exp--; - fp_putmant128(dest, &temp, 32); - } else - fp_putmant128(dest, &temp, 31); - - if (exp >= 0x7fff) { - fp_set_ovrflw(dest); - return dest; - } - dest->exp = exp; - if (exp < 0) { - fp_set_sr(FPSR_EXC_UNFL); - fp_denormalize(dest, -exp); - } - - return dest; -} - -struct fp_ext * -fp_fsglmul(struct fp_ext *dest, struct fp_ext *src) -{ - int exp; - - dprint(PINSTR, "fsglmul\n"); - - fp_dyadic_check(dest, src); - - /* calculate the correct sign now, as it's necessary for infinities */ - dest->sign = src->sign ^ dest->sign; - - /* Handle infinities */ - if (IS_INF(dest)) { - if (IS_ZERO(src)) - fp_set_nan(dest); - return dest; - } - if (IS_INF(src)) { - if (IS_ZERO(dest)) - fp_set_nan(dest); - else - fp_copy_ext(dest, src); - return dest; - } - - /* Of course, as we all know, zero * anything = zero. You may - not have known that it might be a positive or negative - zero... */ - if (IS_ZERO(dest) || IS_ZERO(src)) { - dest->exp = 0; - dest->mant.m64 = 0; - dest->lowmant = 0; - - return dest; - } - - exp = dest->exp + src->exp - 0x3ffe; - - /* do a 32-bit multiply */ - fp_mul64(dest->mant.m32[0], dest->mant.m32[1], - dest->mant.m32[0] & 0xffffff00, - src->mant.m32[0] & 0xffffff00); - - if (exp >= 0x7fff) { - fp_set_ovrflw(dest); - return dest; - } - dest->exp = exp; - if (exp < 0) { - fp_set_sr(FPSR_EXC_UNFL); - fp_denormalize(dest, -exp); - } - - return dest; -} - -struct fp_ext * -fp_fsgldiv(struct fp_ext *dest, struct fp_ext *src) -{ - int exp; - unsigned long quot, rem; - - dprint(PINSTR, "fsgldiv\n"); - - fp_dyadic_check(dest, src); - - /* calculate the correct sign now, as it's necessary for infinities */ - dest->sign = src->sign ^ dest->sign; - - /* Handle infinities */ - if (IS_INF(dest)) { - /* infinity / infinity = NaN (quiet, as always) */ - if (IS_INF(src)) - fp_set_nan(dest); - /* infinity / anything else = infinity (with approprate sign) */ - return dest; - } - if (IS_INF(src)) { - /* anything / infinity = zero (with appropriate sign) */ - dest->exp = 0; - dest->mant.m64 = 0; - dest->lowmant = 0; - - return dest; - } - - /* zeroes */ - if (IS_ZERO(dest)) { - /* zero / zero = NaN */ - if (IS_ZERO(src)) - fp_set_nan(dest); - /* zero / anything else = zero */ - return dest; - } - if (IS_ZERO(src)) { - /* anything / zero = infinity (with appropriate sign) */ - fp_set_sr(FPSR_EXC_DZ); - dest->exp = 0x7fff; - dest->mant.m64 = 0; - - return dest; - } - - exp = dest->exp - src->exp + 0x3fff; - - dest->mant.m32[0] &= 0xffffff00; - src->mant.m32[0] &= 0xffffff00; - - /* do the 32-bit divide */ - if (dest->mant.m32[0] >= src->mant.m32[0]) { - fp_sub64(dest->mant, src->mant); - fp_div64(quot, rem, dest->mant.m32[0], 0, src->mant.m32[0]); - dest->mant.m32[0] = 0x80000000 | (quot >> 1); - dest->mant.m32[1] = (quot & 1) | rem; /* only for rounding */ - } else { - fp_div64(quot, rem, dest->mant.m32[0], 0, src->mant.m32[0]); - dest->mant.m32[0] = quot; - dest->mant.m32[1] = rem; /* only for rounding */ - exp--; - } - - if (exp >= 0x7fff) { - fp_set_ovrflw(dest); - return dest; - } - dest->exp = exp; - if (exp < 0) { - fp_set_sr(FPSR_EXC_UNFL); - fp_denormalize(dest, -exp); - } - - return dest; -} - -/* fp_roundint: Internal rounding function for use by several of these - emulated instructions. - - This one rounds off the fractional part using the rounding mode - specified. */ - -static void fp_roundint(struct fp_ext *dest, int mode) -{ - union fp_mant64 oldmant; - unsigned long mask; - - if (!fp_normalize_ext(dest)) - return; - - /* infinities and zeroes */ - if (IS_INF(dest) || IS_ZERO(dest)) - return; - - /* first truncate the lower bits */ - oldmant = dest->mant; - switch (dest->exp) { - case 0 ... 0x3ffe: - dest->mant.m64 = 0; - break; - case 0x3fff ... 0x401e: - dest->mant.m32[0] &= 0xffffffffU << (0x401e - dest->exp); - dest->mant.m32[1] = 0; - if (oldmant.m64 == dest->mant.m64) - return; - break; - case 0x401f ... 0x403e: - dest->mant.m32[1] &= 0xffffffffU << (0x403e - dest->exp); - if (oldmant.m32[1] == dest->mant.m32[1]) - return; - break; - default: - return; - } - fp_set_sr(FPSR_EXC_INEX2); - - /* We might want to normalize upwards here... however, since - we know that this is only called on the output of fp_fdiv, - or with the input to fp_fint or fp_fintrz, and the inputs - to all these functions are either normal or denormalized - (no subnormals allowed!), there's really no need. - - In the case of fp_fdiv, observe that 0x80000000 / 0xffff = - 0xffff8000, and the same holds for 128-bit / 64-bit. (i.e. the - smallest possible normal dividend and the largest possible normal - divisor will still produce a normal quotient, therefore, (normal - << 64) / normal is normal in all cases) */ - - switch (mode) { - case FPCR_ROUND_RN: - switch (dest->exp) { - case 0 ... 0x3ffd: - return; - case 0x3ffe: - /* As noted above, the input is always normal, so the - guard bit (bit 63) is always set. therefore, the - only case in which we will NOT round to 1.0 is when - the input is exactly 0.5. */ - if (oldmant.m64 == (1ULL << 63)) - return; - break; - case 0x3fff ... 0x401d: - mask = 1 << (0x401d - dest->exp); - if (!(oldmant.m32[0] & mask)) - return; - if (oldmant.m32[0] & (mask << 1)) - break; - if (!(oldmant.m32[0] << (dest->exp - 0x3ffd)) && - !oldmant.m32[1]) - return; - break; - case 0x401e: - if (!(oldmant.m32[1] >= 0)) - return; - if (oldmant.m32[0] & 1) - break; - if (!(oldmant.m32[1] << 1)) - return; - break; - case 0x401f ... 0x403d: - mask = 1 << (0x403d - dest->exp); - if (!(oldmant.m32[1] & mask)) - return; - if (oldmant.m32[1] & (mask << 1)) - break; - if (!(oldmant.m32[1] << (dest->exp - 0x401d))) - return; - break; - default: - return; - } - break; - case FPCR_ROUND_RZ: - return; - default: - if (dest->sign ^ (mode - FPCR_ROUND_RM)) - break; - return; - } - - switch (dest->exp) { - case 0 ... 0x3ffe: - dest->exp = 0x3fff; - dest->mant.m64 = 1ULL << 63; - break; - case 0x3fff ... 0x401e: - mask = 1 << (0x401e - dest->exp); - if (dest->mant.m32[0] += mask) - break; - dest->mant.m32[0] = 0x80000000; - dest->exp++; - break; - case 0x401f ... 0x403e: - mask = 1 << (0x403e - dest->exp); - if (dest->mant.m32[1] += mask) - break; - if (dest->mant.m32[0] += 1) - break; - dest->mant.m32[0] = 0x80000000; - dest->exp++; - break; - } -} - -/* modrem_kernel: Implementation of the FREM and FMOD instructions - (which are exactly the same, except for the rounding used on the - intermediate value) */ - -static struct fp_ext * -modrem_kernel(struct fp_ext *dest, struct fp_ext *src, int mode) -{ - struct fp_ext tmp; - - fp_dyadic_check(dest, src); - - /* Infinities and zeros */ - if (IS_INF(dest) || IS_ZERO(src)) { - fp_set_nan(dest); - return dest; - } - if (IS_ZERO(dest) || IS_INF(src)) - return dest; - - /* FIXME: there is almost certainly a smarter way to do this */ - fp_copy_ext(&tmp, dest); - fp_fdiv(&tmp, src); /* NOTE: src might be modified */ - fp_roundint(&tmp, mode); - fp_fmul(&tmp, src); - fp_fsub(dest, &tmp); - - /* set the quotient byte */ - fp_set_quotient((dest->mant.m64 & 0x7f) | (dest->sign << 7)); - return dest; -} - -/* fp_fmod: Implements the kernel of the FMOD instruction. - - Again, the argument order is backwards. The result, as defined in - the Motorola manuals, is: - - fmod(src,dest) = (dest - (src * floor(dest / src))) */ - -struct fp_ext * -fp_fmod(struct fp_ext *dest, struct fp_ext *src) -{ - dprint(PINSTR, "fmod\n"); - return modrem_kernel(dest, src, FPCR_ROUND_RZ); -} - -/* fp_frem: Implements the kernel of the FREM instruction. - - frem(src,dest) = (dest - (src * round(dest / src))) - */ - -struct fp_ext * -fp_frem(struct fp_ext *dest, struct fp_ext *src) -{ - dprint(PINSTR, "frem\n"); - return modrem_kernel(dest, src, FPCR_ROUND_RN); -} - -struct fp_ext * -fp_fint(struct fp_ext *dest, struct fp_ext *src) -{ - dprint(PINSTR, "fint\n"); - - fp_copy_ext(dest, src); - - fp_roundint(dest, FPDATA->rnd); - - return dest; -} - -struct fp_ext * -fp_fintrz(struct fp_ext *dest, struct fp_ext *src) -{ - dprint(PINSTR, "fintrz\n"); - - fp_copy_ext(dest, src); - - fp_roundint(dest, FPCR_ROUND_RZ); - - return dest; -} - -struct fp_ext * -fp_fscale(struct fp_ext *dest, struct fp_ext *src) -{ - int scale, oldround; - - dprint(PINSTR, "fscale\n"); - - fp_dyadic_check(dest, src); - - /* Infinities */ - if (IS_INF(src)) { - fp_set_nan(dest); - return dest; - } - if (IS_INF(dest)) - return dest; - - /* zeroes */ - if (IS_ZERO(src) || IS_ZERO(dest)) - return dest; - - /* Source exponent out of range */ - if (src->exp >= 0x400c) { - fp_set_ovrflw(dest); - return dest; - } - - /* src must be rounded with round to zero. */ - oldround = FPDATA->rnd; - FPDATA->rnd = FPCR_ROUND_RZ; - scale = fp_conv_ext2long(src); - FPDATA->rnd = oldround; - - /* new exponent */ - scale += dest->exp; - - if (scale >= 0x7fff) { - fp_set_ovrflw(dest); - } else if (scale <= 0) { - fp_set_sr(FPSR_EXC_UNFL); - fp_denormalize(dest, -scale); - } else - dest->exp = scale; - - return dest; -} - diff --git a/ANDROID_3.4.5/arch/m68k/math-emu/fp_arith.h b/ANDROID_3.4.5/arch/m68k/math-emu/fp_arith.h deleted file mode 100644 index 2cc3f846..00000000 --- a/ANDROID_3.4.5/arch/m68k/math-emu/fp_arith.h +++ /dev/null @@ -1,52 +0,0 @@ -/* - - fp_arith.h: floating-point math routines for the Linux-m68k - floating point emulator. - - Copyright (c) 1998 David Huggins-Daines. - - Somewhat based on the AlphaLinux floating point emulator, by David - Mosberger-Tang. - - You may copy, modify, and redistribute this file under the terms of - the GNU General Public License, version 2, or any later version, at - your convenience. - - */ - -#ifndef FP_ARITH_H -#define FP_ARITH_H - -/* easy ones */ -struct fp_ext * -fp_fabs(struct fp_ext *dest, struct fp_ext *src); -struct fp_ext * -fp_fneg(struct fp_ext *dest, struct fp_ext *src); - -/* straightforward arithmetic */ -struct fp_ext * -fp_fadd(struct fp_ext *dest, struct fp_ext *src); -struct fp_ext * -fp_fsub(struct fp_ext *dest, struct fp_ext *src); -struct fp_ext * -fp_fcmp(struct fp_ext *dest, struct fp_ext *src); -struct fp_ext * -fp_ftst(struct fp_ext *dest, struct fp_ext *src); -struct fp_ext * -fp_fmul(struct fp_ext *dest, struct fp_ext *src); -struct fp_ext * -fp_fdiv(struct fp_ext *dest, struct fp_ext *src); - -/* ones that do rounding and integer conversions */ -struct fp_ext * -fp_fmod(struct fp_ext *dest, struct fp_ext *src); -struct fp_ext * -fp_frem(struct fp_ext *dest, struct fp_ext *src); -struct fp_ext * -fp_fint(struct fp_ext *dest, struct fp_ext *src); -struct fp_ext * -fp_fintrz(struct fp_ext *dest, struct fp_ext *src); -struct fp_ext * -fp_fscale(struct fp_ext *dest, struct fp_ext *src); - -#endif /* FP_ARITH__H */ diff --git a/ANDROID_3.4.5/arch/m68k/math-emu/fp_cond.S b/ANDROID_3.4.5/arch/m68k/math-emu/fp_cond.S deleted file mode 100644 index ddae8b1b..00000000 --- a/ANDROID_3.4.5/arch/m68k/math-emu/fp_cond.S +++ /dev/null @@ -1,334 +0,0 @@ -/* - * fp_cond.S - * - * Copyright Roman Zippel, 1997. All rights reserved. - * - * Redistribution and use in source and binary forms, with or without - * modification, are permitted provided that the following conditions - * are met: - * 1. Redistributions of source code must retain the above copyright - * notice, and the entire permission notice in its entirety, - * including the disclaimer of warranties. - * 2. Redistributions in binary form must reproduce the above copyright - * notice, this list of conditions and the following disclaimer in the - * documentation and/or other materials provided with the distribution. - * 3. The name of the author may not be used to endorse or promote - * products derived from this software without specific prior - * written permission. - * - * ALTERNATIVELY, this product may be distributed under the terms of - * the GNU General Public License, in which case the provisions of the GPL are - * required INSTEAD OF the above restrictions. (This clause is - * necessary due to a potential bad interaction between the GPL and - * the restrictions contained in a BSD-style copyright.) - * - * THIS SOFTWARE IS PROVIDED ``AS IS'' AND ANY EXPRESS OR IMPLIED - * WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES - * OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE - * DISCLAIMED. IN NO EVENT SHALL THE AUTHOR BE LIABLE FOR ANY DIRECT, - * INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES - * (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR - * SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) - * HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, - * STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) - * ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED - * OF THE POSSIBILITY OF SUCH DAMAGE. - */ - -#include "fp_emu.h" -#include "fp_decode.h" - - .globl fp_fscc, fp_fbccw, fp_fbccl - -#ifdef FPU_EMU_DEBUG -fp_fnop: - printf PDECODE,"fnop\n" - jra fp_end -#else -#define fp_fnop fp_end -#endif - -fp_fbccw: - tst.w %d2 - jeq fp_fnop - printf PDECODE,"fbccw " - fp_get_pc %a0 - lea (-2,%a0,%d2.w),%a0 - jra 1f - -fp_fbccl: - printf PDECODE,"fbccl " - fp_get_pc %a0 - move.l %d2,%d0 - swap %d0 - fp_get_instr_word %d0,fp_err_ua1 - lea (-2,%a0,%d0.l),%a0 -1: printf PDECODE,"%x",1,%a0 - move.l %d2,%d0 - swap %d0 - jsr fp_compute_cond - tst.l %d0 - jeq 1f - fp_put_pc %a0,1 -1: printf PDECODE,"\n" - jra fp_end - -fp_fdbcc: - printf PDECODE,"fdbcc " - fp_get_pc %a1 | calculate new pc - fp_get_instr_word %d0,fp_err_ua1 - add.w %d0,%a1 - fp_decode_addr_reg - printf PDECODE,"d%d,%x\n",2,%d0,%a1 - swap %d1 | test condition in %d1 - tst.w %d1 - jne 2f - move.l %d0,%d1 - jsr fp_get_data_reg - subq.w #1,%d0 - jcs 1f - fp_put_pc %a1,1 -1: jsr fp_put_data_reg -2: jra fp_end - -| set flags for decode macros for fs<cc> -do_fscc=1 -do_no_pc_mode=1 - -fp_fscc: - printf PDECODE,"fscc " - move.l %d2,%d0 - jsr fp_compute_cond - move.w %d0,%d1 - swap %d1 - - | decode addressing mode - fp_decode_addr_mode - - .long fp_data, fp_fdbcc - .long fp_indirect, fp_postinc - .long fp_predecr, fp_disp16 - .long fp_extmode0, fp_extmode1 - - | addressing mode: data register direct -fp_data: - fp_mode_data_direct - move.w %d0,%d1 | save register nr - jsr fp_get_data_reg - swap %d1 - move.b %d1,%d0 - swap %d1 - jsr fp_put_data_reg - printf PDECODE,"\n" - jra fp_end - -fp_indirect: - fp_mode_addr_indirect - jra fp_do_scc - -fp_postinc: - fp_mode_addr_indirect_postinc - jra fp_do_scc - -fp_predecr: - fp_mode_addr_indirect_predec - jra fp_do_scc - -fp_disp16: - fp_mode_addr_indirect_disp16 - jra fp_do_scc - -fp_extmode0: - fp_mode_addr_indirect_extmode0 - jra fp_do_scc - -fp_extmode1: - bfextu %d2{#13,#3},%d0 - jmp ([0f:w,%pc,%d0*4]) - - .align 4 -0: - .long fp_absolute_short, fp_absolute_long - .long fp_ill, fp_ill | NOTE: jump here to ftrap.x - .long fp_ill, fp_ill - .long fp_ill, fp_ill - -fp_absolute_short: - fp_mode_abs_short - jra fp_do_scc - -fp_absolute_long: - fp_mode_abs_long -| jra fp_do_scc - -fp_do_scc: - swap %d1 - putuser.b %d1,(%a0),fp_err_ua1,%a0 - printf PDECODE,"\n" - jra fp_end - - -#define tst_NAN btst #24,%d1 -#define tst_Z btst #26,%d1 -#define tst_N btst #27,%d1 - -fp_compute_cond: - move.l (FPD_FPSR,FPDATA),%d1 - btst #4,%d0 - jeq 1f - tst_NAN - jeq 1f - bset #15,%d1 - bset #7,%d1 - move.l %d1,(FPD_FPSR,FPDATA) -1: and.w #0xf,%d0 - jmp ([0f:w,%pc,%d0.w*4]) - - .align 4 -0: - .long fp_f , fp_eq , fp_ogt, fp_oge - .long fp_olt, fp_ole, fp_ogl, fp_or - .long fp_un , fp_ueq, fp_ugt, fp_uge - .long fp_ult, fp_ule, fp_ne , fp_t - -fp_f: - moveq #0,%d0 - rts - -fp_eq: - moveq #0,%d0 - tst_Z - jeq 1f - moveq #-1,%d0 -1: rts - -fp_ogt: - moveq #0,%d0 - tst_NAN - jne 1f - tst_Z - jne 1f - tst_N - jne 1f - moveq #-1,%d0 -1: rts - -fp_oge: - moveq #-1,%d0 - tst_Z - jne 2f - tst_NAN - jne 1f - tst_N - jeq 2f -1: moveq #0,%d0 -2: rts - -fp_olt: - moveq #0,%d0 - tst_NAN - jne 1f - tst_Z - jne 1f - tst_N - jeq 1f - moveq #-1,%d0 -1: rts - -fp_ole: - moveq #-1,%d0 - tst_Z - jne 2f - tst_NAN - jne 1f - tst_N - jne 2f -1: moveq #0,%d0 -2: rts - -fp_ogl: - moveq #0,%d0 - tst_NAN - jne 1f - tst_Z - jne 1f - moveq #-1,%d0 -1: rts - -fp_or: - moveq #0,%d0 - tst_NAN - jne 1f - moveq #-1,%d0 -1: rts - -fp_un: - moveq #0,%d0 - tst_NAN - jeq 1f - moveq #-1,%d0 - rts - -fp_ueq: - moveq #-1,%d0 - tst_NAN - jne 1f - tst_Z - jne 1f - moveq #0,%d0 -1: rts - -fp_ugt: - moveq #-1,%d0 - tst_NAN - jne 2f - tst_N - jne 1f - tst_Z - jeq 2f -1: moveq #0,%d0 -2: rts - -fp_uge: - moveq #-1,%d0 - tst_NAN - jne 1f - tst_Z - jne 1f - tst_N - jeq 1f - moveq #0,%d0 -1: rts - -fp_ult: - moveq #-1,%d0 - tst_NAN - jne 2f - tst_Z - jne 1f - tst_N - jne 2f -1: moveq #0,%d0 -2: rts - -fp_ule: - moveq #-1,%d0 - tst_NAN - jne 1f - tst_Z - jne 1f - tst_N - jne 1f - moveq #0,%d0 -1: rts - -fp_ne: - moveq #0,%d0 - tst_Z - jne 1f - moveq #-1,%d0 -1: rts - -fp_t: - moveq #-1,%d0 - rts diff --git a/ANDROID_3.4.5/arch/m68k/math-emu/fp_decode.h b/ANDROID_3.4.5/arch/m68k/math-emu/fp_decode.h deleted file mode 100644 index 759679d9..00000000 --- a/ANDROID_3.4.5/arch/m68k/math-emu/fp_decode.h +++ /dev/null @@ -1,417 +0,0 @@ -/* - * fp_decode.h - * - * Copyright Roman Zippel, 1997. All rights reserved. - * - * Redistribution and use in source and binary forms, with or without - * modification, are permitted provided that the following conditions - * are met: - * 1. Redistributions of source code must retain the above copyright - * notice, and the entire permission notice in its entirety, - * including the disclaimer of warranties. - * 2. Redistributions in binary form must reproduce the above copyright - * notice, this list of conditions and the following disclaimer in the - * documentation and/or other materials provided with the distribution. - * 3. The name of the author may not be used to endorse or promote - * products derived from this software without specific prior - * written permission. - * - * ALTERNATIVELY, this product may be distributed under the terms of - * the GNU General Public License, in which case the provisions of the GPL are - * required INSTEAD OF the above restrictions. (This clause is - * necessary due to a potential bad interaction between the GPL and - * the restrictions contained in a BSD-style copyright.) - * - * THIS SOFTWARE IS PROVIDED ``AS IS'' AND ANY EXPRESS OR IMPLIED - * WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES - * OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE - * DISCLAIMED. IN NO EVENT SHALL THE AUTHOR BE LIABLE FOR ANY DIRECT, - * INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES - * (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR - * SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) - * HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, - * STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) - * ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED - * OF THE POSSIBILITY OF SUCH DAMAGE. - */ - -#ifndef _FP_DECODE_H -#define _FP_DECODE_H - -/* These macros do the dirty work of the instr decoding, several variables - * can be defined in the source file to modify the work of these macros, - * currently the following variables are used: - * ... - * The register usage: - * d0 - will contain source operand for data direct mode, - * otherwise scratch register - * d1 - upper 16bit are reserved for caller - * lower 16bit may contain further arguments, - * is destroyed during decoding - * d2 - contains first two instruction words, - * first word will be used for extension word - * a0 - will point to source/dest operand for any indirect mode - * otherwise scratch register - * a1 - scratch register - * a2 - base addr to the task structure - * - * the current implementation doesn't check for every disallowed - * addressing mode (e.g. pc relative modes as destination), as long - * as it only means a new addressing mode, which should not appear - * in a program and that doesn't crash the emulation, I think it's - * not a problem to allow these modes. - */ - -do_fmovem=0 -do_fmovem_cr=0 -do_no_pc_mode=0 -do_fscc=0 - -| first decoding of the instr type -| this separates the conditional instr -.macro fp_decode_cond_instr_type - bfextu %d2{#8,#2},%d0 - jmp ([0f:w,%pc,%d0*4]) - - .align 4 -0: -| .long "f<op>","fscc/fdbcc" -| .long "fbccw","fbccl" -.endm - -| second decoding of the instr type -| this separates most move instr -.macro fp_decode_move_instr_type - bfextu %d2{#16,#3},%d0 - jmp ([0f:w,%pc,%d0*4]) - - .align 4 -0: -| .long "f<op> fpx,fpx","invalid instr" -| .long "f<op> <ea>,fpx","fmove fpx,<ea>" -| .long "fmovem <ea>,fpcr","fmovem <ea>,fpx" -| .long "fmovem fpcr,<ea>","fmovem fpx,<ea>" -.endm - -| extract the source specifier, specifies -| either source fp register or data format -.macro fp_decode_sourcespec - bfextu %d2{#19,#3},%d0 -.endm - -| decode destination format for fmove reg,ea -.macro fp_decode_dest_format - bfextu %d2{#19,#3},%d0 -.endm - -| decode source register for fmove reg,ea -.macro fp_decode_src_reg - bfextu %d2{#22,#3},%d0 -.endm - -| extract the addressing mode -| it depends on the instr which of the modes is valid -.macro fp_decode_addr_mode - bfextu %d2{#10,#3},%d0 - jmp ([0f:w,%pc,%d0*4]) - - .align 4 -0: -| .long "data register direct","addr register direct" -| .long "addr register indirect" -| .long "addr register indirect postincrement" -| .long "addr register indirect predecrement" -| .long "addr register + index16" -| .long "extension mode1","extension mode2" -.endm - -| extract the register for the addressing mode -.macro fp_decode_addr_reg - bfextu %d2{#13,#3},%d0 -.endm - -| decode the 8bit diplacement from the brief extension word -.macro fp_decode_disp8 - move.b %d2,%d0 - ext.w %d0 -.endm - -| decode the index of the brief/full extension word -.macro fp_decode_index - bfextu %d2{#17,#3},%d0 | get the register nr - btst #15,%d2 | test for data/addr register - jne 1\@f - printf PDECODE,"d%d",1,%d0 - jsr fp_get_data_reg - jra 2\@f -1\@: printf PDECODE,"a%d",1,%d0 - jsr fp_get_addr_reg - move.l %a0,%d0 -2\@: -debug lea "'l'.w,%a0" - btst #11,%d2 | 16/32 bit size? - jne 3\@f -debug lea "'w'.w,%a0" - ext.l %d0 -3\@: printf PDECODE,":%c",1,%a0 - move.w %d2,%d1 | scale factor - rol.w #7,%d1 - and.w #3,%d1 -debug move.l "%d1,-(%sp)" -debug ext.l "%d1" - printf PDECODE,":%d",1,%d1 -debug move.l "(%sp)+,%d1" - lsl.l %d1,%d0 -.endm - -| decode the base displacement size -.macro fp_decode_basedisp - bfextu %d2{#26,#2},%d0 - jmp ([0f:w,%pc,%d0*4]) - - .align 4 -0: -| .long "reserved","null displacement" -| .long "word displacement","long displacement" -.endm - -.macro fp_decode_outerdisp - bfextu %d2{#30,#2},%d0 - jmp ([0f:w,%pc,%d0*4]) - - .align 4 -0: -| .long "no memory indirect action/reserved","null outer displacement" -| .long "word outer displacement","long outer displacement" -.endm - -| get the extension word and test for brief or full extension type -.macro fp_get_test_extword label - fp_get_instr_word %d2,fp_err_ua1 - btst #8,%d2 - jne \label -.endm - - -| test if %pc is the base register for the indirect addr mode -.macro fp_test_basereg_d16 label - btst #20,%d2 - jeq \label -.endm - -| test if %pc is the base register for one of the extended modes -.macro fp_test_basereg_ext label - btst #19,%d2 - jeq \label -.endm - -.macro fp_test_suppr_index label - btst #6,%d2 - jne \label -.endm - - -| addressing mode: data register direct -.macro fp_mode_data_direct - fp_decode_addr_reg - printf PDECODE,"d%d",1,%d0 -.endm - -| addressing mode: address register indirect -.macro fp_mode_addr_indirect - fp_decode_addr_reg - printf PDECODE,"(a%d)",1,%d0 - jsr fp_get_addr_reg -.endm - -| adjust stack for byte moves from/to stack -.macro fp_test_sp_byte_move - .if !do_fmovem - .if do_fscc - move.w #6,%d1 - .endif - cmp.w #7,%d0 - jne 1\@f - .if !do_fscc - cmp.w #6,%d1 - jne 1\@f - .endif - move.w #4,%d1 -1\@: - .endif -.endm - -| addressing mode: address register indirect with postincrement -.macro fp_mode_addr_indirect_postinc - fp_decode_addr_reg - printf PDECODE,"(a%d)+",1,%d0 - fp_test_sp_byte_move - jsr fp_get_addr_reg - move.l %a0,%a1 | save addr - .if do_fmovem - lea (%a0,%d1.w*4),%a0 - .if !do_fmovem_cr - lea (%a0,%d1.w*8),%a0 - .endif - .else - add.w (fp_datasize,%d1.w*2),%a0 - .endif - jsr fp_put_addr_reg - move.l %a1,%a0 -.endm - -| addressing mode: address register indirect with predecrement -.macro fp_mode_addr_indirect_predec - fp_decode_addr_reg - printf PDECODE,"-(a%d)",1,%d0 - fp_test_sp_byte_move - jsr fp_get_addr_reg - .if do_fmovem - .if !do_fmovem_cr - lea (-12,%a0),%a1 | setup to addr of 1st reg to move - neg.w %d1 - lea (%a0,%d1.w*4),%a0 - add.w %d1,%d1 - lea (%a0,%d1.w*4),%a0 - jsr fp_put_addr_reg - move.l %a1,%a0 - .else - neg.w %d1 - lea (%a0,%d1.w*4),%a0 - jsr fp_put_addr_reg - .endif - .else - sub.w (fp_datasize,%d1.w*2),%a0 - jsr fp_put_addr_reg - .endif -.endm - -| addressing mode: address register/programm counter indirect -| with 16bit displacement -.macro fp_mode_addr_indirect_disp16 - .if !do_no_pc_mode - fp_test_basereg_d16 1f - printf PDECODE,"pc" - fp_get_pc %a0 - jra 2f - .endif -1: fp_decode_addr_reg - printf PDECODE,"a%d",1,%d0 - jsr fp_get_addr_reg -2: fp_get_instr_word %a1,fp_err_ua1 - printf PDECODE,"@(%x)",1,%a1 - add.l %a1,%a0 -.endm - -| perform preindex (if I/IS == 0xx and xx != 00) -.macro fp_do_preindex - moveq #3,%d0 - and.w %d2,%d0 - jeq 1f - btst #2,%d2 - jne 1f - printf PDECODE,")@(" - getuser.l (%a1),%a1,fp_err_ua1,%a1 -debug jra "2f" -1: printf PDECODE,"," -2: -.endm - -| perform postindex (if I/IS == 1xx) -.macro fp_do_postindex - btst #2,%d2 - jeq 1f - printf PDECODE,")@(" - getuser.l (%a1),%a1,fp_err_ua1,%a1 -debug jra "2f" -1: printf PDECODE,"," -2: -.endm - -| all other indirect addressing modes will finally end up here -.macro fp_mode_addr_indirect_extmode0 - .if !do_no_pc_mode - fp_test_basereg_ext 1f - printf PDECODE,"pc" - fp_get_pc %a0 - jra 2f - .endif -1: fp_decode_addr_reg - printf PDECODE,"a%d",1,%d0 - jsr fp_get_addr_reg -2: move.l %a0,%a1 - swap %d2 - fp_get_test_extword 3f - | addressing mode: address register/programm counter indirect - | with index and 8bit displacement - fp_decode_disp8 -debug ext.l "%d0" - printf PDECODE,"@(%x,",1,%d0 - add.w %d0,%a1 - fp_decode_index - add.l %d0,%a1 - printf PDECODE,")" - jra 9f -3: | addressing mode: address register/programm counter memory indirect - | with base and/or outer displacement - btst #7,%d2 | base register suppressed? - jeq 1f - printf PDECODE,"!" - sub.l %a1,%a1 -1: printf PDECODE,"@(" - fp_decode_basedisp - - .long fp_ill,1f - .long 2f,3f - -#ifdef FPU_EMU_DEBUG -1: printf PDECODE,"0" | null base displacement - jra 1f -#endif -2: fp_get_instr_word %a0,fp_err_ua1 | 16bit base displacement - printf PDECODE,"%x:w",1,%a0 - jra 4f -3: fp_get_instr_long %a0,fp_err_ua1 | 32bit base displacement - printf PDECODE,"%x:l",1,%a0 -4: add.l %a0,%a1 -1: - fp_do_postindex - fp_test_suppr_index 1f - fp_decode_index - add.l %d0,%a1 -1: fp_do_preindex - - fp_decode_outerdisp - - .long 5f,1f - .long 2f,3f - -#ifdef FPU_EMU_DEBUG -1: printf PDECODE,"0" | null outer displacement - jra 1f -#endif -2: fp_get_instr_word %a0,fp_err_ua1 | 16bit outer displacement - printf PDECODE,"%x:w",1,%a0 - jra 4f -3: fp_get_instr_long %a0,fp_err_ua1 | 32bit outer displacement - printf PDECODE,"%x:l",1,%a0 -4: add.l %a0,%a1 -1: -5: printf PDECODE,")" -9: move.l %a1,%a0 - swap %d2 -.endm - -| get the absolute short address from user space -.macro fp_mode_abs_short - fp_get_instr_word %a0,fp_err_ua1 - printf PDECODE,"%x.w",1,%a0 -.endm - -| get the absolute long address from user space -.macro fp_mode_abs_long - fp_get_instr_long %a0,fp_err_ua1 - printf PDECODE,"%x.l",1,%a0 -.endm - -#endif /* _FP_DECODE_H */ diff --git a/ANDROID_3.4.5/arch/m68k/math-emu/fp_emu.h b/ANDROID_3.4.5/arch/m68k/math-emu/fp_emu.h deleted file mode 100644 index c1ecfef7..00000000 --- a/ANDROID_3.4.5/arch/m68k/math-emu/fp_emu.h +++ /dev/null @@ -1,146 +0,0 @@ -/* - * fp_emu.h - * - * Copyright Roman Zippel, 1997. All rights reserved. - * - * Redistribution and use in source and binary forms, with or without - * modification, are permitted provided that the following conditions - * are met: - * 1. Redistributions of source code must retain the above copyright - * notice, and the entire permission notice in its entirety, - * including the disclaimer of warranties. - * 2. Redistributions in binary form must reproduce the above copyright - * notice, this list of conditions and the following disclaimer in the - * documentation and/or other materials provided with the distribution. - * 3. The name of the author may not be used to endorse or promote - * products derived from this software without specific prior - * written permission. - * - * ALTERNATIVELY, this product may be distributed under the terms of - * the GNU General Public License, in which case the provisions of the GPL are - * required INSTEAD OF the above restrictions. (This clause is - * necessary due to a potential bad interaction between the GPL and - * the restrictions contained in a BSD-style copyright.) - * - * THIS SOFTWARE IS PROVIDED ``AS IS'' AND ANY EXPRESS OR IMPLIED - * WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES - * OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE - * DISCLAIMED. IN NO EVENT SHALL THE AUTHOR BE LIABLE FOR ANY DIRECT, - * INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES - * (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR - * SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) - * HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, - * STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) - * ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED - * OF THE POSSIBILITY OF SUCH DAMAGE. - */ - -#ifndef _FP_EMU_H -#define _FP_EMU_H - -#ifdef __ASSEMBLY__ -#include <asm/asm-offsets.h> -#endif -#include <asm/math-emu.h> - -#ifndef __ASSEMBLY__ - -#define IS_INF(a) ((a)->exp == 0x7fff) -#define IS_ZERO(a) ((a)->mant.m64 == 0) - - -#define fp_set_sr(bit) ({ \ - FPDATA->fpsr |= 1 << (bit); \ -}) - -#define fp_set_quotient(quotient) ({ \ - FPDATA->fpsr &= 0xff00ffff; \ - FPDATA->fpsr |= ((quotient) & 0xff) << 16; \ -}) - -/* linkage for several useful functions */ - -/* Normalize the extended struct, return 0 for a NaN */ -#define fp_normalize_ext(fpreg) ({ \ - register struct fp_ext *reg asm ("a0") = fpreg; \ - register int res asm ("d0"); \ - \ - asm volatile ("jsr fp_conv_ext2ext" \ - : "=d" (res) : "a" (reg) \ - : "a1", "d1", "d2", "memory"); \ - res; \ -}) - -#define fp_copy_ext(dest, src) ({ \ - *dest = *src; \ -}) - -#define fp_monadic_check(dest, src) ({ \ - fp_copy_ext(dest, src); \ - if (!fp_normalize_ext(dest)) \ - return dest; \ -}) - -#define fp_dyadic_check(dest, src) ({ \ - if (!fp_normalize_ext(dest)) \ - return dest; \ - if (!fp_normalize_ext(src)) { \ - fp_copy_ext(dest, src); \ - return dest; \ - } \ -}) - -extern const struct fp_ext fp_QNaN; -extern const struct fp_ext fp_Inf; - -#define fp_set_nan(dest) ({ \ - fp_set_sr(FPSR_EXC_OPERR); \ - *dest = fp_QNaN; \ -}) - -/* TODO check rounding mode? */ -#define fp_set_ovrflw(dest) ({ \ - fp_set_sr(FPSR_EXC_OVFL); \ - dest->exp = 0x7fff; \ - dest->mant.m64 = 0; \ -}) - -#define fp_conv_ext2long(src) ({ \ - register struct fp_ext *__src asm ("a0") = src; \ - register int __res asm ("d0"); \ - \ - asm volatile ("jsr fp_conv_ext2long" \ - : "=d" (__res) : "a" (__src) \ - : "a1", "d1", "d2", "memory"); \ - __res; \ -}) - -#define fp_conv_long2ext(dest, src) ({ \ - register struct fp_ext *__dest asm ("a0") = dest; \ - register int __src asm ("d0") = src; \ - \ - asm volatile ("jsr fp_conv_ext2long" \ - : : "d" (__src), "a" (__dest) \ - : "a1", "d1", "d2", "memory"); \ -}) - -#else /* __ASSEMBLY__ */ - -/* - * set, reset or clear a bit in the fp status register - */ -.macro fp_set_sr bit - bset #(\bit&7),(FPD_FPSR+3-(\bit/8),FPDATA) -.endm - -.macro fp_clr_sr bit - bclr #(\bit&7),(FPD_FPSR+3-(\bit/8),FPDATA) -.endm - -.macro fp_tst_sr bit - btst #(\bit&7),(FPD_FPSR+3-(\bit/8),FPDATA) -.endm - -#endif /* __ASSEMBLY__ */ - -#endif /* _FP_EMU_H */ diff --git a/ANDROID_3.4.5/arch/m68k/math-emu/fp_entry.S b/ANDROID_3.4.5/arch/m68k/math-emu/fp_entry.S deleted file mode 100644 index a3fe1f34..00000000 --- a/ANDROID_3.4.5/arch/m68k/math-emu/fp_entry.S +++ /dev/null @@ -1,324 +0,0 @@ -/* - * fp_emu.S - * - * Copyright Roman Zippel, 1997. All rights reserved. - * - * Redistribution and use in source and binary forms, with or without - * modification, are permitted provided that the following conditions - * are met: - * 1. Redistributions of source code must retain the above copyright - * notice, and the entire permission notice in its entirety, - * including the disclaimer of warranties. - * 2. Redistributions in binary form must reproduce the above copyright - * notice, this list of conditions and the following disclaimer in the - * documentation and/or other materials provided with the distribution. - * 3. The name of the author may not be used to endorse or promote - * products derived from this software without specific prior - * written permission. - * - * ALTERNATIVELY, this product may be distributed under the terms of - * the GNU General Public License, in which case the provisions of the GPL are - * required INSTEAD OF the above restrictions. (This clause is - * necessary due to a potential bad interaction between the GPL and - * the restrictions contained in a BSD-style copyright.) - * - * THIS SOFTWARE IS PROVIDED ``AS IS'' AND ANY EXPRESS OR IMPLIED - * WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES - * OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE - * DISCLAIMED. IN NO EVENT SHALL THE AUTHOR BE LIABLE FOR ANY DIRECT, - * INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES - * (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR - * SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) - * HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, - * STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) - * ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED - * OF THE POSSIBILITY OF SUCH DAMAGE. - */ - -#include <linux/linkage.h> -#include <asm/entry.h> - -#include "fp_emu.h" - - .globl fpu_emu - .globl fp_debugprint - .globl fp_err_ua1,fp_err_ua2 - - .text -fpu_emu: - SAVE_ALL_INT - GET_CURRENT(%d0) - -#if defined(CPU_M68020_OR_M68030) && defined(CPU_M68040_OR_M68060) - tst.l m68k_is040or060 - jeq 1f -#endif -#if defined(CPU_M68040_OR_M68060) - move.l (FPS_PC2,%sp),(FPS_PC,%sp) -#endif -1: - | emulate the instruction - jsr fp_scan - -#if defined(CONFIG_M68060) -#if !defined(CPU_M68060_ONLY) - btst #3,m68k_cputype+3 - jeq 1f -#endif - btst #7,(FPS_SR,%sp) - jne fp_sendtrace060 -#endif -1: - | emulation successful? - tst.l %d0 - jeq ret_from_exception - - | send some signal to program here - - jra ret_from_exception - - | we jump here after an access error while trying to access - | user space, we correct stackpointer and send a SIGSEGV to - | the user process -fp_err_ua2: - addq.l #4,%sp -fp_err_ua1: - addq.l #4,%sp - move.l %a0,-(%sp) - pea LSEGV_MAPERR - pea LSIGSEGV - jsr fpemu_signal - add.w #12,%sp - jra ret_from_exception - -#if defined(CONFIG_M68060) - | send a trace signal if we are debugged - | it does not really belong here, but... -fp_sendtrace060: - move.l (FPS_PC,%sp),-(%sp) - pea LTRAP_TRACE - pea LSIGTRAP - jsr fpemu_signal - add.w #12,%sp - jra ret_from_exception -#endif - - .globl fp_get_data_reg, fp_put_data_reg - .globl fp_get_addr_reg, fp_put_addr_reg - - | Entry points to get/put a register. Some of them can be get/put - | directly, others are on the stack, as we read/write the stack - | directly here, these function may only be called from within - | instruction decoding, otherwise the stack pointer is incorrect - | and the stack gets corrupted. -fp_get_data_reg: - jmp ([0f:w,%pc,%d0.w*4]) - - .align 4 -0: - .long fp_get_d0, fp_get_d1 - .long fp_get_d2, fp_get_d3 - .long fp_get_d4, fp_get_d5 - .long fp_get_d6, fp_get_d7 - -fp_get_d0: - move.l (PT_OFF_D0+8,%sp),%d0 - printf PREGISTER,"{d0->%08x}",1,%d0 - rts - -fp_get_d1: - move.l (PT_OFF_D1+8,%sp),%d0 - printf PREGISTER,"{d1->%08x}",1,%d0 - rts - -fp_get_d2: - move.l (PT_OFF_D2+8,%sp),%d0 - printf PREGISTER,"{d2->%08x}",1,%d0 - rts - -fp_get_d3: - move.l %d3,%d0 - printf PREGISTER,"{d3->%08x}",1,%d0 - rts - -fp_get_d4: - move.l %d4,%d0 - printf PREGISTER,"{d4->%08x}",1,%d0 - rts - -fp_get_d5: - move.l %d5,%d0 - printf PREGISTER,"{d5->%08x}",1,%d0 - rts - -fp_get_d6: - move.l %d6,%d0 - printf PREGISTER,"{d6->%08x}",1,%d0 - rts - -fp_get_d7: - move.l %d7,%d0 - printf PREGISTER,"{d7->%08x}",1,%d0 - rts - -fp_put_data_reg: - jmp ([0f:w,%pc,%d1.w*4]) - - .align 4 -0: - .long fp_put_d0, fp_put_d1 - .long fp_put_d2, fp_put_d3 - .long fp_put_d4, fp_put_d5 - .long fp_put_d6, fp_put_d7 - -fp_put_d0: - printf PREGISTER,"{d0<-%08x}",1,%d0 - move.l %d0,(PT_OFF_D0+8,%sp) - rts - -fp_put_d1: - printf PREGISTER,"{d1<-%08x}",1,%d0 - move.l %d0,(PT_OFF_D1+8,%sp) - rts - -fp_put_d2: - printf PREGISTER,"{d2<-%08x}",1,%d0 - move.l %d0,(PT_OFF_D2+8,%sp) - rts - -fp_put_d3: - printf PREGISTER,"{d3<-%08x}",1,%d0 -| move.l %d0,%d3 - move.l %d0,(PT_OFF_D3+8,%sp) - rts - -fp_put_d4: - printf PREGISTER,"{d4<-%08x}",1,%d0 -| move.l %d0,%d4 - move.l %d0,(PT_OFF_D4+8,%sp) - rts - -fp_put_d5: - printf PREGISTER,"{d5<-%08x}",1,%d0 -| move.l %d0,%d5 - move.l %d0,(PT_OFF_D5+8,%sp) - rts - -fp_put_d6: - printf PREGISTER,"{d6<-%08x}",1,%d0 - move.l %d0,%d6 - rts - -fp_put_d7: - printf PREGISTER,"{d7<-%08x}",1,%d0 - move.l %d0,%d7 - rts - -fp_get_addr_reg: - jmp ([0f:w,%pc,%d0.w*4]) - - .align 4 -0: - .long fp_get_a0, fp_get_a1 - .long fp_get_a2, fp_get_a3 - .long fp_get_a4, fp_get_a5 - .long fp_get_a6, fp_get_a7 - -fp_get_a0: - move.l (PT_OFF_A0+8,%sp),%a0 - printf PREGISTER,"{a0->%08x}",1,%a0 - rts - -fp_get_a1: - move.l (PT_OFF_A1+8,%sp),%a0 - printf PREGISTER,"{a1->%08x}",1,%a0 - rts - -fp_get_a2: - move.l (PT_OFF_A2+8,%sp),%a0 - printf PREGISTER,"{a2->%08x}",1,%a0 - rts - -fp_get_a3: - move.l %a3,%a0 - printf PREGISTER,"{a3->%08x}",1,%a0 - rts - -fp_get_a4: - move.l %a4,%a0 - printf PREGISTER,"{a4->%08x}",1,%a0 - rts - -fp_get_a5: - move.l %a5,%a0 - printf PREGISTER,"{a5->%08x}",1,%a0 - rts - -fp_get_a6: - move.l %a6,%a0 - printf PREGISTER,"{a6->%08x}",1,%a0 - rts - -fp_get_a7: - move.l %usp,%a0 - printf PREGISTER,"{a7->%08x}",1,%a0 - rts - -fp_put_addr_reg: - jmp ([0f:w,%pc,%d0.w*4]) - - .align 4 -0: - .long fp_put_a0, fp_put_a1 - .long fp_put_a2, fp_put_a3 - .long fp_put_a4, fp_put_a5 - .long fp_put_a6, fp_put_a7 - -fp_put_a0: - printf PREGISTER,"{a0<-%08x}",1,%a0 - move.l %a0,(PT_OFF_A0+8,%sp) - rts - -fp_put_a1: - printf PREGISTER,"{a1<-%08x}",1,%a0 - move.l %a0,(PT_OFF_A1+8,%sp) - rts - -fp_put_a2: - printf PREGISTER,"{a2<-%08x}",1,%a0 - move.l %a0,(PT_OFF_A2+8,%sp) - rts - -fp_put_a3: - printf PREGISTER,"{a3<-%08x}",1,%a0 - move.l %a0,%a3 - rts - -fp_put_a4: - printf PREGISTER,"{a4<-%08x}",1,%a0 - move.l %a0,%a4 - rts - -fp_put_a5: - printf PREGISTER,"{a5<-%08x}",1,%a0 - move.l %a0,%a5 - rts - -fp_put_a6: - printf PREGISTER,"{a6<-%08x}",1,%a0 - move.l %a0,%a6 - rts - -fp_put_a7: - printf PREGISTER,"{a7<-%08x}",1,%a0 - move.l %a0,%usp - rts - - .data - .align 4 - -fp_debugprint: -| .long PMDECODE - .long PMINSTR+PMDECODE+PMCONV+PMNORM -| .long PMCONV+PMNORM+PMINSTR -| .long 0 diff --git a/ANDROID_3.4.5/arch/m68k/math-emu/fp_log.c b/ANDROID_3.4.5/arch/m68k/math-emu/fp_log.c deleted file mode 100644 index 3384a524..00000000 --- a/ANDROID_3.4.5/arch/m68k/math-emu/fp_log.c +++ /dev/null @@ -1,219 +0,0 @@ -/* - - fp_trig.c: floating-point math routines for the Linux-m68k - floating point emulator. - - Copyright (c) 1998-1999 David Huggins-Daines / Roman Zippel. - - I hereby give permission, free of charge, to copy, modify, and - redistribute this software, in source or binary form, provided that - the above copyright notice and the following disclaimer are included - in all such copies. - - THIS SOFTWARE IS PROVIDED "AS IS", WITH ABSOLUTELY NO WARRANTY, REAL - OR IMPLIED. - -*/ - -#include "fp_emu.h" - -static const struct fp_ext fp_one = -{ - .exp = 0x3fff, -}; - -extern struct fp_ext *fp_fadd(struct fp_ext *dest, const struct fp_ext *src); -extern struct fp_ext *fp_fdiv(struct fp_ext *dest, const struct fp_ext *src); - -struct fp_ext * -fp_fsqrt(struct fp_ext *dest, struct fp_ext *src) -{ - struct fp_ext tmp, src2; - int i, exp; - - dprint(PINSTR, "fsqrt\n"); - - fp_monadic_check(dest, src); - - if (IS_ZERO(dest)) - return dest; - - if (dest->sign) { - fp_set_nan(dest); - return dest; - } - if (IS_INF(dest)) - return dest; - - /* - * sqrt(m) * 2^(p) , if e = 2*p - * sqrt(m*2^e) = - * sqrt(2*m) * 2^(p) , if e = 2*p + 1 - * - * So we use the last bit of the exponent to decide wether to - * use the m or 2*m. - * - * Since only the fractional part of the mantissa is stored and - * the integer part is assumed to be one, we place a 1 or 2 into - * the fixed point representation. - */ - exp = dest->exp; - dest->exp = 0x3FFF; - if (!(exp & 1)) /* lowest bit of exponent is set */ - dest->exp++; - fp_copy_ext(&src2, dest); - - /* - * The taylor row around a for sqrt(x) is: - * sqrt(x) = sqrt(a) + 1/(2*sqrt(a))*(x-a) + R - * With a=1 this gives: - * sqrt(x) = 1 + 1/2*(x-1) - * = 1/2*(1+x) - */ - fp_fadd(dest, &fp_one); - dest->exp--; /* * 1/2 */ - - /* - * We now apply the newton rule to the function - * f(x) := x^2 - r - * which has a null point on x = sqrt(r). - * - * It gives: - * x' := x - f(x)/f'(x) - * = x - (x^2 -r)/(2*x) - * = x - (x - r/x)/2 - * = (2*x - x + r/x)/2 - * = (x + r/x)/2 - */ - for (i = 0; i < 9; i++) { - fp_copy_ext(&tmp, &src2); - - fp_fdiv(&tmp, dest); - fp_fadd(dest, &tmp); - dest->exp--; - } - - dest->exp += (exp - 0x3FFF) / 2; - - return dest; -} - -struct fp_ext * -fp_fetoxm1(struct fp_ext *dest, struct fp_ext *src) -{ - uprint("fetoxm1\n"); - - fp_monadic_check(dest, src); - - return dest; -} - -struct fp_ext * -fp_fetox(struct fp_ext *dest, struct fp_ext *src) -{ - uprint("fetox\n"); - - fp_monadic_check(dest, src); - - return dest; -} - -struct fp_ext * -fp_ftwotox(struct fp_ext *dest, struct fp_ext *src) -{ - uprint("ftwotox\n"); - - fp_monadic_check(dest, src); - - return dest; -} - -struct fp_ext * -fp_ftentox(struct fp_ext *dest, struct fp_ext *src) -{ - uprint("ftentox\n"); - - fp_monadic_check(dest, src); - - return dest; -} - -struct fp_ext * -fp_flogn(struct fp_ext *dest, struct fp_ext *src) -{ - uprint("flogn\n"); - - fp_monadic_check(dest, src); - - return dest; -} - -struct fp_ext * -fp_flognp1(struct fp_ext *dest, struct fp_ext *src) -{ - uprint("flognp1\n"); - - fp_monadic_check(dest, src); - - return dest; -} - -struct fp_ext * -fp_flog10(struct fp_ext *dest, struct fp_ext *src) -{ - uprint("flog10\n"); - - fp_monadic_check(dest, src); - - return dest; -} - -struct fp_ext * -fp_flog2(struct fp_ext *dest, struct fp_ext *src) -{ - uprint("flog2\n"); - - fp_monadic_check(dest, src); - - return dest; -} - -struct fp_ext * -fp_fgetexp(struct fp_ext *dest, struct fp_ext *src) -{ - dprint(PINSTR, "fgetexp\n"); - - fp_monadic_check(dest, src); - - if (IS_INF(dest)) { - fp_set_nan(dest); - return dest; - } - if (IS_ZERO(dest)) - return dest; - - fp_conv_long2ext(dest, (int)dest->exp - 0x3FFF); - - fp_normalize_ext(dest); - - return dest; -} - -struct fp_ext * -fp_fgetman(struct fp_ext *dest, struct fp_ext *src) -{ - dprint(PINSTR, "fgetman\n"); - - fp_monadic_check(dest, src); - - if (IS_ZERO(dest)) - return dest; - - if (IS_INF(dest)) - return dest; - - dest->exp = 0x3FFF; - - return dest; -} - diff --git a/ANDROID_3.4.5/arch/m68k/math-emu/fp_move.S b/ANDROID_3.4.5/arch/m68k/math-emu/fp_move.S deleted file mode 100644 index 71bdf83b..00000000 --- a/ANDROID_3.4.5/arch/m68k/math-emu/fp_move.S +++ /dev/null @@ -1,244 +0,0 @@ -/* - * fp_move.S - * - * Copyright Roman Zippel, 1997. All rights reserved. - * - * Redistribution and use in source and binary forms, with or without - * modification, are permitted provided that the following conditions - * are met: - * 1. Redistributions of source code must retain the above copyright - * notice, and the entire permission notice in its entirety, - * including the disclaimer of warranties. - * 2. Redistributions in binary form must reproduce the above copyright - * notice, this list of conditions and the following disclaimer in the - * documentation and/or other materials provided with the distribution. - * 3. The name of the author may not be used to endorse or promote - * products derived from this software without specific prior - * written permission. - * - * ALTERNATIVELY, this product may be distributed under the terms of - * the GNU General Public License, in which case the provisions of the GPL are - * required INSTEAD OF the above restrictions. (This clause is - * necessary due to a potential bad interaction between the GPL and - * the restrictions contained in a BSD-style copyright.) - * - * THIS SOFTWARE IS PROVIDED ``AS IS'' AND ANY EXPRESS OR IMPLIED - * WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES - * OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE - * DISCLAIMED. IN NO EVENT SHALL THE AUTHOR BE LIABLE FOR ANY DIRECT, - * INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES - * (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR - * SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) - * HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, - * STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) - * ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED - * OF THE POSSIBILITY OF SUCH DAMAGE. - */ - -#include "fp_emu.h" -#include "fp_decode.h" - -do_no_pc_mode=1 - - .globl fp_fmove_fp2mem - -fp_fmove_fp2mem: - clr.b (2+FPD_FPSR,FPDATA) - fp_decode_dest_format - move.w %d0,%d1 | store data size twice in %d1 - swap %d1 | one can be trashed below - move.w %d0,%d1 -#ifdef FPU_EMU_DEBUG - lea 0f,%a0 - clr.l %d0 - move.b (%a0,%d1.w),%d0 - printf PDECODE,"fmove.%c ",1,%d0 - fp_decode_src_reg - printf PDECODE,"fp%d,",1,%d0 - - .data -0: .byte 'l','s','x','p','w','d','b','p' - .previous -#endif - - | encode addressing mode for dest - fp_decode_addr_mode - - .long fp_data, fp_ill - .long fp_indirect, fp_postinc - .long fp_predecr, fp_disp16 - .long fp_extmode0, fp_extmode1 - - | addressing mode: data register direct -fp_data: - fp_mode_data_direct - move.w %d0,%d1 - fp_decode_src_reg - fp_get_fp_reg - lea (FPD_TEMPFP1,FPDATA),%a1 - move.l (%a0)+,(%a1)+ - move.l (%a0)+,(%a1)+ - move.l (%a0),(%a1) - lea (-8,%a1),%a0 - swap %d1 - move.l %d1,%d2 - printf PDECODE,"\n" - jmp ([0f:w,%pc,%d1.w*4]) - - .align 4 -0: - .long fp_data_long, fp_data_single - .long fp_ill, fp_ill - .long fp_data_word, fp_ill - .long fp_data_byte, fp_ill - -fp_data_byte: - jsr fp_normalize_ext - jsr fp_conv_ext2byte - move.l %d0,%d1 - swap %d2 - move.w %d2,%d0 - jsr fp_get_data_reg - move.b %d1,%d0 - move.w %d2,%d1 - jsr fp_put_data_reg - jra fp_final - -fp_data_word: - jsr fp_normalize_ext - jsr fp_conv_ext2short - move.l %d0,%d1 - swap %d2 - move.w %d2,%d0 - jsr fp_get_data_reg - move.w %d1,%d0 - move.l %d2,%d1 - jsr fp_put_data_reg - jra fp_final - -fp_data_long: - jsr fp_normalize_ext - jsr fp_conv_ext2long - swap %d2 - move.w %d2,%d1 - jsr fp_put_data_reg - jra fp_final - -fp_data_single: - jsr fp_normalize_ext - jsr fp_conv_ext2single - swap %d2 - move.w %d2,%d1 - jsr fp_put_data_reg - jra fp_final - - | addressing mode: address register indirect -fp_indirect: - fp_mode_addr_indirect - jra fp_putdest - - | addressing mode: address register indirect with postincrement -fp_postinc: - fp_mode_addr_indirect_postinc - jra fp_putdest - - | addressing mode: address register indirect with predecrement -fp_predecr: - fp_mode_addr_indirect_predec - jra fp_putdest - - | addressing mode: address register indirect with 16bit displacement -fp_disp16: - fp_mode_addr_indirect_disp16 - jra fp_putdest - -fp_extmode0: - fp_mode_addr_indirect_extmode0 - jra fp_putdest - -fp_extmode1: - fp_decode_addr_reg - jmp ([0f:w,%pc,%d0*4]) - - .align 4 -0: - .long fp_abs_short, fp_abs_long - .long fp_ill, fp_ill - .long fp_ill, fp_ill - .long fp_ill, fp_ill - -fp_abs_short: - fp_mode_abs_short - jra fp_putdest - -fp_abs_long: - fp_mode_abs_long - jra fp_putdest - -fp_putdest: - move.l %a0,%a1 - fp_decode_src_reg - move.l %d1,%d2 | save size - fp_get_fp_reg - printf PDECODE,"\n" - addq.l #8,%a0 - move.l (%a0),-(%sp) - move.l -(%a0),-(%sp) - move.l -(%a0),-(%sp) - move.l %sp,%a0 - jsr fp_normalize_ext - - swap %d2 - jmp ([0f:w,%pc,%d2.w*4]) - - .align 4 -0: - .long fp_format_long, fp_format_single - .long fp_format_extended, fp_format_packed - .long fp_format_word, fp_format_double - .long fp_format_byte, fp_format_packed - -fp_format_long: - jsr fp_conv_ext2long - putuser.l %d0,(%a1),fp_err_ua1,%a1 - jra fp_finish_move - -fp_format_single: - jsr fp_conv_ext2single - putuser.l %d0,(%a1),fp_err_ua1,%a1 - jra fp_finish_move - -fp_format_extended: - move.l (%a0)+,%d0 - lsl.w #1,%d0 - lsl.l #7,%d0 - lsl.l #8,%d0 - putuser.l %d0,(%a1)+,fp_err_ua1,%a1 - move.l (%a0)+,%d0 - putuser.l %d0,(%a1)+,fp_err_ua1,%a1 - move.l (%a0),%d0 - putuser.l %d0,(%a1),fp_err_ua1,%a1 - jra fp_finish_move - -fp_format_packed: - /* not supported yet */ - lea (12,%sp),%sp - jra fp_ill - -fp_format_word: - jsr fp_conv_ext2short - putuser.w %d0,(%a1),fp_err_ua1,%a1 - jra fp_finish_move - -fp_format_double: - jsr fp_conv_ext2double - jra fp_finish_move - -fp_format_byte: - jsr fp_conv_ext2byte - putuser.b %d0,(%a1),fp_err_ua1,%a1 -| jra fp_finish_move - -fp_finish_move: - lea (12,%sp),%sp - jra fp_final diff --git a/ANDROID_3.4.5/arch/m68k/math-emu/fp_movem.S b/ANDROID_3.4.5/arch/m68k/math-emu/fp_movem.S deleted file mode 100644 index 8354d39e..00000000 --- a/ANDROID_3.4.5/arch/m68k/math-emu/fp_movem.S +++ /dev/null @@ -1,368 +0,0 @@ -/* - * fp_movem.S - * - * Copyright Roman Zippel, 1997. All rights reserved. - * - * Redistribution and use in source and binary forms, with or without - * modification, are permitted provided that the following conditions - * are met: - * 1. Redistributions of source code must retain the above copyright - * notice, and the entire permission notice in its entirety, - * including the disclaimer of warranties. - * 2. Redistributions in binary form must reproduce the above copyright - * notice, this list of conditions and the following disclaimer in the - * documentation and/or other materials provided with the distribution. - * 3. The name of the author may not be used to endorse or promote - * products derived from this software without specific prior - * written permission. - * - * ALTERNATIVELY, this product may be distributed under the terms of - * the GNU General Public License, in which case the provisions of the GPL are - * required INSTEAD OF the above restrictions. (This clause is - * necessary due to a potential bad interaction between the GPL and - * the restrictions contained in a BSD-style copyright.) - * - * THIS SOFTWARE IS PROVIDED ``AS IS'' AND ANY EXPRESS OR IMPLIED - * WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES - * OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE - * DISCLAIMED. IN NO EVENT SHALL THE AUTHOR BE LIABLE FOR ANY DIRECT, - * INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES - * (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR - * SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) - * HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, - * STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) - * ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED - * OF THE POSSIBILITY OF SUCH DAMAGE. - */ - -#include "fp_emu.h" -#include "fp_decode.h" - -| set flags for decode macros for fmovem -do_fmovem=1 - - .globl fp_fmovem_fp, fp_fmovem_cr - -| %d1 contains the mask and count of the register list -| for other register usage see fp_decode.h - -fp_fmovem_fp: - printf PDECODE,"fmovem.x " - | get register list and count them - btst #11,%d2 - jne 1f - bfextu %d2{#24,#8},%d0 | static register list - jra 2f -1: bfextu %d2{#25,#3},%d0 | dynamic register list - jsr fp_get_data_reg -2: move.l %d0,%d1 - swap %d1 - jra 2f -1: addq.w #1,%d1 | count the # of registers in -2: lsr.b #1,%d0 | register list and keep it in %d1 - jcs 1b - jne 2b - printf PDECODE,"#%08x",1,%d1 -#ifdef FPU_EMU_DEBUG - btst #12,%d2 - jne 1f - printf PDECODE,"-" | decremental move - jra 2f -1: printf PDECODE,"+" | incremental move -2: btst #13,%d2 - jeq 1f - printf PDECODE,"->" | fpu -> cpu - jra 2f -1: printf PDECODE,"<-" | fpu <- cpu -2: -#endif - - | decode address mode - fp_decode_addr_mode - - .long fp_ill, fp_ill - .long fpr_indirect, fpr_postinc - .long fpr_predecr, fpr_disp16 - .long fpr_extmode0, fpr_extmode1 - - | addressing mode: address register indirect -fpr_indirect: - fp_mode_addr_indirect - jra fpr_do_movem - - | addressing mode: address register indirect with postincrement -fpr_postinc: - fp_mode_addr_indirect_postinc - jra fpr_do_movem - -fpr_predecr: - fp_mode_addr_indirect_predec - jra fpr_do_movem - - | addressing mode: address register/programm counter indirect - | with 16bit displacement -fpr_disp16: - fp_mode_addr_indirect_disp16 - jra fpr_do_movem - -fpr_extmode0: - fp_mode_addr_indirect_extmode0 - jra fpr_do_movem - -fpr_extmode1: - fp_decode_addr_reg - jmp ([0f:w,%pc,%d0*4]) - - .align 4 -0: - .long fpr_absolute_short, fpr_absolute_long - .long fpr_disp16, fpr_extmode0 - .long fp_ill, fp_ill - .long fp_ill, fp_ill - -fpr_absolute_short: - fp_mode_abs_short - jra fpr_do_movem - -fpr_absolute_long: - fp_mode_abs_long -| jra fpr_do_movem - -fpr_do_movem: - swap %d1 | get fpu register list - lea (FPD_FPREG,FPDATA),%a1 - moveq #12,%d0 - btst #12,%d2 - jne 1f - lea (-12,%a1,%d0*8),%a1 - neg.l %d0 -1: btst #13,%d2 - jne 4f - | move register from memory into fpu - jra 3f -1: printf PMOVEM,"(%p>%p)",2,%a0,%a1 - getuser.l (%a0)+,%d2,fp_err_ua1,%a0 - lsr.l #8,%d2 - lsr.l #7,%d2 - lsr.w #1,%d2 - move.l %d2,(%a1)+ - getuser.l (%a0)+,%d2,fp_err_ua1,%a0 - move.l %d2,(%a1)+ - getuser.l (%a0),%d2,fp_err_ua1,%a0 - move.l %d2,(%a1) - subq.l #8,%a0 - subq.l #8,%a1 - add.l %d0,%a0 -2: add.l %d0,%a1 -3: lsl.b #1,%d1 - jcs 1b - jne 2b - jra 5f - | move register from fpu into memory -1: printf PMOVEM,"(%p>%p)",2,%a1,%a0 - move.l (%a1)+,%d2 - lsl.w #1,%d2 - lsl.l #7,%d2 - lsl.l #8,%d2 - putuser.l %d2,(%a0)+,fp_err_ua1,%a0 - move.l (%a1)+,%d2 - putuser.l %d2,(%a0)+,fp_err_ua1,%a0 - move.l (%a1),%d2 - putuser.l %d2,(%a0),fp_err_ua1,%a0 - subq.l #8,%a1 - subq.l #8,%a0 - add.l %d0,%a0 -2: add.l %d0,%a1 -4: lsl.b #1,%d1 - jcs 1b - jne 2b -5: - printf PDECODE,"\n" -#if 0 - lea (FPD_FPREG,FPDATA),%a0 - printf PMOVEM,"fp:" - printx PMOVEM,%a0@(0) - printx PMOVEM,%a0@(12) - printf PMOVEM,"\n " - printx PMOVEM,%a0@(24) - printx PMOVEM,%a0@(36) - printf PMOVEM,"\n " - printx PMOVEM,%a0@(48) - printx PMOVEM,%a0@(60) - printf PMOVEM,"\n " - printx PMOVEM,%a0@(72) - printx PMOVEM,%a0@(84) - printf PMOVEM,"\n" -#endif - jra fp_end - -| set flags for decode macros for fmovem control register -do_fmovem=1 -do_fmovem_cr=1 - -fp_fmovem_cr: - printf PDECODE,"fmovem.cr " - | get register list and count them - bfextu %d2{#19,#3},%d0 - move.l %d0,%d1 - swap %d1 - jra 2f -1: addq.w #1,%d1 -2: lsr.l #1,%d0 - jcs 1b - jne 2b - printf PDECODE,"#%08x",1,%d1 -#ifdef FPU_EMU_DEBUG - btst #13,%d2 - jeq 1f - printf PDECODE,"->" | fpu -> cpu - jra 2f -1: printf PDECODE,"<-" | fpu <- cpu -2: -#endif - - | decode address mode - fp_decode_addr_mode - - .long fpc_data, fpc_addr - .long fpc_indirect, fpc_postinc - .long fpc_predecr, fpc_disp16 - .long fpc_extmode0, fpc_extmode1 - -fpc_data: - fp_mode_data_direct - move.w %d0,%d1 - bfffo %d2{#19,#3},%d0 - sub.w #19,%d0 - lea (FPD_FPCR,FPDATA,%d0.w*4),%a1 - btst #13,%d2 - jne 1f - move.w %d1,%d0 - jsr fp_get_data_reg - move.l %d0,(%a1) - jra fpc_movem_fin -1: move.l (%a1),%d0 - jsr fp_put_data_reg - jra fpc_movem_fin - -fpc_addr: - fp_decode_addr_reg - printf PDECODE,"a%d",1,%d0 - btst #13,%d2 - jne 1f - jsr fp_get_addr_reg - move.l %a0,(FPD_FPIAR,FPDATA) - jra fpc_movem_fin -1: move.l (FPD_FPIAR,FPDATA),%a0 - jsr fp_put_addr_reg - jra fpc_movem_fin - -fpc_indirect: - fp_mode_addr_indirect - jra fpc_do_movem - -fpc_postinc: - fp_mode_addr_indirect_postinc - jra fpc_do_movem - -fpc_predecr: - fp_mode_addr_indirect_predec - jra fpc_do_movem - -fpc_disp16: - fp_mode_addr_indirect_disp16 - jra fpc_do_movem - -fpc_extmode0: - fp_mode_addr_indirect_extmode0 - jra fpc_do_movem - -fpc_extmode1: - fp_decode_addr_reg - jmp ([0f:w,%pc,%d0*4]) - - .align 4 -0: - .long fpc_absolute_short, fpc_absolute_long - .long fpc_disp16, fpc_extmode0 - .long fpc_immediate, fp_ill - .long fp_ill, fp_ill - -fpc_absolute_short: - fp_mode_abs_short - jra fpc_do_movem - -fpc_absolute_long: - fp_mode_abs_long - jra fpc_do_movem - -fpc_immediate: - fp_get_pc %a0 - lea (%a0,%d1.w*4),%a1 - fp_put_pc %a1 - printf PDECODE,"#imm" -| jra fpc_do_movem -#if 0 - swap %d1 - lsl.l #5,%d1 - lea (FPD_FPCR,FPDATA),%a0 - jra 3f -1: move.l %d0,(%a0) -2: addq.l #4,%a0 -3: lsl.b #1,%d1 - jcs 1b - jne 2b - jra fpc_movem_fin -#endif - -fpc_do_movem: - swap %d1 | get fpu register list - lsl.l #5,%d1 - lea (FPD_FPCR,FPDATA),%a1 -1: btst #13,%d2 - jne 4f - - | move register from memory into fpu - jra 3f -1: printf PMOVEM,"(%p>%p)",2,%a0,%a1 - getuser.l (%a0)+,%d0,fp_err_ua1,%a0 - move.l %d0,(%a1) -2: addq.l #4,%a1 -3: lsl.b #1,%d1 - jcs 1b - jne 2b - jra fpc_movem_fin - - | move register from fpu into memory -1: printf PMOVEM,"(%p>%p)",2,%a1,%a0 - move.l (%a1),%d0 - putuser.l %d0,(%a0)+,fp_err_ua1,%a0 -2: addq.l #4,%a1 -4: lsl.b #1,%d1 - jcs 1b - jne 2b - -fpc_movem_fin: - and.l #0x0000fff0,(FPD_FPCR,FPDATA) - and.l #0x0ffffff8,(FPD_FPSR,FPDATA) - move.l (FPD_FPCR,FPDATA),%d0 - lsr.l #4,%d0 - moveq #3,%d1 - and.l %d0,%d1 - move.w %d1,(FPD_RND,FPDATA) - lsr.l #2,%d0 - moveq #3,%d1 - and.l %d0,%d1 - move.w %d1,(FPD_PREC,FPDATA) - printf PDECODE,"\n" -#if 0 - printf PMOVEM,"fpcr : %08x\n",1,FPDATA@(FPD_FPCR) - printf PMOVEM,"fpsr : %08x\n",1,FPDATA@(FPD_FPSR) - printf PMOVEM,"fpiar: %08x\n",1,FPDATA@(FPD_FPIAR) - clr.l %d0 - move.w (FPD_PREC,FPDATA),%d0 - printf PMOVEM,"prec : %04x\n",1,%d0 - move.w (FPD_RND,FPDATA),%d0 - printf PMOVEM,"rnd : %04x\n",1,%d0 -#endif - jra fp_end diff --git a/ANDROID_3.4.5/arch/m68k/math-emu/fp_scan.S b/ANDROID_3.4.5/arch/m68k/math-emu/fp_scan.S deleted file mode 100644 index e4146ed5..00000000 --- a/ANDROID_3.4.5/arch/m68k/math-emu/fp_scan.S +++ /dev/null @@ -1,478 +0,0 @@ -/* - * fp_scan.S - * - * Copyright Roman Zippel, 1997. All rights reserved. - * - * Redistribution and use in source and binary forms, with or without - * modification, are permitted provided that the following conditions - * are met: - * 1. Redistributions of source code must retain the above copyright - * notice, and the entire permission notice in its entirety, - * including the disclaimer of warranties. - * 2. Redistributions in binary form must reproduce the above copyright - * notice, this list of conditions and the following disclaimer in the - * documentation and/or other materials provided with the distribution. - * 3. The name of the author may not be used to endorse or promote - * products derived from this software without specific prior - * written permission. - * - * ALTERNATIVELY, this product may be distributed under the terms of - * the GNU General Public License, in which case the provisions of the GPL are - * required INSTEAD OF the above restrictions. (This clause is - * necessary due to a potential bad interaction between the GPL and - * the restrictions contained in a BSD-style copyright.) - * - * THIS SOFTWARE IS PROVIDED ``AS IS'' AND ANY EXPRESS OR IMPLIED - * WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES - * OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE - * DISCLAIMED. IN NO EVENT SHALL THE AUTHOR BE LIABLE FOR ANY DIRECT, - * INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES - * (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR - * SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) - * HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, - * STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) - * ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED - * OF THE POSSIBILITY OF SUCH DAMAGE. - */ - -#include "fp_emu.h" -#include "fp_decode.h" - - .globl fp_scan, fp_datasize - - .data - -| %d2 - first two instr words -| %d1 - operand size - -/* operand formats are: - - Long = 0, i.e. fmove.l - Single, i.e. fmove.s - Extended, i.e. fmove.x - Packed-BCD, i.e. fmove.p - Word, i.e. fmove.w - Double, i.e. fmove.d -*/ - - .text - -| On entry: -| FPDATA - base of emulated FPU registers - -fp_scan: -| normal fpu instruction? (this excludes fsave/frestore) - fp_get_pc %a0 - printf PDECODE,"%08x: ",1,%a0 - getuser.b (%a0),%d0,fp_err_ua1,%a0 -#if 1 - cmp.b #0xf2,%d0 | cpid = 1 -#else - cmp.b #0xfc,%d0 | cpid = 6 -#endif - jne fp_nonstd -| first two instruction words are kept in %d2 - getuser.l (%a0)+,%d2,fp_err_ua1,%a0 - fp_put_pc %a0 -fp_decode_cond: | separate conditional instr - fp_decode_cond_instr_type - - .long fp_decode_move, fp_fscc - .long fp_fbccw, fp_fbccl - -fp_decode_move: | separate move instr - fp_decode_move_instr_type - - .long fp_fgen_fp, fp_ill - .long fp_fgen_ea, fp_fmove_fp2mem - .long fp_fmovem_cr, fp_fmovem_cr - .long fp_fmovem_fp, fp_fmovem_fp - -| now all arithmetic instr and a few move instr are left -fp_fgen_fp: | source is a fpu register - clr.b (FPD_FPSR+2,FPDATA) | clear the exception byte - fp_decode_sourcespec - printf PDECODE,"f<op>.x fp%d",1,%d0 - fp_get_fp_reg - lea (FPD_TEMPFP1,FPDATA),%a1 | copy src into a temp location - move.l (%a0)+,(%a1)+ - move.l (%a0)+,(%a1)+ - move.l (%a0),(%a1) - lea (-8,%a1),%a0 - jra fp_getdest - -fp_fgen_ea: | source is <ea> - clr.b (FPD_FPSR+2,FPDATA) | clear the exception byte - | sort out fmovecr, keep data size in %d1 - fp_decode_sourcespec - cmp.w #7,%d0 - jeq fp_fmovecr - move.w %d0,%d1 | store data size twice in %d1 - swap %d1 | one can be trashed below - move.w %d0,%d1 -#ifdef FPU_EMU_DEBUG - lea 0f,%a0 - clr.l %d0 - move.b (%a0,%d1.w),%d0 - printf PDECODE,"f<op>.%c ",1,%d0 - - .data -0: .byte 'l','s','x','p','w','d','b',0 - .previous -#endif - -/* - fp_getsource, fp_getdest - - basically, we end up with a pointer to the source operand in - %a1, and a pointer to the destination operand in %a0. both - are, of course, 96-bit extended floating point numbers. -*/ - -fp_getsource: - | decode addressing mode for source - fp_decode_addr_mode - - .long fp_data, fp_ill - .long fp_indirect, fp_postinc - .long fp_predecr, fp_disp16 - .long fp_extmode0, fp_extmode1 - - | addressing mode: data register direct -fp_data: - fp_mode_data_direct - jsr fp_get_data_reg - lea (FPD_TEMPFP1,FPDATA),%a0 - jmp ([0f:w,%pc,%d1.w*4]) - - .align 4 -0: - .long fp_data_long, fp_data_single - .long fp_ill, fp_ill - .long fp_data_word, fp_ill - .long fp_data_byte, fp_ill - - | data types that fit in an integer data register -fp_data_byte: - extb.l %d0 - jra fp_data_long - -fp_data_word: - ext.l %d0 - -fp_data_long: - jsr fp_conv_long2ext - jra fp_getdest - -fp_data_single: - jsr fp_conv_single2ext - jra fp_getdest - - | addressing mode: address register indirect -fp_indirect: - fp_mode_addr_indirect - jra fp_fetchsource - - | addressing mode: address register indirect with postincrement -fp_postinc: - fp_mode_addr_indirect_postinc - jra fp_fetchsource - - | addressing mode: address register indirect with predecrement -fp_predecr: - fp_mode_addr_indirect_predec - jra fp_fetchsource - - | addressing mode: address register/programm counter indirect - | with 16bit displacement -fp_disp16: - fp_mode_addr_indirect_disp16 - jra fp_fetchsource - - | all other indirect addressing modes will finally end up here -fp_extmode0: - fp_mode_addr_indirect_extmode0 - jra fp_fetchsource - -| all pc relative addressing modes and immediate/absolute modes end up here -| the first ones are sent to fp_extmode0 or fp_disp16 -| and only the latter are handled here -fp_extmode1: - fp_decode_addr_reg - jmp ([0f:w,%pc,%d0*4]) - - .align 4 -0: - .long fp_abs_short, fp_abs_long - .long fp_disp16, fp_extmode0 - .long fp_immediate, fp_ill - .long fp_ill, fp_ill - - | addressing mode: absolute short -fp_abs_short: - fp_mode_abs_short - jra fp_fetchsource - - | addressing mode: absolute long -fp_abs_long: - fp_mode_abs_long - jra fp_fetchsource - - | addressing mode: immediate data -fp_immediate: - printf PDECODE,"#" - fp_get_pc %a0 - move.w (fp_datasize,%d1.w*2),%d0 - addq.w #1,%d0 - and.w #-2,%d0 -#ifdef FPU_EMU_DEBUG - movem.l %d0/%d1,-(%sp) - movel %a0,%a1 - clr.l %d1 - jra 2f -1: getuser.b (%a1)+,%d1,fp_err_ua1,%a1 - printf PDECODE,"%02x",1,%d1 -2: dbra %d0,1b - movem.l (%sp)+,%d0/%d1 -#endif - lea (%a0,%d0.w),%a1 - fp_put_pc %a1 -| jra fp_fetchsource - -fp_fetchsource: - move.l %a0,%a1 - swap %d1 - lea (FPD_TEMPFP1,FPDATA),%a0 - jmp ([0f:w,%pc,%d1.w*4]) - - .align 4 -0: .long fp_long, fp_single - .long fp_ext, fp_pack - .long fp_word, fp_double - .long fp_byte, fp_ill - -fp_long: - getuser.l (%a1),%d0,fp_err_ua1,%a1 - jsr fp_conv_long2ext - jra fp_getdest - -fp_single: - getuser.l (%a1),%d0,fp_err_ua1,%a1 - jsr fp_conv_single2ext - jra fp_getdest - -fp_ext: - getuser.l (%a1)+,%d0,fp_err_ua1,%a1 - lsr.l #8,%d0 - lsr.l #7,%d0 - lsr.w #1,%d0 - move.l %d0,(%a0)+ - getuser.l (%a1)+,%d0,fp_err_ua1,%a1 - move.l %d0,(%a0)+ - getuser.l (%a1),%d0,fp_err_ua1,%a1 - move.l %d0,(%a0) - subq.l #8,%a0 - jra fp_getdest - -fp_pack: - /* not supported yet */ - jra fp_ill - -fp_word: - getuser.w (%a1),%d0,fp_err_ua1,%a1 - ext.l %d0 - jsr fp_conv_long2ext - jra fp_getdest - -fp_double: - jsr fp_conv_double2ext - jra fp_getdest - -fp_byte: - getuser.b (%a1),%d0,fp_err_ua1,%a1 - extb.l %d0 - jsr fp_conv_long2ext -| jra fp_getdest - -fp_getdest: - move.l %a0,%a1 - bfextu %d2{#22,#3},%d0 - printf PDECODE,",fp%d\n",1,%d0 - fp_get_fp_reg - movem.l %a0/%a1,-(%sp) - pea fp_finalrounding - bfextu %d2{#25,#7},%d0 - jmp ([0f:w,%pc,%d0*4]) - - .align 4 -0: - .long fp_fmove_mem2fp, fp_fint, fp_fsinh, fp_fintrz - .long fp_fsqrt, fp_ill, fp_flognp1, fp_ill - .long fp_fetoxm1, fp_ftanh, fp_fatan, fp_ill - .long fp_fasin, fp_fatanh, fp_fsin, fp_ftan - .long fp_fetox, fp_ftwotox, fp_ftentox, fp_ill - .long fp_flogn, fp_flog10, fp_flog2, fp_ill - .long fp_fabs, fp_fcosh, fp_fneg, fp_ill - .long fp_facos, fp_fcos, fp_fgetexp, fp_fgetman - .long fp_fdiv, fp_fmod, fp_fadd, fp_fmul - .long fpa_fsgldiv, fp_frem, fp_fscale, fpa_fsglmul - .long fp_fsub, fp_ill, fp_ill, fp_ill - .long fp_ill, fp_ill, fp_ill, fp_ill - .long fp_fsincos0, fp_fsincos1, fp_fsincos2, fp_fsincos3 - .long fp_fsincos4, fp_fsincos5, fp_fsincos6, fp_fsincos7 - .long fp_fcmp, fp_ill, fp_ftst, fp_ill - .long fp_ill, fp_ill, fp_ill, fp_ill - .long fp_fsmove, fp_fssqrt, fp_ill, fp_ill - .long fp_fdmove, fp_fdsqrt, fp_ill, fp_ill - .long fp_ill, fp_ill, fp_ill, fp_ill - .long fp_ill, fp_ill, fp_ill, fp_ill - .long fp_ill, fp_ill, fp_ill, fp_ill - .long fp_ill, fp_ill, fp_ill, fp_ill - .long fp_fsabs, fp_ill, fp_fsneg, fp_ill - .long fp_fdabs, fp_ill, fp_fdneg, fp_ill - .long fp_fsdiv, fp_ill, fp_fsadd, fp_fsmul - .long fp_fddiv, fp_ill, fp_fdadd, fp_fdmul - .long fp_fssub, fp_ill, fp_ill, fp_ill - .long fp_fdsub, fp_ill, fp_ill, fp_ill - .long fp_ill, fp_ill, fp_ill, fp_ill - .long fp_ill, fp_ill, fp_ill, fp_ill - .long fp_ill, fp_ill, fp_ill, fp_ill - .long fp_ill, fp_ill, fp_ill, fp_ill - - | Instructions follow - - | Move an (emulated) ROM constant -fp_fmovecr: - bfextu %d2{#27,#5},%d0 - printf PINSTR,"fp_fmovecr #%d",1,%d0 - move.l %d0,%d1 - add.l %d0,%d0 - add.l %d1,%d0 - lea (fp_constants,%d0*4),%a0 - move.l #0x801cc0ff,%d0 - addq.l #1,%d1 - lsl.l %d1,%d0 - jcc 1f - fp_set_sr FPSR_EXC_INEX2 | INEX2 exception -1: moveq #-128,%d0 | continue with fmove - and.l %d0,%d2 - jra fp_getdest - - .data - .align 4 -fp_constants: - .long 0x00004000,0xc90fdaa2,0x2168c235 | pi - .extend 0,0,0,0,0,0,0,0,0,0 - .long 0x00003ffd,0x9a209a84,0xfbcff798 | log10(2) - .long 0x00004000,0xadf85458,0xa2bb4a9a | e - .long 0x00003fff,0xb8aa3b29,0x5c17f0bc | log2(e) - .long 0x00003ffd,0xde5bd8a9,0x37287195 | log10(e) - .long 0x00000000,0x00000000,0x00000000 | 0.0 - .long 0x00003ffe,0xb17217f7,0xd1cf79ac | 1n(2) - .long 0x00004000,0x935d8ddd,0xaaa8ac17 | 1n(10) - | read this as "1.0 * 2^0" - note the high bit in the mantissa - .long 0x00003fff,0x80000000,0x00000000 | 10^0 - .long 0x00004002,0xa0000000,0x00000000 | 10^1 - .long 0x00004005,0xc8000000,0x00000000 | 10^2 - .long 0x0000400c,0x9c400000,0x00000000 | 10^4 - .long 0x00004019,0xbebc2000,0x00000000 | 10^8 - .long 0x00004034,0x8e1bc9bf,0x04000000 | 10^16 - .long 0x00004069,0x9dc5ada8,0x2b70b59e | 10^32 - .long 0x000040d3,0xc2781f49,0xffcfa6d5 | 10^64 - .long 0x000041a8,0x93ba47c9,0x80e98ce0 | 10^128 - .long 0x00004351,0xaa7eebfb,0x9df9de8e | 10^256 - .long 0x000046a3,0xe319a0ae,0xa60e91c7 | 10^512 - .long 0x00004d48,0xc9767586,0x81750c17 | 10^1024 - .long 0x00005a92,0x9e8b3b5d,0xc53d5de5 | 10^2048 - .long 0x00007525,0xc4605202,0x8a20979b | 10^4096 - .previous - -fp_fmove_mem2fp: - printf PINSTR,"fmove %p,%p\n",2,%a0,%a1 - move.l (%a1)+,(%a0)+ - move.l (%a1)+,(%a0)+ - move.l (%a1),(%a0) - subq.l #8,%a0 - rts - -fpa_fsglmul: - move.l #fp_finalrounding_single_fast,(%sp) - jra fp_fsglmul - -fpa_fsgldiv: - move.l #fp_finalrounding_single_fast,(%sp) - jra fp_fsgldiv - -.macro fp_dosingleprec instr - printf PINSTR,"single " - move.l #fp_finalrounding_single,(%sp) - jra \instr -.endm - -.macro fp_dodoubleprec instr - printf PINSTR,"double " - move.l #fp_finalrounding_double,(%sp) - jra \instr -.endm - -fp_fsmove: - fp_dosingleprec fp_fmove_mem2fp - -fp_fssqrt: - fp_dosingleprec fp_fsqrt - -fp_fdmove: - fp_dodoubleprec fp_fmove_mem2fp - -fp_fdsqrt: - fp_dodoubleprec fp_fsqrt - -fp_fsabs: - fp_dosingleprec fp_fabs - -fp_fsneg: - fp_dosingleprec fp_fneg - -fp_fdabs: - fp_dodoubleprec fp_fabs - -fp_fdneg: - fp_dodoubleprec fp_fneg - -fp_fsdiv: - fp_dosingleprec fp_fdiv - -fp_fsadd: - fp_dosingleprec fp_fadd - -fp_fsmul: - fp_dosingleprec fp_fmul - -fp_fddiv: - fp_dodoubleprec fp_fdiv - -fp_fdadd: - fp_dodoubleprec fp_fadd - -fp_fdmul: - fp_dodoubleprec fp_fmul - -fp_fssub: - fp_dosingleprec fp_fsub - -fp_fdsub: - fp_dodoubleprec fp_fsub - -fp_nonstd: - fp_get_pc %a0 - getuser.l (%a0),%d0,fp_err_ua1,%a0 - printf ,"nonstd ((%08x)=%08x)\n",2,%a0,%d0 - moveq #-1,%d0 - rts - - .data - .align 4 - - | data sizes corresponding to the operand formats -fp_datasize: - .word 4, 4, 12, 12, 2, 8, 1, 0 diff --git a/ANDROID_3.4.5/arch/m68k/math-emu/fp_trig.c b/ANDROID_3.4.5/arch/m68k/math-emu/fp_trig.c deleted file mode 100644 index 6361d078..00000000 --- a/ANDROID_3.4.5/arch/m68k/math-emu/fp_trig.c +++ /dev/null @@ -1,183 +0,0 @@ -/* - - fp_trig.c: floating-point math routines for the Linux-m68k - floating point emulator. - - Copyright (c) 1998-1999 David Huggins-Daines / Roman Zippel. - - I hereby give permission, free of charge, to copy, modify, and - redistribute this software, in source or binary form, provided that - the above copyright notice and the following disclaimer are included - in all such copies. - - THIS SOFTWARE IS PROVIDED "AS IS", WITH ABSOLUTELY NO WARRANTY, REAL - OR IMPLIED. - -*/ - -#include "fp_emu.h" -#include "fp_trig.h" - -struct fp_ext * -fp_fsin(struct fp_ext *dest, struct fp_ext *src) -{ - uprint("fsin\n"); - - fp_monadic_check(dest, src); - - return dest; -} - -struct fp_ext * -fp_fcos(struct fp_ext *dest, struct fp_ext *src) -{ - uprint("fcos\n"); - - fp_monadic_check(dest, src); - - return dest; -} - -struct fp_ext * -fp_ftan(struct fp_ext *dest, struct fp_ext *src) -{ - uprint("ftan\n"); - - fp_monadic_check(dest, src); - - return dest; -} - -struct fp_ext * -fp_fasin(struct fp_ext *dest, struct fp_ext *src) -{ - uprint("fasin\n"); - - fp_monadic_check(dest, src); - - return dest; -} - -struct fp_ext * -fp_facos(struct fp_ext *dest, struct fp_ext *src) -{ - uprint("facos\n"); - - fp_monadic_check(dest, src); - - return dest; -} - -struct fp_ext * -fp_fatan(struct fp_ext *dest, struct fp_ext *src) -{ - uprint("fatan\n"); - - fp_monadic_check(dest, src); - - return dest; -} - -struct fp_ext * -fp_fsinh(struct fp_ext *dest, struct fp_ext *src) -{ - uprint("fsinh\n"); - - fp_monadic_check(dest, src); - - return dest; -} - -struct fp_ext * -fp_fcosh(struct fp_ext *dest, struct fp_ext *src) -{ - uprint("fcosh\n"); - - fp_monadic_check(dest, src); - - return dest; -} - -struct fp_ext * -fp_ftanh(struct fp_ext *dest, struct fp_ext *src) -{ - uprint("ftanh\n"); - - fp_monadic_check(dest, src); - - return dest; -} - -struct fp_ext * -fp_fatanh(struct fp_ext *dest, struct fp_ext *src) -{ - uprint("fatanh\n"); - - fp_monadic_check(dest, src); - - return dest; -} - -struct fp_ext * -fp_fsincos0(struct fp_ext *dest, struct fp_ext *src) -{ - uprint("fsincos0\n"); - - return dest; -} - -struct fp_ext * -fp_fsincos1(struct fp_ext *dest, struct fp_ext *src) -{ - uprint("fsincos1\n"); - - return dest; -} - -struct fp_ext * -fp_fsincos2(struct fp_ext *dest, struct fp_ext *src) -{ - uprint("fsincos2\n"); - - return dest; -} - -struct fp_ext * -fp_fsincos3(struct fp_ext *dest, struct fp_ext *src) -{ - uprint("fsincos3\n"); - - return dest; -} - -struct fp_ext * -fp_fsincos4(struct fp_ext *dest, struct fp_ext *src) -{ - uprint("fsincos4\n"); - - return dest; -} - -struct fp_ext * -fp_fsincos5(struct fp_ext *dest, struct fp_ext *src) -{ - uprint("fsincos5\n"); - - return dest; -} - -struct fp_ext * -fp_fsincos6(struct fp_ext *dest, struct fp_ext *src) -{ - uprint("fsincos6\n"); - - return dest; -} - -struct fp_ext * -fp_fsincos7(struct fp_ext *dest, struct fp_ext *src) -{ - uprint("fsincos7\n"); - - return dest; -} diff --git a/ANDROID_3.4.5/arch/m68k/math-emu/fp_trig.h b/ANDROID_3.4.5/arch/m68k/math-emu/fp_trig.h deleted file mode 100644 index af8b247e..00000000 --- a/ANDROID_3.4.5/arch/m68k/math-emu/fp_trig.h +++ /dev/null @@ -1,32 +0,0 @@ -/* - - fp_trig.h: floating-point math routines for the Linux-m68k - floating point emulator. - - Copyright (c) 1998 David Huggins-Daines. - - I hereby give permission, free of charge, to copy, modify, and - redistribute this software, in source or binary form, provided that - the above copyright notice and the following disclaimer are included - in all such copies. - - THIS SOFTWARE IS PROVIDED "AS IS", WITH ABSOLUTELY NO WARRANTY, REAL - OR IMPLIED. - -*/ - -#ifndef FP_TRIG_H -#define FP_TRIG_H - -#include "fp_emu.h" - -/* floating point trigonometric instructions: - - the arguments to these are in the "internal" extended format, that - is, an "exploded" version of the 96-bit extended fp format used by - the 68881. - - they return a status code, which should end up in %d0, if all goes - well. */ - -#endif /* FP_TRIG__H */ diff --git a/ANDROID_3.4.5/arch/m68k/math-emu/fp_util.S b/ANDROID_3.4.5/arch/m68k/math-emu/fp_util.S deleted file mode 100644 index b093b85f..00000000 --- a/ANDROID_3.4.5/arch/m68k/math-emu/fp_util.S +++ /dev/null @@ -1,1454 +0,0 @@ -/* - * fp_util.S - * - * Copyright Roman Zippel, 1997. All rights reserved. - * - * Redistribution and use in source and binary forms, with or without - * modification, are permitted provided that the following conditions - * are met: - * 1. Redistributions of source code must retain the above copyright - * notice, and the entire permission notice in its entirety, - * including the disclaimer of warranties. - * 2. Redistributions in binary form must reproduce the above copyright - * notice, this list of conditions and the following disclaimer in the - * documentation and/or other materials provided with the distribution. - * 3. The name of the author may not be used to endorse or promote - * products derived from this software without specific prior - * written permission. - * - * ALTERNATIVELY, this product may be distributed under the terms of - * the GNU General Public License, in which case the provisions of the GPL are - * required INSTEAD OF the above restrictions. (This clause is - * necessary due to a potential bad interaction between the GPL and - * the restrictions contained in a BSD-style copyright.) - * - * THIS SOFTWARE IS PROVIDED ``AS IS'' AND ANY EXPRESS OR IMPLIED - * WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES - * OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE - * DISCLAIMED. IN NO EVENT SHALL THE AUTHOR BE LIABLE FOR ANY DIRECT, - * INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES - * (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR - * SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) - * HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, - * STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) - * ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED - * OF THE POSSIBILITY OF SUCH DAMAGE. - */ - -#include "fp_emu.h" - -/* - * Here are lots of conversion and normalization functions mainly - * used by fp_scan.S - * Note that these functions are optimized for "normal" numbers, - * these are handled first and exit as fast as possible, this is - * especially important for fp_normalize_ext/fp_conv_ext2ext, as - * it's called very often. - * The register usage is optimized for fp_scan.S and which register - * is currently at that time unused, be careful if you want change - * something here. %d0 and %d1 is always usable, sometimes %d2 (or - * only the lower half) most function have to return the %a0 - * unmodified, so that the caller can immediately reuse it. - */ - - .globl fp_ill, fp_end - - | exits from fp_scan: - | illegal instruction -fp_ill: - printf ,"fp_illegal\n" - rts - | completed instruction -fp_end: - tst.l (TASK_MM-8,%a2) - jmi 1f - tst.l (TASK_MM-4,%a2) - jmi 1f - tst.l (TASK_MM,%a2) - jpl 2f -1: printf ,"oops:%p,%p,%p\n",3,%a2@(TASK_MM-8),%a2@(TASK_MM-4),%a2@(TASK_MM) -2: clr.l %d0 - rts - - .globl fp_conv_long2ext, fp_conv_single2ext - .globl fp_conv_double2ext, fp_conv_ext2ext - .globl fp_normalize_ext, fp_normalize_double - .globl fp_normalize_single, fp_normalize_single_fast - .globl fp_conv_ext2double, fp_conv_ext2single - .globl fp_conv_ext2long, fp_conv_ext2short - .globl fp_conv_ext2byte - .globl fp_finalrounding_single, fp_finalrounding_single_fast - .globl fp_finalrounding_double - .globl fp_finalrounding, fp_finaltest, fp_final - -/* - * First several conversion functions from a source operand - * into the extended format. Note, that only fp_conv_ext2ext - * normalizes the number and is always called after the other - * conversion functions, which only move the information into - * fp_ext structure. - */ - - | fp_conv_long2ext: - | - | args: %d0 = source (32-bit long) - | %a0 = destination (ptr to struct fp_ext) - -fp_conv_long2ext: - printf PCONV,"l2e: %p -> %p(",2,%d0,%a0 - clr.l %d1 | sign defaults to zero - tst.l %d0 - jeq fp_l2e_zero | is source zero? - jpl 1f | positive? - moveq #1,%d1 - neg.l %d0 -1: swap %d1 - move.w #0x3fff+31,%d1 - move.l %d1,(%a0)+ | set sign / exp - move.l %d0,(%a0)+ | set mantissa - clr.l (%a0) - subq.l #8,%a0 | restore %a0 - printx PCONV,%a0@ - printf PCONV,")\n" - rts - | source is zero -fp_l2e_zero: - clr.l (%a0)+ - clr.l (%a0)+ - clr.l (%a0) - subq.l #8,%a0 - printx PCONV,%a0@ - printf PCONV,")\n" - rts - - | fp_conv_single2ext - | args: %d0 = source (single-precision fp value) - | %a0 = dest (struct fp_ext *) - -fp_conv_single2ext: - printf PCONV,"s2e: %p -> %p(",2,%d0,%a0 - move.l %d0,%d1 - lsl.l #8,%d0 | shift mantissa - lsr.l #8,%d1 | exponent / sign - lsr.l #7,%d1 - lsr.w #8,%d1 - jeq fp_s2e_small | zero / denormal? - cmp.w #0xff,%d1 | NaN / Inf? - jeq fp_s2e_large - bset #31,%d0 | set explizit bit - add.w #0x3fff-0x7f,%d1 | re-bias the exponent. -9: move.l %d1,(%a0)+ | fp_ext.sign, fp_ext.exp - move.l %d0,(%a0)+ | high lword of fp_ext.mant - clr.l (%a0) | low lword = 0 - subq.l #8,%a0 - printx PCONV,%a0@ - printf PCONV,")\n" - rts - | zeros and denormalized -fp_s2e_small: - | exponent is zero, so explizit bit is already zero too - tst.l %d0 - jeq 9b - move.w #0x4000-0x7f,%d1 - jra 9b - | infinities and NAN -fp_s2e_large: - bclr #31,%d0 | clear explizit bit - move.w #0x7fff,%d1 - jra 9b - -fp_conv_double2ext: -#ifdef FPU_EMU_DEBUG - getuser.l %a1@(0),%d0,fp_err_ua2,%a1 - getuser.l %a1@(4),%d1,fp_err_ua2,%a1 - printf PCONV,"d2e: %p%p -> %p(",3,%d0,%d1,%a0 -#endif - getuser.l (%a1)+,%d0,fp_err_ua2,%a1 - move.l %d0,%d1 - lsl.l #8,%d0 | shift high mantissa - lsl.l #3,%d0 - lsr.l #8,%d1 | exponent / sign - lsr.l #7,%d1 - lsr.w #5,%d1 - jeq fp_d2e_small | zero / denormal? - cmp.w #0x7ff,%d1 | NaN / Inf? - jeq fp_d2e_large - bset #31,%d0 | set explizit bit - add.w #0x3fff-0x3ff,%d1 | re-bias the exponent. -9: move.l %d1,(%a0)+ | fp_ext.sign, fp_ext.exp - move.l %d0,(%a0)+ - getuser.l (%a1)+,%d0,fp_err_ua2,%a1 - move.l %d0,%d1 - lsl.l #8,%d0 - lsl.l #3,%d0 - move.l %d0,(%a0) - moveq #21,%d0 - lsr.l %d0,%d1 - or.l %d1,-(%a0) - subq.l #4,%a0 - printx PCONV,%a0@ - printf PCONV,")\n" - rts - | zeros and denormalized -fp_d2e_small: - | exponent is zero, so explizit bit is already zero too - tst.l %d0 - jeq 9b - move.w #0x4000-0x3ff,%d1 - jra 9b - | infinities and NAN -fp_d2e_large: - bclr #31,%d0 | clear explizit bit - move.w #0x7fff,%d1 - jra 9b - - | fp_conv_ext2ext: - | originally used to get longdouble from userspace, now it's - | called before arithmetic operations to make sure the number - | is normalized [maybe rename it?]. - | args: %a0 = dest (struct fp_ext *) - | returns 0 in %d0 for a NaN, otherwise 1 - -fp_conv_ext2ext: - printf PCONV,"e2e: %p(",1,%a0 - printx PCONV,%a0@ - printf PCONV,"), " - move.l (%a0)+,%d0 - cmp.w #0x7fff,%d0 | Inf / NaN? - jeq fp_e2e_large - move.l (%a0),%d0 - jpl fp_e2e_small | zero / denorm? - | The high bit is set, so normalization is irrelevant. -fp_e2e_checkround: - subq.l #4,%a0 -#ifdef CONFIG_M68KFPU_EMU_EXTRAPREC - move.b (%a0),%d0 - jne fp_e2e_round -#endif - printf PCONV,"%p(",1,%a0 - printx PCONV,%a0@ - printf PCONV,")\n" - moveq #1,%d0 - rts -#ifdef CONFIG_M68KFPU_EMU_EXTRAPREC -fp_e2e_round: - fp_set_sr FPSR_EXC_INEX2 - clr.b (%a0) - move.w (FPD_RND,FPDATA),%d2 - jne fp_e2e_roundother | %d2 == 0, round to nearest - tst.b %d0 | test guard bit - jpl 9f | zero is closer - btst #0,(11,%a0) | test lsb bit - jne fp_e2e_doroundup | round to infinity - lsl.b #1,%d0 | check low bits - jeq 9f | round to zero -fp_e2e_doroundup: - addq.l #1,(8,%a0) - jcc 9f - addq.l #1,(4,%a0) - jcc 9f - move.w #0x8000,(4,%a0) - addq.w #1,(2,%a0) -9: printf PNORM,"%p(",1,%a0 - printx PNORM,%a0@ - printf PNORM,")\n" - rts -fp_e2e_roundother: - subq.w #2,%d2 - jcs 9b | %d2 < 2, round to zero - jhi 1f | %d2 > 2, round to +infinity - tst.b (1,%a0) | to -inf - jne fp_e2e_doroundup | negative, round to infinity - jra 9b | positive, round to zero -1: tst.b (1,%a0) | to +inf - jeq fp_e2e_doroundup | positive, round to infinity - jra 9b | negative, round to zero -#endif - | zeros and subnormals: - | try to normalize these anyway. -fp_e2e_small: - jne fp_e2e_small1 | high lword zero? - move.l (4,%a0),%d0 - jne fp_e2e_small2 -#ifdef CONFIG_M68KFPU_EMU_EXTRAPREC - clr.l %d0 - move.b (-4,%a0),%d0 - jne fp_e2e_small3 -#endif - | Genuine zero. - clr.w -(%a0) - subq.l #2,%a0 - printf PNORM,"%p(",1,%a0 - printx PNORM,%a0@ - printf PNORM,")\n" - moveq #1,%d0 - rts - | definitely subnormal, need to shift all 64 bits -fp_e2e_small1: - bfffo %d0{#0,#32},%d1 - move.w -(%a0),%d2 - sub.w %d1,%d2 - jcc 1f - | Pathologically small, denormalize. - add.w %d2,%d1 - clr.w %d2 -1: move.w %d2,(%a0)+ - move.w %d1,%d2 - jeq fp_e2e_checkround - | fancy 64-bit double-shift begins here - lsl.l %d2,%d0 - move.l %d0,(%a0)+ - move.l (%a0),%d0 - move.l %d0,%d1 - lsl.l %d2,%d0 - move.l %d0,(%a0) - neg.w %d2 - and.w #0x1f,%d2 - lsr.l %d2,%d1 - or.l %d1,-(%a0) -#ifdef CONFIG_M68KFPU_EMU_EXTRAPREC -fp_e2e_extra1: - clr.l %d0 - move.b (-4,%a0),%d0 - neg.w %d2 - add.w #24,%d2 - jcc 1f - clr.b (-4,%a0) - lsl.l %d2,%d0 - or.l %d0,(4,%a0) - jra fp_e2e_checkround -1: addq.w #8,%d2 - lsl.l %d2,%d0 - move.b %d0,(-4,%a0) - lsr.l #8,%d0 - or.l %d0,(4,%a0) -#endif - jra fp_e2e_checkround - | pathologically small subnormal -fp_e2e_small2: - bfffo %d0{#0,#32},%d1 - add.w #32,%d1 - move.w -(%a0),%d2 - sub.w %d1,%d2 - jcc 1f - | Beyond pathologically small, denormalize. - add.w %d2,%d1 - clr.w %d2 -1: move.w %d2,(%a0)+ - ext.l %d1 - jeq fp_e2e_checkround - clr.l (4,%a0) - sub.w #32,%d2 - jcs 1f - lsl.l %d1,%d0 | lower lword needs only to be shifted - move.l %d0,(%a0) | into the higher lword -#ifdef CONFIG_M68KFPU_EMU_EXTRAPREC - clr.l %d0 - move.b (-4,%a0),%d0 - clr.b (-4,%a0) - neg.w %d1 - add.w #32,%d1 - bfins %d0,(%a0){%d1,#8} -#endif - jra fp_e2e_checkround -1: neg.w %d1 | lower lword is splitted between - bfins %d0,(%a0){%d1,#32} | higher and lower lword -#ifndef CONFIG_M68KFPU_EMU_EXTRAPREC - jra fp_e2e_checkround -#else - move.w %d1,%d2 - jra fp_e2e_extra1 - | These are extremely small numbers, that will mostly end up as zero - | anyway, so this is only important for correct rounding. -fp_e2e_small3: - bfffo %d0{#24,#8},%d1 - add.w #40,%d1 - move.w -(%a0),%d2 - sub.w %d1,%d2 - jcc 1f - | Pathologically small, denormalize. - add.w %d2,%d1 - clr.w %d2 -1: move.w %d2,(%a0)+ - ext.l %d1 - jeq fp_e2e_checkround - cmp.w #8,%d1 - jcs 2f -1: clr.b (-4,%a0) - sub.w #64,%d1 - jcs 1f - add.w #24,%d1 - lsl.l %d1,%d0 - move.l %d0,(%a0) - jra fp_e2e_checkround -1: neg.w %d1 - bfins %d0,(%a0){%d1,#8} - jra fp_e2e_checkround -2: lsl.l %d1,%d0 - move.b %d0,(-4,%a0) - lsr.l #8,%d0 - move.b %d0,(7,%a0) - jra fp_e2e_checkround -#endif -1: move.l %d0,%d1 | lower lword is splitted between - lsl.l %d2,%d0 | higher and lower lword - move.l %d0,(%a0) - move.l %d1,%d0 - neg.w %d2 - add.w #32,%d2 - lsr.l %d2,%d0 - move.l %d0,-(%a0) - jra fp_e2e_checkround - | Infinities and NaNs -fp_e2e_large: - move.l (%a0)+,%d0 - jne 3f -1: tst.l (%a0) - jne 4f - moveq #1,%d0 -2: subq.l #8,%a0 - printf PCONV,"%p(",1,%a0 - printx PCONV,%a0@ - printf PCONV,")\n" - rts - | we have maybe a NaN, shift off the highest bit -3: lsl.l #1,%d0 - jeq 1b - | we have a NaN, clear the return value -4: clrl %d0 - jra 2b - - -/* - * Normalization functions. Call these on the output of general - * FP operators, and before any conversion into the destination - * formats. fp_normalize_ext has always to be called first, the - * following conversion functions expect an already normalized - * number. - */ - - | fp_normalize_ext: - | normalize an extended in extended (unpacked) format, basically - | it does the same as fp_conv_ext2ext, additionally it also does - | the necessary postprocessing checks. - | args: %a0 (struct fp_ext *) - | NOTE: it does _not_ modify %a0/%a1 and the upper word of %d2 - -fp_normalize_ext: - printf PNORM,"ne: %p(",1,%a0 - printx PNORM,%a0@ - printf PNORM,"), " - move.l (%a0)+,%d0 - cmp.w #0x7fff,%d0 | Inf / NaN? - jeq fp_ne_large - move.l (%a0),%d0 - jpl fp_ne_small | zero / denorm? - | The high bit is set, so normalization is irrelevant. -fp_ne_checkround: - subq.l #4,%a0 -#ifdef CONFIG_M68KFPU_EMU_EXTRAPREC - move.b (%a0),%d0 - jne fp_ne_round -#endif - printf PNORM,"%p(",1,%a0 - printx PNORM,%a0@ - printf PNORM,")\n" - rts -#ifdef CONFIG_M68KFPU_EMU_EXTRAPREC -fp_ne_round: - fp_set_sr FPSR_EXC_INEX2 - clr.b (%a0) - move.w (FPD_RND,FPDATA),%d2 - jne fp_ne_roundother | %d2 == 0, round to nearest - tst.b %d0 | test guard bit - jpl 9f | zero is closer - btst #0,(11,%a0) | test lsb bit - jne fp_ne_doroundup | round to infinity - lsl.b #1,%d0 | check low bits - jeq 9f | round to zero -fp_ne_doroundup: - addq.l #1,(8,%a0) - jcc 9f - addq.l #1,(4,%a0) - jcc 9f - addq.w #1,(2,%a0) - move.w #0x8000,(4,%a0) -9: printf PNORM,"%p(",1,%a0 - printx PNORM,%a0@ - printf PNORM,")\n" - rts -fp_ne_roundother: - subq.w #2,%d2 - jcs 9b | %d2 < 2, round to zero - jhi 1f | %d2 > 2, round to +infinity - tst.b (1,%a0) | to -inf - jne fp_ne_doroundup | negative, round to infinity - jra 9b | positive, round to zero -1: tst.b (1,%a0) | to +inf - jeq fp_ne_doroundup | positive, round to infinity - jra 9b | negative, round to zero -#endif - | Zeros and subnormal numbers - | These are probably merely subnormal, rather than "denormalized" - | numbers, so we will try to make them normal again. -fp_ne_small: - jne fp_ne_small1 | high lword zero? - move.l (4,%a0),%d0 - jne fp_ne_small2 -#ifdef CONFIG_M68KFPU_EMU_EXTRAPREC - clr.l %d0 - move.b (-4,%a0),%d0 - jne fp_ne_small3 -#endif - | Genuine zero. - clr.w -(%a0) - subq.l #2,%a0 - printf PNORM,"%p(",1,%a0 - printx PNORM,%a0@ - printf PNORM,")\n" - rts - | Subnormal. -fp_ne_small1: - bfffo %d0{#0,#32},%d1 - move.w -(%a0),%d2 - sub.w %d1,%d2 - jcc 1f - | Pathologically small, denormalize. - add.w %d2,%d1 - clr.w %d2 - fp_set_sr FPSR_EXC_UNFL -1: move.w %d2,(%a0)+ - move.w %d1,%d2 - jeq fp_ne_checkround - | This is exactly the same 64-bit double shift as seen above. - lsl.l %d2,%d0 - move.l %d0,(%a0)+ - move.l (%a0),%d0 - move.l %d0,%d1 - lsl.l %d2,%d0 - move.l %d0,(%a0) - neg.w %d2 - and.w #0x1f,%d2 - lsr.l %d2,%d1 - or.l %d1,-(%a0) -#ifdef CONFIG_M68KFPU_EMU_EXTRAPREC -fp_ne_extra1: - clr.l %d0 - move.b (-4,%a0),%d0 - neg.w %d2 - add.w #24,%d2 - jcc 1f - clr.b (-4,%a0) - lsl.l %d2,%d0 - or.l %d0,(4,%a0) - jra fp_ne_checkround -1: addq.w #8,%d2 - lsl.l %d2,%d0 - move.b %d0,(-4,%a0) - lsr.l #8,%d0 - or.l %d0,(4,%a0) -#endif - jra fp_ne_checkround - | May or may not be subnormal, if so, only 32 bits to shift. -fp_ne_small2: - bfffo %d0{#0,#32},%d1 - add.w #32,%d1 - move.w -(%a0),%d2 - sub.w %d1,%d2 - jcc 1f - | Beyond pathologically small, denormalize. - add.w %d2,%d1 - clr.w %d2 - fp_set_sr FPSR_EXC_UNFL -1: move.w %d2,(%a0)+ - ext.l %d1 - jeq fp_ne_checkround - clr.l (4,%a0) - sub.w #32,%d1 - jcs 1f - lsl.l %d1,%d0 | lower lword needs only to be shifted - move.l %d0,(%a0) | into the higher lword -#ifdef CONFIG_M68KFPU_EMU_EXTRAPREC - clr.l %d0 - move.b (-4,%a0),%d0 - clr.b (-4,%a0) - neg.w %d1 - add.w #32,%d1 - bfins %d0,(%a0){%d1,#8} -#endif - jra fp_ne_checkround -1: neg.w %d1 | lower lword is splitted between - bfins %d0,(%a0){%d1,#32} | higher and lower lword -#ifndef CONFIG_M68KFPU_EMU_EXTRAPREC - jra fp_ne_checkround -#else - move.w %d1,%d2 - jra fp_ne_extra1 - | These are extremely small numbers, that will mostly end up as zero - | anyway, so this is only important for correct rounding. -fp_ne_small3: - bfffo %d0{#24,#8},%d1 - add.w #40,%d1 - move.w -(%a0),%d2 - sub.w %d1,%d2 - jcc 1f - | Pathologically small, denormalize. - add.w %d2,%d1 - clr.w %d2 -1: move.w %d2,(%a0)+ - ext.l %d1 - jeq fp_ne_checkround - cmp.w #8,%d1 - jcs 2f -1: clr.b (-4,%a0) - sub.w #64,%d1 - jcs 1f - add.w #24,%d1 - lsl.l %d1,%d0 - move.l %d0,(%a0) - jra fp_ne_checkround -1: neg.w %d1 - bfins %d0,(%a0){%d1,#8} - jra fp_ne_checkround -2: lsl.l %d1,%d0 - move.b %d0,(-4,%a0) - lsr.l #8,%d0 - move.b %d0,(7,%a0) - jra fp_ne_checkround -#endif - | Infinities and NaNs, again, same as above. -fp_ne_large: - move.l (%a0)+,%d0 - jne 3f -1: tst.l (%a0) - jne 4f -2: subq.l #8,%a0 - printf PNORM,"%p(",1,%a0 - printx PNORM,%a0@ - printf PNORM,")\n" - rts - | we have maybe a NaN, shift off the highest bit -3: move.l %d0,%d1 - lsl.l #1,%d1 - jne 4f - clr.l (-4,%a0) - jra 1b - | we have a NaN, test if it is signaling -4: bset #30,%d0 - jne 2b - fp_set_sr FPSR_EXC_SNAN - move.l %d0,(-4,%a0) - jra 2b - - | these next two do rounding as per the IEEE standard. - | values for the rounding modes appear to be: - | 0: Round to nearest - | 1: Round to zero - | 2: Round to -Infinity - | 3: Round to +Infinity - | both functions expect that fp_normalize was already - | called (and extended argument is already normalized - | as far as possible), these are used if there is different - | rounding precision is selected and before converting - | into single/double - - | fp_normalize_double: - | normalize an extended with double (52-bit) precision - | args: %a0 (struct fp_ext *) - -fp_normalize_double: - printf PNORM,"nd: %p(",1,%a0 - printx PNORM,%a0@ - printf PNORM,"), " - move.l (%a0)+,%d2 - tst.w %d2 - jeq fp_nd_zero | zero / denormalized - cmp.w #0x7fff,%d2 - jeq fp_nd_huge | NaN / infinitive. - sub.w #0x4000-0x3ff,%d2 | will the exponent fit? - jcs fp_nd_small | too small. - cmp.w #0x7fe,%d2 - jcc fp_nd_large | too big. - addq.l #4,%a0 - move.l (%a0),%d0 | low lword of mantissa - | now, round off the low 11 bits. -fp_nd_round: - moveq #21,%d1 - lsl.l %d1,%d0 | keep 11 low bits. - jne fp_nd_checkround | Are they non-zero? - | nothing to do here -9: subq.l #8,%a0 - printf PNORM,"%p(",1,%a0 - printx PNORM,%a0@ - printf PNORM,")\n" - rts - | Be careful with the X bit! It contains the lsb - | from the shift above, it is needed for round to nearest. -fp_nd_checkround: - fp_set_sr FPSR_EXC_INEX2 | INEX2 bit - and.w #0xf800,(2,%a0) | clear bits 0-10 - move.w (FPD_RND,FPDATA),%d2 | rounding mode - jne 2f | %d2 == 0, round to nearest - tst.l %d0 | test guard bit - jpl 9b | zero is closer - | here we test the X bit by adding it to %d2 - clr.w %d2 | first set z bit, addx only clears it - addx.w %d2,%d2 | test lsb bit - | IEEE754-specified "round to even" behaviour. If the guard - | bit is set, then the number is odd, so rounding works like - | in grade-school arithmetic (i.e. 1.5 rounds to 2.0) - | Otherwise, an equal distance rounds towards zero, so as not - | to produce an odd number. This is strange, but it is what - | the standard says. - jne fp_nd_doroundup | round to infinity - lsl.l #1,%d0 | check low bits - jeq 9b | round to zero -fp_nd_doroundup: - | round (the mantissa, that is) towards infinity - add.l #0x800,(%a0) - jcc 9b | no overflow, good. - addq.l #1,-(%a0) | extend to high lword - jcc 1f | no overflow, good. - | Yow! we have managed to overflow the mantissa. Since this - | only happens when %d1 was 0xfffff800, it is now zero, so - | reset the high bit, and increment the exponent. - move.w #0x8000,(%a0) - addq.w #1,-(%a0) - cmp.w #0x43ff,(%a0)+ | exponent now overflown? - jeq fp_nd_large | yes, so make it infinity. -1: subq.l #4,%a0 - printf PNORM,"%p(",1,%a0 - printx PNORM,%a0@ - printf PNORM,")\n" - rts -2: subq.w #2,%d2 - jcs 9b | %d2 < 2, round to zero - jhi 3f | %d2 > 2, round to +infinity - | Round to +Inf or -Inf. High word of %d2 contains the - | sign of the number, by the way. - swap %d2 | to -inf - tst.b %d2 - jne fp_nd_doroundup | negative, round to infinity - jra 9b | positive, round to zero -3: swap %d2 | to +inf - tst.b %d2 - jeq fp_nd_doroundup | positive, round to infinity - jra 9b | negative, round to zero - | Exponent underflow. Try to make a denormal, and set it to - | the smallest possible fraction if this fails. -fp_nd_small: - fp_set_sr FPSR_EXC_UNFL | set UNFL bit - move.w #0x3c01,(-2,%a0) | 2**-1022 - neg.w %d2 | degree of underflow - cmp.w #32,%d2 | single or double shift? - jcc 1f - | Again, another 64-bit double shift. - move.l (%a0),%d0 - move.l %d0,%d1 - lsr.l %d2,%d0 - move.l %d0,(%a0)+ - move.l (%a0),%d0 - lsr.l %d2,%d0 - neg.w %d2 - add.w #32,%d2 - lsl.l %d2,%d1 - or.l %d1,%d0 - move.l (%a0),%d1 - move.l %d0,(%a0) - | Check to see if we shifted off any significant bits - lsl.l %d2,%d1 - jeq fp_nd_round | Nope, round. - bset #0,%d0 | Yes, so set the "sticky bit". - jra fp_nd_round | Now, round. - | Another 64-bit single shift and store -1: sub.w #32,%d2 - cmp.w #32,%d2 | Do we really need to shift? - jcc 2f | No, the number is too small. - move.l (%a0),%d0 - clr.l (%a0)+ - move.l %d0,%d1 - lsr.l %d2,%d0 - neg.w %d2 - add.w #32,%d2 - | Again, check to see if we shifted off any significant bits. - tst.l (%a0) - jeq 1f - bset #0,%d0 | Sticky bit. -1: move.l %d0,(%a0) - lsl.l %d2,%d1 - jeq fp_nd_round - bset #0,%d0 - jra fp_nd_round - | Sorry, the number is just too small. -2: clr.l (%a0)+ - clr.l (%a0) - moveq #1,%d0 | Smallest possible fraction, - jra fp_nd_round | round as desired. - | zero and denormalized -fp_nd_zero: - tst.l (%a0)+ - jne 1f - tst.l (%a0) - jne 1f - subq.l #8,%a0 - printf PNORM,"%p(",1,%a0 - printx PNORM,%a0@ - printf PNORM,")\n" - rts | zero. nothing to do. - | These are not merely subnormal numbers, but true denormals, - | i.e. pathologically small (exponent is 2**-16383) numbers. - | It is clearly impossible for even a normal extended number - | with that exponent to fit into double precision, so just - | write these ones off as "too darn small". -1: fp_set_sr FPSR_EXC_UNFL | Set UNFL bit - clr.l (%a0) - clr.l -(%a0) - move.w #0x3c01,-(%a0) | i.e. 2**-1022 - addq.l #6,%a0 - moveq #1,%d0 - jra fp_nd_round | round. - | Exponent overflow. Just call it infinity. -fp_nd_large: - move.w #0x7ff,%d0 - and.w (6,%a0),%d0 - jeq 1f - fp_set_sr FPSR_EXC_INEX2 -1: fp_set_sr FPSR_EXC_OVFL - move.w (FPD_RND,FPDATA),%d2 - jne 3f | %d2 = 0 round to nearest -1: move.w #0x7fff,(-2,%a0) - clr.l (%a0)+ - clr.l (%a0) -2: subq.l #8,%a0 - printf PNORM,"%p(",1,%a0 - printx PNORM,%a0@ - printf PNORM,")\n" - rts -3: subq.w #2,%d2 - jcs 5f | %d2 < 2, round to zero - jhi 4f | %d2 > 2, round to +infinity - tst.b (-3,%a0) | to -inf - jne 1b - jra 5f -4: tst.b (-3,%a0) | to +inf - jeq 1b -5: move.w #0x43fe,(-2,%a0) - moveq #-1,%d0 - move.l %d0,(%a0)+ - move.w #0xf800,%d0 - move.l %d0,(%a0) - jra 2b - | Infinities or NaNs -fp_nd_huge: - subq.l #4,%a0 - printf PNORM,"%p(",1,%a0 - printx PNORM,%a0@ - printf PNORM,")\n" - rts - - | fp_normalize_single: - | normalize an extended with single (23-bit) precision - | args: %a0 (struct fp_ext *) - -fp_normalize_single: - printf PNORM,"ns: %p(",1,%a0 - printx PNORM,%a0@ - printf PNORM,") " - addq.l #2,%a0 - move.w (%a0)+,%d2 - jeq fp_ns_zero | zero / denormalized - cmp.w #0x7fff,%d2 - jeq fp_ns_huge | NaN / infinitive. - sub.w #0x4000-0x7f,%d2 | will the exponent fit? - jcs fp_ns_small | too small. - cmp.w #0xfe,%d2 - jcc fp_ns_large | too big. - move.l (%a0)+,%d0 | get high lword of mantissa -fp_ns_round: - tst.l (%a0) | check the low lword - jeq 1f - | Set a sticky bit if it is non-zero. This should only - | affect the rounding in what would otherwise be equal- - | distance situations, which is what we want it to do. - bset #0,%d0 -1: clr.l (%a0) | zap it from memory. - | now, round off the low 8 bits of the hi lword. - tst.b %d0 | 8 low bits. - jne fp_ns_checkround | Are they non-zero? - | nothing to do here - subq.l #8,%a0 - printf PNORM,"%p(",1,%a0 - printx PNORM,%a0@ - printf PNORM,")\n" - rts -fp_ns_checkround: - fp_set_sr FPSR_EXC_INEX2 | INEX2 bit - clr.b -(%a0) | clear low byte of high lword - subq.l #3,%a0 - move.w (FPD_RND,FPDATA),%d2 | rounding mode - jne 2f | %d2 == 0, round to nearest - tst.b %d0 | test guard bit - jpl 9f | zero is closer - btst #8,%d0 | test lsb bit - | round to even behaviour, see above. - jne fp_ns_doroundup | round to infinity - lsl.b #1,%d0 | check low bits - jeq 9f | round to zero -fp_ns_doroundup: - | round (the mantissa, that is) towards infinity - add.l #0x100,(%a0) - jcc 9f | no overflow, good. - | Overflow. This means that the %d1 was 0xffffff00, so it - | is now zero. We will set the mantissa to reflect this, and - | increment the exponent (checking for overflow there too) - move.w #0x8000,(%a0) - addq.w #1,-(%a0) - cmp.w #0x407f,(%a0)+ | exponent now overflown? - jeq fp_ns_large | yes, so make it infinity. -9: subq.l #4,%a0 - printf PNORM,"%p(",1,%a0 - printx PNORM,%a0@ - printf PNORM,")\n" - rts - | check nondefault rounding modes -2: subq.w #2,%d2 - jcs 9b | %d2 < 2, round to zero - jhi 3f | %d2 > 2, round to +infinity - tst.b (-3,%a0) | to -inf - jne fp_ns_doroundup | negative, round to infinity - jra 9b | positive, round to zero -3: tst.b (-3,%a0) | to +inf - jeq fp_ns_doroundup | positive, round to infinity - jra 9b | negative, round to zero - | Exponent underflow. Try to make a denormal, and set it to - | the smallest possible fraction if this fails. -fp_ns_small: - fp_set_sr FPSR_EXC_UNFL | set UNFL bit - move.w #0x3f81,(-2,%a0) | 2**-126 - neg.w %d2 | degree of underflow - cmp.w #32,%d2 | single or double shift? - jcc 2f - | a 32-bit shift. - move.l (%a0),%d0 - move.l %d0,%d1 - lsr.l %d2,%d0 - move.l %d0,(%a0)+ - | Check to see if we shifted off any significant bits. - neg.w %d2 - add.w #32,%d2 - lsl.l %d2,%d1 - jeq 1f - bset #0,%d0 | Sticky bit. - | Check the lower lword -1: tst.l (%a0) - jeq fp_ns_round - clr (%a0) - bset #0,%d0 | Sticky bit. - jra fp_ns_round - | Sorry, the number is just too small. -2: clr.l (%a0)+ - clr.l (%a0) - moveq #1,%d0 | Smallest possible fraction, - jra fp_ns_round | round as desired. - | Exponent overflow. Just call it infinity. -fp_ns_large: - tst.b (3,%a0) - jeq 1f - fp_set_sr FPSR_EXC_INEX2 -1: fp_set_sr FPSR_EXC_OVFL - move.w (FPD_RND,FPDATA),%d2 - jne 3f | %d2 = 0 round to nearest -1: move.w #0x7fff,(-2,%a0) - clr.l (%a0)+ - clr.l (%a0) -2: subq.l #8,%a0 - printf PNORM,"%p(",1,%a0 - printx PNORM,%a0@ - printf PNORM,")\n" - rts -3: subq.w #2,%d2 - jcs 5f | %d2 < 2, round to zero - jhi 4f | %d2 > 2, round to +infinity - tst.b (-3,%a0) | to -inf - jne 1b - jra 5f -4: tst.b (-3,%a0) | to +inf - jeq 1b -5: move.w #0x407e,(-2,%a0) - move.l #0xffffff00,(%a0)+ - clr.l (%a0) - jra 2b - | zero and denormalized -fp_ns_zero: - tst.l (%a0)+ - jne 1f - tst.l (%a0) - jne 1f - subq.l #8,%a0 - printf PNORM,"%p(",1,%a0 - printx PNORM,%a0@ - printf PNORM,")\n" - rts | zero. nothing to do. - | These are not merely subnormal numbers, but true denormals, - | i.e. pathologically small (exponent is 2**-16383) numbers. - | It is clearly impossible for even a normal extended number - | with that exponent to fit into single precision, so just - | write these ones off as "too darn small". -1: fp_set_sr FPSR_EXC_UNFL | Set UNFL bit - clr.l (%a0) - clr.l -(%a0) - move.w #0x3f81,-(%a0) | i.e. 2**-126 - addq.l #6,%a0 - moveq #1,%d0 - jra fp_ns_round | round. - | Infinities or NaNs -fp_ns_huge: - subq.l #4,%a0 - printf PNORM,"%p(",1,%a0 - printx PNORM,%a0@ - printf PNORM,")\n" - rts - - | fp_normalize_single_fast: - | normalize an extended with single (23-bit) precision - | this is only used by fsgldiv/fsgdlmul, where the - | operand is not completly normalized. - | args: %a0 (struct fp_ext *) - -fp_normalize_single_fast: - printf PNORM,"nsf: %p(",1,%a0 - printx PNORM,%a0@ - printf PNORM,") " - addq.l #2,%a0 - move.w (%a0)+,%d2 - cmp.w #0x7fff,%d2 - jeq fp_nsf_huge | NaN / infinitive. - move.l (%a0)+,%d0 | get high lword of mantissa -fp_nsf_round: - tst.l (%a0) | check the low lword - jeq 1f - | Set a sticky bit if it is non-zero. This should only - | affect the rounding in what would otherwise be equal- - | distance situations, which is what we want it to do. - bset #0,%d0 -1: clr.l (%a0) | zap it from memory. - | now, round off the low 8 bits of the hi lword. - tst.b %d0 | 8 low bits. - jne fp_nsf_checkround | Are they non-zero? - | nothing to do here - subq.l #8,%a0 - printf PNORM,"%p(",1,%a0 - printx PNORM,%a0@ - printf PNORM,")\n" - rts -fp_nsf_checkround: - fp_set_sr FPSR_EXC_INEX2 | INEX2 bit - clr.b -(%a0) | clear low byte of high lword - subq.l #3,%a0 - move.w (FPD_RND,FPDATA),%d2 | rounding mode - jne 2f | %d2 == 0, round to nearest - tst.b %d0 | test guard bit - jpl 9f | zero is closer - btst #8,%d0 | test lsb bit - | round to even behaviour, see above. - jne fp_nsf_doroundup | round to infinity - lsl.b #1,%d0 | check low bits - jeq 9f | round to zero -fp_nsf_doroundup: - | round (the mantissa, that is) towards infinity - add.l #0x100,(%a0) - jcc 9f | no overflow, good. - | Overflow. This means that the %d1 was 0xffffff00, so it - | is now zero. We will set the mantissa to reflect this, and - | increment the exponent (checking for overflow there too) - move.w #0x8000,(%a0) - addq.w #1,-(%a0) - cmp.w #0x407f,(%a0)+ | exponent now overflown? - jeq fp_nsf_large | yes, so make it infinity. -9: subq.l #4,%a0 - printf PNORM,"%p(",1,%a0 - printx PNORM,%a0@ - printf PNORM,")\n" - rts - | check nondefault rounding modes -2: subq.w #2,%d2 - jcs 9b | %d2 < 2, round to zero - jhi 3f | %d2 > 2, round to +infinity - tst.b (-3,%a0) | to -inf - jne fp_nsf_doroundup | negative, round to infinity - jra 9b | positive, round to zero -3: tst.b (-3,%a0) | to +inf - jeq fp_nsf_doroundup | positive, round to infinity - jra 9b | negative, round to zero - | Exponent overflow. Just call it infinity. -fp_nsf_large: - tst.b (3,%a0) - jeq 1f - fp_set_sr FPSR_EXC_INEX2 -1: fp_set_sr FPSR_EXC_OVFL - move.w (FPD_RND,FPDATA),%d2 - jne 3f | %d2 = 0 round to nearest -1: move.w #0x7fff,(-2,%a0) - clr.l (%a0)+ - clr.l (%a0) -2: subq.l #8,%a0 - printf PNORM,"%p(",1,%a0 - printx PNORM,%a0@ - printf PNORM,")\n" - rts -3: subq.w #2,%d2 - jcs 5f | %d2 < 2, round to zero - jhi 4f | %d2 > 2, round to +infinity - tst.b (-3,%a0) | to -inf - jne 1b - jra 5f -4: tst.b (-3,%a0) | to +inf - jeq 1b -5: move.w #0x407e,(-2,%a0) - move.l #0xffffff00,(%a0)+ - clr.l (%a0) - jra 2b - | Infinities or NaNs -fp_nsf_huge: - subq.l #4,%a0 - printf PNORM,"%p(",1,%a0 - printx PNORM,%a0@ - printf PNORM,")\n" - rts - - | conv_ext2int (macro): - | Generates a subroutine that converts an extended value to an - | integer of a given size, again, with the appropriate type of - | rounding. - - | Macro arguments: - | s: size, as given in an assembly instruction. - | b: number of bits in that size. - - | Subroutine arguments: - | %a0: source (struct fp_ext *) - - | Returns the integer in %d0 (like it should) - -.macro conv_ext2int s,b - .set inf,(1<<(\b-1))-1 | i.e. MAXINT - printf PCONV,"e2i%d: %p(",2,#\b,%a0 - printx PCONV,%a0@ - printf PCONV,") " - addq.l #2,%a0 - move.w (%a0)+,%d2 | exponent - jeq fp_e2i_zero\b | zero / denorm (== 0, here) - cmp.w #0x7fff,%d2 - jeq fp_e2i_huge\b | Inf / NaN - sub.w #0x3ffe,%d2 - jcs fp_e2i_small\b - cmp.w #\b,%d2 - jhi fp_e2i_large\b - move.l (%a0),%d0 - move.l %d0,%d1 - lsl.l %d2,%d1 - jne fp_e2i_round\b - tst.l (4,%a0) - jne fp_e2i_round\b - neg.w %d2 - add.w #32,%d2 - lsr.l %d2,%d0 -9: tst.w (-4,%a0) - jne 1f - tst.\s %d0 - jmi fp_e2i_large\b - printf PCONV,"-> %p\n",1,%d0 - rts -1: neg.\s %d0 - jeq 1f - jpl fp_e2i_large\b -1: printf PCONV,"-> %p\n",1,%d0 - rts -fp_e2i_round\b: - fp_set_sr FPSR_EXC_INEX2 | INEX2 bit - neg.w %d2 - add.w #32,%d2 - .if \b>16 - jeq 5f - .endif - lsr.l %d2,%d0 - move.w (FPD_RND,FPDATA),%d2 | rounding mode - jne 2f | %d2 == 0, round to nearest - tst.l %d1 | test guard bit - jpl 9b | zero is closer - btst %d2,%d0 | test lsb bit (%d2 still 0) - jne fp_e2i_doroundup\b - lsl.l #1,%d1 | check low bits - jne fp_e2i_doroundup\b - tst.l (4,%a0) - jeq 9b -fp_e2i_doroundup\b: - addq.l #1,%d0 - jra 9b - | check nondefault rounding modes -2: subq.w #2,%d2 - jcs 9b | %d2 < 2, round to zero - jhi 3f | %d2 > 2, round to +infinity - tst.w (-4,%a0) | to -inf - jne fp_e2i_doroundup\b | negative, round to infinity - jra 9b | positive, round to zero -3: tst.w (-4,%a0) | to +inf - jeq fp_e2i_doroundup\b | positive, round to infinity - jra 9b | negative, round to zero - | we are only want -2**127 get correctly rounded here, - | since the guard bit is in the lower lword. - | everything else ends up anyway as overflow. - .if \b>16 -5: move.w (FPD_RND,FPDATA),%d2 | rounding mode - jne 2b | %d2 == 0, round to nearest - move.l (4,%a0),%d1 | test guard bit - jpl 9b | zero is closer - lsl.l #1,%d1 | check low bits - jne fp_e2i_doroundup\b - jra 9b - .endif -fp_e2i_zero\b: - clr.l %d0 - tst.l (%a0)+ - jne 1f - tst.l (%a0) - jeq 3f -1: subq.l #4,%a0 - fp_clr_sr FPSR_EXC_UNFL | fp_normalize_ext has set this bit -fp_e2i_small\b: - fp_set_sr FPSR_EXC_INEX2 - clr.l %d0 - move.w (FPD_RND,FPDATA),%d2 | rounding mode - subq.w #2,%d2 - jcs 3f | %d2 < 2, round to nearest/zero - jhi 2f | %d2 > 2, round to +infinity - tst.w (-4,%a0) | to -inf - jeq 3f - subq.\s #1,%d0 - jra 3f -2: tst.w (-4,%a0) | to +inf - jne 3f - addq.\s #1,%d0 -3: printf PCONV,"-> %p\n",1,%d0 - rts -fp_e2i_large\b: - fp_set_sr FPSR_EXC_OPERR - move.\s #inf,%d0 - tst.w (-4,%a0) - jeq 1f - addq.\s #1,%d0 -1: printf PCONV,"-> %p\n",1,%d0 - rts -fp_e2i_huge\b: - move.\s (%a0),%d0 - tst.l (%a0) - jne 1f - tst.l (%a0) - jeq fp_e2i_large\b - | fp_normalize_ext has set this bit already - | and made the number nonsignaling -1: fp_tst_sr FPSR_EXC_SNAN - jne 1f - fp_set_sr FPSR_EXC_OPERR -1: printf PCONV,"-> %p\n",1,%d0 - rts -.endm - -fp_conv_ext2long: - conv_ext2int l,32 - -fp_conv_ext2short: - conv_ext2int w,16 - -fp_conv_ext2byte: - conv_ext2int b,8 - -fp_conv_ext2double: - jsr fp_normalize_double - printf PCONV,"e2d: %p(",1,%a0 - printx PCONV,%a0@ - printf PCONV,"), " - move.l (%a0)+,%d2 - cmp.w #0x7fff,%d2 - jne 1f - move.w #0x7ff,%d2 - move.l (%a0)+,%d0 - jra 2f -1: sub.w #0x3fff-0x3ff,%d2 - move.l (%a0)+,%d0 - jmi 2f - clr.w %d2 -2: lsl.w #5,%d2 - lsl.l #7,%d2 - lsl.l #8,%d2 - move.l %d0,%d1 - lsl.l #1,%d0 - lsr.l #4,%d0 - lsr.l #8,%d0 - or.l %d2,%d0 - putuser.l %d0,(%a1)+,fp_err_ua2,%a1 - moveq #21,%d0 - lsl.l %d0,%d1 - move.l (%a0),%d0 - lsr.l #4,%d0 - lsr.l #7,%d0 - or.l %d1,%d0 - putuser.l %d0,(%a1),fp_err_ua2,%a1 -#ifdef FPU_EMU_DEBUG - getuser.l %a1@(-4),%d0,fp_err_ua2,%a1 - getuser.l %a1@(0),%d1,fp_err_ua2,%a1 - printf PCONV,"%p(%08x%08x)\n",3,%a1,%d0,%d1 -#endif - rts - -fp_conv_ext2single: - jsr fp_normalize_single - printf PCONV,"e2s: %p(",1,%a0 - printx PCONV,%a0@ - printf PCONV,"), " - move.l (%a0)+,%d1 - cmp.w #0x7fff,%d1 - jne 1f - move.w #0xff,%d1 - move.l (%a0)+,%d0 - jra 2f -1: sub.w #0x3fff-0x7f,%d1 - move.l (%a0)+,%d0 - jmi 2f - clr.w %d1 -2: lsl.w #8,%d1 - lsl.l #7,%d1 - lsl.l #8,%d1 - bclr #31,%d0 - lsr.l #8,%d0 - or.l %d1,%d0 - printf PCONV,"%08x\n",1,%d0 - rts - - | special return addresses for instr that - | encode the rounding precision in the opcode - | (e.g. fsmove,fdmove) - -fp_finalrounding_single: - addq.l #8,%sp - jsr fp_normalize_ext - jsr fp_normalize_single - jra fp_finaltest - -fp_finalrounding_single_fast: - addq.l #8,%sp - jsr fp_normalize_ext - jsr fp_normalize_single_fast - jra fp_finaltest - -fp_finalrounding_double: - addq.l #8,%sp - jsr fp_normalize_ext - jsr fp_normalize_double - jra fp_finaltest - - | fp_finaltest: - | set the emulated status register based on the outcome of an - | emulated instruction. - -fp_finalrounding: - addq.l #8,%sp -| printf ,"f: %p\n",1,%a0 - jsr fp_normalize_ext - move.w (FPD_PREC,FPDATA),%d0 - subq.w #1,%d0 - jcs fp_finaltest - jne 1f - jsr fp_normalize_single - jra 2f -1: jsr fp_normalize_double -2:| printf ,"f: %p\n",1,%a0 -fp_finaltest: - | First, we do some of the obvious tests for the exception - | status byte and condition code bytes of fp_sr here, so that - | they do not have to be handled individually by every - | emulated instruction. - clr.l %d0 - addq.l #1,%a0 - tst.b (%a0)+ | sign - jeq 1f - bset #FPSR_CC_NEG-24,%d0 | N bit -1: cmp.w #0x7fff,(%a0)+ | exponent - jeq 2f - | test for zero - moveq #FPSR_CC_Z-24,%d1 - tst.l (%a0)+ - jne 9f - tst.l (%a0) - jne 9f - jra 8f - | infinitiv and NAN -2: moveq #FPSR_CC_NAN-24,%d1 - move.l (%a0)+,%d2 - lsl.l #1,%d2 | ignore high bit - jne 8f - tst.l (%a0) - jne 8f - moveq #FPSR_CC_INF-24,%d1 -8: bset %d1,%d0 -9: move.b %d0,(FPD_FPSR+0,FPDATA) | set condition test result - | move instructions enter here - | Here, we test things in the exception status byte, and set - | other things in the accrued exception byte accordingly. - | Emulated instructions can set various things in the former, - | as defined in fp_emu.h. -fp_final: - move.l (FPD_FPSR,FPDATA),%d0 -#if 0 - btst #FPSR_EXC_SNAN,%d0 | EXC_SNAN - jne 1f - btst #FPSR_EXC_OPERR,%d0 | EXC_OPERR - jeq 2f -1: bset #FPSR_AEXC_IOP,%d0 | set IOP bit -2: btst #FPSR_EXC_OVFL,%d0 | EXC_OVFL - jeq 1f - bset #FPSR_AEXC_OVFL,%d0 | set OVFL bit -1: btst #FPSR_EXC_UNFL,%d0 | EXC_UNFL - jeq 1f - btst #FPSR_EXC_INEX2,%d0 | EXC_INEX2 - jeq 1f - bset #FPSR_AEXC_UNFL,%d0 | set UNFL bit -1: btst #FPSR_EXC_DZ,%d0 | EXC_INEX1 - jeq 1f - bset #FPSR_AEXC_DZ,%d0 | set DZ bit -1: btst #FPSR_EXC_OVFL,%d0 | EXC_OVFL - jne 1f - btst #FPSR_EXC_INEX2,%d0 | EXC_INEX2 - jne 1f - btst #FPSR_EXC_INEX1,%d0 | EXC_INEX1 - jeq 2f -1: bset #FPSR_AEXC_INEX,%d0 | set INEX bit -2: move.l %d0,(FPD_FPSR,FPDATA) -#else - | same as above, greatly optimized, but untested (yet) - move.l %d0,%d2 - lsr.l #5,%d0 - move.l %d0,%d1 - lsr.l #4,%d1 - or.l %d0,%d1 - and.b #0x08,%d1 - move.l %d2,%d0 - lsr.l #6,%d0 - or.l %d1,%d0 - move.l %d2,%d1 - lsr.l #4,%d1 - or.b #0xdf,%d1 - and.b %d1,%d0 - move.l %d2,%d1 - lsr.l #7,%d1 - and.b #0x80,%d1 - or.b %d1,%d0 - and.b #0xf8,%d0 - or.b %d0,%d2 - move.l %d2,(FPD_FPSR,FPDATA) -#endif - move.b (FPD_FPSR+2,FPDATA),%d0 - and.b (FPD_FPCR+2,FPDATA),%d0 - jeq 1f - printf ,"send signal!!!\n" -1: jra fp_end diff --git a/ANDROID_3.4.5/arch/m68k/math-emu/multi_arith.h b/ANDROID_3.4.5/arch/m68k/math-emu/multi_arith.h deleted file mode 100644 index 4b5eb3d8..00000000 --- a/ANDROID_3.4.5/arch/m68k/math-emu/multi_arith.h +++ /dev/null @@ -1,289 +0,0 @@ -/* multi_arith.h: multi-precision integer arithmetic functions, needed - to do extended-precision floating point. - - (c) 1998 David Huggins-Daines. - - Somewhat based on arch/alpha/math-emu/ieee-math.c, which is (c) - David Mosberger-Tang. - - You may copy, modify, and redistribute this file under the terms of - the GNU General Public License, version 2, or any later version, at - your convenience. */ - -/* Note: - - These are not general multi-precision math routines. Rather, they - implement the subset of integer arithmetic that we need in order to - multiply, divide, and normalize 128-bit unsigned mantissae. */ - -#ifndef MULTI_ARITH_H -#define MULTI_ARITH_H - -static inline void fp_denormalize(struct fp_ext *reg, unsigned int cnt) -{ - reg->exp += cnt; - - switch (cnt) { - case 0 ... 8: - reg->lowmant = reg->mant.m32[1] << (8 - cnt); - reg->mant.m32[1] = (reg->mant.m32[1] >> cnt) | - (reg->mant.m32[0] << (32 - cnt)); - reg->mant.m32[0] = reg->mant.m32[0] >> cnt; - break; - case 9 ... 32: - reg->lowmant = reg->mant.m32[1] >> (cnt - 8); - if (reg->mant.m32[1] << (40 - cnt)) - reg->lowmant |= 1; - reg->mant.m32[1] = (reg->mant.m32[1] >> cnt) | - (reg->mant.m32[0] << (32 - cnt)); - reg->mant.m32[0] = reg->mant.m32[0] >> cnt; - break; - case 33 ... 39: - asm volatile ("bfextu %1{%2,#8},%0" : "=d" (reg->lowmant) - : "m" (reg->mant.m32[0]), "d" (64 - cnt)); - if (reg->mant.m32[1] << (40 - cnt)) - reg->lowmant |= 1; - reg->mant.m32[1] = reg->mant.m32[0] >> (cnt - 32); - reg->mant.m32[0] = 0; - break; - case 40 ... 71: - reg->lowmant = reg->mant.m32[0] >> (cnt - 40); - if ((reg->mant.m32[0] << (72 - cnt)) || reg->mant.m32[1]) - reg->lowmant |= 1; - reg->mant.m32[1] = reg->mant.m32[0] >> (cnt - 32); - reg->mant.m32[0] = 0; - break; - default: - reg->lowmant = reg->mant.m32[0] || reg->mant.m32[1]; - reg->mant.m32[0] = 0; - reg->mant.m32[1] = 0; - break; - } -} - -static inline int fp_overnormalize(struct fp_ext *reg) -{ - int shift; - - if (reg->mant.m32[0]) { - asm ("bfffo %1{#0,#32},%0" : "=d" (shift) : "dm" (reg->mant.m32[0])); - reg->mant.m32[0] = (reg->mant.m32[0] << shift) | (reg->mant.m32[1] >> (32 - shift)); - reg->mant.m32[1] = (reg->mant.m32[1] << shift); - } else { - asm ("bfffo %1{#0,#32},%0" : "=d" (shift) : "dm" (reg->mant.m32[1])); - reg->mant.m32[0] = (reg->mant.m32[1] << shift); - reg->mant.m32[1] = 0; - shift += 32; - } - - return shift; -} - -static inline int fp_addmant(struct fp_ext *dest, struct fp_ext *src) -{ - int carry; - - /* we assume here, gcc only insert move and a clr instr */ - asm volatile ("add.b %1,%0" : "=d,g" (dest->lowmant) - : "g,d" (src->lowmant), "0,0" (dest->lowmant)); - asm volatile ("addx.l %1,%0" : "=d" (dest->mant.m32[1]) - : "d" (src->mant.m32[1]), "0" (dest->mant.m32[1])); - asm volatile ("addx.l %1,%0" : "=d" (dest->mant.m32[0]) - : "d" (src->mant.m32[0]), "0" (dest->mant.m32[0])); - asm volatile ("addx.l %0,%0" : "=d" (carry) : "0" (0)); - - return carry; -} - -static inline int fp_addcarry(struct fp_ext *reg) -{ - if (++reg->exp == 0x7fff) { - if (reg->mant.m64) - fp_set_sr(FPSR_EXC_INEX2); - reg->mant.m64 = 0; - fp_set_sr(FPSR_EXC_OVFL); - return 0; - } - reg->lowmant = (reg->mant.m32[1] << 7) | (reg->lowmant ? 1 : 0); - reg->mant.m32[1] = (reg->mant.m32[1] >> 1) | - (reg->mant.m32[0] << 31); - reg->mant.m32[0] = (reg->mant.m32[0] >> 1) | 0x80000000; - - return 1; -} - -static inline void fp_submant(struct fp_ext *dest, struct fp_ext *src1, - struct fp_ext *src2) -{ - /* we assume here, gcc only insert move and a clr instr */ - asm volatile ("sub.b %1,%0" : "=d,g" (dest->lowmant) - : "g,d" (src2->lowmant), "0,0" (src1->lowmant)); - asm volatile ("subx.l %1,%0" : "=d" (dest->mant.m32[1]) - : "d" (src2->mant.m32[1]), "0" (src1->mant.m32[1])); - asm volatile ("subx.l %1,%0" : "=d" (dest->mant.m32[0]) - : "d" (src2->mant.m32[0]), "0" (src1->mant.m32[0])); -} - -#define fp_mul64(desth, destl, src1, src2) ({ \ - asm ("mulu.l %2,%1:%0" : "=d" (destl), "=d" (desth) \ - : "dm" (src1), "0" (src2)); \ -}) -#define fp_div64(quot, rem, srch, srcl, div) \ - asm ("divu.l %2,%1:%0" : "=d" (quot), "=d" (rem) \ - : "dm" (div), "1" (srch), "0" (srcl)) -#define fp_add64(dest1, dest2, src1, src2) ({ \ - asm ("add.l %1,%0" : "=d,dm" (dest2) \ - : "dm,d" (src2), "0,0" (dest2)); \ - asm ("addx.l %1,%0" : "=d" (dest1) \ - : "d" (src1), "0" (dest1)); \ -}) -#define fp_addx96(dest, src) ({ \ - /* we assume here, gcc only insert move and a clr instr */ \ - asm volatile ("add.l %1,%0" : "=d,g" (dest->m32[2]) \ - : "g,d" (temp.m32[1]), "0,0" (dest->m32[2])); \ - asm volatile ("addx.l %1,%0" : "=d" (dest->m32[1]) \ - : "d" (temp.m32[0]), "0" (dest->m32[1])); \ - asm volatile ("addx.l %1,%0" : "=d" (dest->m32[0]) \ - : "d" (0), "0" (dest->m32[0])); \ -}) -#define fp_sub64(dest, src) ({ \ - asm ("sub.l %1,%0" : "=d,dm" (dest.m32[1]) \ - : "dm,d" (src.m32[1]), "0,0" (dest.m32[1])); \ - asm ("subx.l %1,%0" : "=d" (dest.m32[0]) \ - : "d" (src.m32[0]), "0" (dest.m32[0])); \ -}) -#define fp_sub96c(dest, srch, srcm, srcl) ({ \ - char carry; \ - asm ("sub.l %1,%0" : "=d,dm" (dest.m32[2]) \ - : "dm,d" (srcl), "0,0" (dest.m32[2])); \ - asm ("subx.l %1,%0" : "=d" (dest.m32[1]) \ - : "d" (srcm), "0" (dest.m32[1])); \ - asm ("subx.l %2,%1; scs %0" : "=d" (carry), "=d" (dest.m32[0]) \ - : "d" (srch), "1" (dest.m32[0])); \ - carry; \ -}) - -static inline void fp_multiplymant(union fp_mant128 *dest, struct fp_ext *src1, - struct fp_ext *src2) -{ - union fp_mant64 temp; - - fp_mul64(dest->m32[0], dest->m32[1], src1->mant.m32[0], src2->mant.m32[0]); - fp_mul64(dest->m32[2], dest->m32[3], src1->mant.m32[1], src2->mant.m32[1]); - - fp_mul64(temp.m32[0], temp.m32[1], src1->mant.m32[0], src2->mant.m32[1]); - fp_addx96(dest, temp); - - fp_mul64(temp.m32[0], temp.m32[1], src1->mant.m32[1], src2->mant.m32[0]); - fp_addx96(dest, temp); -} - -static inline void fp_dividemant(union fp_mant128 *dest, struct fp_ext *src, - struct fp_ext *div) -{ - union fp_mant128 tmp; - union fp_mant64 tmp64; - unsigned long *mantp = dest->m32; - unsigned long fix, rem, first, dummy; - int i; - - /* the algorithm below requires dest to be smaller than div, - but both have the high bit set */ - if (src->mant.m64 >= div->mant.m64) { - fp_sub64(src->mant, div->mant); - *mantp = 1; - } else - *mantp = 0; - mantp++; - - /* basic idea behind this algorithm: we can't divide two 64bit numbers - (AB/CD) directly, but we can calculate AB/C0, but this means this - quotient is off by C0/CD, so we have to multiply the first result - to fix the result, after that we have nearly the correct result - and only a few corrections are needed. */ - - /* C0/CD can be precalculated, but it's an 64bit division again, but - we can make it a bit easier, by dividing first through C so we get - 10/1D and now only a single shift and the value fits into 32bit. */ - fix = 0x80000000; - dummy = div->mant.m32[1] / div->mant.m32[0] + 1; - dummy = (dummy >> 1) | fix; - fp_div64(fix, dummy, fix, 0, dummy); - fix--; - - for (i = 0; i < 3; i++, mantp++) { - if (src->mant.m32[0] == div->mant.m32[0]) { - fp_div64(first, rem, 0, src->mant.m32[1], div->mant.m32[0]); - - fp_mul64(*mantp, dummy, first, fix); - *mantp += fix; - } else { - fp_div64(first, rem, src->mant.m32[0], src->mant.m32[1], div->mant.m32[0]); - - fp_mul64(*mantp, dummy, first, fix); - } - - fp_mul64(tmp.m32[0], tmp.m32[1], div->mant.m32[0], first - *mantp); - fp_add64(tmp.m32[0], tmp.m32[1], 0, rem); - tmp.m32[2] = 0; - - fp_mul64(tmp64.m32[0], tmp64.m32[1], *mantp, div->mant.m32[1]); - fp_sub96c(tmp, 0, tmp64.m32[0], tmp64.m32[1]); - - src->mant.m32[0] = tmp.m32[1]; - src->mant.m32[1] = tmp.m32[2]; - - while (!fp_sub96c(tmp, 0, div->mant.m32[0], div->mant.m32[1])) { - src->mant.m32[0] = tmp.m32[1]; - src->mant.m32[1] = tmp.m32[2]; - *mantp += 1; - } - } -} - -static inline void fp_putmant128(struct fp_ext *dest, union fp_mant128 *src, - int shift) -{ - unsigned long tmp; - - switch (shift) { - case 0: - dest->mant.m64 = src->m64[0]; - dest->lowmant = src->m32[2] >> 24; - if (src->m32[3] || (src->m32[2] << 8)) - dest->lowmant |= 1; - break; - case 1: - asm volatile ("lsl.l #1,%0" - : "=d" (tmp) : "0" (src->m32[2])); - asm volatile ("roxl.l #1,%0" - : "=d" (dest->mant.m32[1]) : "0" (src->m32[1])); - asm volatile ("roxl.l #1,%0" - : "=d" (dest->mant.m32[0]) : "0" (src->m32[0])); - dest->lowmant = tmp >> 24; - if (src->m32[3] || (tmp << 8)) - dest->lowmant |= 1; - break; - case 31: - asm volatile ("lsr.l #1,%1; roxr.l #1,%0" - : "=d" (dest->mant.m32[0]) - : "d" (src->m32[0]), "0" (src->m32[1])); - asm volatile ("roxr.l #1,%0" - : "=d" (dest->mant.m32[1]) : "0" (src->m32[2])); - asm volatile ("roxr.l #1,%0" - : "=d" (tmp) : "0" (src->m32[3])); - dest->lowmant = tmp >> 24; - if (src->m32[3] << 7) - dest->lowmant |= 1; - break; - case 32: - dest->mant.m32[0] = src->m32[1]; - dest->mant.m32[1] = src->m32[2]; - dest->lowmant = src->m32[3] >> 24; - if (src->m32[3] << 8) - dest->lowmant |= 1; - break; - } -} - -#endif /* MULTI_ARITH_H */ |