From 87a27b3337a64d58b3a39fdc7f9e02a460bddcbc Mon Sep 17 00:00:00 2001
From: Ulrich Drepper <drepper@redhat.com>
Date: Sun, 13 Feb 2000 19:40:31 +0000
Subject: Update.

	* sysdeps/i386/fpu/s_cosf.S: Domain of opcode is not large enough
	so test for overflow and handle it.
	* sysdeps/i386/fpu/s_sinf.S: Likewise.
	* sysdeps/i386/fpu/s_tanf.S: Likewise.
	Patch by Miloslav Trmac <mitr@volny.cz> (PR libc/1563).
---
 sysdeps/i386/fpu/s_cosf.S | 15 ++++++++++++++-
 sysdeps/i386/fpu/s_sinf.S | 15 ++++++++++++++-
 sysdeps/i386/fpu/s_tanf.S | 15 ++++++++++++++-
 3 files changed, 42 insertions(+), 3 deletions(-)

(limited to 'sysdeps/i386/fpu')

diff --git a/sysdeps/i386/fpu/s_cosf.S b/sysdeps/i386/fpu/s_cosf.S
index d8e8090..21f87aa 100644
--- a/sysdeps/i386/fpu/s_cosf.S
+++ b/sysdeps/i386/fpu/s_cosf.S
@@ -7,10 +7,23 @@
 
 RCSID("$NetBSD: s_cosf.S,v 1.3 1995/05/08 23:55:16 jtc Exp $")
 
-/* A float's domain isn't large enough to require argument reduction. */
 ENTRY(__cosf)
 	flds	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 (__cosf)
 weak_alias (__cosf, cosf)
diff --git a/sysdeps/i386/fpu/s_sinf.S b/sysdeps/i386/fpu/s_sinf.S
index a010d60..5ca45f5 100644
--- a/sysdeps/i386/fpu/s_sinf.S
+++ b/sysdeps/i386/fpu/s_sinf.S
@@ -7,10 +7,23 @@
 
 RCSID("$NetBSD: s_sinf.S,v 1.3 1995/05/09 00:27:53 jtc Exp $")
 
-/* A float's domain isn't large enough to require argument reduction. */
 ENTRY(__sinf)
 	flds	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 (__sinf)
 weak_alias (__sinf, sinf)
diff --git a/sysdeps/i386/fpu/s_tanf.S b/sysdeps/i386/fpu/s_tanf.S
index 7a75091..355dff9 100644
--- a/sysdeps/i386/fpu/s_tanf.S
+++ b/sysdeps/i386/fpu/s_tanf.S
@@ -7,10 +7,23 @@
 
 RCSID("$NetBSD: s_tanf.S,v 1.3 1995/05/09 00:31:09 jtc Exp $")
 
-/* A float's domain isn't large enough to require argument reduction. */
 ENTRY(__tanf)
 	flds	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 (__tanf)
-- 
cgit v1.1