Implement fcosh, fsinh, ftanh, fatanh, facos, fasin,
fatan, fsincos, fcos, fsin, ftan.
This is the last series to implement FPU instructions for 680x0...
As previously, all the floatx80 functions are copied from "Previous",
the NeXT Computer Emulator, and written by Andreas Grabher.
I did not improve or clean up the code, it's a simple port of
Andreas' work from "Previous".
I think this series will need some improvemants in the future:
for intance sin and cos functions can be merged in sincos
function (and we have unmodified "adjn" variable in floatx80_cos()
and floatx80_sin() that can be removed otherwise).
"Previous" code I have ported to QEMU can be found in (r844):
http://svn.code.sf.net/p/previous/code/trunk/src/softfloat/
"Previous" code is a C implementation of m68040 assembly language functions
found in ($NetBSD: copyright.s,v 1.2 1994/10/26 07:48:57 cgd Exp)
https://ftp.netbsd.org/pub/NetBSD/NetBSD-current/src/sys/arch/m68k/fpsp/
All the original work is under SoftFloat-2a license, and additional work
under BSD license or GPL-v2-or-later license.
All softfloat new functions are added in target/m68k as they are derived
from m68k code.
I have compared results of these instructions from a real m68040 and from
QEMU, and they match (sincos differs [1] because in QEMU we compute it as
sin and cos, and on m68040 sin and cos results differ also with
sincos results. It looks like a rounding problem, perhaps a bug in m68040
FPU? or in my test program?)
I know this will need more work, but for the moment I only would
like to provide these new instructions to the maintream of QEMU.
[1] Example:
sincos -6.125952 -> (0.156586, 0.987664)
QEMU:
sincos c0010000c407cd1182c2bb27 -> 3ffc0000a05812beea449a4d
3ffe0000fcd791d65887d19a
sin c0010000c407cd1182c2bb27 -> 3ffc0000a05812beea449a4d
cos c0010000c407cd1182c2bb27 -> 3ffe0000fcd791d65887d19a
m68040:
sincos c0010000c407cd1182c2bb27 -> 3ffc0000a05812beea449a4d
3ffe0000fcd791d65887d199 <<--
sin c0010000c407cd1182c2bb27 -> 3ffc0000a05812beea449a4d
cos c0010000c407cd1182c2bb27 -> 3ffe0000fcd791d65887d19a
test program is basically:
register long double a, sin, cos, pi;
long double step;
asm("fmovecr #0, %0" : "=f" (pi));
step = pi / 1024;
for (a = -2*pi; a <= 2*pi; a+= step) {
asm("fsincos.x %2, %0, %1" : "=f" (cos), "=f" (sin) : "f" (a));
print_result(a, sin, cos);
asm("fsin.x %1, %0" : "=f" (sin) : "f" (a));
asm("fcos.x %1, %0" : "=f" (cos) : "f" (a));
print_result(a, sin, cos);
}
Laurent Vivier (11):
target/m68k: implement ftan
target/m68k: implement fsin
target/m68k: implement fcos
target/m68k: implement fsincos
target/m68k: implement fatan
target/m68k: implement fasin
target/m68k: implement facos
target/m68k: implement fatanh
target/m68k: implement ftanh
target/m68k: implement fsinh
target/m68k: implement fcosh
target/m68k/fpu_helper.c | 61 ++
target/m68k/helper.h | 11 +
target/m68k/softfloat.c | 1637 +++++++++++++++++++++++++++++++++++
target/m68k/softfloat.h | 11 +
target/m68k/softfloat_fpsp_tables.h | 267 ++++++
target/m68k/translate.c | 38 +
6 files changed, 2025 insertions(+)
--
2.14.3