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