aboutsummaryrefslogtreecommitdiffstats
path: root/emulators
diff options
context:
space:
mode:
authornork <nork@FreeBSD.org>2005-05-03 12:02:46 +0800
committernork <nork@FreeBSD.org>2005-05-03 12:02:46 +0800
commit13f88574e0d34999b1baba8acb242ba03e33daf3 (patch)
tree17663f31066c22885a25339a80c9964f47a911f6 /emulators
parent3ebf5b276eb09d35a4915ea69d8166edc3d990c6 (diff)
downloadfreebsd-ports-gnome-13f88574e0d34999b1baba8acb242ba03e33daf3.tar.gz
freebsd-ports-gnome-13f88574e0d34999b1baba8acb242ba03e33daf3.tar.zst
freebsd-ports-gnome-13f88574e0d34999b1baba8acb242ba03e33daf3.zip
o Oops, I forgot to commit 2 patches.
o In some environment, qemu missing 2 patches might be able to compiled. So I bump PORTREVISION. Pointed out by: Jung-uk Kim <jkim@niksun.com> Juergen Lock <nox@jelal.kn-bremen.de> (maintainer) Mark Treacy <mark.treacy@gmail.com> pointyhat via kris
Diffstat (limited to 'emulators')
-rw-r--r--emulators/qemu-devel/Makefile1
-rw-r--r--emulators/qemu-devel/files/patch-fbsd141
-rw-r--r--emulators/qemu-devel/files/patch-libmath2524
-rw-r--r--emulators/qemu/Makefile1
-rw-r--r--emulators/qemu/files/patch-fbsd141
-rw-r--r--emulators/qemu/files/patch-libmath2524
6 files changed, 5332 insertions, 0 deletions
diff --git a/emulators/qemu-devel/Makefile b/emulators/qemu-devel/Makefile
index 8d31b2da7e38..6d1a7c1a4c54 100644
--- a/emulators/qemu-devel/Makefile
+++ b/emulators/qemu-devel/Makefile
@@ -7,6 +7,7 @@
PORTNAME= qemu
PORTVERSION= 0.7.0
+POTREVISION= 1
CATEGORIES= emulators
MASTER_SITES= http://www.qemu.org/
diff --git a/emulators/qemu-devel/files/patch-fbsd b/emulators/qemu-devel/files/patch-fbsd
new file mode 100644
index 000000000000..dab1f93fc8a0
--- /dev/null
+++ b/emulators/qemu-devel/files/patch-fbsd
@@ -0,0 +1,141 @@
+Index: qemu/Makefile
+@@ -13,7 +13,7 @@
+ endif
+ DOCS=qemu-doc.html qemu-tech.html qemu.1 qemu-img.1
+
+-all: dyngen$(EXESUF) $(TOOLS) $(DOCS)
++all: bsd/libmath.a dyngen$(EXESUF) $(TOOLS) $(DOCS)
+ for d in $(TARGET_DIRS); do \
+ $(MAKE) -C $$d $@ || exit 1 ; \
+ done
+@@ -25,6 +25,9 @@
+ endif
+ endif
+
++bsd/libmath.a:
++ ( cd bsd ; $(BSD_MAKE) CC=$(CC) )
++
+ qemu-img$(EXESUF): qemu-img.c block.c block-cow.c block-qcow.c aes.c block-vmdk.c block-cloop.c block-dmg.c block-bochs.c block-vpc.c
+ $(CC) -DQEMU_TOOL $(CFLAGS) $(LDFLAGS) $(DEFINES) -o $@ $^ -lz $(LIBS)
+
+@@ -33,6 +36,7 @@
+
+ clean:
+ # avoid old build problems by removing potentially incorrect old files
++ ( cd bsd ; $(BSD_MAKE) clean )
+ rm -f config.mak config.h op-i386.h opc-i386.h gen-op-i386.h op-arm.h opc-arm.h gen-op-arm.h
+ rm -f *.o *.a $(TOOLS) dyngen$(EXESUF) TAGS *.pod *~ */*~
+ $(MAKE) -C tests clean
+Index: qemu/Makefile.target
+@@ -391,8 +391,8 @@
+ VL_LDFLAGS+=-Wl,-G0 -Wl,-T,$(SRC_PATH)/ia64.ld
+ endif
+
+-$(QEMU_SYSTEM): $(VL_OBJS) libqemu.a
+- $(CC) $(VL_LDFLAGS) -o $@ $^ $(LIBS) $(SDL_LIBS) $(COCOA_LIBS) $(VL_LIBS)
++$(QEMU_SYSTEM): $(VL_OBJS) libqemu.a ../bsd/libmath.a
++ $(CC) $(VL_LDFLAGS) -o $@ $^ $(LIBS) $(SDL_LIBS) $(COCOA_LIBS) $(VL_LIBS) ../bsd/libmath.a
+
+ cocoa.o: cocoa.m
+ $(CC) $(CFLAGS) $(DEFINES) -c -o $@ $<
+Index: qemu/fpu/softfloat-native.c
+@@ -2,11 +2,15 @@
+ context is supported */
+ #include "softfloat.h"
+ #include <math.h>
++#if defined(__FreeBSD__) && __FreeBSD_version < 500000
++#include <ieeefp.h>
++#endif
+
+ void set_float_rounding_mode(int val STATUS_PARAM)
+ {
+ STATUS(float_rounding_mode) = val;
+-#if defined(_BSD) && !defined(__APPLE__)
++#if defined(_BSD) && !defined(__APPLE__) && \
++ (defined(__FreeBSD__) && __FreeBSD_version < 500000)
+ fpsetround(val);
+ #elif defined(__arm__)
+ /* nothing to do */
+@@ -22,7 +26,7 @@
+ }
+ #endif
+
+-#if defined(_BSD)
++#if defined(_BSD) && !defined(__FreeBSD__)
+ #define lrint(d) ((int32_t)rint(d))
+ #define llrint(d) ((int64_t)rint(d))
+ #endif
+Index: qemu/fpu/softfloat-native.h
+@@ -1,7 +1,17 @@
+ /* Native implementation of soft float functions */
+ #include <math.h>
+-#if defined(_BSD) && !defined(__APPLE__)
++#if defined(_BSD) && !defined(__APPLE__) && \
++ (!defined(__FreeBSD__) || __FreeBSD_version < 500000)
+ #include <ieeefp.h>
++#if defined(__FreeBSD__)
++#define isgreater(x, y) __builtin_isgreater((x), (y))
++#define isgreaterequal(x, y) __builtin_isgreaterequal((x), (y))
++#define isless(x, y) __builtin_isless((x), (y))
++#define islessequal(x, y) __builtin_islessequal((x), (y))
++#define islessgreater(x, y) __builtin_islessgreater((x), (y))
++#define isunordered(x, y) __builtin_isunordered((x), (y))
++long double fabsl(long double x);
++#endif
+ #else
+ #include <fenv.h>
+ #endif
+@@ -33,12 +43,13 @@
+ /*----------------------------------------------------------------------------
+ | Software IEC/IEEE floating-point rounding mode.
+ *----------------------------------------------------------------------------*/
+-#if defined(_BSD) && !defined(__APPLE__)
++#if defined(_BSD) && !defined(__APPLE__) && \
++ (!defined(__FreeBSD__) || __FreeBSD_version < 500000)
+ enum {
+ float_round_nearest_even = FP_RN,
+- float_round_down = FE_RM,
+- float_round_up = FE_RP,
+- float_round_to_zero = FE_RZ
++ float_round_down = FP_RM,
++ float_round_up = FP_RP,
++ float_round_to_zero = FP_RZ
+ };
+ #elif defined(__arm__)
+ enum {
+Index: qemu/fpu/softfloat.h
+@@ -84,7 +84,8 @@
+ #define FLOAT128
+ #else
+ /* native float support */
+-#if (defined(__i386__) || defined(__x86_64__)) && !defined(_BSD)
++#if (defined(__i386__) || defined(__x86_64__)) && \
++ (!defined(_BSD) || defined(__FreeBSD__))
+ #define FLOATX80
+ #endif
+ #endif /* !CONFIG_SOFTFLOAT */
+Index: qemu/target-ppc/op_helper.c
+@@ -303,6 +303,13 @@
+ FT0 = sqrt(FT0);
+ }
+
++#ifndef isnormal
++#define isnormal(x) \
++ ((sizeof (x) == sizeof (float)) ? __isnormalf(x) \
++ : (sizeof (x) == sizeof (double)) ? __isnormal(x) \
++ : __isnormall(x))
++#endif
++
+ void do_fres (void)
+ {
+ union {
+Index: qemu/x86_64.ld
+@@ -2,7 +2,7 @@
+ OUTPUT_FORMAT("elf64-x86-64", "elf64-x86-64", "elf64-x86-64")
+ OUTPUT_ARCH(i386:x86-64)
+ ENTRY(_start)
+-SEARCH_DIR("/lib64"); SEARCH_DIR("/usr/lib64"); SEARCH_DIR("/usr/local/lib64");
++SEARCH_DIR("/lib"); SEARCH_DIR("/usr/lib"); SEARCH_DIR("/usr/local/lib");
+ SECTIONS
+ {
+ /* Read-only sections, merged into text segment: */
diff --git a/emulators/qemu-devel/files/patch-libmath b/emulators/qemu-devel/files/patch-libmath
new file mode 100644
index 000000000000..d0c4cbaa363c
--- /dev/null
+++ b/emulators/qemu-devel/files/patch-libmath
@@ -0,0 +1,2524 @@
+diff -Nru qemu-0.7.0/bsd.orig/Makefile qemu-0.7.0/bsd/Makefile
+--- qemu-0.7.0/bsd.orig/Makefile Wed Dec 31 19:00:00 1969
++++ qemu-0.7.0/bsd/Makefile Fri Apr 29 02:11:27 2005
+@@ -0,0 +1,32 @@
++SRCS= ${MACHINE_ARCH}/e_atan2l.c \
++ ${MACHINE_ARCH}/e_logl.S \
++ ${MACHINE_ARCH}/e_powl.S \
++ ${MACHINE_ARCH}/e_remainderl.S \
++ ${MACHINE_ARCH}/e_sqrtl.c \
++ ${MACHINE_ARCH}/s_ceill.S \
++ ${MACHINE_ARCH}/s_cosl.S \
++ ${MACHINE_ARCH}/s_floorl.S \
++ ${MACHINE_ARCH}/s_isnormal.c \
++ ${MACHINE_ARCH}/s_llrint.S \
++ ${MACHINE_ARCH}/s_llrintf.S \
++ ${MACHINE_ARCH}/s_llrintl.S \
++ ${MACHINE_ARCH}/s_lrint.S \
++ ${MACHINE_ARCH}/s_lrintf.S \
++ ${MACHINE_ARCH}/s_lrintl.S \
++ ${MACHINE_ARCH}/s_rintl.c \
++ ${MACHINE_ARCH}/s_round.c \
++ ${MACHINE_ARCH}/s_sinl.S \
++ ${MACHINE_ARCH}/s_tanl.S
++
++OBJS= ${SRCS:R:S/$/.o/}
++
++CFLAGS= -O2 -Wall -I.
++
++all: libmath.a
++
++libmath.a: ${OBJS}
++ rm -f $@
++ ${AR} rcs $@ ${OBJS:T}
++
++clean:
++ rm -f ${OBJS:T} libmath.a
+diff -Nru qemu-0.7.0/bsd.orig/amd64/e_atan2l.c qemu-0.7.0/bsd/amd64/e_atan2l.c
+--- qemu-0.7.0/bsd.orig/amd64/e_atan2l.c Wed Dec 31 19:00:00 1969
++++ qemu-0.7.0/bsd/amd64/e_atan2l.c Fri Apr 29 02:11:27 2005
+@@ -0,0 +1,20 @@
++/*
++ * Written by J.T. Conklin <jtc@netbsd.org>.
++ * Public domain.
++ *
++ * Adapted for `long double' by Ulrich Drepper <drepper@cygnus.com>.
++ */
++
++#include <sysdep.h>
++
++long double
++__ieee754_atan2l (long double y, long double x)
++{
++ long double res;
++
++ asm ("fpatan" : "=t" (res) : "u" (y), "0" (x) : "st(1)");
++
++ return res;
++}
++
++weak_alias(__ieee754_atan2l, atan2l)
+diff -Nru qemu-0.7.0/bsd.orig/amd64/e_logl.S qemu-0.7.0/bsd/amd64/e_logl.S
+--- qemu-0.7.0/bsd.orig/amd64/e_logl.S Wed Dec 31 19:00:00 1969
++++ qemu-0.7.0/bsd/amd64/e_logl.S Fri Apr 29 02:11:27 2005
+@@ -0,0 +1,59 @@
++/*
++ * Written by J.T. Conklin <jtc@netbsd.org>.
++ * Public domain.
++ *
++ * Adapted for `long double' by Ulrich Drepper <drepper@cygnus.com>.
++ * Adapted for x86-64 by Andreas Jaeger <aj@suse.de>.
++ */
++
++#include <machine/asm.h>
++#include <sysdep.h>
++
++RCSID("$NetBSD: $")
++
++
++#ifdef __ELF__
++ .section .rodata
++#else
++ .text
++#endif
++ .align ALIGNARG(4)
++ ASM_TYPE_DIRECTIVE(one,@object)
++one: .double 1.0
++ ASM_SIZE_DIRECTIVE(one)
++ /* It is not important that this constant is precise. It is only
++ a value which is known to be on the safe side for using the
++ fyl2xp1 instruction. */
++ ASM_TYPE_DIRECTIVE(limit,@object)
++limit: .double 0.29
++ ASM_SIZE_DIRECTIVE(limit)
++
++
++#ifdef PIC
++#define MO(op) op##(%rip)
++#else
++#define MO(op) op
++#endif
++
++ .text
++ENTRY(__ieee754_logl)
++ fldln2 // log(2)
++ fldt 8(%rsp) // x : log(2)
++ fld %st // x : x : log(2)
++ fsubl MO(one) // x-1 : x : log(2)
++ fld %st // x-1 : x-1 : x : log(2)
++ fabs // |x-1| : x-1 : x : log(2)
++ fcompl MO(limit) // x-1 : x : log(2)
++ fnstsw // x-1 : x : log(2)
++ andb $0x45, %ah
++ jz 2f
++ fstp %st(1) // x-1 : log(2)
++ fyl2xp1 // log(x)
++ ret
++
++2: fstp %st(0) // x : log(2)
++ fyl2x // log(x)
++ ret
++END (__ieee754_logl)
++
++weak_alias(__ieee754_logl,logl)
+diff -Nru qemu-0.7.0/bsd.orig/amd64/e_powl.S qemu-0.7.0/bsd/amd64/e_powl.S
+--- qemu-0.7.0/bsd.orig/amd64/e_powl.S Wed Dec 31 19:00:00 1969
++++ qemu-0.7.0/bsd/amd64/e_powl.S Fri Apr 29 02:11:27 2005
+@@ -0,0 +1,341 @@
++/* ix87 specific implementation of pow function.
++ Copyright (C) 1996, 1997, 1998, 1999, 2001, 2004 Free Software Foundation, Inc.
++ This file is part of the GNU C Library.
++ Contributed by Ulrich Drepper <drepper@cygnus.com>, 1996.
++
++ The GNU C Library is free software; you can redistribute it and/or
++ modify it under the terms of the GNU Lesser General Public
++ License as published by the Free Software Foundation; either
++ version 2.1 of the License, or (at your option) any later version.
++
++ The GNU C Library is distributed in the hope that it will be useful,
++ but WITHOUT ANY WARRANTY; without even the implied warranty of
++ MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
++ Lesser General Public License for more details.
++
++ You should have received a copy of the GNU Lesser General Public
++ License along with the GNU C Library; if not, write to the Free
++ Software Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA
++ 02111-1307 USA. */
++
++#include <machine/asm.h>
++#include <sysdep.h>
++
++#ifdef __ELF__
++ .section .rodata
++#else
++ .text
++#endif
++
++ .align ALIGNARG(4)
++ ASM_TYPE_DIRECTIVE(infinity,@object)
++inf_zero:
++infinity:
++ .byte 0, 0, 0, 0, 0, 0, 0xf0, 0x7f
++ ASM_SIZE_DIRECTIVE(infinity)
++ ASM_TYPE_DIRECTIVE(zero,@object)
++zero: .double 0.0
++ ASM_SIZE_DIRECTIVE(zero)
++ ASM_TYPE_DIRECTIVE(minf_mzero,@object)
++minf_mzero:
++minfinity:
++ .byte 0, 0, 0, 0, 0, 0, 0xf0, 0xff
++mzero:
++ .byte 0, 0, 0, 0, 0, 0, 0, 0x80
++ ASM_SIZE_DIRECTIVE(minf_mzero)
++ ASM_TYPE_DIRECTIVE(one,@object)
++one: .double 1.0
++ ASM_SIZE_DIRECTIVE(one)
++ ASM_TYPE_DIRECTIVE(limit,@object)
++limit: .double 0.29
++ ASM_SIZE_DIRECTIVE(limit)
++ ASM_TYPE_DIRECTIVE(p63,@object)
++p63:
++ .byte 0, 0, 0, 0, 0, 0, 0xe0, 0x43
++ ASM_SIZE_DIRECTIVE(p63)
++
++#ifdef PIC
++#define MO(op) op##(%rip)
++#else
++#define MO(op) op
++#endif
++
++ .text
++ENTRY(__ieee754_powl)
++ fldt 24(%rsp) // y
++ fxam
++
++
++ fnstsw
++ movb %ah, %dl
++ andb $0x45, %ah
++ cmpb $0x40, %ah // is y == 0 ?
++ je 11f
++
++ cmpb $0x05, %ah // is y == ±inf ?
++ je 12f
++
++ cmpb $0x01, %ah // is y == NaN ?
++ je 30f
++
++ fldt 8(%rsp) // x : y
++
++ fxam
++ fnstsw
++ movb %ah, %dh
++ andb $0x45, %ah
++ cmpb $0x40, %ah
++ je 20f // x is ±0
++
++ cmpb $0x05, %ah
++ je 15f // x is ±inf
++
++ fxch // y : x
++
++ /* fistpll raises invalid exception for |y| >= 1L<<63. */
++ fldl MO(p63) // 1L<<63 : y : x
++ fld %st(1) // y : 1L<<63 : y : x
++ fabs // |y| : 1L<<63 : y : x
++ fcomip %st(1), %st // 1L<<63 : y : x
++ fstp %st(0) // y : x
++ jnc 2f
++
++ /* First see whether `y' is a natural number. In this case we
++ can use a more precise algorithm. */
++ fld %st // y : y : x
++ fistpll -8(%rsp) // y : x
++ fildll -8(%rsp) // int(y) : y : x
++ fucomip %st(1),%st // y : x
++ jne 2f
++
++ /* OK, we have an integer value for y. */
++ mov -8(%rsp),%eax
++ mov -4(%rsp),%edx
++ orl $0, %edx
++ fstp %st(0) // x
++ jns 4f // y >= 0, jump
++ fdivrl MO(one) // 1/x (now referred to as x)
++ negl %eax
++ adcl $0, %edx
++ negl %edx
++4: fldl MO(one) // 1 : x
++ fxch
++
++6: shrdl $1, %edx, %eax
++ jnc 5f
++ fxch
++ fmul %st(1) // x : ST*x
++ fxch
++5: fmul %st(0), %st // x*x : ST*x
++ shrl $1, %edx
++ movl %eax, %ecx
++ orl %edx, %ecx
++ jnz 6b
++ fstp %st(0) // ST*x
++ ret
++
++ /* y is ±NAN */
++30: fldt 8(%rsp) // x : y
++ fldl MO(one) // 1.0 : x : y
++ fucomip %st(1),%st // x : y
++ je 31f
++ fxch // y : x
++31: fstp %st(1)
++ ret
++
++ .align ALIGNARG(4)
++2: /* y is a real number. */
++ fxch // x : y
++ fldl MO(one) // 1.0 : x : y
++ fld %st(1) // x : 1.0 : x : y
++ fsub %st(1) // x-1 : 1.0 : x : y
++ fabs // |x-1| : 1.0 : x : y
++ fcompl MO(limit) // 1.0 : x : y
++ fnstsw
++ fxch // x : 1.0 : y
++ test $4500,%eax
++ jz 7f
++ fsub %st(1) // x-1 : 1.0 : y
++ fyl2xp1 // log2(x) : y
++ jmp 8f
++
++7: fyl2x // log2(x) : y
++8: fmul %st(1) // y*log2(x) : y
++ fxam
++ fnstsw
++ andb $0x45, %ah
++ cmpb $0x05, %ah // is y*log2(x) == ±inf ?
++ je 28f
++ fst %st(1) // y*log2(x) : y*log2(x)
++ frndint // int(y*log2(x)) : y*log2(x)
++ fsubr %st, %st(1) // int(y*log2(x)) : fract(y*log2(x))
++ fxch // fract(y*log2(x)) : int(y*log2(x))
++ f2xm1 // 2^fract(y*log2(x))-1 : int(y*log2(x))
++ faddl MO(one) // 2^fract(y*log2(x)) : int(y*log2(x))
++ fscale // 2^fract(y*log2(x))*2^int(y*log2(x)) : int(y*log2(x))
++ fstp %st(1) // 2^fract(y*log2(x))*2^int(y*log2(x))
++ ret
++
++28: fstp %st(1) // y*log2(x)
++ fldl MO(one) // 1 : y*log2(x)
++ fscale // 2^(y*log2(x)) : y*log2(x)
++ fstp %st(1) // 2^(y*log2(x))
++ ret
++
++ // pow(x,±0) = 1
++ .align ALIGNARG(4)
++11: fstp %st(0) // pop y
++ fldl MO(one)
++ ret
++
++ // y == ±inf
++ .align ALIGNARG(4)
++12: fstp %st(0) // pop y
++ fldt 8(%rsp) // x
++ fabs
++ fcompl MO(one) // < 1, == 1, or > 1
++ fnstsw
++ andb $0x45, %ah
++ cmpb $0x45, %ah
++ je 13f // jump if x is NaN
++
++ cmpb $0x40, %ah
++ je 14f // jump if |x| == 1
++
++ shlb $1, %ah
++ xorb %ah, %dl
++ andl $2, %edx
++#ifdef PIC
++ lea inf_zero(%rip),%rcx
++ fldl (%rcx, %rdx, 4)
++#else
++ fldl inf_zero(,%rdx, 4)
++#endif
++ ret
++
++ .align ALIGNARG(4)
++14: fldl MO(one)
++ ret
++
++ .align ALIGNARG(4)
++13: fldt 8(%rsp) // load x == NaN
++ ret
++
++ .align ALIGNARG(4)
++ // x is ±inf
++15: fstp %st(0) // y
++ testb $2, %dh
++ jz 16f // jump if x == +inf
++
++ // We must find out whether y is an odd integer.
++ fld %st // y : y
++ fistpll -8(%rsp) // y
++ fildll -8(%rsp) // int(y) : y
++ fucomip %st(1),%st
++ ffreep %st // <empty>
++ jne 17f
++
++ // OK, the value is an integer, but is it odd?
++ mov -8(%rsp), %eax
++ mov -4(%rsp), %edx
++ andb $1, %al
++ jz 18f // jump if not odd
++ // It's an odd integer.
++ shrl $31, %edx
++#ifdef PIC
++ lea minf_mzero(%rip),%rcx
++ fldl (%rcx, %rdx, 8)
++#else
++ fldl minf_mzero(,%rdx, 8)
++#endif
++ ret
++
++ .align ALIGNARG(4)
++16: fcompl MO(zero)
++ fnstsw
++ shrl $5, %eax
++ andl $8, %eax
++#ifdef PIC
++ lea inf_zero(%rip),%rcx
++ fldl (%rcx, %rax, 1)
++#else
++ fldl inf_zero(,%rax, 1)
++#endif
++ ret
++
++ .align ALIGNARG(4)
++17: shll $30, %edx // sign bit for y in right position
++18: shrl $31, %edx
++#ifdef PIC
++ lea inf_zero(%rip),%rcx
++ fldl (%rcx, %rdx, 8)
++#else
++ fldl inf_zero(,%rdx, 8)
++#endif
++ ret
++
++ .align ALIGNARG(4)
++ // x is ±0
++20: fstp %st(0) // y
++ testb $2, %dl
++ jz 21f // y > 0
++
++ // x is ±0 and y is < 0. We must find out whether y is an odd integer.
++ testb $2, %dh
++ jz 25f
++
++ fld %st // y : y
++ fistpll -8(%rsp) // y
++ fildll -8(%rsp) // int(y) : y
++ fucomip %st(1),%st
++ ffreep %st // <empty>
++ jne 26f
++
++ // OK, the value is an integer, but is it odd?
++ mov -8(%rsp),%eax
++ mov -4(%rsp),%edx
++ andb $1, %al
++ jz 27f // jump if not odd
++ // It's an odd integer.
++ // Raise divide-by-zero exception and get minus infinity value.
++ fldl MO(one)
++ fdivl MO(zero)
++ fchs
++ ret
++
++25: fstp %st(0)
++26:
++27: // Raise divide-by-zero exception and get infinity value.
++ fldl MO(one)
++ fdivl MO(zero)
++ ret
++
++ .align ALIGNARG(4)
++ // x is ±0 and y is > 0. We must find out whether y is an odd integer.
++21: testb $2, %dh
++ jz 22f
++
++ fld %st // y : y
++ fistpll -8(%rsp) // y
++ fildll -8(%rsp) // int(y) : y
++ fucomip %st(1),%st
++ ffreep %st // <empty>
++ jne 23f
++
++ // OK, the value is an integer, but is it odd?
++ mov -8(%rsp),%eax
++ mov -4(%rsp),%edx
++ andb $1, %al
++ jz 24f // jump if not odd
++ // It's an odd integer.
++ fldl MO(mzero)
++ ret
++
++22: fstp %st(0)
++23:
++24: fldl MO(zero)
++ ret
++
++END(__ieee754_powl)
++
++weak_alias(__ieee754_powl,powl)
+diff -Nru qemu-0.7.0/bsd.orig/amd64/e_remainderl.S qemu-0.7.0/bsd/amd64/e_remainderl.S
+--- qemu-0.7.0/bsd.orig/amd64/e_remainderl.S Wed Dec 31 19:00:00 1969
++++ qemu-0.7.0/bsd/amd64/e_remainderl.S Fri Apr 29 02:11:27 2005
+@@ -0,0 +1,23 @@
++/*
++ * Written by J.T. Conklin <jtc@netbsd.org>.
++ * Public domain.
++ *
++ * Adapted for `long double' by Ulrich Drepper <drepper@cygnus.com>.
++ * Adapted for x86-64 by Andreas Jaeger <aj@suse.de>.
++ */
++
++#include <machine/asm.h>
++#include <sysdep.h>
++
++ENTRY(__ieee754_remainderl)
++ fldt 24(%rsp)
++ fldt 8(%rsp)
++1: fprem1
++ fstsw %ax
++ testl $0x400,%eax
++ jnz 1b
++ fstp %st(1)
++ ret
++END (__ieee754_remainderl)
++
++weak_alias(__ieee754_remainderl,remainderl)
+diff -Nru qemu-0.7.0/bsd.orig/amd64/e_sqrtl.c qemu-0.7.0/bsd/amd64/e_sqrtl.c
+--- qemu-0.7.0/bsd.orig/amd64/e_sqrtl.c Wed Dec 31 19:00:00 1969
++++ qemu-0.7.0/bsd/amd64/e_sqrtl.c Fri Apr 29 02:11:27 2005
+@@ -0,0 +1,20 @@
++/*
++ * Written by J.T. Conklin <jtc@netbsd.org>.
++ * Public domain.
++ *
++ * Adapted for `long double' by Ulrich Drepper <drepper@cygnus.com>.
++ */
++
++#include <sysdep.h>
++
++long double
++__ieee754_sqrtl (long double x)
++{
++ long double res;
++
++ asm ("fsqrt" : "=t" (res) : "0" (x));
++
++ return res;
++}
++
++weak_alias(__ieee754_sqrtl,sqrtl)
+diff -Nru qemu-0.7.0/bsd.orig/amd64/s_ceill.S qemu-0.7.0/bsd/amd64/s_ceill.S
+--- qemu-0.7.0/bsd.orig/amd64/s_ceill.S Wed Dec 31 19:00:00 1969
++++ qemu-0.7.0/bsd/amd64/s_ceill.S Fri Apr 29 02:22:18 2005
+@@ -0,0 +1,246 @@
++/*
++ * ====================================================
++ * Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved.
++ *
++ * Developed at SunPro, a Sun Microsystems, Inc. business.
++ * Permission to use, copy, modify, and distribute this
++ * software is freely granted, provided that this notice
++ * is preserved.
++ * ====================================================
++ *
++ * From: @(#)s_ceil.c 5.1 93/09/24
++ */
++/* XXX: generated from src/lib/msun/src/s_ceill.c */
++
++#include <machine/asm.h>
++
++__FBSDID("$FreeBSD: /tmp/pcvs/ports/emulators/qemu-devel/files/Attic/patch-libmath,v 1.1 2005-05-03 04:02:46 nork Exp $")
++
++ .file "s_ceill.c"
++ .section .rodata.cst8,"aM",@progbits,8
++ .p2align 3
++.LC0:
++ .long 2281731484
++ .long 2117592124
++ .text
++ .p2align 4,,15
++.globl ceill
++ .type ceill, @function
++ceill:
++.LFB17:
++ pushq %rbp
++.LCFI0:
++ pushq %rbx
++.LCFI1:
++ fldt 24(%rsp)
++ movq $0, -16(%rsp)
++ fld %st(0)
++ fstpt -40(%rsp)
++ movl -32(%rsp), %edi
++ movq -40(%rsp), %rsi
++ movl %edi, -16(%rsp)
++ movl -16(%rsp), %r11d
++ movq %rsi, -24(%rsp)
++ movl %r11d, %r10d
++ andl $32767, %r10d
++ leal -16383(%r10), %r8d
++ cmpl $30, %r8d
++ jg .L2
++ testl %r8d, %r8d
++ js .L38
++ movl -20(%rsp), %r9d
++ leal 1(%r8), %ecx
++ mov -24(%rsp), %eax
++ movl $4294967295, %ebp
++ movq %rsi, -40(%rsp)
++ movl %edi, -32(%rsp)
++ shrq %cl, %rbp
++ mov %r9d, %ebx
++ movq %rbx, %rdx
++ andq %rbp, %rdx
++ orq %rax, %rdx
++ fldt -40(%rsp)
++ je .L42
++ ffreep %st(0)
++ testb $-128, -15(%rsp)
++ jne .L12
++ movl $31, %ecx
++ movl $1, %eax
++ subl %r8d, %ecx
++ salq %cl, %rax
++ addl %eax, %r9d
++ mov %r9d, %eax
++ cmpq %rbx, %rax
++ jae .L32
++ leal 1(%r10), %edx
++ movl %r11d, %eax
++ orl $-2147483648, %r9d
++ andw $-32768, %ax
++ andw $32767, %dx
++ orl %edx, %eax
++ movw %ax, -16(%rsp)
++.L32:
++ movl %r9d, -20(%rsp)
++.L12:
++ faddl .LC0(%rip)
++ fldz
++ fxch %st(1)
++ fucomip %st(1), %st
++ fstp %st(0)
++ jbe .L31
++ movl %ebp, %eax
++ movl $0, -24(%rsp)
++ notl %eax
++ andl %eax, %r9d
++ movl %r9d, -20(%rsp)
++ .p2align 4,,7
++.L31:
++ movq -24(%rsp), %rsi
++ movl -16(%rsp), %edi
++ movq %rsi, -40(%rsp)
++ movl %edi, -32(%rsp)
++ fldt -40(%rsp)
++ popq %rbx
++ popq %rbp
++ ret
++ .p2align 4,,7
++.L2:
++ cmpl $62, %r8d
++ jle .L45
++.L44:
++ ffreep %st(0)
++.L17:
++ movq %rsi, -40(%rsp)
++ movl %edi, -32(%rsp)
++ fldt -40(%rsp)
++ popq %rbx
++ popq %rbp
++ ret
++ .p2align 4,,7
++.L45:
++ movl -24(%rsp), %edx
++ leal 1(%r8), %ecx
++ movq $-1, %rbx
++ movq %rsi, -40(%rsp)
++ movl %edi, -32(%rsp)
++ shrq %cl, %rbx
++ mov %edx, %r9d
++ testq %rbx, %r9
++ fldt -40(%rsp)
++ je .L42
++ ffreep %st(0)
++ testb $-128, -15(%rsp)
++ jne .L20
++ cmpl $31, %r8d
++ je .L36
++ movl $63, %ecx
++ movl $1, %eax
++ subl %r8d, %ecx
++ salq %cl, %rax
++ leal (%rdx,%rax), %eax
++ movl %eax, -24(%rsp)
++ mov %eax, %eax
++ cmpq %r9, %rax
++ jae .L20
++.L36:
++ movl -20(%rsp), %eax
++ leal 1(%rax), %ecx
++ cmpl %eax, %ecx
++ jae .L34
++ leal 1(%r10), %edx
++ movl %r11d, %eax
++ orl $-2147483648, %ecx
++ andw $-32768, %ax
++ andw $32767, %dx
++ orl %edx, %eax
++ movw %ax, -16(%rsp)
++.L34:
++ movl %ecx, -20(%rsp)
++ .p2align 4,,7
++.L20:
++ faddl .LC0(%rip)
++ fldz
++ fxch %st(1)
++ fucomip %st(1), %st
++ fstp %st(0)
++ jbe .L31
++ movl %ebx, %eax
++ notl %eax
++ andl %eax, -24(%rsp)
++ jmp .L31
++ .p2align 4,,7
++.L42:
++ fstp %st(1)
++ popq %rbx
++ popq %rbp
++ ret
++ .p2align 4,,7
++.L38:
++ fldl .LC0(%rip)
++ faddp %st, %st(1)
++ fldz
++ fxch %st(1)
++ fucomip %st(1), %st
++ jbe .L44
++ testl %r10d, %r10d
++ jle .L39
++.L7:
++ movabsq $-9223372036854775808, %rsi
++ movl $16383, %edi
++ testb $-128, -15(%rsp)
++ movq %rsi, -40(%rsp)
++ movl %edi, -32(%rsp)
++ fldt -40(%rsp)
++ fcmovne %st(1), %st
++ fstp %st(1)
++ fstpt -40(%rsp)
++ movq -40(%rsp), %rsi
++ movl -32(%rsp), %edi
++ movq %rsi, -24(%rsp)
++ movl %edi, -16(%rsp)
++ jmp .L17
++.L39:
++ movl -24(%rsp), %eax
++ orl -20(%rsp), %eax
++ je .L44
++ jmp .L7
++.LFE17:
++ .size ceill, .-ceill
++ .section .eh_frame,"a",@progbits
++.Lframe1:
++ .long .LECIE1-.LSCIE1
++.LSCIE1:
++ .long 0x0
++ .byte 0x1
++ .string ""
++ .uleb128 0x1
++ .sleb128 -8
++ .byte 0x10
++ .byte 0xc
++ .uleb128 0x7
++ .uleb128 0x8
++ .byte 0x90
++ .uleb128 0x1
++ .p2align 3
++.LECIE1:
++.LSFDE1:
++ .long .LEFDE1-.LASFDE1
++.LASFDE1:
++ .long .LASFDE1-.Lframe1
++ .quad .LFB17
++ .quad .LFE17-.LFB17
++ .byte 0x4
++ .long .LCFI0-.LFB17
++ .byte 0xe
++ .uleb128 0x10
++ .byte 0x4
++ .long .LCFI1-.LCFI0
++ .byte 0xe
++ .uleb128 0x18
++ .byte 0x83
++ .uleb128 0x3
++ .byte 0x86
++ .uleb128 0x2
++ .p2align 3
++.LEFDE1:
++ .ident "GCC: (GNU) 3.4.4 [FreeBSD] 20050421"
+diff -Nru qemu-0.7.0/bsd.orig/amd64/s_cosl.S qemu-0.7.0/bsd/amd64/s_cosl.S
+--- qemu-0.7.0/bsd.orig/amd64/s_cosl.S Wed Dec 31 19:00:00 1969
++++ qemu-0.7.0/bsd/amd64/s_cosl.S Fri Apr 29 02:11:27 2005
+@@ -0,0 +1,33 @@
++/*
++ * Written by J.T. Conklin <jtc@netbsd.org>.
++ * Public domain.
++ *
++ * Adapted for `long double' by Ulrich Drepper <drepper@cygnus.com>.
++ * Adapted for x86-64 by Andreas Jaeger <aj@suse.de>.
++ */
++
++#include <machine/asm.h>
++#include <sysdep.h>
++
++RCSID("$NetBSD: $")
++
++ENTRY(__cosl)
++ fldt 8(%rsp)
++ fcos
++ fnstsw %ax
++ testl $0x400,%eax
++ jnz 1f
++ ret
++ .align ALIGNARG(4)
++1: fldpi
++ fadd %st(0)
++ fxch %st(1)
++2: fprem1
++ fnstsw %ax
++ testl $0x400,%eax
++ jnz 2b
++ fstp %st(1)
++ fcos
++ ret
++END (__cosl)
++weak_alias (__cosl, cosl)
+diff -Nru qemu-0.7.0/bsd.orig/amd64/s_floorl.S qemu-0.7.0/bsd/amd64/s_floorl.S
+--- qemu-0.7.0/bsd.orig/amd64/s_floorl.S Wed Dec 31 19:00:00 1969
++++ qemu-0.7.0/bsd/amd64/s_floorl.S Fri Apr 29 02:24:32 2005
+@@ -0,0 +1,247 @@
++/*
++ * ====================================================
++ * Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved.
++ *
++ * Developed at SunPro, a Sun Microsystems, Inc. business.
++ * Permission to use, copy, modify, and distribute this
++ * software is freely granted, provided that this notice
++ * is preserved.
++ * ====================================================
++ *
++ * From: @(#)s_floor.c 5.1 93/09/24
++ */
++/* XXX: generated from src/lib/msun/src/s_floorl.c */
++
++#include <machine/asm.h>
++
++__FBSDID("$FreeBSD: /tmp/pcvs/ports/emulators/qemu-devel/files/Attic/patch-libmath,v 1.1 2005-05-03 04:02:46 nork Exp $")
++
++ .file "s_floorl.c"
++ .section .rodata.cst8,"aM",@progbits,8
++ .p2align 3
++.LC0:
++ .long 2281731484
++ .long 2117592124
++ .section .rodata.cst4,"aM",@progbits,4
++ .p2align 2
++.LC2:
++ .long 3212836864
++ .text
++ .p2align 4,,15
++.globl floorl
++ .type floorl, @function
++floorl:
++.LFB17:
++ pushq %rbp
++.LCFI0:
++ pushq %rbx
++.LCFI1:
++ fldt 24(%rsp)
++ movq $0, -16(%rsp)
++ fld %st(0)
++ fstpt -40(%rsp)
++ movl -32(%rsp), %edi
++ movq -40(%rsp), %rsi
++ movl %edi, -16(%rsp)
++ movl -16(%rsp), %r11d
++ movq %rsi, -24(%rsp)
++ movl %r11d, %r10d
++ andl $32767, %r10d
++ leal -16383(%r10), %r8d
++ cmpl $30, %r8d
++ jg .L2
++ testl %r8d, %r8d
++ js .L38
++ movl -20(%rsp), %r9d
++ leal 1(%r8), %ecx
++ mov -24(%rsp), %eax
++ movl $4294967295, %ebp
++ movq %rsi, -40(%rsp)
++ movl %edi, -32(%rsp)
++ shrq %cl, %rbp
++ mov %r9d, %ebx
++ movq %rbx, %rdx
++ andq %rbp, %rdx
++ orq %rax, %rdx
++ fldt -40(%rsp)
++ je .L42
++ ffreep %st(0)
++ testb $-128, -15(%rsp)
++ je .L12
++ movl $31, %ecx
++ movl $1, %eax
++ subl %r8d, %ecx
++ salq %cl, %rax
++ addl %eax, %r9d
++ mov %r9d, %eax
++ cmpq %rbx, %rax
++ jae .L32
++ leal 1(%r10), %edx
++ movl %r11d, %eax
++ orl $-2147483648, %r9d
++ andw $-32768, %ax
++ andw $32767, %dx
++ orl %edx, %eax
++ movw %ax, -16(%rsp)
++.L32:
++ movl %r9d, -20(%rsp)
++.L12:
++ faddl .LC0(%rip)
++ fldz
++ fxch %st(1)
++ fucomip %st(1), %st
++ fstp %st(0)
++ jbe .L31
++ movl %ebp, %eax
++ movl $0, -24(%rsp)
++ notl %eax
++ andl %eax, %r9d
++ movl %r9d, -20(%rsp)
++ .p2align 4,,7
++.L31:
++ movq -24(%rsp), %rsi
++ movl -16(%rsp), %edi
++ movq %rsi, -40(%rsp)
++ movl %edi, -32(%rsp)
++ fldt -40(%rsp)
++ popq %rbx
++ popq %rbp
++ ret
++ .p2align 4,,7
++.L2:
++ cmpl $62, %r8d
++ jle .L45
++.L44:
++ ffreep %st(0)
++.L17:
++ movq %rsi, -40(%rsp)
++ movl %edi, -32(%rsp)
++ fldt -40(%rsp)
++ popq %rbx
++ popq %rbp
++ ret
++ .p2align 4,,7
++.L45:
++ movl -24(%rsp), %edx
++ leal 1(%r8), %ecx
++ movq $-1, %rbx
++ movq %rsi, -40(%rsp)
++ movl %edi, -32(%rsp)
++ shrq %cl, %rbx
++ mov %edx, %r9d
++ testq %rbx, %r9
++ fldt -40(%rsp)
++ je .L42
++ ffreep %st(0)
++ testb $-128, -15(%rsp)
++ je .L20
++ cmpl $31, %r8d
++ je .L36
++ movl $63, %ecx
++ movl $1, %eax
++ subl %r8d, %ecx
++ salq %cl, %rax
++ leal (%rdx,%rax), %eax
++ movl %eax, -24(%rsp)
++ mov %eax, %eax
++ cmpq %r9, %rax
++ jae .L20
++.L36:
++ movl -20(%rsp), %eax
++ leal 1(%rax), %ecx
++ cmpl %eax, %ecx
++ jae .L34
++ leal 1(%r10), %edx
++ movl %r11d, %eax
++ orl $-2147483648, %ecx
++ andw $-32768, %ax
++ andw $32767, %dx
++ orl %edx, %eax
++ movw %ax, -16(%rsp)
++.L34:
++ movl %ecx, -20(%rsp)
++ .p2align 4,,7
++.L20:
++ faddl .LC0(%rip)
++ fldz
++ fxch %st(1)
++ fucomip %st(1), %st
++ fstp %st(0)
++ jbe .L31
++ movl %ebx, %eax
++ notl %eax
++ andl %eax, -24(%rsp)
++ jmp .L31
++ .p2align 4,,7
++.L42:
++ fstp %st(1)
++ popq %rbx
++ popq %rbp
++ ret
++ .p2align 4,,7
++.L38:
++ fldl .LC0(%rip)
++ faddp %st, %st(1)
++ fldz
++ fxch %st(1)
++ fucomip %st(1), %st
++ jbe .L44
++ testl %r10d, %r10d
++ jle .L39
++.L7:
++ testb $-128, -15(%rsp)
++ je .L9
++ ffreep %st(0)
++ flds .LC2(%rip)
++.L9:
++ fstpt -40(%rsp)
++ movq -40(%rsp), %rsi
++ movl -32(%rsp), %edi
++ movq %rsi, -24(%rsp)
++ movl %edi, -16(%rsp)
++ jmp .L17
++.L39:
++ movl -24(%rsp), %eax
++ orl -20(%rsp), %eax
++ je .L44
++ jmp .L7
++.LFE17:
++ .size floorl, .-floorl
++ .section .eh_frame,"a",@progbits
++.Lframe1:
++ .long .LECIE1-.LSCIE1
++.LSCIE1:
++ .long 0x0
++ .byte 0x1
++ .string ""
++ .uleb128 0x1
++ .sleb128 -8
++ .byte 0x10
++ .byte 0xc
++ .uleb128 0x7
++ .uleb128 0x8
++ .byte 0x90
++ .uleb128 0x1
++ .p2align 3
++.LECIE1:
++.LSFDE1:
++ .long .LEFDE1-.LASFDE1
++.LASFDE1:
++ .long .LASFDE1-.Lframe1
++ .quad .LFB17
++ .quad .LFE17-.LFB17
++ .byte 0x4
++ .long .LCFI0-.LFB17
++ .byte 0xe
++ .uleb128 0x10
++ .byte 0x4
++ .long .LCFI1-.LCFI0
++ .byte 0xe
++ .uleb128 0x18
++ .byte 0x83
++ .uleb128 0x3
++ .byte 0x86
++ .uleb128 0x2
++ .p2align 3
++.LEFDE1:
++ .ident "GCC: (GNU) 3.4.4 [FreeBSD] 20050421"
+diff -Nru qemu-0.7.0/bsd.orig/amd64/s_isnormal.c qemu-0.7.0/bsd/amd64/s_isnormal.c
+--- qemu-0.7.0/bsd.orig/amd64/s_isnormal.c Wed Dec 31 19:00:00 1969
++++ qemu-0.7.0/bsd/amd64/s_isnormal.c Fri Apr 29 02:11:27 2005
+@@ -0,0 +1,85 @@
++/*-
++ * Copyright (c) 2003 Mike Barcroft <mike@FreeBSD.org>
++ * Copyright (c) 2002-2004 David Schultz <das@FreeBSD.ORG>
++ * 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, this list of conditions and the following disclaimer.
++ * 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.
++ *
++ * THIS SOFTWARE IS PROVIDED BY THE AUTHOR AND CONTRIBUTORS ``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 OR CONTRIBUTORS 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.
++ *
++ * $FreeBSD: /tmp/pcvs/ports/emulators/qemu-devel/files/Attic/patch-libmath,v 1.1 2005-05-03 04:02:46 nork Exp $
++ */
++
++union IEEEf2bits {
++ float f;
++ struct {
++ unsigned int man :23;
++ unsigned int exp :8;
++ unsigned int sign :1;
++ } bits;
++};
++
++union IEEEd2bits {
++ double d;
++ struct {
++ unsigned int manl :32;
++ unsigned int manh :20;
++ unsigned int exp :11;
++ unsigned int sign :1;
++ } bits;
++};
++
++union IEEEl2bits {
++ long double e;
++ struct {
++ unsigned int manl :32;
++ unsigned int manh :32;
++ unsigned int exp :15;
++ unsigned int sign :1;
++ unsigned int junk :16;
++ } bits;
++};
++
++int
++__isnormal(double d)
++{
++ union IEEEd2bits u;
++
++ u.d = d;
++ return (u.bits.exp != 0 && u.bits.exp != 2047);
++}
++
++int
++__isnormalf(float f)
++{
++ union IEEEf2bits u;
++
++ u.f = f;
++ return (u.bits.exp != 0 && u.bits.exp != 255);
++}
++
++int
++__isnormall(long double e)
++{
++ union IEEEl2bits u;
++
++ u.e = e;
++ return (u.bits.exp != 0 && u.bits.exp != 32767);
++}
+diff -Nru qemu-0.7.0/bsd.orig/amd64/s_llrint.S qemu-0.7.0/bsd/amd64/s_llrint.S
+--- qemu-0.7.0/bsd.orig/amd64/s_llrint.S Wed Dec 31 19:00:00 1969
++++ qemu-0.7.0/bsd/amd64/s_llrint.S Fri Apr 29 02:11:27 2005
+@@ -0,0 +1,6 @@
++#include <machine/asm.h>
++__FBSDID("$FreeBSD: /tmp/pcvs/ports/emulators/qemu-devel/files/Attic/patch-libmath,v 1.1 2005-05-03 04:02:46 nork Exp $")
++
++/* sizeof(long) == sizeof(long long) */
++#define fn llrint
++#include "s_lrint.S"
+diff -Nru qemu-0.7.0/bsd.orig/amd64/s_llrintf.S qemu-0.7.0/bsd/amd64/s_llrintf.S
+--- qemu-0.7.0/bsd.orig/amd64/s_llrintf.S Wed Dec 31 19:00:00 1969
++++ qemu-0.7.0/bsd/amd64/s_llrintf.S Fri Apr 29 02:11:27 2005
+@@ -0,0 +1,6 @@
++#include <machine/asm.h>
++__FBSDID("$FreeBSD: /tmp/pcvs/ports/emulators/qemu-devel/files/Attic/patch-libmath,v 1.1 2005-05-03 04:02:46 nork Exp $")
++
++/* sizeof(long) == sizeof(long long) */
++#define fn llrintf
++#include "s_lrintf.S"
+diff -Nru qemu-0.7.0/bsd.orig/amd64/s_llrintl.S qemu-0.7.0/bsd/amd64/s_llrintl.S
+--- qemu-0.7.0/bsd.orig/amd64/s_llrintl.S Wed Dec 31 19:00:00 1969
++++ qemu-0.7.0/bsd/amd64/s_llrintl.S Fri Apr 29 02:11:27 2005
+@@ -0,0 +1,35 @@
++/* Round argument to nearest integral value according to current rounding
++ direction.
++ Copyright (C) 1997, 2002 Free Software Foundation, Inc.
++ This file is part of the GNU C Library.
++
++ The GNU C Library is free software; you can redistribute it and/or
++ modify it under the terms of the GNU Lesser General Public
++ License as published by the Free Software Foundation; either
++ version 2.1 of the License, or (at your option) any later version.
++
++ The GNU C Library is distributed in the hope that it will be useful,
++ but WITHOUT ANY WARRANTY; without even the implied warranty of
++ MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
++ Lesser General Public License for more details.
++
++ You should have received a copy of the GNU Lesser General Public
++ License along with the GNU C Library; if not, write to the Free
++ Software Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA
++ 02111-1307 USA. */
++
++#include <machine/asm.h>
++#include <sysdep.h>
++
++ .text
++ENTRY(__llrintl)
++ fldt 8(%rsp)
++ fistpll -8(%rsp)
++ fwait
++ movq -8(%rsp),%rax
++ ret
++END(__llrintl)
++weak_alias (__llrintl, llrintl)
++strong_alias (__llrintl, __lrintl)
++weak_alias (__llrintl, lrintl)
++
+diff -Nru qemu-0.7.0/bsd.orig/amd64/s_lrint.S qemu-0.7.0/bsd/amd64/s_lrint.S
+--- qemu-0.7.0/bsd.orig/amd64/s_lrint.S Wed Dec 31 19:00:00 1969
++++ qemu-0.7.0/bsd/amd64/s_lrint.S Fri Apr 29 02:11:27 2005
+@@ -0,0 +1,36 @@
++/*-
++ * Copyright (c) 2005 David Schultz <das@FreeBSD.ORG>
++ * 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, this list of conditions and the following disclaimer.
++ * 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.
++ *
++ * THIS SOFTWARE IS PROVIDED BY THE AUTHOR AND CONTRIBUTORS ``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 OR CONTRIBUTORS 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 <machine/asm.h>
++
++#ifndef fn
++__FBSDID("$FreeBSD: /tmp/pcvs/ports/emulators/qemu-devel/files/Attic/patch-libmath,v 1.1 2005-05-03 04:02:46 nork Exp $")
++#define fn lrint
++#endif
++
++ENTRY(fn)
++ cvtsd2si %xmm0, %rax
++ ret
+diff -Nru qemu-0.7.0/bsd.orig/amd64/s_lrintf.S qemu-0.7.0/bsd/amd64/s_lrintf.S
+--- qemu-0.7.0/bsd.orig/amd64/s_lrintf.S Wed Dec 31 19:00:00 1969
++++ qemu-0.7.0/bsd/amd64/s_lrintf.S Fri Apr 29 02:11:27 2005
+@@ -0,0 +1,36 @@
++/*-
++ * Copyright (c) 2005 David Schultz <das@FreeBSD.ORG>
++ * 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, this list of conditions and the following disclaimer.
++ * 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.
++ *
++ * THIS SOFTWARE IS PROVIDED BY THE AUTHOR AND CONTRIBUTORS ``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 OR CONTRIBUTORS 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 <machine/asm.h>
++
++#ifndef fn
++__FBSDID("$FreeBSD: /tmp/pcvs/ports/emulators/qemu-devel/files/Attic/patch-libmath,v 1.1 2005-05-03 04:02:46 nork Exp $")
++#define fn lrintf
++#endif
++
++ENTRY(fn)
++ cvtss2si %xmm0, %rax
++ ret
+diff -Nru qemu-0.7.0/bsd.orig/amd64/s_lrintl.S qemu-0.7.0/bsd/amd64/s_lrintl.S
+--- qemu-0.7.0/bsd.orig/amd64/s_lrintl.S Wed Dec 31 19:00:00 1969
++++ qemu-0.7.0/bsd/amd64/s_lrintl.S Fri Apr 29 02:11:27 2005
+@@ -0,0 +1 @@
++/* Not needed, see s_llrintl.S. */
+diff -Nru qemu-0.7.0/bsd.orig/amd64/s_rintl.c qemu-0.7.0/bsd/amd64/s_rintl.c
+--- qemu-0.7.0/bsd.orig/amd64/s_rintl.c Wed Dec 31 19:00:00 1969
++++ qemu-0.7.0/bsd/amd64/s_rintl.c Fri Apr 29 02:11:27 2005
+@@ -0,0 +1,18 @@
++/*
++ * Written by J.T. Conklin <jtc@netbsd.org>.
++ * Changes for long double by Ulrich Drepper <drepper@cygnus.com>
++ * Public domain.
++ */
++
++#include <sysdep.h>
++
++long double
++__rintl (long double x)
++{
++ long double res;
++
++ asm ("frndint" : "=t" (res) : "0" (x));
++ return res;
++}
++
++weak_alias (__rintl, rintl)
+diff -Nru qemu-0.7.0/bsd.orig/amd64/s_round.c qemu-0.7.0/bsd/amd64/s_round.c
+--- qemu-0.7.0/bsd.orig/amd64/s_round.c Wed Dec 31 19:00:00 1969
++++ qemu-0.7.0/bsd/amd64/s_round.c Fri Apr 29 02:11:27 2005
+@@ -0,0 +1,51 @@
++/*-
++ * Copyright (c) 2003, Steven G. Kargl
++ * 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 unmodified, this list of conditions, and the following
++ * disclaimer.
++ * 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.
++ *
++ * THIS SOFTWARE IS PROVIDED BY THE AUTHOR ``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 <sys/cdefs.h>
++__FBSDID("$FreeBSD: /tmp/pcvs/ports/emulators/qemu-devel/files/Attic/patch-libmath,v 1.1 2005-05-03 04:02:46 nork Exp $");
++
++#include <math.h>
++
++double
++round(double x)
++{
++ double t;
++
++ if (!isfinite(x))
++ return (x);
++
++ if (x >= 0.0) {
++ t = ceil(x);
++ if (t - x > 0.5)
++ t -= 1.0;
++ return (t);
++ } else {
++ t = ceil(-x);
++ if (t + x > 0.5)
++ t -= 1.0;
++ return (-t);
++ }
++}
+diff -Nru qemu-0.7.0/bsd.orig/amd64/s_sinl.S qemu-0.7.0/bsd/amd64/s_sinl.S
+--- qemu-0.7.0/bsd.orig/amd64/s_sinl.S Wed Dec 31 19:00:00 1969
++++ qemu-0.7.0/bsd/amd64/s_sinl.S Fri Apr 29 02:11:27 2005
+@@ -0,0 +1,31 @@
++/*
++ * Written by J.T. Conklin <jtc@netbsd.org>.
++ * Public domain.
++ *
++ * Adapted for `long double' by Ulrich Drepper <drepper@cygnus.com>.
++ * Adapted for x86-64 by Andreas Jaeger <aj@suse.de>.
++ */
++
++#include <machine/asm.h>
++#include <sysdep.h>
++
++ENTRY(__sinl)
++ fldt 8(%rsp)
++ fsin
++ fnstsw %ax
++ testl $0x400,%eax
++ jnz 1f
++ ret
++ .align ALIGNARG(4)
++1: fldpi
++ fadd %st(0)
++ fxch %st(1)
++2: fprem1
++ fnstsw %ax
++ testl $0x400,%eax
++ jnz 2b
++ fstp %st(1)
++ fsin
++ ret
++END (__sinl)
++weak_alias (__sinl, sinl)
+diff -Nru qemu-0.7.0/bsd.orig/amd64/s_tanl.S qemu-0.7.0/bsd/amd64/s_tanl.S
+--- qemu-0.7.0/bsd.orig/amd64/s_tanl.S Wed Dec 31 19:00:00 1969
++++ qemu-0.7.0/bsd/amd64/s_tanl.S Fri Apr 29 02:11:27 2005
+@@ -0,0 +1,34 @@
++/*
++ * Written by J.T. Conklin <jtc@netbsd.org>.
++ * Public domain.
++ *
++ * Adapted for `long double' by Ulrich Drepper <drepper@cygnus.com>.
++ * Adapted for x86-64 by Andreas Jaeger <aj@suse.de>.
++ */
++
++#include <machine/asm.h>
++#include <sysdep.h>
++
++RCSID("$NetBSD: $")
++
++ENTRY(__tanl)
++ fldt 8(%rsp)
++ fptan
++ fnstsw %ax
++ testl $0x400,%eax
++ jnz 1f
++ fstp %st(0)
++ ret
++1: fldpi
++ fadd %st(0)
++ fxch %st(1)
++2: fprem1
++ fstsw %ax
++ testl $0x400,%eax
++ jnz 2b
++ fstp %st(1)
++ fptan
++ fstp %st(0)
++ ret
++END (__tanl)
++weak_alias (__tanl, tanl)
+diff -Nru qemu-0.7.0/bsd.orig/i386/e_atan2l.c qemu-0.7.0/bsd/i386/e_atan2l.c
+--- qemu-0.7.0/bsd.orig/i386/e_atan2l.c Wed Dec 31 19:00:00 1969
++++ qemu-0.7.0/bsd/i386/e_atan2l.c Fri Apr 29 02:11:27 2005
+@@ -0,0 +1,20 @@
++/*
++ * Written by J.T. Conklin <jtc@netbsd.org>.
++ * Public domain.
++ *
++ * Adapted for `long double' by Ulrich Drepper <drepper@cygnus.com>.
++ */
++
++#include <sysdep.h>
++
++long double
++__ieee754_atan2l (long double y, long double x)
++{
++ long double res;
++
++ asm ("fpatan" : "=t" (res) : "u" (y), "0" (x) : "st(1)");
++
++ return res;
++}
++
++weak_alias(__ieee754_atan2l, atan2l)
+diff -Nru qemu-0.7.0/bsd.orig/i386/e_logl.S qemu-0.7.0/bsd/i386/e_logl.S
+--- qemu-0.7.0/bsd.orig/i386/e_logl.S Wed Dec 31 19:00:00 1969
++++ qemu-0.7.0/bsd/i386/e_logl.S Fri Apr 29 02:11:27 2005
+@@ -0,0 +1,63 @@
++/*
++ * Written by J.T. Conklin <jtc@netbsd.org>.
++ * Public domain.
++ *
++ * Adapted for `long double' by Ulrich Drepper <drepper@cygnus.com>.
++ */
++
++#include <machine/asm.h>
++#include <sysdep.h>
++
++RCSID("$NetBSD: $")
++
++
++#ifdef __ELF__
++ .section .rodata
++#else
++ .text
++#endif
++ .align ALIGNARG(4)
++ ASM_TYPE_DIRECTIVE(one,@object)
++one: .double 1.0
++ ASM_SIZE_DIRECTIVE(one)
++ /* It is not important that this constant is precise. It is only
++ a value which is known to be on the safe side for using the
++ fyl2xp1 instruction. */
++ ASM_TYPE_DIRECTIVE(limit,@object)
++limit: .double 0.29
++ ASM_SIZE_DIRECTIVE(limit)
++
++
++#ifdef PIC
++#define MO(op) op##@GOTOFF(%edx)
++#else
++#define MO(op) op
++#endif
++
++ .text
++ENTRY(__ieee754_logl)
++ fldln2 // log(2)
++ fldt 4(%esp) // x : log(2)
++#ifdef PIC
++ call 1f
++1: popl %edx
++ addl $_GLOBAL_OFFSET_TABLE_+[.-1b], %edx
++#endif
++ fld %st // x : x : log(2)
++ fsubl MO(one) // x-1 : x : log(2)
++ fld %st // x-1 : x-1 : x : log(2)
++ fabs // |x-1| : x-1 : x : log(2)
++ fcompl MO(limit) // x-1 : x : log(2)
++ fnstsw // x-1 : x : log(2)
++ andb $0x45, %ah
++ jz 2f
++ fstp %st(1) // x-1 : log(2)
++ fyl2xp1 // log(x)
++ ret
++
++2: fstp %st(0) // x : log(2)
++ fyl2x // log(x)
++ ret
++END (__ieee754_logl)
++
++weak_alias(__ieee754_logl,logl)
+diff -Nru qemu-0.7.0/bsd.orig/i386/e_powl.S qemu-0.7.0/bsd/i386/e_powl.S
+--- qemu-0.7.0/bsd.orig/i386/e_powl.S Wed Dec 31 19:00:00 1969
++++ qemu-0.7.0/bsd/i386/e_powl.S Fri Apr 29 02:11:27 2005
+@@ -0,0 +1,341 @@
++/* ix87 specific implementation of pow function.
++ Copyright (C) 1996, 1997, 1998, 1999, 2001, 2004
++ Free Software Foundation, Inc.
++ This file is part of the GNU C Library.
++ Contributed by Ulrich Drepper <drepper@cygnus.com>, 1996.
++
++ The GNU C Library is free software; you can redistribute it and/or
++ modify it under the terms of the GNU Lesser General Public
++ License as published by the Free Software Foundation; either
++ version 2.1 of the License, or (at your option) any later version.
++
++ The GNU C Library is distributed in the hope that it will be useful,
++ but WITHOUT ANY WARRANTY; without even the implied warranty of
++ MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
++ Lesser General Public License for more details.
++
++ You should have received a copy of the GNU Lesser General Public
++ License along with the GNU C Library; if not, write to the Free
++ Software Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA
++ 02111-1307 USA. */
++
++#include <machine/asm.h>
++#include <sysdep.h>
++
++#ifdef __ELF__
++ .section .rodata
++#else
++ .text
++#endif
++
++ .align ALIGNARG(4)
++ ASM_TYPE_DIRECTIVE(infinity,@object)
++inf_zero:
++infinity:
++ .byte 0, 0, 0, 0, 0, 0, 0xf0, 0x7f
++ ASM_SIZE_DIRECTIVE(infinity)
++ ASM_TYPE_DIRECTIVE(zero,@object)
++zero: .double 0.0
++ ASM_SIZE_DIRECTIVE(zero)
++ ASM_TYPE_DIRECTIVE(minf_mzero,@object)
++minf_mzero:
++minfinity:
++ .byte 0, 0, 0, 0, 0, 0, 0xf0, 0xff
++mzero:
++ .byte 0, 0, 0, 0, 0, 0, 0, 0x80
++ ASM_SIZE_DIRECTIVE(minf_mzero)
++ ASM_TYPE_DIRECTIVE(one,@object)
++one: .double 1.0
++ ASM_SIZE_DIRECTIVE(one)
++ ASM_TYPE_DIRECTIVE(limit,@object)
++limit: .double 0.29
++ ASM_SIZE_DIRECTIVE(limit)
++ ASM_TYPE_DIRECTIVE(p63,@object)
++p63: .byte 0, 0, 0, 0, 0, 0, 0xe0, 0x43
++ ASM_SIZE_DIRECTIVE(p63)
++
++#ifdef PIC
++#define MO(op) op##@GOTOFF(%ecx)
++#define MOX(op,x,f) op##@GOTOFF(%ecx,x,f)
++#else
++#define MO(op) op
++#define MOX(op,x,f) op(,x,f)
++#endif
++
++ .text
++ENTRY(__ieee754_powl)
++ fldt 16(%esp) // y
++ fxam
++
++#ifdef PIC
++ call 1f
++1: popl %ecx
++ addl $_GLOBAL_OFFSET_TABLE_+[.-1b], %ecx
++#endif
++
++ fnstsw
++ movb %ah, %dl
++ andb $0x45, %ah
++ cmpb $0x40, %ah // is y == 0 ?
++ je 11f
++
++ cmpb $0x05, %ah // is y == ±inf ?
++ je 12f
++
++ cmpb $0x01, %ah // is y == NaN ?
++ je 30f
++
++ fldt 4(%esp) // x : y
++
++ subl $8,%esp
++
++ fxam
++ fnstsw
++ movb %ah, %dh
++ andb $0x45, %ah
++ cmpb $0x40, %ah
++ je 20f // x is ±0
++
++ cmpb $0x05, %ah
++ je 15f // x is ±inf
++
++ fxch // y : x
++
++ /* fistpll raises invalid exception for |y| >= 1L<<63. */
++ fld %st // y : y : x
++ fabs // |y| : y : x
++ fcompl MO(p63) // y : x
++ fnstsw
++ sahf
++ jnc 2f
++
++ /* First see whether `y' is a natural number. In this case we
++ can use a more precise algorithm. */
++ fld %st // y : y : x
++ fistpll (%esp) // y : x
++ fildll (%esp) // int(y) : y : x
++ fucomp %st(1) // y : x
++ fnstsw
++ sahf
++ jne 2f
++
++ /* OK, we have an integer value for y. */
++ popl %eax
++ popl %edx
++ orl $0, %edx
++ fstp %st(0) // x
++ jns 4f // y >= 0, jump
++ fdivrl MO(one) // 1/x (now referred to as x)
++ negl %eax
++ adcl $0, %edx
++ negl %edx
++4: fldl MO(one) // 1 : x
++ fxch
++
++6: shrdl $1, %edx, %eax
++ jnc 5f
++ fxch
++ fmul %st(1) // x : ST*x
++ fxch
++5: fmul %st(0), %st // x*x : ST*x
++ shrl $1, %edx
++ movl %eax, %ecx
++ orl %edx, %ecx
++ jnz 6b
++ fstp %st(0) // ST*x
++ ret
++
++ /* y is ±NAN */
++30: fldt 4(%esp) // x : y
++ fldl MO(one) // 1.0 : x : y
++ fucomp %st(1) // x : y
++ fnstsw
++ sahf
++ je 31f
++ fxch // y : x
++31: fstp %st(1)
++ ret
++
++ .align ALIGNARG(4)
++2: /* y is a real number. */
++ fxch // x : y
++ fldl MO(one) // 1.0 : x : y
++ fld %st(1) // x : 1.0 : x : y
++ fsub %st(1) // x-1 : 1.0 : x : y
++ fabs // |x-1| : 1.0 : x : y
++ fcompl MO(limit) // 1.0 : x : y
++ fnstsw
++ fxch // x : 1.0 : y
++ sahf
++ ja 7f
++ fsub %st(1) // x-1 : 1.0 : y
++ fyl2xp1 // log2(x) : y
++ jmp 8f
++
++7: fyl2x // log2(x) : y
++8: fmul %st(1) // y*log2(x) : y
++ fxam
++ fnstsw
++ andb $0x45, %ah
++ cmpb $0x05, %ah // is y*log2(x) == ±inf ?
++ je 28f
++ fst %st(1) // y*log2(x) : y*log2(x)
++ frndint // int(y*log2(x)) : y*log2(x)
++ fsubr %st, %st(1) // int(y*log2(x)) : fract(y*log2(x))
++ fxch // fract(y*log2(x)) : int(y*log2(x))
++ f2xm1 // 2^fract(y*log2(x))-1 : int(y*log2(x))
++ faddl MO(one) // 2^fract(y*log2(x)) : int(y*log2(x))
++ fscale // 2^fract(y*log2(x))*2^int(y*log2(x)) : int(y*log2(x))
++ addl $8, %esp
++ fstp %st(1) // 2^fract(y*log2(x))*2^int(y*log2(x))
++ ret
++
++28: fstp %st(1) // y*log2(x)
++ fldl MO(one) // 1 : y*log2(x)
++ fscale // 2^(y*log2(x)) : y*log2(x)
++ addl $8, %esp
++ fstp %st(1) // 2^(y*log2(x))
++ ret
++
++ // pow(x,±0) = 1
++ .align ALIGNARG(4)
++11: fstp %st(0) // pop y
++ fldl MO(one)
++ ret
++
++ // y == ±inf
++ .align ALIGNARG(4)
++12: fstp %st(0) // pop y
++ fldt 4(%esp) // x
++ fabs
++ fcompl MO(one) // < 1, == 1, or > 1
++ fnstsw
++ andb $0x45, %ah
++ cmpb $0x45, %ah
++ je 13f // jump if x is NaN
++
++ cmpb $0x40, %ah
++ je 14f // jump if |x| == 1
++
++ shlb $1, %ah
++ xorb %ah, %dl
++ andl $2, %edx
++ fldl MOX(inf_zero, %edx, 4)
++ ret
++
++ .align ALIGNARG(4)
++14: fldl MO(one)
++ ret
++
++ .align ALIGNARG(4)
++13: fldt 4(%esp) // load x == NaN
++ ret
++
++ .align ALIGNARG(4)
++ // x is ±inf
++15: fstp %st(0) // y
++ testb $2, %dh
++ jz 16f // jump if x == +inf
++
++ // We must find out whether y is an odd integer.
++ fld %st // y : y
++ fistpll (%esp) // y
++ fildll (%esp) // int(y) : y
++ fucompp // <empty>
++ fnstsw
++ sahf
++ jne 17f
++
++ // OK, the value is an integer, but is it odd?
++ popl %eax
++ popl %edx
++ andb $1, %al
++ jz 18f // jump if not odd
++ // It's an odd integer.
++ shrl $31, %edx
++ fldl MOX(minf_mzero, %edx, 8)
++ ret
++
++ .align ALIGNARG(4)
++16: fcompl MO(zero)
++ addl $8, %esp
++ fnstsw
++ shrl $5, %eax
++ andl $8, %eax
++ fldl MOX(inf_zero, %eax, 1)
++ ret
++
++ .align ALIGNARG(4)
++17: shll $30, %edx // sign bit for y in right position
++ addl $8, %esp
++18: shrl $31, %edx
++ fldl MOX(inf_zero, %edx, 8)
++ ret
++
++ .align ALIGNARG(4)
++ // x is ±0
++20: fstp %st(0) // y
++ testb $2, %dl
++ jz 21f // y > 0
++
++ // x is ±0 and y is < 0. We must find out whether y is an odd integer.
++ testb $2, %dh
++ jz 25f
++
++ fld %st // y : y
++ fistpll (%esp) // y
++ fildll (%esp) // int(y) : y
++ fucompp // <empty>
++ fnstsw
++ sahf
++ jne 26f
++
++ // OK, the value is an integer, but is it odd?
++ popl %eax
++ popl %edx
++ andb $1, %al
++ jz 27f // jump if not odd
++ // It's an odd integer.
++ // Raise divide-by-zero exception and get minus infinity value.
++ fldl MO(one)
++ fdivl MO(zero)
++ fchs
++ ret
++
++25: fstp %st(0)
++26: addl $8, %esp
++27: // Raise divide-by-zero exception and get infinity value.
++ fldl MO(one)
++ fdivl MO(zero)
++ ret
++
++ .align ALIGNARG(4)
++ // x is ±0 and y is > 0. We must find out whether y is an odd integer.
++21: testb $2, %dh
++ jz 22f
++
++ fld %st // y : y
++ fistpll (%esp) // y
++ fildll (%esp) // int(y) : y
++ fucompp // <empty>
++ fnstsw
++ sahf
++ jne 23f
++
++ // OK, the value is an integer, but is it odd?
++ popl %eax
++ popl %edx
++ andb $1, %al
++ jz 24f // jump if not odd
++ // It's an odd integer.
++ fldl MO(mzero)
++ ret
++
++22: fstp %st(0)
++23: addl $8, %esp // Don't use 2 x pop
++24: fldl MO(zero)
++ ret
++
++END(__ieee754_powl)
++
++weak_alias(__ieee754_powl,powl)
+diff -Nru qemu-0.7.0/bsd.orig/i386/e_remainderl.S qemu-0.7.0/bsd/i386/e_remainderl.S
+--- qemu-0.7.0/bsd.orig/i386/e_remainderl.S Wed Dec 31 19:00:00 1969
++++ qemu-0.7.0/bsd/i386/e_remainderl.S Fri Apr 29 02:11:27 2005
+@@ -0,0 +1,24 @@
++/*
++ * Written by J.T. Conklin <jtc@netbsd.org>.
++ * Public domain.
++ *
++ * Adapted for `long double' by Ulrich Drepper <drepper@cygnus.com>.
++ */
++
++#include <machine/asm.h>
++#include <sysdep.h>
++
++RCSID("$NetBSD: $")
++
++ENTRY(__ieee754_remainderl)
++ fldt 16(%esp)
++ fldt 4(%esp)
++1: fprem1
++ fstsw %ax
++ sahf
++ jp 1b
++ fstp %st(1)
++ ret
++END (__ieee754_remainderl)
++
++weak_alias(__ieee754_remainderl,remainderl)
+diff -Nru qemu-0.7.0/bsd.orig/i386/e_sqrtl.c qemu-0.7.0/bsd/i386/e_sqrtl.c
+--- qemu-0.7.0/bsd.orig/i386/e_sqrtl.c Wed Dec 31 19:00:00 1969
++++ qemu-0.7.0/bsd/i386/e_sqrtl.c Fri Apr 29 02:11:27 2005
+@@ -0,0 +1,20 @@
++/*
++ * Written by J.T. Conklin <jtc@netbsd.org>.
++ * Public domain.
++ *
++ * Adapted for `long double' by Ulrich Drepper <drepper@cygnus.com>.
++ */
++
++#include <sysdep.h>
++
++long double
++__ieee754_sqrtl (long double x)
++{
++ long double res;
++
++ asm ("fsqrt" : "=t" (res) : "0" (x));
++
++ return res;
++}
++
++weak_alias(__ieee754_sqrtl,sqrtl)
+diff -Nru qemu-0.7.0/bsd.orig/i386/s_ceill.S qemu-0.7.0/bsd/i386/s_ceill.S
+--- qemu-0.7.0/bsd.orig/i386/s_ceill.S Wed Dec 31 19:00:00 1969
++++ qemu-0.7.0/bsd/i386/s_ceill.S Fri Apr 29 02:11:27 2005
+@@ -0,0 +1,27 @@
++/*
++ * Based on code written by J.T. Conklin <jtc@NetBSD.org>.
++ * Public domain.
++ */
++
++#include <machine/asm.h>
++RCSID("$FreeBSD: /tmp/pcvs/ports/emulators/qemu-devel/files/Attic/patch-libmath,v 1.1 2005-05-03 04:02:46 nork Exp $")
++
++ENTRY(ceill)
++ pushl %ebp
++ movl %esp,%ebp
++ subl $8,%esp
++
++ fstcw -4(%ebp) /* store fpu control word */
++ movw -4(%ebp),%dx
++ orw $0x0800,%dx /* round towards +oo */
++ andw $0xfbff,%dx
++ movw %dx,-8(%ebp)
++ fldcw -8(%ebp) /* load modfied control word */
++
++ fldt 8(%ebp) /* round */
++ frndint
++
++ fldcw -4(%ebp) /* restore original control word */
++
++ leave
++ ret
+diff -Nru qemu-0.7.0/bsd.orig/i386/s_cosl.S qemu-0.7.0/bsd/i386/s_cosl.S
+--- qemu-0.7.0/bsd.orig/i386/s_cosl.S Wed Dec 31 19:00:00 1969
++++ qemu-0.7.0/bsd/i386/s_cosl.S Fri Apr 29 02:11:27 2005
+@@ -0,0 +1,32 @@
++/*
++ * Written by J.T. Conklin <jtc@netbsd.org>.
++ * Public domain.
++ *
++ * Adapted for `long double' by Ulrich Drepper <drepper@cygnus.com>.
++ */
++
++#include <machine/asm.h>
++#include <sysdep.h>
++
++RCSID("$NetBSD: $")
++
++ENTRY(__cosl)
++ fldt 4(%esp)
++ fcos
++ fnstsw %ax
++ testl $0x400,%eax
++ jnz 1f
++ ret
++ .align ALIGNARG(4)
++1: fldpi
++ fadd %st(0)
++ fxch %st(1)
++2: fprem1
++ fnstsw %ax
++ testl $0x400,%eax
++ jnz 2b
++ fstp %st(1)
++ fcos
++ ret
++END (__cosl)
++weak_alias (__cosl, cosl)
+diff -Nru qemu-0.7.0/bsd.orig/i386/s_floorl.S qemu-0.7.0/bsd/i386/s_floorl.S
+--- qemu-0.7.0/bsd.orig/i386/s_floorl.S Wed Dec 31 19:00:00 1969
++++ qemu-0.7.0/bsd/i386/s_floorl.S Fri Apr 29 02:11:27 2005
+@@ -0,0 +1,27 @@
++/*
++ * Based on code written by J.T. Conklin <jtc@NetBSD.org>.
++ * Public domain.
++ */
++
++#include <machine/asm.h>
++RCSID("$FreeBSD: /tmp/pcvs/ports/emulators/qemu-devel/files/Attic/patch-libmath,v 1.1 2005-05-03 04:02:46 nork Exp $")
++
++ENTRY(floorl)
++ pushl %ebp
++ movl %esp,%ebp
++ subl $8,%esp
++
++ fstcw -4(%ebp) /* store fpu control word */
++ movw -4(%ebp),%dx
++ orw $0x0400,%dx /* round towards -oo */
++ andw $0xf7ff,%dx
++ movw %dx,-8(%ebp)
++ fldcw -8(%ebp) /* load modfied control word */
++
++ fldt 8(%ebp) /* round */
++ frndint
++
++ fldcw -4(%ebp) /* restore original control word */
++
++ leave
++ ret
+diff -Nru qemu-0.7.0/bsd.orig/i386/s_isnormal.c qemu-0.7.0/bsd/i386/s_isnormal.c
+--- qemu-0.7.0/bsd.orig/i386/s_isnormal.c Wed Dec 31 19:00:00 1969
++++ qemu-0.7.0/bsd/i386/s_isnormal.c Fri Apr 29 02:11:27 2005
+@@ -0,0 +1,85 @@
++/*-
++ * Copyright (c) 2003 Mike Barcroft <mike@FreeBSD.org>
++ * Copyright (c) 2002-2004 David Schultz <das@FreeBSD.ORG>
++ * 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, this list of conditions and the following disclaimer.
++ * 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.
++ *
++ * THIS SOFTWARE IS PROVIDED BY THE AUTHOR AND CONTRIBUTORS ``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 OR CONTRIBUTORS 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.
++ *
++ * $FreeBSD: /tmp/pcvs/ports/emulators/qemu-devel/files/Attic/patch-libmath,v 1.1 2005-05-03 04:02:46 nork Exp $
++ */
++
++union IEEEf2bits {
++ float f;
++ struct {
++ unsigned int man :23;
++ unsigned int exp :8;
++ unsigned int sign :1;
++ } bits;
++};
++
++union IEEEd2bits {
++ double d;
++ struct {
++ unsigned int manl :32;
++ unsigned int manh :20;
++ unsigned int exp :11;
++ unsigned int sign :1;
++ } bits;
++};
++
++union IEEEl2bits {
++ long double e;
++ struct {
++ unsigned int manl :32;
++ unsigned int manh :32;
++ unsigned int exp :15;
++ unsigned int sign :1;
++ unsigned int junk :16;
++ } bits;
++};
++
++int
++__isnormal(double d)
++{
++ union IEEEd2bits u;
++
++ u.d = d;
++ return (u.bits.exp != 0 && u.bits.exp != 2047);
++}
++
++int
++__isnormalf(float f)
++{
++ union IEEEf2bits u;
++
++ u.f = f;
++ return (u.bits.exp != 0 && u.bits.exp != 255);
++}
++
++int
++__isnormall(long double e)
++{
++ union IEEEl2bits u;
++
++ u.e = e;
++ return (u.bits.exp != 0 && u.bits.exp != 32767);
++}
+diff -Nru qemu-0.7.0/bsd.orig/i386/s_llrint.S qemu-0.7.0/bsd/i386/s_llrint.S
+--- qemu-0.7.0/bsd.orig/i386/s_llrint.S Wed Dec 31 19:00:00 1969
++++ qemu-0.7.0/bsd/i386/s_llrint.S Fri Apr 29 02:11:27 2005
+@@ -0,0 +1,36 @@
++/*-
++ * Copyright (c) 2005 David Schultz <das@FreeBSD.ORG>
++ * 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, this list of conditions and the following disclaimer.
++ * 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.
++ *
++ * THIS SOFTWARE IS PROVIDED BY THE AUTHOR AND CONTRIBUTORS ``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 OR CONTRIBUTORS 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 <machine/asm.h>
++RCSID("$FreeBSD: /tmp/pcvs/ports/emulators/qemu-devel/files/Attic/patch-libmath,v 1.1 2005-05-03 04:02:46 nork Exp $");
++
++ENTRY(llrint)
++ fldl 4(%esp)
++ subl $8,%esp
++ fistpll (%esp)
++ popl %eax
++ popl %edx
++ ret
+diff -Nru qemu-0.7.0/bsd.orig/i386/s_llrintf.S qemu-0.7.0/bsd/i386/s_llrintf.S
+--- qemu-0.7.0/bsd.orig/i386/s_llrintf.S Wed Dec 31 19:00:00 1969
++++ qemu-0.7.0/bsd/i386/s_llrintf.S Fri Apr 29 02:11:27 2005
+@@ -0,0 +1,36 @@
++/*-
++ * Copyright (c) 2005 David Schultz <das@FreeBSD.ORG>
++ * 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, this list of conditions and the following disclaimer.
++ * 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.
++ *
++ * THIS SOFTWARE IS PROVIDED BY THE AUTHOR AND CONTRIBUTORS ``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 OR CONTRIBUTORS 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 <machine/asm.h>
++RCSID("$FreeBSD: /tmp/pcvs/ports/emulators/qemu-devel/files/Attic/patch-libmath,v 1.1 2005-05-03 04:02:46 nork Exp $")
++
++ENTRY(llrintf)
++ flds 4(%esp)
++ subl $8,%esp
++ fistpll (%esp)
++ popl %eax
++ popl %edx
++ ret
+diff -Nru qemu-0.7.0/bsd.orig/i386/s_llrintl.S qemu-0.7.0/bsd/i386/s_llrintl.S
+--- qemu-0.7.0/bsd.orig/i386/s_llrintl.S Wed Dec 31 19:00:00 1969
++++ qemu-0.7.0/bsd/i386/s_llrintl.S Fri Apr 29 02:11:27 2005
+@@ -0,0 +1,35 @@
++/* Round argument to nearest integral value according to current rounding
++ direction.
++ Copyright (C) 1997 Free Software Foundation, Inc.
++ This file is part of the GNU C Library.
++ Contributed by Ulrich Drepper <drepper@cygnus.com>, 1997.
++
++ The GNU C Library is free software; you can redistribute it and/or
++ modify it under the terms of the GNU Lesser General Public
++ License as published by the Free Software Foundation; either
++ version 2.1 of the License, or (at your option) any later version.
++
++ The GNU C Library is distributed in the hope that it will be useful,
++ but WITHOUT ANY WARRANTY; without even the implied warranty of
++ MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
++ Lesser General Public License for more details.
++
++ You should have received a copy of the GNU Lesser General Public
++ License along with the GNU C Library; if not, write to the Free
++ Software Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA
++ 02111-1307 USA. */
++
++#include <machine/asm.h>
++#include <sysdep.h>
++
++ .text
++ENTRY(__llrintl)
++ fldt 4(%esp)
++ subl $8, %esp
++ fistpll (%esp)
++ fwait
++ popl %eax
++ popl %edx
++ ret
++END(__llrintl)
++weak_alias (__llrintl, llrintl)
+diff -Nru qemu-0.7.0/bsd.orig/i386/s_lrint.S qemu-0.7.0/bsd/i386/s_lrint.S
+--- qemu-0.7.0/bsd.orig/i386/s_lrint.S Wed Dec 31 19:00:00 1969
++++ qemu-0.7.0/bsd/i386/s_lrint.S Fri Apr 29 02:11:27 2005
+@@ -0,0 +1,35 @@
++/*-
++ * Copyright (c) 2005 David Schultz <das@FreeBSD.ORG>
++ * 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, this list of conditions and the following disclaimer.
++ * 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.
++ *
++ * THIS SOFTWARE IS PROVIDED BY THE AUTHOR AND CONTRIBUTORS ``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 OR CONTRIBUTORS 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 <machine/asm.h>
++RCSID("$FreeBSD: /tmp/pcvs/ports/emulators/qemu-devel/files/Attic/patch-libmath,v 1.1 2005-05-03 04:02:46 nork Exp $");
++
++ENTRY(lrint)
++ fldl 4(%esp)
++ subl $4,%esp
++ fistpl (%esp)
++ popl %eax
++ ret
+diff -Nru qemu-0.7.0/bsd.orig/i386/s_lrintf.S qemu-0.7.0/bsd/i386/s_lrintf.S
+--- qemu-0.7.0/bsd.orig/i386/s_lrintf.S Wed Dec 31 19:00:00 1969
++++ qemu-0.7.0/bsd/i386/s_lrintf.S Fri Apr 29 02:11:27 2005
+@@ -0,0 +1,35 @@
++/*-
++ * Copyright (c) 2005 David Schultz <das@FreeBSD.ORG>
++ * 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, this list of conditions and the following disclaimer.
++ * 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.
++ *
++ * THIS SOFTWARE IS PROVIDED BY THE AUTHOR AND CONTRIBUTORS ``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 OR CONTRIBUTORS 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 <machine/asm.h>
++RCSID("$FreeBSD: /tmp/pcvs/ports/emulators/qemu-devel/files/Attic/patch-libmath,v 1.1 2005-05-03 04:02:46 nork Exp $")
++
++ENTRY(lrintf)
++ flds 4(%esp)
++ subl $4,%esp
++ fistpl (%esp)
++ popl %eax
++ ret
+diff -Nru qemu-0.7.0/bsd.orig/i386/s_lrintl.S qemu-0.7.0/bsd/i386/s_lrintl.S
+--- qemu-0.7.0/bsd.orig/i386/s_lrintl.S Wed Dec 31 19:00:00 1969
++++ qemu-0.7.0/bsd/i386/s_lrintl.S Fri Apr 29 02:11:27 2005
+@@ -0,0 +1,34 @@
++/* Round argument to nearest integral value according to current rounding
++ direction.
++ Copyright (C) 1997 Free Software Foundation, Inc.
++ This file is part of the GNU C Library.
++ Contributed by Ulrich Drepper <drepper@cygnus.com>, 1997.
++
++ The GNU C Library is free software; you can redistribute it and/or
++ modify it under the terms of the GNU Lesser General Public
++ License as published by the Free Software Foundation; either
++ version 2.1 of the License, or (at your option) any later version.
++
++ The GNU C Library is distributed in the hope that it will be useful,
++ but WITHOUT ANY WARRANTY; without even the implied warranty of
++ MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
++ Lesser General Public License for more details.
++
++ You should have received a copy of the GNU Lesser General Public
++ License along with the GNU C Library; if not, write to the Free
++ Software Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA
++ 02111-1307 USA. */
++
++#include <machine/asm.h>
++#include <sysdep.h>
++
++ .text
++ENTRY(__lrintl)
++ fldt 4(%esp)
++ subl $4, %esp
++ fistpl (%esp)
++ fwait
++ popl %eax
++ ret
++END(__lrintl)
++weak_alias (__lrintl, lrintl)
+diff -Nru qemu-0.7.0/bsd.orig/i386/s_rintl.c qemu-0.7.0/bsd/i386/s_rintl.c
+--- qemu-0.7.0/bsd.orig/i386/s_rintl.c Wed Dec 31 19:00:00 1969
++++ qemu-0.7.0/bsd/i386/s_rintl.c Fri Apr 29 02:11:27 2005
+@@ -0,0 +1,18 @@
++/*
++ * Written by J.T. Conklin <jtc@netbsd.org>.
++ * Changes for long double by Ulrich Drepper <drepper@cygnus.com>
++ * Public domain.
++ */
++
++#include <sysdep.h>
++
++long double
++__rintl (long double x)
++{
++ long double res;
++
++ asm ("frndint" : "=t" (res) : "0" (x));
++ return res;
++}
++
++weak_alias (__rintl, rintl)
+diff -Nru qemu-0.7.0/bsd.orig/i386/s_round.c qemu-0.7.0/bsd/i386/s_round.c
+--- qemu-0.7.0/bsd.orig/i386/s_round.c Wed Dec 31 19:00:00 1969
++++ qemu-0.7.0/bsd/i386/s_round.c Fri Apr 29 02:11:27 2005
+@@ -0,0 +1,51 @@
++/*-
++ * Copyright (c) 2003, Steven G. Kargl
++ * 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 unmodified, this list of conditions, and the following
++ * disclaimer.
++ * 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.
++ *
++ * THIS SOFTWARE IS PROVIDED BY THE AUTHOR ``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 <sys/cdefs.h>
++__FBSDID("$FreeBSD: /tmp/pcvs/ports/emulators/qemu-devel/files/Attic/patch-libmath,v 1.1 2005-05-03 04:02:46 nork Exp $");
++
++#include <math.h>
++
++double
++round(double x)
++{
++ double t;
++
++ if (!isfinite(x))
++ return (x);
++
++ if (x >= 0.0) {
++ t = ceil(x);
++ if (t - x > 0.5)
++ t -= 1.0;
++ return (t);
++ } else {
++ t = ceil(-x);
++ if (t + x > 0.5)
++ t -= 1.0;
++ return (-t);
++ }
++}
+diff -Nru qemu-0.7.0/bsd.orig/i386/s_sinl.S qemu-0.7.0/bsd/i386/s_sinl.S
+--- qemu-0.7.0/bsd.orig/i386/s_sinl.S Wed Dec 31 19:00:00 1969
++++ qemu-0.7.0/bsd/i386/s_sinl.S Fri Apr 29 02:11:27 2005
+@@ -0,0 +1,32 @@
++/*
++ * Written by J.T. Conklin <jtc@netbsd.org>.
++ * Public domain.
++ *
++ * Adapted for `long double' by Ulrich Drepper <drepper@cygnus.com>.
++ */
++
++#include <machine/asm.h>
++#include <sysdep.h>
++
++RCSID("$NetBSD: $")
++
++ENTRY(__sinl)
++ fldt 4(%esp)
++ fsin
++ fnstsw %ax
++ testl $0x400,%eax
++ jnz 1f
++ ret
++ .align ALIGNARG(4)
++1: fldpi
++ fadd %st(0)
++ fxch %st(1)
++2: fprem1
++ fnstsw %ax
++ testl $0x400,%eax
++ jnz 2b
++ fstp %st(1)
++ fsin
++ ret
++END (__sinl)
++weak_alias (__sinl, sinl)
+diff -Nru qemu-0.7.0/bsd.orig/i386/s_tanl.S qemu-0.7.0/bsd/i386/s_tanl.S
+--- qemu-0.7.0/bsd.orig/i386/s_tanl.S Wed Dec 31 19:00:00 1969
++++ qemu-0.7.0/bsd/i386/s_tanl.S Fri Apr 29 02:11:27 2005
+@@ -0,0 +1,33 @@
++/*
++ * Written by J.T. Conklin <jtc@netbsd.org>.
++ * Public domain.
++ *
++ * Adapted for `long double' by Ulrich Drepper <drepper@cygnus.com>.
++ */
++
++#include <machine/asm.h>
++#include <sysdep.h>
++
++RCSID("$NetBSD: $")
++
++ENTRY(__tanl)
++ fldt 4(%esp)
++ fptan
++ fnstsw %ax
++ testl $0x400,%eax
++ jnz 1f
++ fstp %st(0)
++ ret
++1: fldpi
++ fadd %st(0)
++ fxch %st(1)
++2: fprem1
++ fstsw %ax
++ testl $0x400,%eax
++ jnz 2b
++ fstp %st(1)
++ fptan
++ fstp %st(0)
++ ret
++END (__tanl)
++weak_alias (__tanl, tanl)
+diff -Nru qemu-0.7.0/bsd.orig/sysdep.h qemu-0.7.0/bsd/sysdep.h
+--- qemu-0.7.0/bsd.orig/sysdep.h Wed Dec 31 19:00:00 1969
++++ qemu-0.7.0/bsd/sysdep.h Fri Apr 29 02:11:27 2005
+@@ -0,0 +1,20 @@
++#ifndef _QEMU_BSD_SYSDEP_H_
++#define _QEMU_BSD_SYSDEP_H_
++
++#include <sys/cdefs.h>
++
++#define HAVE_ELF
++
++#ifdef __ASSEMBLER__
++#define ALIGNARG(log2) 1<<log2
++#define ASM_TYPE_DIRECTIVE(name,typearg) .type name,typearg;
++#define ASM_SIZE_DIRECTIVE(name) .size name,.-name;
++#define END(x)
++#define strong_alias(sym,alias) .set alias,sym;
++#define weak_alias(sym,alias) .weak alias; .equ alias,sym;
++#else
++#define strong_alias(sym,alias) __strong_reference(sym,alias);
++#define weak_alias(sym,alias) __weak_reference(sym,alias);
++#endif
++
++#endif
diff --git a/emulators/qemu/Makefile b/emulators/qemu/Makefile
index 8d31b2da7e38..6d1a7c1a4c54 100644
--- a/emulators/qemu/Makefile
+++ b/emulators/qemu/Makefile
@@ -7,6 +7,7 @@
PORTNAME= qemu
PORTVERSION= 0.7.0
+POTREVISION= 1
CATEGORIES= emulators
MASTER_SITES= http://www.qemu.org/
diff --git a/emulators/qemu/files/patch-fbsd b/emulators/qemu/files/patch-fbsd
new file mode 100644
index 000000000000..dab1f93fc8a0
--- /dev/null
+++ b/emulators/qemu/files/patch-fbsd
@@ -0,0 +1,141 @@
+Index: qemu/Makefile
+@@ -13,7 +13,7 @@
+ endif
+ DOCS=qemu-doc.html qemu-tech.html qemu.1 qemu-img.1
+
+-all: dyngen$(EXESUF) $(TOOLS) $(DOCS)
++all: bsd/libmath.a dyngen$(EXESUF) $(TOOLS) $(DOCS)
+ for d in $(TARGET_DIRS); do \
+ $(MAKE) -C $$d $@ || exit 1 ; \
+ done
+@@ -25,6 +25,9 @@
+ endif
+ endif
+
++bsd/libmath.a:
++ ( cd bsd ; $(BSD_MAKE) CC=$(CC) )
++
+ qemu-img$(EXESUF): qemu-img.c block.c block-cow.c block-qcow.c aes.c block-vmdk.c block-cloop.c block-dmg.c block-bochs.c block-vpc.c
+ $(CC) -DQEMU_TOOL $(CFLAGS) $(LDFLAGS) $(DEFINES) -o $@ $^ -lz $(LIBS)
+
+@@ -33,6 +36,7 @@
+
+ clean:
+ # avoid old build problems by removing potentially incorrect old files
++ ( cd bsd ; $(BSD_MAKE) clean )
+ rm -f config.mak config.h op-i386.h opc-i386.h gen-op-i386.h op-arm.h opc-arm.h gen-op-arm.h
+ rm -f *.o *.a $(TOOLS) dyngen$(EXESUF) TAGS *.pod *~ */*~
+ $(MAKE) -C tests clean
+Index: qemu/Makefile.target
+@@ -391,8 +391,8 @@
+ VL_LDFLAGS+=-Wl,-G0 -Wl,-T,$(SRC_PATH)/ia64.ld
+ endif
+
+-$(QEMU_SYSTEM): $(VL_OBJS) libqemu.a
+- $(CC) $(VL_LDFLAGS) -o $@ $^ $(LIBS) $(SDL_LIBS) $(COCOA_LIBS) $(VL_LIBS)
++$(QEMU_SYSTEM): $(VL_OBJS) libqemu.a ../bsd/libmath.a
++ $(CC) $(VL_LDFLAGS) -o $@ $^ $(LIBS) $(SDL_LIBS) $(COCOA_LIBS) $(VL_LIBS) ../bsd/libmath.a
+
+ cocoa.o: cocoa.m
+ $(CC) $(CFLAGS) $(DEFINES) -c -o $@ $<
+Index: qemu/fpu/softfloat-native.c
+@@ -2,11 +2,15 @@
+ context is supported */
+ #include "softfloat.h"
+ #include <math.h>
++#if defined(__FreeBSD__) && __FreeBSD_version < 500000
++#include <ieeefp.h>
++#endif
+
+ void set_float_rounding_mode(int val STATUS_PARAM)
+ {
+ STATUS(float_rounding_mode) = val;
+-#if defined(_BSD) && !defined(__APPLE__)
++#if defined(_BSD) && !defined(__APPLE__) && \
++ (defined(__FreeBSD__) && __FreeBSD_version < 500000)
+ fpsetround(val);
+ #elif defined(__arm__)
+ /* nothing to do */
+@@ -22,7 +26,7 @@
+ }
+ #endif
+
+-#if defined(_BSD)
++#if defined(_BSD) && !defined(__FreeBSD__)
+ #define lrint(d) ((int32_t)rint(d))
+ #define llrint(d) ((int64_t)rint(d))
+ #endif
+Index: qemu/fpu/softfloat-native.h
+@@ -1,7 +1,17 @@
+ /* Native implementation of soft float functions */
+ #include <math.h>
+-#if defined(_BSD) && !defined(__APPLE__)
++#if defined(_BSD) && !defined(__APPLE__) && \
++ (!defined(__FreeBSD__) || __FreeBSD_version < 500000)
+ #include <ieeefp.h>
++#if defined(__FreeBSD__)
++#define isgreater(x, y) __builtin_isgreater((x), (y))
++#define isgreaterequal(x, y) __builtin_isgreaterequal((x), (y))
++#define isless(x, y) __builtin_isless((x), (y))
++#define islessequal(x, y) __builtin_islessequal((x), (y))
++#define islessgreater(x, y) __builtin_islessgreater((x), (y))
++#define isunordered(x, y) __builtin_isunordered((x), (y))
++long double fabsl(long double x);
++#endif
+ #else
+ #include <fenv.h>
+ #endif
+@@ -33,12 +43,13 @@
+ /*----------------------------------------------------------------------------
+ | Software IEC/IEEE floating-point rounding mode.
+ *----------------------------------------------------------------------------*/
+-#if defined(_BSD) && !defined(__APPLE__)
++#if defined(_BSD) && !defined(__APPLE__) && \
++ (!defined(__FreeBSD__) || __FreeBSD_version < 500000)
+ enum {
+ float_round_nearest_even = FP_RN,
+- float_round_down = FE_RM,
+- float_round_up = FE_RP,
+- float_round_to_zero = FE_RZ
++ float_round_down = FP_RM,
++ float_round_up = FP_RP,
++ float_round_to_zero = FP_RZ
+ };
+ #elif defined(__arm__)
+ enum {
+Index: qemu/fpu/softfloat.h
+@@ -84,7 +84,8 @@
+ #define FLOAT128
+ #else
+ /* native float support */
+-#if (defined(__i386__) || defined(__x86_64__)) && !defined(_BSD)
++#if (defined(__i386__) || defined(__x86_64__)) && \
++ (!defined(_BSD) || defined(__FreeBSD__))
+ #define FLOATX80
+ #endif
+ #endif /* !CONFIG_SOFTFLOAT */
+Index: qemu/target-ppc/op_helper.c
+@@ -303,6 +303,13 @@
+ FT0 = sqrt(FT0);
+ }
+
++#ifndef isnormal
++#define isnormal(x) \
++ ((sizeof (x) == sizeof (float)) ? __isnormalf(x) \
++ : (sizeof (x) == sizeof (double)) ? __isnormal(x) \
++ : __isnormall(x))
++#endif
++
+ void do_fres (void)
+ {
+ union {
+Index: qemu/x86_64.ld
+@@ -2,7 +2,7 @@
+ OUTPUT_FORMAT("elf64-x86-64", "elf64-x86-64", "elf64-x86-64")
+ OUTPUT_ARCH(i386:x86-64)
+ ENTRY(_start)
+-SEARCH_DIR("/lib64"); SEARCH_DIR("/usr/lib64"); SEARCH_DIR("/usr/local/lib64");
++SEARCH_DIR("/lib"); SEARCH_DIR("/usr/lib"); SEARCH_DIR("/usr/local/lib");
+ SECTIONS
+ {
+ /* Read-only sections, merged into text segment: */
diff --git a/emulators/qemu/files/patch-libmath b/emulators/qemu/files/patch-libmath
new file mode 100644
index 000000000000..08ba92fdb6c8
--- /dev/null
+++ b/emulators/qemu/files/patch-libmath
@@ -0,0 +1,2524 @@
+diff -Nru qemu-0.7.0/bsd.orig/Makefile qemu-0.7.0/bsd/Makefile
+--- qemu-0.7.0/bsd.orig/Makefile Wed Dec 31 19:00:00 1969
++++ qemu-0.7.0/bsd/Makefile Fri Apr 29 02:11:27 2005
+@@ -0,0 +1,32 @@
++SRCS= ${MACHINE_ARCH}/e_atan2l.c \
++ ${MACHINE_ARCH}/e_logl.S \
++ ${MACHINE_ARCH}/e_powl.S \
++ ${MACHINE_ARCH}/e_remainderl.S \
++ ${MACHINE_ARCH}/e_sqrtl.c \
++ ${MACHINE_ARCH}/s_ceill.S \
++ ${MACHINE_ARCH}/s_cosl.S \
++ ${MACHINE_ARCH}/s_floorl.S \
++ ${MACHINE_ARCH}/s_isnormal.c \
++ ${MACHINE_ARCH}/s_llrint.S \
++ ${MACHINE_ARCH}/s_llrintf.S \
++ ${MACHINE_ARCH}/s_llrintl.S \
++ ${MACHINE_ARCH}/s_lrint.S \
++ ${MACHINE_ARCH}/s_lrintf.S \
++ ${MACHINE_ARCH}/s_lrintl.S \
++ ${MACHINE_ARCH}/s_rintl.c \
++ ${MACHINE_ARCH}/s_round.c \
++ ${MACHINE_ARCH}/s_sinl.S \
++ ${MACHINE_ARCH}/s_tanl.S
++
++OBJS= ${SRCS:R:S/$/.o/}
++
++CFLAGS= -O2 -Wall -I.
++
++all: libmath.a
++
++libmath.a: ${OBJS}
++ rm -f $@
++ ${AR} rcs $@ ${OBJS:T}
++
++clean:
++ rm -f ${OBJS:T} libmath.a
+diff -Nru qemu-0.7.0/bsd.orig/amd64/e_atan2l.c qemu-0.7.0/bsd/amd64/e_atan2l.c
+--- qemu-0.7.0/bsd.orig/amd64/e_atan2l.c Wed Dec 31 19:00:00 1969
++++ qemu-0.7.0/bsd/amd64/e_atan2l.c Fri Apr 29 02:11:27 2005
+@@ -0,0 +1,20 @@
++/*
++ * Written by J.T. Conklin <jtc@netbsd.org>.
++ * Public domain.
++ *
++ * Adapted for `long double' by Ulrich Drepper <drepper@cygnus.com>.
++ */
++
++#include <sysdep.h>
++
++long double
++__ieee754_atan2l (long double y, long double x)
++{
++ long double res;
++
++ asm ("fpatan" : "=t" (res) : "u" (y), "0" (x) : "st(1)");
++
++ return res;
++}
++
++weak_alias(__ieee754_atan2l, atan2l)
+diff -Nru qemu-0.7.0/bsd.orig/amd64/e_logl.S qemu-0.7.0/bsd/amd64/e_logl.S
+--- qemu-0.7.0/bsd.orig/amd64/e_logl.S Wed Dec 31 19:00:00 1969
++++ qemu-0.7.0/bsd/amd64/e_logl.S Fri Apr 29 02:11:27 2005
+@@ -0,0 +1,59 @@
++/*
++ * Written by J.T. Conklin <jtc@netbsd.org>.
++ * Public domain.
++ *
++ * Adapted for `long double' by Ulrich Drepper <drepper@cygnus.com>.
++ * Adapted for x86-64 by Andreas Jaeger <aj@suse.de>.
++ */
++
++#include <machine/asm.h>
++#include <sysdep.h>
++
++RCSID("$NetBSD: $")
++
++
++#ifdef __ELF__
++ .section .rodata
++#else
++ .text
++#endif
++ .align ALIGNARG(4)
++ ASM_TYPE_DIRECTIVE(one,@object)
++one: .double 1.0
++ ASM_SIZE_DIRECTIVE(one)
++ /* It is not important that this constant is precise. It is only
++ a value which is known to be on the safe side for using the
++ fyl2xp1 instruction. */
++ ASM_TYPE_DIRECTIVE(limit,@object)
++limit: .double 0.29
++ ASM_SIZE_DIRECTIVE(limit)
++
++
++#ifdef PIC
++#define MO(op) op##(%rip)
++#else
++#define MO(op) op
++#endif
++
++ .text
++ENTRY(__ieee754_logl)
++ fldln2 // log(2)
++ fldt 8(%rsp) // x : log(2)
++ fld %st // x : x : log(2)
++ fsubl MO(one) // x-1 : x : log(2)
++ fld %st // x-1 : x-1 : x : log(2)
++ fabs // |x-1| : x-1 : x : log(2)
++ fcompl MO(limit) // x-1 : x : log(2)
++ fnstsw // x-1 : x : log(2)
++ andb $0x45, %ah
++ jz 2f
++ fstp %st(1) // x-1 : log(2)
++ fyl2xp1 // log(x)
++ ret
++
++2: fstp %st(0) // x : log(2)
++ fyl2x // log(x)
++ ret
++END (__ieee754_logl)
++
++weak_alias(__ieee754_logl,logl)
+diff -Nru qemu-0.7.0/bsd.orig/amd64/e_powl.S qemu-0.7.0/bsd/amd64/e_powl.S
+--- qemu-0.7.0/bsd.orig/amd64/e_powl.S Wed Dec 31 19:00:00 1969
++++ qemu-0.7.0/bsd/amd64/e_powl.S Fri Apr 29 02:11:27 2005
+@@ -0,0 +1,341 @@
++/* ix87 specific implementation of pow function.
++ Copyright (C) 1996, 1997, 1998, 1999, 2001, 2004 Free Software Foundation, Inc.
++ This file is part of the GNU C Library.
++ Contributed by Ulrich Drepper <drepper@cygnus.com>, 1996.
++
++ The GNU C Library is free software; you can redistribute it and/or
++ modify it under the terms of the GNU Lesser General Public
++ License as published by the Free Software Foundation; either
++ version 2.1 of the License, or (at your option) any later version.
++
++ The GNU C Library is distributed in the hope that it will be useful,
++ but WITHOUT ANY WARRANTY; without even the implied warranty of
++ MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
++ Lesser General Public License for more details.
++
++ You should have received a copy of the GNU Lesser General Public
++ License along with the GNU C Library; if not, write to the Free
++ Software Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA
++ 02111-1307 USA. */
++
++#include <machine/asm.h>
++#include <sysdep.h>
++
++#ifdef __ELF__
++ .section .rodata
++#else
++ .text
++#endif
++
++ .align ALIGNARG(4)
++ ASM_TYPE_DIRECTIVE(infinity,@object)
++inf_zero:
++infinity:
++ .byte 0, 0, 0, 0, 0, 0, 0xf0, 0x7f
++ ASM_SIZE_DIRECTIVE(infinity)
++ ASM_TYPE_DIRECTIVE(zero,@object)
++zero: .double 0.0
++ ASM_SIZE_DIRECTIVE(zero)
++ ASM_TYPE_DIRECTIVE(minf_mzero,@object)
++minf_mzero:
++minfinity:
++ .byte 0, 0, 0, 0, 0, 0, 0xf0, 0xff
++mzero:
++ .byte 0, 0, 0, 0, 0, 0, 0, 0x80
++ ASM_SIZE_DIRECTIVE(minf_mzero)
++ ASM_TYPE_DIRECTIVE(one,@object)
++one: .double 1.0
++ ASM_SIZE_DIRECTIVE(one)
++ ASM_TYPE_DIRECTIVE(limit,@object)
++limit: .double 0.29
++ ASM_SIZE_DIRECTIVE(limit)
++ ASM_TYPE_DIRECTIVE(p63,@object)
++p63:
++ .byte 0, 0, 0, 0, 0, 0, 0xe0, 0x43
++ ASM_SIZE_DIRECTIVE(p63)
++
++#ifdef PIC
++#define MO(op) op##(%rip)
++#else
++#define MO(op) op
++#endif
++
++ .text
++ENTRY(__ieee754_powl)
++ fldt 24(%rsp) // y
++ fxam
++
++
++ fnstsw
++ movb %ah, %dl
++ andb $0x45, %ah
++ cmpb $0x40, %ah // is y == 0 ?
++ je 11f
++
++ cmpb $0x05, %ah // is y == ±inf ?
++ je 12f
++
++ cmpb $0x01, %ah // is y == NaN ?
++ je 30f
++
++ fldt 8(%rsp) // x : y
++
++ fxam
++ fnstsw
++ movb %ah, %dh
++ andb $0x45, %ah
++ cmpb $0x40, %ah
++ je 20f // x is ±0
++
++ cmpb $0x05, %ah
++ je 15f // x is ±inf
++
++ fxch // y : x
++
++ /* fistpll raises invalid exception for |y| >= 1L<<63. */
++ fldl MO(p63) // 1L<<63 : y : x
++ fld %st(1) // y : 1L<<63 : y : x
++ fabs // |y| : 1L<<63 : y : x
++ fcomip %st(1), %st // 1L<<63 : y : x
++ fstp %st(0) // y : x
++ jnc 2f
++
++ /* First see whether `y' is a natural number. In this case we
++ can use a more precise algorithm. */
++ fld %st // y : y : x
++ fistpll -8(%rsp) // y : x
++ fildll -8(%rsp) // int(y) : y : x
++ fucomip %st(1),%st // y : x
++ jne 2f
++
++ /* OK, we have an integer value for y. */
++ mov -8(%rsp),%eax
++ mov -4(%rsp),%edx
++ orl $0, %edx
++ fstp %st(0) // x
++ jns 4f // y >= 0, jump
++ fdivrl MO(one) // 1/x (now referred to as x)
++ negl %eax
++ adcl $0, %edx
++ negl %edx
++4: fldl MO(one) // 1 : x
++ fxch
++
++6: shrdl $1, %edx, %eax
++ jnc 5f
++ fxch
++ fmul %st(1) // x : ST*x
++ fxch
++5: fmul %st(0), %st // x*x : ST*x
++ shrl $1, %edx
++ movl %eax, %ecx
++ orl %edx, %ecx
++ jnz 6b
++ fstp %st(0) // ST*x
++ ret
++
++ /* y is ±NAN */
++30: fldt 8(%rsp) // x : y
++ fldl MO(one) // 1.0 : x : y
++ fucomip %st(1),%st // x : y
++ je 31f
++ fxch // y : x
++31: fstp %st(1)
++ ret
++
++ .align ALIGNARG(4)
++2: /* y is a real number. */
++ fxch // x : y
++ fldl MO(one) // 1.0 : x : y
++ fld %st(1) // x : 1.0 : x : y
++ fsub %st(1) // x-1 : 1.0 : x : y
++ fabs // |x-1| : 1.0 : x : y
++ fcompl MO(limit) // 1.0 : x : y
++ fnstsw
++ fxch // x : 1.0 : y
++ test $4500,%eax
++ jz 7f
++ fsub %st(1) // x-1 : 1.0 : y
++ fyl2xp1 // log2(x) : y
++ jmp 8f
++
++7: fyl2x // log2(x) : y
++8: fmul %st(1) // y*log2(x) : y
++ fxam
++ fnstsw
++ andb $0x45, %ah
++ cmpb $0x05, %ah // is y*log2(x) == ±inf ?
++ je 28f
++ fst %st(1) // y*log2(x) : y*log2(x)
++ frndint // int(y*log2(x)) : y*log2(x)
++ fsubr %st, %st(1) // int(y*log2(x)) : fract(y*log2(x))
++ fxch // fract(y*log2(x)) : int(y*log2(x))
++ f2xm1 // 2^fract(y*log2(x))-1 : int(y*log2(x))
++ faddl MO(one) // 2^fract(y*log2(x)) : int(y*log2(x))
++ fscale // 2^fract(y*log2(x))*2^int(y*log2(x)) : int(y*log2(x))
++ fstp %st(1) // 2^fract(y*log2(x))*2^int(y*log2(x))
++ ret
++
++28: fstp %st(1) // y*log2(x)
++ fldl MO(one) // 1 : y*log2(x)
++ fscale // 2^(y*log2(x)) : y*log2(x)
++ fstp %st(1) // 2^(y*log2(x))
++ ret
++
++ // pow(x,±0) = 1
++ .align ALIGNARG(4)
++11: fstp %st(0) // pop y
++ fldl MO(one)
++ ret
++
++ // y == ±inf
++ .align ALIGNARG(4)
++12: fstp %st(0) // pop y
++ fldt 8(%rsp) // x
++ fabs
++ fcompl MO(one) // < 1, == 1, or > 1
++ fnstsw
++ andb $0x45, %ah
++ cmpb $0x45, %ah
++ je 13f // jump if x is NaN
++
++ cmpb $0x40, %ah
++ je 14f // jump if |x| == 1
++
++ shlb $1, %ah
++ xorb %ah, %dl
++ andl $2, %edx
++#ifdef PIC
++ lea inf_zero(%rip),%rcx
++ fldl (%rcx, %rdx, 4)
++#else
++ fldl inf_zero(,%rdx, 4)
++#endif
++ ret
++
++ .align ALIGNARG(4)
++14: fldl MO(one)
++ ret
++
++ .align ALIGNARG(4)
++13: fldt 8(%rsp) // load x == NaN
++ ret
++
++ .align ALIGNARG(4)
++ // x is ±inf
++15: fstp %st(0) // y
++ testb $2, %dh
++ jz 16f // jump if x == +inf
++
++ // We must find out whether y is an odd integer.
++ fld %st // y : y
++ fistpll -8(%rsp) // y
++ fildll -8(%rsp) // int(y) : y
++ fucomip %st(1),%st
++ ffreep %st // <empty>
++ jne 17f
++
++ // OK, the value is an integer, but is it odd?
++ mov -8(%rsp), %eax
++ mov -4(%rsp), %edx
++ andb $1, %al
++ jz 18f // jump if not odd
++ // It's an odd integer.
++ shrl $31, %edx
++#ifdef PIC
++ lea minf_mzero(%rip),%rcx
++ fldl (%rcx, %rdx, 8)
++#else
++ fldl minf_mzero(,%rdx, 8)
++#endif
++ ret
++
++ .align ALIGNARG(4)
++16: fcompl MO(zero)
++ fnstsw
++ shrl $5, %eax
++ andl $8, %eax
++#ifdef PIC
++ lea inf_zero(%rip),%rcx
++ fldl (%rcx, %rax, 1)
++#else
++ fldl inf_zero(,%rax, 1)
++#endif
++ ret
++
++ .align ALIGNARG(4)
++17: shll $30, %edx // sign bit for y in right position
++18: shrl $31, %edx
++#ifdef PIC
++ lea inf_zero(%rip),%rcx
++ fldl (%rcx, %rdx, 8)
++#else
++ fldl inf_zero(,%rdx, 8)
++#endif
++ ret
++
++ .align ALIGNARG(4)
++ // x is ±0
++20: fstp %st(0) // y
++ testb $2, %dl
++ jz 21f // y > 0
++
++ // x is ±0 and y is < 0. We must find out whether y is an odd integer.
++ testb $2, %dh
++ jz 25f
++
++ fld %st // y : y
++ fistpll -8(%rsp) // y
++ fildll -8(%rsp) // int(y) : y
++ fucomip %st(1),%st
++ ffreep %st // <empty>
++ jne 26f
++
++ // OK, the value is an integer, but is it odd?
++ mov -8(%rsp),%eax
++ mov -4(%rsp),%edx
++ andb $1, %al
++ jz 27f // jump if not odd
++ // It's an odd integer.
++ // Raise divide-by-zero exception and get minus infinity value.
++ fldl MO(one)
++ fdivl MO(zero)
++ fchs
++ ret
++
++25: fstp %st(0)
++26:
++27: // Raise divide-by-zero exception and get infinity value.
++ fldl MO(one)
++ fdivl MO(zero)
++ ret
++
++ .align ALIGNARG(4)
++ // x is ±0 and y is > 0. We must find out whether y is an odd integer.
++21: testb $2, %dh
++ jz 22f
++
++ fld %st // y : y
++ fistpll -8(%rsp) // y
++ fildll -8(%rsp) // int(y) : y
++ fucomip %st(1),%st
++ ffreep %st // <empty>
++ jne 23f
++
++ // OK, the value is an integer, but is it odd?
++ mov -8(%rsp),%eax
++ mov -4(%rsp),%edx
++ andb $1, %al
++ jz 24f // jump if not odd
++ // It's an odd integer.
++ fldl MO(mzero)
++ ret
++
++22: fstp %st(0)
++23:
++24: fldl MO(zero)
++ ret
++
++END(__ieee754_powl)
++
++weak_alias(__ieee754_powl,powl)
+diff -Nru qemu-0.7.0/bsd.orig/amd64/e_remainderl.S qemu-0.7.0/bsd/amd64/e_remainderl.S
+--- qemu-0.7.0/bsd.orig/amd64/e_remainderl.S Wed Dec 31 19:00:00 1969
++++ qemu-0.7.0/bsd/amd64/e_remainderl.S Fri Apr 29 02:11:27 2005
+@@ -0,0 +1,23 @@
++/*
++ * Written by J.T. Conklin <jtc@netbsd.org>.
++ * Public domain.
++ *
++ * Adapted for `long double' by Ulrich Drepper <drepper@cygnus.com>.
++ * Adapted for x86-64 by Andreas Jaeger <aj@suse.de>.
++ */
++
++#include <machine/asm.h>
++#include <sysdep.h>
++
++ENTRY(__ieee754_remainderl)
++ fldt 24(%rsp)
++ fldt 8(%rsp)
++1: fprem1
++ fstsw %ax
++ testl $0x400,%eax
++ jnz 1b
++ fstp %st(1)
++ ret
++END (__ieee754_remainderl)
++
++weak_alias(__ieee754_remainderl,remainderl)
+diff -Nru qemu-0.7.0/bsd.orig/amd64/e_sqrtl.c qemu-0.7.0/bsd/amd64/e_sqrtl.c
+--- qemu-0.7.0/bsd.orig/amd64/e_sqrtl.c Wed Dec 31 19:00:00 1969
++++ qemu-0.7.0/bsd/amd64/e_sqrtl.c Fri Apr 29 02:11:27 2005
+@@ -0,0 +1,20 @@
++/*
++ * Written by J.T. Conklin <jtc@netbsd.org>.
++ * Public domain.
++ *
++ * Adapted for `long double' by Ulrich Drepper <drepper@cygnus.com>.
++ */
++
++#include <sysdep.h>
++
++long double
++__ieee754_sqrtl (long double x)
++{
++ long double res;
++
++ asm ("fsqrt" : "=t" (res) : "0" (x));
++
++ return res;
++}
++
++weak_alias(__ieee754_sqrtl,sqrtl)
+diff -Nru qemu-0.7.0/bsd.orig/amd64/s_ceill.S qemu-0.7.0/bsd/amd64/s_ceill.S
+--- qemu-0.7.0/bsd.orig/amd64/s_ceill.S Wed Dec 31 19:00:00 1969
++++ qemu-0.7.0/bsd/amd64/s_ceill.S Fri Apr 29 02:22:18 2005
+@@ -0,0 +1,246 @@
++/*
++ * ====================================================
++ * Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved.
++ *
++ * Developed at SunPro, a Sun Microsystems, Inc. business.
++ * Permission to use, copy, modify, and distribute this
++ * software is freely granted, provided that this notice
++ * is preserved.
++ * ====================================================
++ *
++ * From: @(#)s_ceil.c 5.1 93/09/24
++ */
++/* XXX: generated from src/lib/msun/src/s_ceill.c */
++
++#include <machine/asm.h>
++
++__FBSDID("$FreeBSD: /tmp/pcvs/ports/emulators/qemu/files/patch-libmath,v 1.1 2005-05-03 04:02:46 nork Exp $")
++
++ .file "s_ceill.c"
++ .section .rodata.cst8,"aM",@progbits,8
++ .p2align 3
++.LC0:
++ .long 2281731484
++ .long 2117592124
++ .text
++ .p2align 4,,15
++.globl ceill
++ .type ceill, @function
++ceill:
++.LFB17:
++ pushq %rbp
++.LCFI0:
++ pushq %rbx
++.LCFI1:
++ fldt 24(%rsp)
++ movq $0, -16(%rsp)
++ fld %st(0)
++ fstpt -40(%rsp)
++ movl -32(%rsp), %edi
++ movq -40(%rsp), %rsi
++ movl %edi, -16(%rsp)
++ movl -16(%rsp), %r11d
++ movq %rsi, -24(%rsp)
++ movl %r11d, %r10d
++ andl $32767, %r10d
++ leal -16383(%r10), %r8d
++ cmpl $30, %r8d
++ jg .L2
++ testl %r8d, %r8d
++ js .L38
++ movl -20(%rsp), %r9d
++ leal 1(%r8), %ecx
++ mov -24(%rsp), %eax
++ movl $4294967295, %ebp
++ movq %rsi, -40(%rsp)
++ movl %edi, -32(%rsp)
++ shrq %cl, %rbp
++ mov %r9d, %ebx
++ movq %rbx, %rdx
++ andq %rbp, %rdx
++ orq %rax, %rdx
++ fldt -40(%rsp)
++ je .L42
++ ffreep %st(0)
++ testb $-128, -15(%rsp)
++ jne .L12
++ movl $31, %ecx
++ movl $1, %eax
++ subl %r8d, %ecx
++ salq %cl, %rax
++ addl %eax, %r9d
++ mov %r9d, %eax
++ cmpq %rbx, %rax
++ jae .L32
++ leal 1(%r10), %edx
++ movl %r11d, %eax
++ orl $-2147483648, %r9d
++ andw $-32768, %ax
++ andw $32767, %dx
++ orl %edx, %eax
++ movw %ax, -16(%rsp)
++.L32:
++ movl %r9d, -20(%rsp)
++.L12:
++ faddl .LC0(%rip)
++ fldz
++ fxch %st(1)
++ fucomip %st(1), %st
++ fstp %st(0)
++ jbe .L31
++ movl %ebp, %eax
++ movl $0, -24(%rsp)
++ notl %eax
++ andl %eax, %r9d
++ movl %r9d, -20(%rsp)
++ .p2align 4,,7
++.L31:
++ movq -24(%rsp), %rsi
++ movl -16(%rsp), %edi
++ movq %rsi, -40(%rsp)
++ movl %edi, -32(%rsp)
++ fldt -40(%rsp)
++ popq %rbx
++ popq %rbp
++ ret
++ .p2align 4,,7
++.L2:
++ cmpl $62, %r8d
++ jle .L45
++.L44:
++ ffreep %st(0)
++.L17:
++ movq %rsi, -40(%rsp)
++ movl %edi, -32(%rsp)
++ fldt -40(%rsp)
++ popq %rbx
++ popq %rbp
++ ret
++ .p2align 4,,7
++.L45:
++ movl -24(%rsp), %edx
++ leal 1(%r8), %ecx
++ movq $-1, %rbx
++ movq %rsi, -40(%rsp)
++ movl %edi, -32(%rsp)
++ shrq %cl, %rbx
++ mov %edx, %r9d
++ testq %rbx, %r9
++ fldt -40(%rsp)
++ je .L42
++ ffreep %st(0)
++ testb $-128, -15(%rsp)
++ jne .L20
++ cmpl $31, %r8d
++ je .L36
++ movl $63, %ecx
++ movl $1, %eax
++ subl %r8d, %ecx
++ salq %cl, %rax
++ leal (%rdx,%rax), %eax
++ movl %eax, -24(%rsp)
++ mov %eax, %eax
++ cmpq %r9, %rax
++ jae .L20
++.L36:
++ movl -20(%rsp), %eax
++ leal 1(%rax), %ecx
++ cmpl %eax, %ecx
++ jae .L34
++ leal 1(%r10), %edx
++ movl %r11d, %eax
++ orl $-2147483648, %ecx
++ andw $-32768, %ax
++ andw $32767, %dx
++ orl %edx, %eax
++ movw %ax, -16(%rsp)
++.L34:
++ movl %ecx, -20(%rsp)
++ .p2align 4,,7
++.L20:
++ faddl .LC0(%rip)
++ fldz
++ fxch %st(1)
++ fucomip %st(1), %st
++ fstp %st(0)
++ jbe .L31
++ movl %ebx, %eax
++ notl %eax
++ andl %eax, -24(%rsp)
++ jmp .L31
++ .p2align 4,,7
++.L42:
++ fstp %st(1)
++ popq %rbx
++ popq %rbp
++ ret
++ .p2align 4,,7
++.L38:
++ fldl .LC0(%rip)
++ faddp %st, %st(1)
++ fldz
++ fxch %st(1)
++ fucomip %st(1), %st
++ jbe .L44
++ testl %r10d, %r10d
++ jle .L39
++.L7:
++ movabsq $-9223372036854775808, %rsi
++ movl $16383, %edi
++ testb $-128, -15(%rsp)
++ movq %rsi, -40(%rsp)
++ movl %edi, -32(%rsp)
++ fldt -40(%rsp)
++ fcmovne %st(1), %st
++ fstp %st(1)
++ fstpt -40(%rsp)
++ movq -40(%rsp), %rsi
++ movl -32(%rsp), %edi
++ movq %rsi, -24(%rsp)
++ movl %edi, -16(%rsp)
++ jmp .L17
++.L39:
++ movl -24(%rsp), %eax
++ orl -20(%rsp), %eax
++ je .L44
++ jmp .L7
++.LFE17:
++ .size ceill, .-ceill
++ .section .eh_frame,"a",@progbits
++.Lframe1:
++ .long .LECIE1-.LSCIE1
++.LSCIE1:
++ .long 0x0
++ .byte 0x1
++ .string ""
++ .uleb128 0x1
++ .sleb128 -8
++ .byte 0x10
++ .byte 0xc
++ .uleb128 0x7
++ .uleb128 0x8
++ .byte 0x90
++ .uleb128 0x1
++ .p2align 3
++.LECIE1:
++.LSFDE1:
++ .long .LEFDE1-.LASFDE1
++.LASFDE1:
++ .long .LASFDE1-.Lframe1
++ .quad .LFB17
++ .quad .LFE17-.LFB17
++ .byte 0x4
++ .long .LCFI0-.LFB17
++ .byte 0xe
++ .uleb128 0x10
++ .byte 0x4
++ .long .LCFI1-.LCFI0
++ .byte 0xe
++ .uleb128 0x18
++ .byte 0x83
++ .uleb128 0x3
++ .byte 0x86
++ .uleb128 0x2
++ .p2align 3
++.LEFDE1:
++ .ident "GCC: (GNU) 3.4.4 [FreeBSD] 20050421"
+diff -Nru qemu-0.7.0/bsd.orig/amd64/s_cosl.S qemu-0.7.0/bsd/amd64/s_cosl.S
+--- qemu-0.7.0/bsd.orig/amd64/s_cosl.S Wed Dec 31 19:00:00 1969
++++ qemu-0.7.0/bsd/amd64/s_cosl.S Fri Apr 29 02:11:27 2005
+@@ -0,0 +1,33 @@
++/*
++ * Written by J.T. Conklin <jtc@netbsd.org>.
++ * Public domain.
++ *
++ * Adapted for `long double' by Ulrich Drepper <drepper@cygnus.com>.
++ * Adapted for x86-64 by Andreas Jaeger <aj@suse.de>.
++ */
++
++#include <machine/asm.h>
++#include <sysdep.h>
++
++RCSID("$NetBSD: $")
++
++ENTRY(__cosl)
++ fldt 8(%rsp)
++ fcos
++ fnstsw %ax
++ testl $0x400,%eax
++ jnz 1f
++ ret
++ .align ALIGNARG(4)
++1: fldpi
++ fadd %st(0)
++ fxch %st(1)
++2: fprem1
++ fnstsw %ax
++ testl $0x400,%eax
++ jnz 2b
++ fstp %st(1)
++ fcos
++ ret
++END (__cosl)
++weak_alias (__cosl, cosl)
+diff -Nru qemu-0.7.0/bsd.orig/amd64/s_floorl.S qemu-0.7.0/bsd/amd64/s_floorl.S
+--- qemu-0.7.0/bsd.orig/amd64/s_floorl.S Wed Dec 31 19:00:00 1969
++++ qemu-0.7.0/bsd/amd64/s_floorl.S Fri Apr 29 02:24:32 2005
+@@ -0,0 +1,247 @@
++/*
++ * ====================================================
++ * Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved.
++ *
++ * Developed at SunPro, a Sun Microsystems, Inc. business.
++ * Permission to use, copy, modify, and distribute this
++ * software is freely granted, provided that this notice
++ * is preserved.
++ * ====================================================
++ *
++ * From: @(#)s_floor.c 5.1 93/09/24
++ */
++/* XXX: generated from src/lib/msun/src/s_floorl.c */
++
++#include <machine/asm.h>
++
++__FBSDID("$FreeBSD: /tmp/pcvs/ports/emulators/qemu/files/patch-libmath,v 1.1 2005-05-03 04:02:46 nork Exp $")
++
++ .file "s_floorl.c"
++ .section .rodata.cst8,"aM",@progbits,8
++ .p2align 3
++.LC0:
++ .long 2281731484
++ .long 2117592124
++ .section .rodata.cst4,"aM",@progbits,4
++ .p2align 2
++.LC2:
++ .long 3212836864
++ .text
++ .p2align 4,,15
++.globl floorl
++ .type floorl, @function
++floorl:
++.LFB17:
++ pushq %rbp
++.LCFI0:
++ pushq %rbx
++.LCFI1:
++ fldt 24(%rsp)
++ movq $0, -16(%rsp)
++ fld %st(0)
++ fstpt -40(%rsp)
++ movl -32(%rsp), %edi
++ movq -40(%rsp), %rsi
++ movl %edi, -16(%rsp)
++ movl -16(%rsp), %r11d
++ movq %rsi, -24(%rsp)
++ movl %r11d, %r10d
++ andl $32767, %r10d
++ leal -16383(%r10), %r8d
++ cmpl $30, %r8d
++ jg .L2
++ testl %r8d, %r8d
++ js .L38
++ movl -20(%rsp), %r9d
++ leal 1(%r8), %ecx
++ mov -24(%rsp), %eax
++ movl $4294967295, %ebp
++ movq %rsi, -40(%rsp)
++ movl %edi, -32(%rsp)
++ shrq %cl, %rbp
++ mov %r9d, %ebx
++ movq %rbx, %rdx
++ andq %rbp, %rdx
++ orq %rax, %rdx
++ fldt -40(%rsp)
++ je .L42
++ ffreep %st(0)
++ testb $-128, -15(%rsp)
++ je .L12
++ movl $31, %ecx
++ movl $1, %eax
++ subl %r8d, %ecx
++ salq %cl, %rax
++ addl %eax, %r9d
++ mov %r9d, %eax
++ cmpq %rbx, %rax
++ jae .L32
++ leal 1(%r10), %edx
++ movl %r11d, %eax
++ orl $-2147483648, %r9d
++ andw $-32768, %ax
++ andw $32767, %dx
++ orl %edx, %eax
++ movw %ax, -16(%rsp)
++.L32:
++ movl %r9d, -20(%rsp)
++.L12:
++ faddl .LC0(%rip)
++ fldz
++ fxch %st(1)
++ fucomip %st(1), %st
++ fstp %st(0)
++ jbe .L31
++ movl %ebp, %eax
++ movl $0, -24(%rsp)
++ notl %eax
++ andl %eax, %r9d
++ movl %r9d, -20(%rsp)
++ .p2align 4,,7
++.L31:
++ movq -24(%rsp), %rsi
++ movl -16(%rsp), %edi
++ movq %rsi, -40(%rsp)
++ movl %edi, -32(%rsp)
++ fldt -40(%rsp)
++ popq %rbx
++ popq %rbp
++ ret
++ .p2align 4,,7
++.L2:
++ cmpl $62, %r8d
++ jle .L45
++.L44:
++ ffreep %st(0)
++.L17:
++ movq %rsi, -40(%rsp)
++ movl %edi, -32(%rsp)
++ fldt -40(%rsp)
++ popq %rbx
++ popq %rbp
++ ret
++ .p2align 4,,7
++.L45:
++ movl -24(%rsp), %edx
++ leal 1(%r8), %ecx
++ movq $-1, %rbx
++ movq %rsi, -40(%rsp)
++ movl %edi, -32(%rsp)
++ shrq %cl, %rbx
++ mov %edx, %r9d
++ testq %rbx, %r9
++ fldt -40(%rsp)
++ je .L42
++ ffreep %st(0)
++ testb $-128, -15(%rsp)
++ je .L20
++ cmpl $31, %r8d
++ je .L36
++ movl $63, %ecx
++ movl $1, %eax
++ subl %r8d, %ecx
++ salq %cl, %rax
++ leal (%rdx,%rax), %eax
++ movl %eax, -24(%rsp)
++ mov %eax, %eax
++ cmpq %r9, %rax
++ jae .L20
++.L36:
++ movl -20(%rsp), %eax
++ leal 1(%rax), %ecx
++ cmpl %eax, %ecx
++ jae .L34
++ leal 1(%r10), %edx
++ movl %r11d, %eax
++ orl $-2147483648, %ecx
++ andw $-32768, %ax
++ andw $32767, %dx
++ orl %edx, %eax
++ movw %ax, -16(%rsp)
++.L34:
++ movl %ecx, -20(%rsp)
++ .p2align 4,,7
++.L20:
++ faddl .LC0(%rip)
++ fldz
++ fxch %st(1)
++ fucomip %st(1), %st
++ fstp %st(0)
++ jbe .L31
++ movl %ebx, %eax
++ notl %eax
++ andl %eax, -24(%rsp)
++ jmp .L31
++ .p2align 4,,7
++.L42:
++ fstp %st(1)
++ popq %rbx
++ popq %rbp
++ ret
++ .p2align 4,,7
++.L38:
++ fldl .LC0(%rip)
++ faddp %st, %st(1)
++ fldz
++ fxch %st(1)
++ fucomip %st(1), %st
++ jbe .L44
++ testl %r10d, %r10d
++ jle .L39
++.L7:
++ testb $-128, -15(%rsp)
++ je .L9
++ ffreep %st(0)
++ flds .LC2(%rip)
++.L9:
++ fstpt -40(%rsp)
++ movq -40(%rsp), %rsi
++ movl -32(%rsp), %edi
++ movq %rsi, -24(%rsp)
++ movl %edi, -16(%rsp)
++ jmp .L17
++.L39:
++ movl -24(%rsp), %eax
++ orl -20(%rsp), %eax
++ je .L44
++ jmp .L7
++.LFE17:
++ .size floorl, .-floorl
++ .section .eh_frame,"a",@progbits
++.Lframe1:
++ .long .LECIE1-.LSCIE1
++.LSCIE1:
++ .long 0x0
++ .byte 0x1
++ .string ""
++ .uleb128 0x1
++ .sleb128 -8
++ .byte 0x10
++ .byte 0xc
++ .uleb128 0x7
++ .uleb128 0x8
++ .byte 0x90
++ .uleb128 0x1
++ .p2align 3
++.LECIE1:
++.LSFDE1:
++ .long .LEFDE1-.LASFDE1
++.LASFDE1:
++ .long .LASFDE1-.Lframe1
++ .quad .LFB17
++ .quad .LFE17-.LFB17
++ .byte 0x4
++ .long .LCFI0-.LFB17
++ .byte 0xe
++ .uleb128 0x10
++ .byte 0x4
++ .long .LCFI1-.LCFI0
++ .byte 0xe
++ .uleb128 0x18
++ .byte 0x83
++ .uleb128 0x3
++ .byte 0x86
++ .uleb128 0x2
++ .p2align 3
++.LEFDE1:
++ .ident "GCC: (GNU) 3.4.4 [FreeBSD] 20050421"
+diff -Nru qemu-0.7.0/bsd.orig/amd64/s_isnormal.c qemu-0.7.0/bsd/amd64/s_isnormal.c
+--- qemu-0.7.0/bsd.orig/amd64/s_isnormal.c Wed Dec 31 19:00:00 1969
++++ qemu-0.7.0/bsd/amd64/s_isnormal.c Fri Apr 29 02:11:27 2005
+@@ -0,0 +1,85 @@
++/*-
++ * Copyright (c) 2003 Mike Barcroft <mike@FreeBSD.org>
++ * Copyright (c) 2002-2004 David Schultz <das@FreeBSD.ORG>
++ * 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, this list of conditions and the following disclaimer.
++ * 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.
++ *
++ * THIS SOFTWARE IS PROVIDED BY THE AUTHOR AND CONTRIBUTORS ``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 OR CONTRIBUTORS 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.
++ *
++ * $FreeBSD: /tmp/pcvs/ports/emulators/qemu/files/patch-libmath,v 1.1 2005-05-03 04:02:46 nork Exp $
++ */
++
++union IEEEf2bits {
++ float f;
++ struct {
++ unsigned int man :23;
++ unsigned int exp :8;
++ unsigned int sign :1;
++ } bits;
++};
++
++union IEEEd2bits {
++ double d;
++ struct {
++ unsigned int manl :32;
++ unsigned int manh :20;
++ unsigned int exp :11;
++ unsigned int sign :1;
++ } bits;
++};
++
++union IEEEl2bits {
++ long double e;
++ struct {
++ unsigned int manl :32;
++ unsigned int manh :32;
++ unsigned int exp :15;
++ unsigned int sign :1;
++ unsigned int junk :16;
++ } bits;
++};
++
++int
++__isnormal(double d)
++{
++ union IEEEd2bits u;
++
++ u.d = d;
++ return (u.bits.exp != 0 && u.bits.exp != 2047);
++}
++
++int
++__isnormalf(float f)
++{
++ union IEEEf2bits u;
++
++ u.f = f;
++ return (u.bits.exp != 0 && u.bits.exp != 255);
++}
++
++int
++__isnormall(long double e)
++{
++ union IEEEl2bits u;
++
++ u.e = e;
++ return (u.bits.exp != 0 && u.bits.exp != 32767);
++}
+diff -Nru qemu-0.7.0/bsd.orig/amd64/s_llrint.S qemu-0.7.0/bsd/amd64/s_llrint.S
+--- qemu-0.7.0/bsd.orig/amd64/s_llrint.S Wed Dec 31 19:00:00 1969
++++ qemu-0.7.0/bsd/amd64/s_llrint.S Fri Apr 29 02:11:27 2005
+@@ -0,0 +1,6 @@
++#include <machine/asm.h>
++__FBSDID("$FreeBSD: /tmp/pcvs/ports/emulators/qemu/files/patch-libmath,v 1.1 2005-05-03 04:02:46 nork Exp $")
++
++/* sizeof(long) == sizeof(long long) */
++#define fn llrint
++#include "s_lrint.S"
+diff -Nru qemu-0.7.0/bsd.orig/amd64/s_llrintf.S qemu-0.7.0/bsd/amd64/s_llrintf.S
+--- qemu-0.7.0/bsd.orig/amd64/s_llrintf.S Wed Dec 31 19:00:00 1969
++++ qemu-0.7.0/bsd/amd64/s_llrintf.S Fri Apr 29 02:11:27 2005
+@@ -0,0 +1,6 @@
++#include <machine/asm.h>
++__FBSDID("$FreeBSD: /tmp/pcvs/ports/emulators/qemu/files/patch-libmath,v 1.1 2005-05-03 04:02:46 nork Exp $")
++
++/* sizeof(long) == sizeof(long long) */
++#define fn llrintf
++#include "s_lrintf.S"
+diff -Nru qemu-0.7.0/bsd.orig/amd64/s_llrintl.S qemu-0.7.0/bsd/amd64/s_llrintl.S
+--- qemu-0.7.0/bsd.orig/amd64/s_llrintl.S Wed Dec 31 19:00:00 1969
++++ qemu-0.7.0/bsd/amd64/s_llrintl.S Fri Apr 29 02:11:27 2005
+@@ -0,0 +1,35 @@
++/* Round argument to nearest integral value according to current rounding
++ direction.
++ Copyright (C) 1997, 2002 Free Software Foundation, Inc.
++ This file is part of the GNU C Library.
++
++ The GNU C Library is free software; you can redistribute it and/or
++ modify it under the terms of the GNU Lesser General Public
++ License as published by the Free Software Foundation; either
++ version 2.1 of the License, or (at your option) any later version.
++
++ The GNU C Library is distributed in the hope that it will be useful,
++ but WITHOUT ANY WARRANTY; without even the implied warranty of
++ MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
++ Lesser General Public License for more details.
++
++ You should have received a copy of the GNU Lesser General Public
++ License along with the GNU C Library; if not, write to the Free
++ Software Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA
++ 02111-1307 USA. */
++
++#include <machine/asm.h>
++#include <sysdep.h>
++
++ .text
++ENTRY(__llrintl)
++ fldt 8(%rsp)
++ fistpll -8(%rsp)
++ fwait
++ movq -8(%rsp),%rax
++ ret
++END(__llrintl)
++weak_alias (__llrintl, llrintl)
++strong_alias (__llrintl, __lrintl)
++weak_alias (__llrintl, lrintl)
++
+diff -Nru qemu-0.7.0/bsd.orig/amd64/s_lrint.S qemu-0.7.0/bsd/amd64/s_lrint.S
+--- qemu-0.7.0/bsd.orig/amd64/s_lrint.S Wed Dec 31 19:00:00 1969
++++ qemu-0.7.0/bsd/amd64/s_lrint.S Fri Apr 29 02:11:27 2005
+@@ -0,0 +1,36 @@
++/*-
++ * Copyright (c) 2005 David Schultz <das@FreeBSD.ORG>
++ * 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, this list of conditions and the following disclaimer.
++ * 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.
++ *
++ * THIS SOFTWARE IS PROVIDED BY THE AUTHOR AND CONTRIBUTORS ``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 OR CONTRIBUTORS 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 <machine/asm.h>
++
++#ifndef fn
++__FBSDID("$FreeBSD: /tmp/pcvs/ports/emulators/qemu/files/patch-libmath,v 1.1 2005-05-03 04:02:46 nork Exp $")
++#define fn lrint
++#endif
++
++ENTRY(fn)
++ cvtsd2si %xmm0, %rax
++ ret
+diff -Nru qemu-0.7.0/bsd.orig/amd64/s_lrintf.S qemu-0.7.0/bsd/amd64/s_lrintf.S
+--- qemu-0.7.0/bsd.orig/amd64/s_lrintf.S Wed Dec 31 19:00:00 1969
++++ qemu-0.7.0/bsd/amd64/s_lrintf.S Fri Apr 29 02:11:27 2005
+@@ -0,0 +1,36 @@
++/*-
++ * Copyright (c) 2005 David Schultz <das@FreeBSD.ORG>
++ * 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, this list of conditions and the following disclaimer.
++ * 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.
++ *
++ * THIS SOFTWARE IS PROVIDED BY THE AUTHOR AND CONTRIBUTORS ``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 OR CONTRIBUTORS 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 <machine/asm.h>
++
++#ifndef fn
++__FBSDID("$FreeBSD: /tmp/pcvs/ports/emulators/qemu/files/patch-libmath,v 1.1 2005-05-03 04:02:46 nork Exp $")
++#define fn lrintf
++#endif
++
++ENTRY(fn)
++ cvtss2si %xmm0, %rax
++ ret
+diff -Nru qemu-0.7.0/bsd.orig/amd64/s_lrintl.S qemu-0.7.0/bsd/amd64/s_lrintl.S
+--- qemu-0.7.0/bsd.orig/amd64/s_lrintl.S Wed Dec 31 19:00:00 1969
++++ qemu-0.7.0/bsd/amd64/s_lrintl.S Fri Apr 29 02:11:27 2005
+@@ -0,0 +1 @@
++/* Not needed, see s_llrintl.S. */
+diff -Nru qemu-0.7.0/bsd.orig/amd64/s_rintl.c qemu-0.7.0/bsd/amd64/s_rintl.c
+--- qemu-0.7.0/bsd.orig/amd64/s_rintl.c Wed Dec 31 19:00:00 1969
++++ qemu-0.7.0/bsd/amd64/s_rintl.c Fri Apr 29 02:11:27 2005
+@@ -0,0 +1,18 @@
++/*
++ * Written by J.T. Conklin <jtc@netbsd.org>.
++ * Changes for long double by Ulrich Drepper <drepper@cygnus.com>
++ * Public domain.
++ */
++
++#include <sysdep.h>
++
++long double
++__rintl (long double x)
++{
++ long double res;
++
++ asm ("frndint" : "=t" (res) : "0" (x));
++ return res;
++}
++
++weak_alias (__rintl, rintl)
+diff -Nru qemu-0.7.0/bsd.orig/amd64/s_round.c qemu-0.7.0/bsd/amd64/s_round.c
+--- qemu-0.7.0/bsd.orig/amd64/s_round.c Wed Dec 31 19:00:00 1969
++++ qemu-0.7.0/bsd/amd64/s_round.c Fri Apr 29 02:11:27 2005
+@@ -0,0 +1,51 @@
++/*-
++ * Copyright (c) 2003, Steven G. Kargl
++ * 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 unmodified, this list of conditions, and the following
++ * disclaimer.
++ * 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.
++ *
++ * THIS SOFTWARE IS PROVIDED BY THE AUTHOR ``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 <sys/cdefs.h>
++__FBSDID("$FreeBSD: /tmp/pcvs/ports/emulators/qemu/files/patch-libmath,v 1.1 2005-05-03 04:02:46 nork Exp $");
++
++#include <math.h>
++
++double
++round(double x)
++{
++ double t;
++
++ if (!isfinite(x))
++ return (x);
++
++ if (x >= 0.0) {
++ t = ceil(x);
++ if (t - x > 0.5)
++ t -= 1.0;
++ return (t);
++ } else {
++ t = ceil(-x);
++ if (t + x > 0.5)
++ t -= 1.0;
++ return (-t);
++ }
++}
+diff -Nru qemu-0.7.0/bsd.orig/amd64/s_sinl.S qemu-0.7.0/bsd/amd64/s_sinl.S
+--- qemu-0.7.0/bsd.orig/amd64/s_sinl.S Wed Dec 31 19:00:00 1969
++++ qemu-0.7.0/bsd/amd64/s_sinl.S Fri Apr 29 02:11:27 2005
+@@ -0,0 +1,31 @@
++/*
++ * Written by J.T. Conklin <jtc@netbsd.org>.
++ * Public domain.
++ *
++ * Adapted for `long double' by Ulrich Drepper <drepper@cygnus.com>.
++ * Adapted for x86-64 by Andreas Jaeger <aj@suse.de>.
++ */
++
++#include <machine/asm.h>
++#include <sysdep.h>
++
++ENTRY(__sinl)
++ fldt 8(%rsp)
++ fsin
++ fnstsw %ax
++ testl $0x400,%eax
++ jnz 1f
++ ret
++ .align ALIGNARG(4)
++1: fldpi
++ fadd %st(0)
++ fxch %st(1)
++2: fprem1
++ fnstsw %ax
++ testl $0x400,%eax
++ jnz 2b
++ fstp %st(1)
++ fsin
++ ret
++END (__sinl)
++weak_alias (__sinl, sinl)
+diff -Nru qemu-0.7.0/bsd.orig/amd64/s_tanl.S qemu-0.7.0/bsd/amd64/s_tanl.S
+--- qemu-0.7.0/bsd.orig/amd64/s_tanl.S Wed Dec 31 19:00:00 1969
++++ qemu-0.7.0/bsd/amd64/s_tanl.S Fri Apr 29 02:11:27 2005
+@@ -0,0 +1,34 @@
++/*
++ * Written by J.T. Conklin <jtc@netbsd.org>.
++ * Public domain.
++ *
++ * Adapted for `long double' by Ulrich Drepper <drepper@cygnus.com>.
++ * Adapted for x86-64 by Andreas Jaeger <aj@suse.de>.
++ */
++
++#include <machine/asm.h>
++#include <sysdep.h>
++
++RCSID("$NetBSD: $")
++
++ENTRY(__tanl)
++ fldt 8(%rsp)
++ fptan
++ fnstsw %ax
++ testl $0x400,%eax
++ jnz 1f
++ fstp %st(0)
++ ret
++1: fldpi
++ fadd %st(0)
++ fxch %st(1)
++2: fprem1
++ fstsw %ax
++ testl $0x400,%eax
++ jnz 2b
++ fstp %st(1)
++ fptan
++ fstp %st(0)
++ ret
++END (__tanl)
++weak_alias (__tanl, tanl)
+diff -Nru qemu-0.7.0/bsd.orig/i386/e_atan2l.c qemu-0.7.0/bsd/i386/e_atan2l.c
+--- qemu-0.7.0/bsd.orig/i386/e_atan2l.c Wed Dec 31 19:00:00 1969
++++ qemu-0.7.0/bsd/i386/e_atan2l.c Fri Apr 29 02:11:27 2005
+@@ -0,0 +1,20 @@
++/*
++ * Written by J.T. Conklin <jtc@netbsd.org>.
++ * Public domain.
++ *
++ * Adapted for `long double' by Ulrich Drepper <drepper@cygnus.com>.
++ */
++
++#include <sysdep.h>
++
++long double
++__ieee754_atan2l (long double y, long double x)
++{
++ long double res;
++
++ asm ("fpatan" : "=t" (res) : "u" (y), "0" (x) : "st(1)");
++
++ return res;
++}
++
++weak_alias(__ieee754_atan2l, atan2l)
+diff -Nru qemu-0.7.0/bsd.orig/i386/e_logl.S qemu-0.7.0/bsd/i386/e_logl.S
+--- qemu-0.7.0/bsd.orig/i386/e_logl.S Wed Dec 31 19:00:00 1969
++++ qemu-0.7.0/bsd/i386/e_logl.S Fri Apr 29 02:11:27 2005
+@@ -0,0 +1,63 @@
++/*
++ * Written by J.T. Conklin <jtc@netbsd.org>.
++ * Public domain.
++ *
++ * Adapted for `long double' by Ulrich Drepper <drepper@cygnus.com>.
++ */
++
++#include <machine/asm.h>
++#include <sysdep.h>
++
++RCSID("$NetBSD: $")
++
++
++#ifdef __ELF__
++ .section .rodata
++#else
++ .text
++#endif
++ .align ALIGNARG(4)
++ ASM_TYPE_DIRECTIVE(one,@object)
++one: .double 1.0
++ ASM_SIZE_DIRECTIVE(one)
++ /* It is not important that this constant is precise. It is only
++ a value which is known to be on the safe side for using the
++ fyl2xp1 instruction. */
++ ASM_TYPE_DIRECTIVE(limit,@object)
++limit: .double 0.29
++ ASM_SIZE_DIRECTIVE(limit)
++
++
++#ifdef PIC
++#define MO(op) op##@GOTOFF(%edx)
++#else
++#define MO(op) op
++#endif
++
++ .text
++ENTRY(__ieee754_logl)
++ fldln2 // log(2)
++ fldt 4(%esp) // x : log(2)
++#ifdef PIC
++ call 1f
++1: popl %edx
++ addl $_GLOBAL_OFFSET_TABLE_+[.-1b], %edx
++#endif
++ fld %st // x : x : log(2)
++ fsubl MO(one) // x-1 : x : log(2)
++ fld %st // x-1 : x-1 : x : log(2)
++ fabs // |x-1| : x-1 : x : log(2)
++ fcompl MO(limit) // x-1 : x : log(2)
++ fnstsw // x-1 : x : log(2)
++ andb $0x45, %ah
++ jz 2f
++ fstp %st(1) // x-1 : log(2)
++ fyl2xp1 // log(x)
++ ret
++
++2: fstp %st(0) // x : log(2)
++ fyl2x // log(x)
++ ret
++END (__ieee754_logl)
++
++weak_alias(__ieee754_logl,logl)
+diff -Nru qemu-0.7.0/bsd.orig/i386/e_powl.S qemu-0.7.0/bsd/i386/e_powl.S
+--- qemu-0.7.0/bsd.orig/i386/e_powl.S Wed Dec 31 19:00:00 1969
++++ qemu-0.7.0/bsd/i386/e_powl.S Fri Apr 29 02:11:27 2005
+@@ -0,0 +1,341 @@
++/* ix87 specific implementation of pow function.
++ Copyright (C) 1996, 1997, 1998, 1999, 2001, 2004
++ Free Software Foundation, Inc.
++ This file is part of the GNU C Library.
++ Contributed by Ulrich Drepper <drepper@cygnus.com>, 1996.
++
++ The GNU C Library is free software; you can redistribute it and/or
++ modify it under the terms of the GNU Lesser General Public
++ License as published by the Free Software Foundation; either
++ version 2.1 of the License, or (at your option) any later version.
++
++ The GNU C Library is distributed in the hope that it will be useful,
++ but WITHOUT ANY WARRANTY; without even the implied warranty of
++ MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
++ Lesser General Public License for more details.
++
++ You should have received a copy of the GNU Lesser General Public
++ License along with the GNU C Library; if not, write to the Free
++ Software Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA
++ 02111-1307 USA. */
++
++#include <machine/asm.h>
++#include <sysdep.h>
++
++#ifdef __ELF__
++ .section .rodata
++#else
++ .text
++#endif
++
++ .align ALIGNARG(4)
++ ASM_TYPE_DIRECTIVE(infinity,@object)
++inf_zero:
++infinity:
++ .byte 0, 0, 0, 0, 0, 0, 0xf0, 0x7f
++ ASM_SIZE_DIRECTIVE(infinity)
++ ASM_TYPE_DIRECTIVE(zero,@object)
++zero: .double 0.0
++ ASM_SIZE_DIRECTIVE(zero)
++ ASM_TYPE_DIRECTIVE(minf_mzero,@object)
++minf_mzero:
++minfinity:
++ .byte 0, 0, 0, 0, 0, 0, 0xf0, 0xff
++mzero:
++ .byte 0, 0, 0, 0, 0, 0, 0, 0x80
++ ASM_SIZE_DIRECTIVE(minf_mzero)
++ ASM_TYPE_DIRECTIVE(one,@object)
++one: .double 1.0
++ ASM_SIZE_DIRECTIVE(one)
++ ASM_TYPE_DIRECTIVE(limit,@object)
++limit: .double 0.29
++ ASM_SIZE_DIRECTIVE(limit)
++ ASM_TYPE_DIRECTIVE(p63,@object)
++p63: .byte 0, 0, 0, 0, 0, 0, 0xe0, 0x43
++ ASM_SIZE_DIRECTIVE(p63)
++
++#ifdef PIC
++#define MO(op) op##@GOTOFF(%ecx)
++#define MOX(op,x,f) op##@GOTOFF(%ecx,x,f)
++#else
++#define MO(op) op
++#define MOX(op,x,f) op(,x,f)
++#endif
++
++ .text
++ENTRY(__ieee754_powl)
++ fldt 16(%esp) // y
++ fxam
++
++#ifdef PIC
++ call 1f
++1: popl %ecx
++ addl $_GLOBAL_OFFSET_TABLE_+[.-1b], %ecx
++#endif
++
++ fnstsw
++ movb %ah, %dl
++ andb $0x45, %ah
++ cmpb $0x40, %ah // is y == 0 ?
++ je 11f
++
++ cmpb $0x05, %ah // is y == ±inf ?
++ je 12f
++
++ cmpb $0x01, %ah // is y == NaN ?
++ je 30f
++
++ fldt 4(%esp) // x : y
++
++ subl $8,%esp
++
++ fxam
++ fnstsw
++ movb %ah, %dh
++ andb $0x45, %ah
++ cmpb $0x40, %ah
++ je 20f // x is ±0
++
++ cmpb $0x05, %ah
++ je 15f // x is ±inf
++
++ fxch // y : x
++
++ /* fistpll raises invalid exception for |y| >= 1L<<63. */
++ fld %st // y : y : x
++ fabs // |y| : y : x
++ fcompl MO(p63) // y : x
++ fnstsw
++ sahf
++ jnc 2f
++
++ /* First see whether `y' is a natural number. In this case we
++ can use a more precise algorithm. */
++ fld %st // y : y : x
++ fistpll (%esp) // y : x
++ fildll (%esp) // int(y) : y : x
++ fucomp %st(1) // y : x
++ fnstsw
++ sahf
++ jne 2f
++
++ /* OK, we have an integer value for y. */
++ popl %eax
++ popl %edx
++ orl $0, %edx
++ fstp %st(0) // x
++ jns 4f // y >= 0, jump
++ fdivrl MO(one) // 1/x (now referred to as x)
++ negl %eax
++ adcl $0, %edx
++ negl %edx
++4: fldl MO(one) // 1 : x
++ fxch
++
++6: shrdl $1, %edx, %eax
++ jnc 5f
++ fxch
++ fmul %st(1) // x : ST*x
++ fxch
++5: fmul %st(0), %st // x*x : ST*x
++ shrl $1, %edx
++ movl %eax, %ecx
++ orl %edx, %ecx
++ jnz 6b
++ fstp %st(0) // ST*x
++ ret
++
++ /* y is ±NAN */
++30: fldt 4(%esp) // x : y
++ fldl MO(one) // 1.0 : x : y
++ fucomp %st(1) // x : y
++ fnstsw
++ sahf
++ je 31f
++ fxch // y : x
++31: fstp %st(1)
++ ret
++
++ .align ALIGNARG(4)
++2: /* y is a real number. */
++ fxch // x : y
++ fldl MO(one) // 1.0 : x : y
++ fld %st(1) // x : 1.0 : x : y
++ fsub %st(1) // x-1 : 1.0 : x : y
++ fabs // |x-1| : 1.0 : x : y
++ fcompl MO(limit) // 1.0 : x : y
++ fnstsw
++ fxch // x : 1.0 : y
++ sahf
++ ja 7f
++ fsub %st(1) // x-1 : 1.0 : y
++ fyl2xp1 // log2(x) : y
++ jmp 8f
++
++7: fyl2x // log2(x) : y
++8: fmul %st(1) // y*log2(x) : y
++ fxam
++ fnstsw
++ andb $0x45, %ah
++ cmpb $0x05, %ah // is y*log2(x) == ±inf ?
++ je 28f
++ fst %st(1) // y*log2(x) : y*log2(x)
++ frndint // int(y*log2(x)) : y*log2(x)
++ fsubr %st, %st(1) // int(y*log2(x)) : fract(y*log2(x))
++ fxch // fract(y*log2(x)) : int(y*log2(x))
++ f2xm1 // 2^fract(y*log2(x))-1 : int(y*log2(x))
++ faddl MO(one) // 2^fract(y*log2(x)) : int(y*log2(x))
++ fscale // 2^fract(y*log2(x))*2^int(y*log2(x)) : int(y*log2(x))
++ addl $8, %esp
++ fstp %st(1) // 2^fract(y*log2(x))*2^int(y*log2(x))
++ ret
++
++28: fstp %st(1) // y*log2(x)
++ fldl MO(one) // 1 : y*log2(x)
++ fscale // 2^(y*log2(x)) : y*log2(x)
++ addl $8, %esp
++ fstp %st(1) // 2^(y*log2(x))
++ ret
++
++ // pow(x,±0) = 1
++ .align ALIGNARG(4)
++11: fstp %st(0) // pop y
++ fldl MO(one)
++ ret
++
++ // y == ±inf
++ .align ALIGNARG(4)
++12: fstp %st(0) // pop y
++ fldt 4(%esp) // x
++ fabs
++ fcompl MO(one) // < 1, == 1, or > 1
++ fnstsw
++ andb $0x45, %ah
++ cmpb $0x45, %ah
++ je 13f // jump if x is NaN
++
++ cmpb $0x40, %ah
++ je 14f // jump if |x| == 1
++
++ shlb $1, %ah
++ xorb %ah, %dl
++ andl $2, %edx
++ fldl MOX(inf_zero, %edx, 4)
++ ret
++
++ .align ALIGNARG(4)
++14: fldl MO(one)
++ ret
++
++ .align ALIGNARG(4)
++13: fldt 4(%esp) // load x == NaN
++ ret
++
++ .align ALIGNARG(4)
++ // x is ±inf
++15: fstp %st(0) // y
++ testb $2, %dh
++ jz 16f // jump if x == +inf
++
++ // We must find out whether y is an odd integer.
++ fld %st // y : y
++ fistpll (%esp) // y
++ fildll (%esp) // int(y) : y
++ fucompp // <empty>
++ fnstsw
++ sahf
++ jne 17f
++
++ // OK, the value is an integer, but is it odd?
++ popl %eax
++ popl %edx
++ andb $1, %al
++ jz 18f // jump if not odd
++ // It's an odd integer.
++ shrl $31, %edx
++ fldl MOX(minf_mzero, %edx, 8)
++ ret
++
++ .align ALIGNARG(4)
++16: fcompl MO(zero)
++ addl $8, %esp
++ fnstsw
++ shrl $5, %eax
++ andl $8, %eax
++ fldl MOX(inf_zero, %eax, 1)
++ ret
++
++ .align ALIGNARG(4)
++17: shll $30, %edx // sign bit for y in right position
++ addl $8, %esp
++18: shrl $31, %edx
++ fldl MOX(inf_zero, %edx, 8)
++ ret
++
++ .align ALIGNARG(4)
++ // x is ±0
++20: fstp %st(0) // y
++ testb $2, %dl
++ jz 21f // y > 0
++
++ // x is ±0 and y is < 0. We must find out whether y is an odd integer.
++ testb $2, %dh
++ jz 25f
++
++ fld %st // y : y
++ fistpll (%esp) // y
++ fildll (%esp) // int(y) : y
++ fucompp // <empty>
++ fnstsw
++ sahf
++ jne 26f
++
++ // OK, the value is an integer, but is it odd?
++ popl %eax
++ popl %edx
++ andb $1, %al
++ jz 27f // jump if not odd
++ // It's an odd integer.
++ // Raise divide-by-zero exception and get minus infinity value.
++ fldl MO(one)
++ fdivl MO(zero)
++ fchs
++ ret
++
++25: fstp %st(0)
++26: addl $8, %esp
++27: // Raise divide-by-zero exception and get infinity value.
++ fldl MO(one)
++ fdivl MO(zero)
++ ret
++
++ .align ALIGNARG(4)
++ // x is ±0 and y is > 0. We must find out whether y is an odd integer.
++21: testb $2, %dh
++ jz 22f
++
++ fld %st // y : y
++ fistpll (%esp) // y
++ fildll (%esp) // int(y) : y
++ fucompp // <empty>
++ fnstsw
++ sahf
++ jne 23f
++
++ // OK, the value is an integer, but is it odd?
++ popl %eax
++ popl %edx
++ andb $1, %al
++ jz 24f // jump if not odd
++ // It's an odd integer.
++ fldl MO(mzero)
++ ret
++
++22: fstp %st(0)
++23: addl $8, %esp // Don't use 2 x pop
++24: fldl MO(zero)
++ ret
++
++END(__ieee754_powl)
++
++weak_alias(__ieee754_powl,powl)
+diff -Nru qemu-0.7.0/bsd.orig/i386/e_remainderl.S qemu-0.7.0/bsd/i386/e_remainderl.S
+--- qemu-0.7.0/bsd.orig/i386/e_remainderl.S Wed Dec 31 19:00:00 1969
++++ qemu-0.7.0/bsd/i386/e_remainderl.S Fri Apr 29 02:11:27 2005
+@@ -0,0 +1,24 @@
++/*
++ * Written by J.T. Conklin <jtc@netbsd.org>.
++ * Public domain.
++ *
++ * Adapted for `long double' by Ulrich Drepper <drepper@cygnus.com>.
++ */
++
++#include <machine/asm.h>
++#include <sysdep.h>
++
++RCSID("$NetBSD: $")
++
++ENTRY(__ieee754_remainderl)
++ fldt 16(%esp)
++ fldt 4(%esp)
++1: fprem1
++ fstsw %ax
++ sahf
++ jp 1b
++ fstp %st(1)
++ ret
++END (__ieee754_remainderl)
++
++weak_alias(__ieee754_remainderl,remainderl)
+diff -Nru qemu-0.7.0/bsd.orig/i386/e_sqrtl.c qemu-0.7.0/bsd/i386/e_sqrtl.c
+--- qemu-0.7.0/bsd.orig/i386/e_sqrtl.c Wed Dec 31 19:00:00 1969
++++ qemu-0.7.0/bsd/i386/e_sqrtl.c Fri Apr 29 02:11:27 2005
+@@ -0,0 +1,20 @@
++/*
++ * Written by J.T. Conklin <jtc@netbsd.org>.
++ * Public domain.
++ *
++ * Adapted for `long double' by Ulrich Drepper <drepper@cygnus.com>.
++ */
++
++#include <sysdep.h>
++
++long double
++__ieee754_sqrtl (long double x)
++{
++ long double res;
++
++ asm ("fsqrt" : "=t" (res) : "0" (x));
++
++ return res;
++}
++
++weak_alias(__ieee754_sqrtl,sqrtl)
+diff -Nru qemu-0.7.0/bsd.orig/i386/s_ceill.S qemu-0.7.0/bsd/i386/s_ceill.S
+--- qemu-0.7.0/bsd.orig/i386/s_ceill.S Wed Dec 31 19:00:00 1969
++++ qemu-0.7.0/bsd/i386/s_ceill.S Fri Apr 29 02:11:27 2005
+@@ -0,0 +1,27 @@
++/*
++ * Based on code written by J.T. Conklin <jtc@NetBSD.org>.
++ * Public domain.
++ */
++
++#include <machine/asm.h>
++RCSID("$FreeBSD: /tmp/pcvs/ports/emulators/qemu/files/patch-libmath,v 1.1 2005-05-03 04:02:46 nork Exp $")
++
++ENTRY(ceill)
++ pushl %ebp
++ movl %esp,%ebp
++ subl $8,%esp
++
++ fstcw -4(%ebp) /* store fpu control word */
++ movw -4(%ebp),%dx
++ orw $0x0800,%dx /* round towards +oo */
++ andw $0xfbff,%dx
++ movw %dx,-8(%ebp)
++ fldcw -8(%ebp) /* load modfied control word */
++
++ fldt 8(%ebp) /* round */
++ frndint
++
++ fldcw -4(%ebp) /* restore original control word */
++
++ leave
++ ret
+diff -Nru qemu-0.7.0/bsd.orig/i386/s_cosl.S qemu-0.7.0/bsd/i386/s_cosl.S
+--- qemu-0.7.0/bsd.orig/i386/s_cosl.S Wed Dec 31 19:00:00 1969
++++ qemu-0.7.0/bsd/i386/s_cosl.S Fri Apr 29 02:11:27 2005
+@@ -0,0 +1,32 @@
++/*
++ * Written by J.T. Conklin <jtc@netbsd.org>.
++ * Public domain.
++ *
++ * Adapted for `long double' by Ulrich Drepper <drepper@cygnus.com>.
++ */
++
++#include <machine/asm.h>
++#include <sysdep.h>
++
++RCSID("$NetBSD: $")
++
++ENTRY(__cosl)
++ fldt 4(%esp)
++ fcos
++ fnstsw %ax
++ testl $0x400,%eax
++ jnz 1f
++ ret
++ .align ALIGNARG(4)
++1: fldpi
++ fadd %st(0)
++ fxch %st(1)
++2: fprem1
++ fnstsw %ax
++ testl $0x400,%eax
++ jnz 2b
++ fstp %st(1)
++ fcos
++ ret
++END (__cosl)
++weak_alias (__cosl, cosl)
+diff -Nru qemu-0.7.0/bsd.orig/i386/s_floorl.S qemu-0.7.0/bsd/i386/s_floorl.S
+--- qemu-0.7.0/bsd.orig/i386/s_floorl.S Wed Dec 31 19:00:00 1969
++++ qemu-0.7.0/bsd/i386/s_floorl.S Fri Apr 29 02:11:27 2005
+@@ -0,0 +1,27 @@
++/*
++ * Based on code written by J.T. Conklin <jtc@NetBSD.org>.
++ * Public domain.
++ */
++
++#include <machine/asm.h>
++RCSID("$FreeBSD: /tmp/pcvs/ports/emulators/qemu/files/patch-libmath,v 1.1 2005-05-03 04:02:46 nork Exp $")
++
++ENTRY(floorl)
++ pushl %ebp
++ movl %esp,%ebp
++ subl $8,%esp
++
++ fstcw -4(%ebp) /* store fpu control word */
++ movw -4(%ebp),%dx
++ orw $0x0400,%dx /* round towards -oo */
++ andw $0xf7ff,%dx
++ movw %dx,-8(%ebp)
++ fldcw -8(%ebp) /* load modfied control word */
++
++ fldt 8(%ebp) /* round */
++ frndint
++
++ fldcw -4(%ebp) /* restore original control word */
++
++ leave
++ ret
+diff -Nru qemu-0.7.0/bsd.orig/i386/s_isnormal.c qemu-0.7.0/bsd/i386/s_isnormal.c
+--- qemu-0.7.0/bsd.orig/i386/s_isnormal.c Wed Dec 31 19:00:00 1969
++++ qemu-0.7.0/bsd/i386/s_isnormal.c Fri Apr 29 02:11:27 2005
+@@ -0,0 +1,85 @@
++/*-
++ * Copyright (c) 2003 Mike Barcroft <mike@FreeBSD.org>
++ * Copyright (c) 2002-2004 David Schultz <das@FreeBSD.ORG>
++ * 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, this list of conditions and the following disclaimer.
++ * 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.
++ *
++ * THIS SOFTWARE IS PROVIDED BY THE AUTHOR AND CONTRIBUTORS ``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 OR CONTRIBUTORS 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.
++ *
++ * $FreeBSD: /tmp/pcvs/ports/emulators/qemu/files/patch-libmath,v 1.1 2005-05-03 04:02:46 nork Exp $
++ */
++
++union IEEEf2bits {
++ float f;
++ struct {
++ unsigned int man :23;
++ unsigned int exp :8;
++ unsigned int sign :1;
++ } bits;
++};
++
++union IEEEd2bits {
++ double d;
++ struct {
++ unsigned int manl :32;
++ unsigned int manh :20;
++ unsigned int exp :11;
++ unsigned int sign :1;
++ } bits;
++};
++
++union IEEEl2bits {
++ long double e;
++ struct {
++ unsigned int manl :32;
++ unsigned int manh :32;
++ unsigned int exp :15;
++ unsigned int sign :1;
++ unsigned int junk :16;
++ } bits;
++};
++
++int
++__isnormal(double d)
++{
++ union IEEEd2bits u;
++
++ u.d = d;
++ return (u.bits.exp != 0 && u.bits.exp != 2047);
++}
++
++int
++__isnormalf(float f)
++{
++ union IEEEf2bits u;
++
++ u.f = f;
++ return (u.bits.exp != 0 && u.bits.exp != 255);
++}
++
++int
++__isnormall(long double e)
++{
++ union IEEEl2bits u;
++
++ u.e = e;
++ return (u.bits.exp != 0 && u.bits.exp != 32767);
++}
+diff -Nru qemu-0.7.0/bsd.orig/i386/s_llrint.S qemu-0.7.0/bsd/i386/s_llrint.S
+--- qemu-0.7.0/bsd.orig/i386/s_llrint.S Wed Dec 31 19:00:00 1969
++++ qemu-0.7.0/bsd/i386/s_llrint.S Fri Apr 29 02:11:27 2005
+@@ -0,0 +1,36 @@
++/*-
++ * Copyright (c) 2005 David Schultz <das@FreeBSD.ORG>
++ * 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, this list of conditions and the following disclaimer.
++ * 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.
++ *
++ * THIS SOFTWARE IS PROVIDED BY THE AUTHOR AND CONTRIBUTORS ``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 OR CONTRIBUTORS 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 <machine/asm.h>
++RCSID("$FreeBSD: /tmp/pcvs/ports/emulators/qemu/files/patch-libmath,v 1.1 2005-05-03 04:02:46 nork Exp $");
++
++ENTRY(llrint)
++ fldl 4(%esp)
++ subl $8,%esp
++ fistpll (%esp)
++ popl %eax
++ popl %edx
++ ret
+diff -Nru qemu-0.7.0/bsd.orig/i386/s_llrintf.S qemu-0.7.0/bsd/i386/s_llrintf.S
+--- qemu-0.7.0/bsd.orig/i386/s_llrintf.S Wed Dec 31 19:00:00 1969
++++ qemu-0.7.0/bsd/i386/s_llrintf.S Fri Apr 29 02:11:27 2005
+@@ -0,0 +1,36 @@
++/*-
++ * Copyright (c) 2005 David Schultz <das@FreeBSD.ORG>
++ * 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, this list of conditions and the following disclaimer.
++ * 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.
++ *
++ * THIS SOFTWARE IS PROVIDED BY THE AUTHOR AND CONTRIBUTORS ``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 OR CONTRIBUTORS 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 <machine/asm.h>
++RCSID("$FreeBSD: /tmp/pcvs/ports/emulators/qemu/files/patch-libmath,v 1.1 2005-05-03 04:02:46 nork Exp $")
++
++ENTRY(llrintf)
++ flds 4(%esp)
++ subl $8,%esp
++ fistpll (%esp)
++ popl %eax
++ popl %edx
++ ret
+diff -Nru qemu-0.7.0/bsd.orig/i386/s_llrintl.S qemu-0.7.0/bsd/i386/s_llrintl.S
+--- qemu-0.7.0/bsd.orig/i386/s_llrintl.S Wed Dec 31 19:00:00 1969
++++ qemu-0.7.0/bsd/i386/s_llrintl.S Fri Apr 29 02:11:27 2005
+@@ -0,0 +1,35 @@
++/* Round argument to nearest integral value according to current rounding
++ direction.
++ Copyright (C) 1997 Free Software Foundation, Inc.
++ This file is part of the GNU C Library.
++ Contributed by Ulrich Drepper <drepper@cygnus.com>, 1997.
++
++ The GNU C Library is free software; you can redistribute it and/or
++ modify it under the terms of the GNU Lesser General Public
++ License as published by the Free Software Foundation; either
++ version 2.1 of the License, or (at your option) any later version.
++
++ The GNU C Library is distributed in the hope that it will be useful,
++ but WITHOUT ANY WARRANTY; without even the implied warranty of
++ MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
++ Lesser General Public License for more details.
++
++ You should have received a copy of the GNU Lesser General Public
++ License along with the GNU C Library; if not, write to the Free
++ Software Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA
++ 02111-1307 USA. */
++
++#include <machine/asm.h>
++#include <sysdep.h>
++
++ .text
++ENTRY(__llrintl)
++ fldt 4(%esp)
++ subl $8, %esp
++ fistpll (%esp)
++ fwait
++ popl %eax
++ popl %edx
++ ret
++END(__llrintl)
++weak_alias (__llrintl, llrintl)
+diff -Nru qemu-0.7.0/bsd.orig/i386/s_lrint.S qemu-0.7.0/bsd/i386/s_lrint.S
+--- qemu-0.7.0/bsd.orig/i386/s_lrint.S Wed Dec 31 19:00:00 1969
++++ qemu-0.7.0/bsd/i386/s_lrint.S Fri Apr 29 02:11:27 2005
+@@ -0,0 +1,35 @@
++/*-
++ * Copyright (c) 2005 David Schultz <das@FreeBSD.ORG>
++ * 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, this list of conditions and the following disclaimer.
++ * 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.
++ *
++ * THIS SOFTWARE IS PROVIDED BY THE AUTHOR AND CONTRIBUTORS ``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 OR CONTRIBUTORS 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 <machine/asm.h>
++RCSID("$FreeBSD: /tmp/pcvs/ports/emulators/qemu/files/patch-libmath,v 1.1 2005-05-03 04:02:46 nork Exp $");
++
++ENTRY(lrint)
++ fldl 4(%esp)
++ subl $4,%esp
++ fistpl (%esp)
++ popl %eax
++ ret
+diff -Nru qemu-0.7.0/bsd.orig/i386/s_lrintf.S qemu-0.7.0/bsd/i386/s_lrintf.S
+--- qemu-0.7.0/bsd.orig/i386/s_lrintf.S Wed Dec 31 19:00:00 1969
++++ qemu-0.7.0/bsd/i386/s_lrintf.S Fri Apr 29 02:11:27 2005
+@@ -0,0 +1,35 @@
++/*-
++ * Copyright (c) 2005 David Schultz <das@FreeBSD.ORG>
++ * 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, this list of conditions and the following disclaimer.
++ * 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.
++ *
++ * THIS SOFTWARE IS PROVIDED BY THE AUTHOR AND CONTRIBUTORS ``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 OR CONTRIBUTORS 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 <machine/asm.h>
++RCSID("$FreeBSD: /tmp/pcvs/ports/emulators/qemu/files/patch-libmath,v 1.1 2005-05-03 04:02:46 nork Exp $")
++
++ENTRY(lrintf)
++ flds 4(%esp)
++ subl $4,%esp
++ fistpl (%esp)
++ popl %eax
++ ret
+diff -Nru qemu-0.7.0/bsd.orig/i386/s_lrintl.S qemu-0.7.0/bsd/i386/s_lrintl.S
+--- qemu-0.7.0/bsd.orig/i386/s_lrintl.S Wed Dec 31 19:00:00 1969
++++ qemu-0.7.0/bsd/i386/s_lrintl.S Fri Apr 29 02:11:27 2005
+@@ -0,0 +1,34 @@
++/* Round argument to nearest integral value according to current rounding
++ direction.
++ Copyright (C) 1997 Free Software Foundation, Inc.
++ This file is part of the GNU C Library.
++ Contributed by Ulrich Drepper <drepper@cygnus.com>, 1997.
++
++ The GNU C Library is free software; you can redistribute it and/or
++ modify it under the terms of the GNU Lesser General Public
++ License as published by the Free Software Foundation; either
++ version 2.1 of the License, or (at your option) any later version.
++
++ The GNU C Library is distributed in the hope that it will be useful,
++ but WITHOUT ANY WARRANTY; without even the implied warranty of
++ MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
++ Lesser General Public License for more details.
++
++ You should have received a copy of the GNU Lesser General Public
++ License along with the GNU C Library; if not, write to the Free
++ Software Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA
++ 02111-1307 USA. */
++
++#include <machine/asm.h>
++#include <sysdep.h>
++
++ .text
++ENTRY(__lrintl)
++ fldt 4(%esp)
++ subl $4, %esp
++ fistpl (%esp)
++ fwait
++ popl %eax
++ ret
++END(__lrintl)
++weak_alias (__lrintl, lrintl)
+diff -Nru qemu-0.7.0/bsd.orig/i386/s_rintl.c qemu-0.7.0/bsd/i386/s_rintl.c
+--- qemu-0.7.0/bsd.orig/i386/s_rintl.c Wed Dec 31 19:00:00 1969
++++ qemu-0.7.0/bsd/i386/s_rintl.c Fri Apr 29 02:11:27 2005
+@@ -0,0 +1,18 @@
++/*
++ * Written by J.T. Conklin <jtc@netbsd.org>.
++ * Changes for long double by Ulrich Drepper <drepper@cygnus.com>
++ * Public domain.
++ */
++
++#include <sysdep.h>
++
++long double
++__rintl (long double x)
++{
++ long double res;
++
++ asm ("frndint" : "=t" (res) : "0" (x));
++ return res;
++}
++
++weak_alias (__rintl, rintl)
+diff -Nru qemu-0.7.0/bsd.orig/i386/s_round.c qemu-0.7.0/bsd/i386/s_round.c
+--- qemu-0.7.0/bsd.orig/i386/s_round.c Wed Dec 31 19:00:00 1969
++++ qemu-0.7.0/bsd/i386/s_round.c Fri Apr 29 02:11:27 2005
+@@ -0,0 +1,51 @@
++/*-
++ * Copyright (c) 2003, Steven G. Kargl
++ * 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 unmodified, this list of conditions, and the following
++ * disclaimer.
++ * 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.
++ *
++ * THIS SOFTWARE IS PROVIDED BY THE AUTHOR ``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 <sys/cdefs.h>
++__FBSDID("$FreeBSD: /tmp/pcvs/ports/emulators/qemu/files/patch-libmath,v 1.1 2005-05-03 04:02:46 nork Exp $");
++
++#include <math.h>
++
++double
++round(double x)
++{
++ double t;
++
++ if (!isfinite(x))
++ return (x);
++
++ if (x >= 0.0) {
++ t = ceil(x);
++ if (t - x > 0.5)
++ t -= 1.0;
++ return (t);
++ } else {
++ t = ceil(-x);
++ if (t + x > 0.5)
++ t -= 1.0;
++ return (-t);
++ }
++}
+diff -Nru qemu-0.7.0/bsd.orig/i386/s_sinl.S qemu-0.7.0/bsd/i386/s_sinl.S
+--- qemu-0.7.0/bsd.orig/i386/s_sinl.S Wed Dec 31 19:00:00 1969
++++ qemu-0.7.0/bsd/i386/s_sinl.S Fri Apr 29 02:11:27 2005
+@@ -0,0 +1,32 @@
++/*
++ * Written by J.T. Conklin <jtc@netbsd.org>.
++ * Public domain.
++ *
++ * Adapted for `long double' by Ulrich Drepper <drepper@cygnus.com>.
++ */
++
++#include <machine/asm.h>
++#include <sysdep.h>
++
++RCSID("$NetBSD: $")
++
++ENTRY(__sinl)
++ fldt 4(%esp)
++ fsin
++ fnstsw %ax
++ testl $0x400,%eax
++ jnz 1f
++ ret
++ .align ALIGNARG(4)
++1: fldpi
++ fadd %st(0)
++ fxch %st(1)
++2: fprem1
++ fnstsw %ax
++ testl $0x400,%eax
++ jnz 2b
++ fstp %st(1)
++ fsin
++ ret
++END (__sinl)
++weak_alias (__sinl, sinl)
+diff -Nru qemu-0.7.0/bsd.orig/i386/s_tanl.S qemu-0.7.0/bsd/i386/s_tanl.S
+--- qemu-0.7.0/bsd.orig/i386/s_tanl.S Wed Dec 31 19:00:00 1969
++++ qemu-0.7.0/bsd/i386/s_tanl.S Fri Apr 29 02:11:27 2005
+@@ -0,0 +1,33 @@
++/*
++ * Written by J.T. Conklin <jtc@netbsd.org>.
++ * Public domain.
++ *
++ * Adapted for `long double' by Ulrich Drepper <drepper@cygnus.com>.
++ */
++
++#include <machine/asm.h>
++#include <sysdep.h>
++
++RCSID("$NetBSD: $")
++
++ENTRY(__tanl)
++ fldt 4(%esp)
++ fptan
++ fnstsw %ax
++ testl $0x400,%eax
++ jnz 1f
++ fstp %st(0)
++ ret
++1: fldpi
++ fadd %st(0)
++ fxch %st(1)
++2: fprem1
++ fstsw %ax
++ testl $0x400,%eax
++ jnz 2b
++ fstp %st(1)
++ fptan
++ fstp %st(0)
++ ret
++END (__tanl)
++weak_alias (__tanl, tanl)
+diff -Nru qemu-0.7.0/bsd.orig/sysdep.h qemu-0.7.0/bsd/sysdep.h
+--- qemu-0.7.0/bsd.orig/sysdep.h Wed Dec 31 19:00:00 1969
++++ qemu-0.7.0/bsd/sysdep.h Fri Apr 29 02:11:27 2005
+@@ -0,0 +1,20 @@
++#ifndef _QEMU_BSD_SYSDEP_H_
++#define _QEMU_BSD_SYSDEP_H_
++
++#include <sys/cdefs.h>
++
++#define HAVE_ELF
++
++#ifdef __ASSEMBLER__
++#define ALIGNARG(log2) 1<<log2
++#define ASM_TYPE_DIRECTIVE(name,typearg) .type name,typearg;
++#define ASM_SIZE_DIRECTIVE(name) .size name,.-name;
++#define END(x)
++#define strong_alias(sym,alias) .set alias,sym;
++#define weak_alias(sym,alias) .weak alias; .equ alias,sym;
++#else
++#define strong_alias(sym,alias) __strong_reference(sym,alias);
++#define weak_alias(sym,alias) __weak_reference(sym,alias);
++#endif
++
++#endif