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] |
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] |