This is the mail archive of the libc-hacker@sources.redhat.com mailing list for the glibc project.

Note that libc-hacker is a closed list. You may look at the archives of this list, but subscription and posting are not open.


Index Nav: [Date Index] [Subject Index] [Author Index] [Thread Index]
Message Nav: [Date Prev] [Date Next] [Thread Prev] [Thread Next]

First patch for i386/fpu function



Here's a small patch for some easy functions.

Ok to commit?

Andreas

2001-05-06  Andreas Jaeger  <aj@suse.de>

	* sysdeps/i386/fpu/e_acosl.c: New, rewrite of e_acos.S.
	* sysdeps/i386/fpu/s_atanl.c: New, rewrite of e_atanl.S.
	* sysdeps/i386/fpu/e_sqrtl.c: New, rewrite of e_sqrt.S.
	* sysdeps/i386/fpu/e_acosl.S: Removed.
	* sysdeps/i386/fpu/s_atanl.S: Removed.
	* sysdeps/i386/fpu/e_sqrtl.S: Removed.

============================================================
Index: sysdeps/i386/fpu/e_acosl.c
--- sysdeps/i386/fpu/e_acosl.c	created
+++ sysdeps/i386/fpu/e_acosl.c	Sun May  6 13:39:44 2001	1.1
@@ -0,0 +1,25 @@
+/*
+ * Written by J.T. Conklin <jtc@netbsd.org>.
+ * Public domain.
+ *
+ * Adapted for `long double' by Ulrich Drepper <drepper@cygnus.com>.
+ */
+
+#include <math_private.h>
+
+long double
+__ieee754_acosl (long double x)
+{
+  long double res;
+
+  /* acosl = atanl (sqrtl(1 - x^2) / x) */
+  asm (	"fld	%%st\n"
+	"fmul	%%st(0)\n"		/* x^2 */
+	"fld1\n"
+	"fsubp\n"			/* 1 - x^2 */
+	"fsqrt\n"			/* sqrtl (1 - x^2) */
+	"fxch	%%st(1)\n"
+	"fpatan"
+	: "=t" (res) : "0" (x) : "st(1)");
+  return res;
+}
============================================================
Index: sysdeps/i386/fpu/s_atanl.c
--- sysdeps/i386/fpu/s_atanl.c	created
+++ sysdeps/i386/fpu/s_atanl.c	Sun May  6 13:36:01 2001	1.1
@@ -0,0 +1,22 @@
+/*
+ * Written by J.T. Conklin <jtc@netbsd.org>.
+ * Public domain.
+ *
+ * Adapted for `long double' by Ulrich Drepper <drepper@cygnus.com>.
+ */
+
+#include <math_private.h>
+
+long double
+__atanl (long double x)
+{
+  long double res;
+
+  asm ("fld1\n"
+       "fpatan"
+       : "=t" (res) : "0" (x));
+
+  return res;
+}
+
+weak_alias (__atanl, atanl)
============================================================
Index: sysdeps/i386/fpu/e_sqrtl.c
--- sysdeps/i386/fpu/e_sqrtl.c	created
+++ sysdeps/i386/fpu/e_sqrtl.c	Sun May  6 13:35:46 2001	1.1
@@ -0,0 +1,18 @@
+/*
+ * Written by J.T. Conklin <jtc@netbsd.org>.
+ * Public domain.
+ *
+ * Adapted for `long double' by Ulrich Drepper <drepper@cygnus.com>.
+ */
+
+#include <math_private.h>
+
+long double
+__ieee754_sqrtl (long double x)
+{
+  long double res;
+
+  asm ("fsqrt" : "=t" (res) : "0" (x));
+
+  return res;
+}

============================================================
Index: sysdeps/i386/fpu/e_acosl.S
--- sysdeps/i386/fpu/e_acosl.S	Sun May  6 13:42:34 2001	1.1
+++ sysdeps/i386/fpu/e_acosl.S	removed
@@ -1,22 +0,0 @@
-/*
- * Written by J.T. Conklin <jtc@netbsd.org>.
- * Public domain.
- *
- * Adapted for `long double' by Ulrich Drepper <drepper@cygnus.com>.
- */
-
-#include <machine/asm.h>
-
-
-/* acosl = atanl (sqrtl(1 - x^2) / x) */
-ENTRY(__ieee754_acosl)
-	fldt	4(%esp)			/* x */
-	fld	%st
-	fmul	%st(0)			/* x^2 */
-	fld1
-	fsubp				/* 1 - x^2 */
-	fsqrt				/* sqrtl (1 - x^2) */
-	fxch	%st(1)
-	fpatan
-	ret
-END (__ieee754_acosl)
============================================================
Index: sysdeps/i386/fpu/s_atanl.S
--- sysdeps/i386/fpu/s_atanl.S	Sun May  6 13:42:34 2001	1.1
+++ sysdeps/i386/fpu/s_atanl.S	removed
@@ -1,18 +0,0 @@
-/*
- * Written by J.T. Conklin <jtc@netbsd.org>.
- * Public domain.
- *
- * Adapted for `long double' by Ulrich Drepper <drepper@cygnus.com>.
- */
-
-#include <machine/asm.h>
-
-RCSID("$NetBSD: $")
-
-ENTRY(__atanl)
-	fldt	4(%esp)
-	fld1
-	fpatan
-	ret
-END (__atanl)
-weak_alias (__atanl, atanl)
============================================================
Index: sysdeps/i386/fpu/e_sqrtl.S
--- sysdeps/i386/fpu/e_sqrtl.S	Sun May  6 13:42:34 2001	1.1
+++ sysdeps/i386/fpu/e_sqrtl.S	removed
@@ -1,16 +0,0 @@
-/*
- * Written by J.T. Conklin <jtc@netbsd.org>.
- * Public domain.
- *
- * Adapted for `long double' by Ulrich Drepper <drepper@cygnus.com>.
- */
-
-#include <machine/asm.h>
-
-RCSID("$NetBSD: $")
-
-ENTRY(__ieee754_sqrtl)
-	fldt	4(%esp)
-	fsqrt
-	ret
-END (__ieee754_sqrtl)

-- 
 Andreas Jaeger
  SuSE Labs aj@suse.de
   private aj@arthur.inka.de
    http://www.suse.de/~aj


Index Nav: [Date Index] [Subject Index] [Author Index] [Thread Index]
Message Nav: [Date Prev] [Date Next] [Thread Prev] [Thread Next]