mirror of
https://sourceware.org/git/glibc.git
synced 2024-10-01 03:27:48 +00:00
9a0a462ceb
1997-09-11 04:36 Ulrich Drepper <drepper@cygnus.com> * db2/db_int.h: Define __set_errno if not yet available. * db2/btree/bt_rec.c: Use __set_errno to set errno value. * db2/clib/getlong.c: Likewise. * db2/db185/db185.c: Likewise. * db2/db185/db185_int.h: Likewise. * db2/dbm/dbm.c: Likewise. * db2/lock/lock_deadlock.c: Likewise. * db2/log/log_archive.c: Likewise. * elf/dl-profile.c: Implement mcount function. * gmon/gmon.c: Use __profil not profil because of namespace pollution. * gmon/mcount.c: Remove BSD kernel code. Use compare&swap instruction if possible to change state variable. Optimize frompc folding. * gmon/sys/gmon.h (struct gmonparam): Change state field to long int. * sysdeps/i386/i486/atomicity.h: New file. * sysdeps/stub/atomicity.h: New file. * sysdeps/mach/hurd/profil.c: Define function as __profil and make profil weak alias. * sysdeps/posix/profil.c: Likewise. * string/bits/string2.h: New file. * include/bits/string2.h: New file. * string/Makefile (routines): Add mempcpy. (tests): Add inl-tester. Remove _D__NO_STRING_INLINES from CFLAGS-* variables. * sysdeps/generic/mempcpy.c: New file. * sysdeps/generic/memccpy.c: Undef function name to enable definition as macro. * sysdeps/generic/memchr.c: Likewise. * sysdeps/generic/memcmp.c: Likewise. * sysdeps/generic/memmem.c: Likewise. * sysdeps/generic/memmove.c: Likewise. * sysdeps/generic/strcat.c: Likewise. * sysdeps/generic/strchr.c: Likewise. * sysdeps/generic/strcmp.c: Likewise. * sysdeps/generic/strcpy.c: Likewise. * sysdeps/generic/strcspn.c: Likewise. * sysdeps/generic/strlen.c: Likewise. * sysdeps/generic/strncat.c: Likewise. * sysdeps/generic/strncmp.c: Likewise. * sysdeps/generic/strncpy.c: Likewise. * sysdeps/generic/strpbrk.c: Likewise. * sysdeps/generic/strrchr.c: Likewise. * sysdeps/generic/strsep.c: Likewise. * sysdeps/generic/strspn.c: Likewise. * sysdeps/generic/strstr.c: Likewise. * sysdeps/generic/strtok.c: Likewise. * sysdeps/generic/strtok_r.c: Likewise. * sysdeps/i386/memset.c: Likewise. * sysdeps/i386/bits/string.h: Correct a few types and constraints. * sysdeps/i386/i486/bits/string.h: Heavy rewrites and optimizations. * string/stratcliff.c: Undefine __USE_STRING_INLINES. * string/tst-strlen.c: Likewise. * string/string.h: Add prototype for mempcpy. Include bits/string2.h header always if optimizing. * intl/dcgettext.c: Don't unconditionally define stpcpy, only if not yet defined. * intl/l10nflist.c: Likewise. * string/tester.c: Add copyright and make little cleanups. * inet/test_ifindex.c: Change type of ni variable to unsigned int. * locale/programs/ld-ctype.c (struct locale_ctype_t): Change type of fields map_collection_max and map_collection_act to size_t. * nss/libnss_files.map: Group entries. * posix/unistd.h: Add prototype for __setpgid and __profil. * sysdeps/generic/crypt.h: Declare __crypt_r. * sysdeps/i386/bits/select.h: Fix fatal bugs, use correct casts now. * sysdeps/i386/fpu/bits/mathinline.h (isgreater, isgreaterequal, isless, islessequal, islessgreater, isunordered): Optimize a bit. * sysdeps/stub/ftruncate.c: Include missing header for prototype. * sysdeps/stub/getdents.c: Likewise. * sysdeps/stub/reboot.c: Likewise. * sysdeps/stub/swapon.c: Likewise. * sysdeps/stub/syscall.c: Likewise. * sysdeps/stub/ualarm.c: Likewise. * sysdeps/stub/usleep.c: Likewise. * sysdeps/unix/sysv/linux/if_index.c: Don't compile or use opensock if SIOGIFINDEX and SIOGIFNAME are not defined. * sysdeps/unix/sysv/linux/net/if.h: Add IFF_PORTSEL and IFF_AUTOMEDIA according to recent kernel changes. 1997-09-10 21:32 Klaus Espenlaub <kespenla@student.informatik.uni-ulm.de> * Makeconfig: Use $(have-initfini) instead of $(elf) to figure out the installed name of the startup code. (common-generated): Add version.mk. * Makefile (distclean-1): Add glibcbug. * Makerules: Replace -lgcc by $(gnulib). * catgets/Makefile (generated): Add xmalloc.o. * csu/Makefile (generated): Replace align.h and end.h by defs.h to match the generated file. * manual/Makefile (mostlyclean): Add stub-manual and stamp.o. (realclean): Changed to remove chapters-incl[12]. * po/Makefile (realclean): New rule to remove the generated .mo files. * time/Makefile: Only include zonefile dependencies if $(no_deps) is not true to avoid make clean failure when directory time doesn't exist yet. (generated): Add tzselect. * stdio/fgets.c (fgets): Add casts to reduce gcc warning noise. * stdio/internals.c (flushbuf): Likewise. * stdio/linewrap.c (lwupdate): Likewise. * stdio/memstream.c (enlarge_buffer): Likewise. * stdio-common/vfscanf.c (_IO_vfscanf): Likewise. * time/tzset.c (compute_change): Likewise. * misc/init-misc.c (__init_misc): Only declare static if HAVE_GNU_LD is defined. * sysdeps/posix/pipestream.c (FUNC): Change to generate ANSI C style functions. * sysdeps/stub/init-posix.c: Likewise. * sysdeps/stub/profil.c: Likewise. * munch-tmpl.c (__libc_init): Convert to ANSI C style declaration to reduce gcc warning noise. * stdio/glue.c (_filbuf, _flsbuf): Likewise. * stdio/obstream.c (grow, seek, input, init_obstream): Likewise. * stdio/vasprintf.c (enlarge_buffer): Likewise. * sysdeps/generic/sysd-stdio.c (__stdio_read, __stdio_write, __stdio_seek, __stdio_close, __stdio_fileno, __stdio_open, __stdio_reopen): Likewise. * sysdeps/posix/defs.c (_cleanup): Likewise. * time/offtime.c (__offtime): Add cast. * posix/getopt.c: Don't use text_set_element if not defined. * configure.in: Provide a check for underscores before user labels that works even when the compiler used for building doesn't work (like when there is no C library). Use the old way if the compiler works. 1997-09-10 05:08 David S. Miller <davem@caip.rutgers.edu> * sysdeps/unix/sysv/linux/sparc/bits/ioctls.h: The TC* ioctls use 'T' not 't' on SparcLinux. * sysdeps/unix/sysv/linux/sparc/bits/termios.h: tcflag_t is 32 bits. * sysdeps/unix/sysv/linux/sparc/sparc64/longjmp.S: Add aliases for _longjmp and siglongjmp. 1997-09-09 Andreas Schwab <schwab@issan.informatik.uni-dortmund.de> * libio/stdio.h: Add format attributes to the extra printf and scanf like functions. * stdio/stdio.h: Likewise. 1997-09-09 Andreas Schwab <schwab@issan.informatik.uni-dortmund.de> * nis/nss_nisplus/nisplus-network.c (_nss_nisplus_getnetbyaddr_r): Print tablename_val, not tablename_len. * nis/nss_nisplus/nisplus-ethers.c (_nss_nisplus_getntohost_r): Use sprintf instead of sprintf, the string always fits. * nis/nss_nisplus/nisplus-hosts.c (_nss_nisplus_gethostbyaddr_r): Likewise. * nis/nss_nisplus/nisplus-network.c (_nss_nisplus_getnetbyaddr_r): Likewise. * nis/nss_nisplus/nisplus-proto.c (_nss_nisplus_getprotobynumber_r): Likewise. * nis/nss_nisplus/nisplus-rpc.c (_nss_nisplus_getrpcbynumber_r): Likewise. * nis/nss_nisplus/nisplus-service.c (_nss_nisplus_getservbynumber_r): Likewise. * nis/nss_nisplus/nisplus-alias.c (_nss_create_tablename): Use __stpcpy, __stpncpy and __strdup instead of public names. * nis/nss_nisplus/nisplus-ethers.c (_nss_create_tablename): Likewise. * nis/nss_nisplus/nisplus-grp.c (_nss_create_tablename): Likewise. * nis/nss_nisplus/nisplus-hosts.c (_nss_create_tablename): Likewise. * nis/nss_nisplus/nisplus-netgrp.c (_nss_nisplus_parse_netgroup): Likewise. * nis/nss_nisplus/nisplus-network.c (_nss_nisplus_parse_netent): Likewise. (_nss_create_tablename): Likewise. * nis/nss_nisplus/nisplus-proto.c (_nss_nisplus_parse_protoent): Likewise. (_nss_create_tablename): Likewise. * nis/nss_nisplus/nisplus-pwd.c (_nss_create_tablename): Likewise. * nis/nss_nisplus/nisplus-rpc.c (_nss_nisplus_parse_rpcent): Likewise. (_nss_create_tablename): Likewise. * nis/nss_nisplus/nisplus-service.c (_nss_nisplus_parse_servent): Likewise. (_nss_create_tablename): Likewise. * nis/nss_nisplus/nisplus-spwd.c (_nss_create_tablename): Likewise. * libc.map: Export __stpcpy and __strdup. 1997-09-09 Andreas Schwab <schwab@issan.informatik.uni-dortmund.de> * math/Makefile (CFLAGS-test-float.c, CFLAGS-test-double.c, CFLAGS-test-ldouble.c): Pass -ffloat-store to avoid excessive precision. 1997-09-09 Andreas Schwab <schwab@issan.informatik.uni-dortmund.de> * include/rpc/auth_des.h: New file. 1997-09-09 Paul Eggert <eggert@twinsun.com> * time/mktime.c (__mktime_internal): Declare sec_requested even if !LEAP_SECONDS_POSSIBLE, since it's needed at the end when checking for time_t overflow. 1997-09-09 22:11 Ulrich Drepper <drepper@cygnus.com> * sysdeps/posix/getcwd.c: Correct test for too small buffer. Reported by Erik Troan <ewt@redhat.com>. * elf/dl-close.c: Include <bits/libc-lock.h>, not <libc-lock.h>. * elf/dl-open.c: Likewise. 1997-09-07 17:09 Richard Henderson <rth@cygnus.com> * sysdeps/alpha/Makefile: Kill setjmp_aux. * sysdeps/alpha/bits/setjmp.h: Rewrite in terms of an array. * sysdeps/alpha/__longjmp.c: Remove. * sysdeps/alpha/setjmp_aux.c: Remove. * sysdeps/alpha/__longjmp.S: New file. * sysdeps/alpha/bsd-_setjmp.S: Stub out. * sysdeps/alpha/bsd-setjmp.S: Likewise. * sysdeps/alpha/setjmp.S: Do the work; don't call __setjmp_aux. Move _setjmp and setjmp from bsd-*.S. 1997-09-06 20:20 Ulrich Drepper <drepper@cygnus.com> * include/rpc/auth.h: New file. * include/rpc/auth_unix.h: New file. 1997-09-06 Paul Eggert <eggert@twinsun.com> Fix gmtime so that it reports leap seconds when TZ indicates that leap seconds are desired. * time/gmtime.c (<stddef.h>): Remove unnecessary include. (gmtime): Put after gmtime_r, to help the compiler inline. (__tz_convert): New decl. (gmtime_r): Use __tz_convert instead of __offtime, so that leap seconds are handled correctly. * time/localtime.c (<errno.h>, <libc-lock.h>): Remove includes that are now unnecessary. (__tzset_internal, __tz_compute, __tzfile_compute, __use_tzfile, __tzset_lock): Remove extern decls that are now unnecessary. (localtime_internal): Moved to __tz_convert in tzset.c. so that localtime and gmtime can both use it easily. (localtime): Put after localtime_r, to help the compiler inline. (localtime_r): Use __tz_convert instead of localtime_internal. * time/strftime.c (__tz_compute): Remove unused (and now incorrect) decl. * time/tzfile.c (__tzfile_compute): New arg USE_LOCALTIME. * time/tzset.c (<errno.h>): Include. (_tmbuf): New decl. (__tzfile_compute): New function. (tz_compute): Renamed from __tz_compute. No longer extern. Remove redundant call to tzset_internal. (tzset_internal): Renamed from __tzset_internal. No longer extern. (tzset_lock): Renamed from __tzset_lock. No longer extern. (__tz_convert): New function, containing functionality of old localtime_internal function, plus locking and optional UTC. 1997-09-06 Paul Eggert <eggert@twinsun.com> * time/tzfile.c (__tzfile_read): Don't read a file if TZ is the empty string, just use UTC without leap seconds. This is for compatibility with the Olson code. 1997-09-06 Paul Eggert <eggert@twinsun.com> * time/tzset.c (__tzname_max): Lock tz data structures before invoking tzset_internal. * time/tzfile.c: Define compute_tzname_max statically. 1997-09-07 10:57 Thorsten Kukuk <kukuk@vt.uni-paderborn.de> * nis/nis_call.c: Remove not longer necessary HAVE_SECURE_RPC ifdefs. * nis/nis_intern.h: Likewise. * nis/nss_nis/nis-publickey.c: Likewise. * nis/nss_nisplus/nisplus-publickey.c: Likewise. * nis/ypclnt.c: Likewise. * sunrpc/auth_des.c: Don't dereference NULL pointer, initialize ad->ad_timediff. * sunrpc/auth_none.c: Don't define our own prototypes, use the one from the header files. * sunrpc/auth_unix.c: Likewise. * sunrpc/clnt_raw.c: Likewise. * sunrpc/clnt_tcp.c: Likewise. * sunrpc/rpc_cmsg.c: Likewise. * sunrpc/key_call.c: Fix signal handling. * sunrpc/openchild.c: Don't use /bin/sh to start /usr/etc/keyenvoy, or we will get a deadlock with NIS+. * sunrpc/rpc/auth.h: Add prototype for xdr_opaque_auth, don't define HAVE_SECURE_RPC. 1997-09-07 15:51 Ulrich Drepper <drepper@cygnus.com> * sysdeps/i386/bits/select.h [__GNUC__] (__FD_ZERO, __FD_SET, __FD_CLR, __FD_ISSET): Use correct casts to address array correctly. Reported by urbanw@cs.umu.se. 1997-09-07 05:07 Ulrich Drepper <drepper@cygnus.com> * elf/dl-close.c: Include <bits/libc-lock.h>, not <libc-lock.h>. * elf/dl-open.c: Likewise. * sysdeps/i386/memset.c: Undefine memset in case the header with the optimized functions is included. Patches by NIIBE Yutaka <gniibe@mri.co.jp>. * sysdeps/i386/bits/string.h [__PIC__] (strcspn, strspn, strpbrk, strsep): Use register for second parameter. * sysdeps/i386/i486/bits/string.h: Likewise. Reported by NIIBE Yutaka <gniibe@mri.co.jp>. 1997-09-03 09:48 Geoff Keating <geoffk@ozemail.com.au> * math/libm-test.c: Change various tolerances to match what the tested routines can actually provide. * math/Makefile: Add new tests. * math/atest-sincos.c: New file. * math/atest-exp.c: New file. * csu/Makefile: Give initfini.s and initfiniS.s their own CFLAGS-* macros so they can be overridden. * sysdeps/powerpc/Makefile [subdir=csu]: Override flags for initfiniS.s to use -fpic instead of -fPIC, because the sed script breaks otherwise. * sysdeps/powerpc/Makefile [build-shared]: Use -fpic not -fPIC for efficiency. * sysdeps/powerpc/dl-machine.h (ELF_MACHINE_RUNTIME_TRAMPOLINE): Don't use register 0, to let _mcount be in a shared object. * sysdeps/powerpc/dl-machine.h: Use full sentences in comments. Generally clean up. Suppress some code we don't need when relocating ld.so. * sysdeps/powerpc/test-arith.c: Change loop indices to size_t when appropriate to suppress gcc warning. * resolv/res_send.c: Suppress warning. * sunrpc/xdr_sizeof.c: Suppress warning. * FAQ: Add ppc-linux. * manual/maint.texi: Add ppc-linux. Explain that gcc can't build it yet. * sysdeps/unix/sysv/linux/powerpc/profil-counter.h: Correct for current kernels. 1997-08-15 07:45 Geoff Keating <geoffk@ozemail.com.au> * stdlib/fmtmsg.c: Use two parameters for __libc_once_define. * sysdeps/i386/machine-gmon.h: Correct typo. * sysdeps/unix/sysv/linux/powerpc/bits/mman.h: Change to match kernel. * sysdeps/generic/dl-sysdep.c: Add hook for bizzare PPC argument hack. * sysdeps/unix/sysv/linux/powerpc/dl-sysdep.c: Rewrite to use sysdeps/linux/dl-sysdep.c. * sysdeps/powerpc/Makefile [subdir=gmon]: Compile ppc-mcount. * sysdeps/powerpc/machine-gmon.h: Use ppc-mcount. * sysdeps/powerpc/ppc-mcount: New file. The following are mostly changes to allow profiling: * sysdeps/powerpc/add_n.S: Added. * sysdeps/powerpc/add_n.s: Removed. * sysdeps/powerpc/addmul_1.S: Added. * sysdeps/powerpc/addmul_1.s: Removed. * sysdeps/powerpc/bsd-_setjmp.S: Use JUMPTARGET macro. * sysdeps/powerpc/bsd-setjmp.S: Use JUMPTARGET macro. * sysdeps/powerpc/lshift.S: Added. * sysdeps/powerpc/lshift.s: Removed. * sysdeps/powerpc/memset.S: Added. * sysdeps/powerpc/memset.s: Removed. * sysdeps/powerpc/mul_1.S: Added. * sysdeps/powerpc/mul_1.s: Removed. * sysdeps/powerpc/rshift.S: Added. * sysdeps/powerpc/rshift.s: Removed. * sysdeps/powerpc/s_copysign.S: Use ENTRY, END, weak_alias macros. * sysdeps/powerpc/s_fabs.S: Use ENTRY, END, weak_alias macros. * sysdeps/powerpc/setjmp.S: Use JUMPTARGET macro. * sysdeps/powerpc/strchr.S: Added. * sysdeps/powerpc/strchr.s: Removed. * sysdeps/powerpc/strcmp.S: Added. * sysdeps/powerpc/strcmp.s: Removed. * sysdeps/powerpc/strlen.S: Added. * sysdeps/powerpc/strlen.s: Removed. * sysdeps/powerpc/sub_n.S: Added. * sysdeps/powerpc/sub_n.s: Removed. * sysdeps/powerpc/submul_1.S: Added. * sysdeps/powerpc/submul_1.s: Removed. * sysdeps/unix/sysv/linux/powerpc/_exit.S: Removed. * sysdeps/unix/sysv/linux/powerpc/brk.S: Added. * sysdeps/unix/sysv/linux/powerpc/brk.c: Removed. * sysdeps/unix/sysv/linux/powerpc/clone.S: Use new macros. Fix various bugs. Document that it isn't tested. * sysdeps/unix/sysv/linux/powerpc/sigreturn.S: Make look like sysdeps/unix/_exit.S. * sysdeps/unix/sysv/linux/powerpc/socket.S: Use new macros. * sysdeps/unix/sysv/linux/powerpc/syscall.S: Use new macros. * sysdeps/unix/sysv/linux/powerpc/sysdep.h: Define some new macros to make assembler (possibly) more portable, allow profiling, etc.
547 lines
18 KiB
C
547 lines
18 KiB
C
/* Inline math functions for i387.
|
|
Copyright (C) 1995, 1996, 1997 Free Software Foundation, Inc.
|
|
This file is part of the GNU C Library.
|
|
Contributed by John C. Bowman <bowman@math.ualberta.ca>, 1995.
|
|
|
|
The GNU C Library is free software; you can redistribute it and/or
|
|
modify it under the terms of the GNU Library General Public License as
|
|
published by the Free Software Foundation; either version 2 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
|
|
Library General Public License for more details.
|
|
|
|
You should have received a copy of the GNU Library General Public
|
|
License along with the GNU C Library; see the file COPYING.LIB. If not,
|
|
write to the Free Software Foundation, Inc., 59 Temple Place - Suite 330,
|
|
Boston, MA 02111-1307, USA. */
|
|
|
|
#ifndef _BITS_MATHINLINE_H
|
|
#define _BITS_MATHINLINE_H 1
|
|
|
|
|
|
#if defined __USE_ISOC9X && defined __GNUC__ && __GNUC__ >= 2
|
|
/* ISO C 9X defines some macros to perform unordered comparisons. The
|
|
ix87 FPU supports this with special opcodes and we should use them.
|
|
These must not be inline functions since we have to be able to handle
|
|
all floating-point types. */
|
|
# define isgreater(x, y) \
|
|
({ register char __result; \
|
|
__asm__ ("fucompp; fnstsw; testb $0x45, %%ah; setz %%al;" \
|
|
: "=a" (__result) : "u" (y), "t" (x) : "cc", "st", "st(1)"); \
|
|
__result; })
|
|
|
|
# define isgreaterequal(x, y) \
|
|
({ register char __result; \
|
|
__asm__ ("fucompp; fnstsw; testb $0x05, %%ah; setz %%al;" \
|
|
: "=a" (__result) : "u" (y), "t" (x) : "cc", "st", "st(1)"); \
|
|
__result; })
|
|
|
|
# define isless(x, y) \
|
|
({ register char __result; \
|
|
__asm__ ("fucompp; fnstsw; xorb $0x01, %%ah; testb $0x45, %%ah;" \
|
|
"setz %%al" \
|
|
: "=a" (__result) : "u" (y), "t" (x) : "cc", "st", "st(1)"); \
|
|
__result; })
|
|
|
|
# define islessequal(x, y) \
|
|
({ register char __result; \
|
|
__asm__ ("fucompp; fnstsw; xorb $0x01, %%ah; testb $0x05, %%ah;" \
|
|
"setz %%al" \
|
|
: "=a" (__result) : "u" (y), "t" (x) : "cc", "st", "st(1)"); \
|
|
__result; })
|
|
|
|
# define islessgreater(x, y) \
|
|
({ register char __result; \
|
|
__asm__ ("fucompp; fnstsw; testb $0x44, %%ah; setz %%al;" \
|
|
: "=a" (__result) : "u" (y), "t" (x) : "cc", "st", "st(1)"); \
|
|
__result; })
|
|
|
|
# define isunordered(x, y) \
|
|
({ register char __result; \
|
|
__asm__ ("fucompp; fnstsw; sahf; setp %%al" \
|
|
: "=a" (__result) : "u" (y), "t" (x) : "cc", "st", "st(1)"); \
|
|
__result; })
|
|
#endif
|
|
|
|
|
|
#ifdef __GNUC__
|
|
#if !defined __NO_MATH_INLINES && defined __OPTIMIZE__
|
|
|
|
/* The gcc, version 2.7 or below, has problems with all this inlining
|
|
code. So disable it for this version of the compiler. */
|
|
#if __GNUC__ > 2 || (__GNUC__ == 2 && __GNUC_MINOR__ > 7)
|
|
|
|
#ifdef __cplusplus
|
|
# define __MATH_INLINE __inline
|
|
#else
|
|
# define __MATH_INLINE extern __inline
|
|
#endif
|
|
|
|
/* A macro to define float, double, and long double versions of various
|
|
math functions for the ix87 FPU. FUNC is the function name (which will
|
|
be suffixed with f and l for the float and long double version,
|
|
respectively). OP is the name of the FPU operation. */
|
|
|
|
#if defined __USE_MISC || defined __USE_ISOC9X
|
|
# define __inline_mathop(func, op) \
|
|
__inline_mathop_ (double, func, op) \
|
|
__inline_mathop_ (float, __CONCAT(func,f), op) \
|
|
__inline_mathop_ (long double, __CONCAT(func,l), op)
|
|
#else
|
|
# define __inline_mathop(func, op) \
|
|
__inline_mathop_ (double, func, op)
|
|
#endif
|
|
|
|
#define __inline_mathop_(float_type, func, op) \
|
|
__inline_mathop_decl_ (float_type, func, op, "0" (__x))
|
|
|
|
|
|
#if defined __USE_MISC || defined __USE_ISOC9X
|
|
# define __inline_mathop_decl(func, op, params...) \
|
|
__inline_mathop_decl_ (double, func, op, params) \
|
|
__inline_mathop_decl_ (float, __CONCAT(func,f), op, params) \
|
|
__inline_mathop_decl_ (long double, __CONCAT(func,l), op, params)
|
|
#else
|
|
# define __inline_mathop_decl(func, op, params...) \
|
|
__inline_mathop_decl_ (double, func, op, params)
|
|
#endif
|
|
|
|
#define __inline_mathop_decl_(float_type, func, op, params...) \
|
|
__MATH_INLINE float_type func (float_type); \
|
|
__MATH_INLINE float_type func (float_type __x) \
|
|
{ \
|
|
register float_type __result; \
|
|
__asm __volatile__ (op : "=t" (__result) : params); \
|
|
return __result; \
|
|
}
|
|
|
|
|
|
#if defined __USE_MISC || defined __USE_ISOC9X
|
|
# define __inline_mathcode(func, arg, code) \
|
|
__inline_mathcode_ (double, func, arg, code) \
|
|
__inline_mathcode_ (float, __CONCAT(func,f), arg, code) \
|
|
__inline_mathcode_ (long double, __CONCAT(func,l), arg, code)
|
|
#else
|
|
# define __inline_mathcode(func, arg, code) \
|
|
__inline_mathcode_ (double, func, arg, code)
|
|
#endif
|
|
|
|
#define __inline_mathcode_(float_type, func, arg, code) \
|
|
__MATH_INLINE float_type func (float_type); \
|
|
__MATH_INLINE float_type func (float_type arg) \
|
|
{ \
|
|
code; \
|
|
}
|
|
|
|
|
|
#if defined __USE_MISC || defined __USE_ISOC9X
|
|
# define __inline_mathcode2(func, arg1, arg2, code) \
|
|
__inline_mathcode2_ (double, func, arg1, arg2, code) \
|
|
__inline_mathcode2_ (float, __CONCAT(func,f), arg1, arg2, code) \
|
|
__inline_mathcode2_ (long double, __CONCAT(func,l), arg1, arg2, code)
|
|
#else
|
|
# define __inline_mathcode2(func, arg1, arg2, code) \
|
|
__inline_mathcode2_ (double, func, arg1, arg2, code)
|
|
#endif
|
|
|
|
#define __inline_mathcode2_(float_type, func, arg1, arg2, code) \
|
|
__MATH_INLINE float_type func (float_type, float_type); \
|
|
__MATH_INLINE float_type func (float_type arg1, float_type arg2) \
|
|
{ \
|
|
code; \
|
|
}
|
|
|
|
|
|
/* Miscellaneous functions */
|
|
|
|
__inline_mathcode (__sgn, __x, \
|
|
return __x == 0.0 ? 0.0 : (__x > 0.0 ? 1.0 : -1.0))
|
|
|
|
__inline_mathcode (__pow2, __x, \
|
|
register long double __value; \
|
|
register long double __exponent; \
|
|
long int __p = (long int) __x; \
|
|
if (__x == (long double) __p) \
|
|
{ \
|
|
__asm __volatile__ \
|
|
("fscale" \
|
|
: "=t" (__value) : "0" (1.0), "u" (__x)); \
|
|
return __value; \
|
|
} \
|
|
__asm __volatile__ \
|
|
("fldl %%st(0)\n\t" \
|
|
"frndint # int(x)\n\t" \
|
|
"fxch\n\t" \
|
|
"fsub %%st(1) # fract(x)\n\t" \
|
|
"f2xm1 # 2^(fract(x)) - 1\n\t" \
|
|
: "=t" (__value), "=u" (__exponent) : "0" (__x)); \
|
|
__value += 1.0; \
|
|
__asm __volatile__ \
|
|
("fscale" \
|
|
: "=t" (__value) : "0" (__value), "u" (__exponent)); \
|
|
return __value)
|
|
|
|
#define __sincos_code \
|
|
register long double __cosr; \
|
|
register long double __sinr; \
|
|
__asm __volatile__ \
|
|
("fsincos\n\t" \
|
|
"fnstsw %%ax\n\t" \
|
|
"testl $0x400, %%eax\n\t" \
|
|
"jz 1f\n\t" \
|
|
"fldpi\n\t" \
|
|
"fadd %%st(0)\n\t" \
|
|
"fxch %%st(1)\n\t" \
|
|
"2: fprem1\n\t" \
|
|
"fnstsw %%ax\n\t" \
|
|
"testl $0x400, %%eax\n\t" \
|
|
"jnz 2b\n\t" \
|
|
"fstp %%st(1)\n\t" \
|
|
"fsincos\n\t" \
|
|
"1:" \
|
|
: "=t" (__cosr), "=u" (__sinr) : "0" (__x)); \
|
|
*__sinx = __sinr; \
|
|
*__cosx = __cosr
|
|
|
|
__MATH_INLINE void __sincos (double __x, double *__sinx, double *__cosx);
|
|
__MATH_INLINE void
|
|
__sincos (double __x, double *__sinx, double *__cosx)
|
|
{
|
|
__sincos_code;
|
|
}
|
|
|
|
__MATH_INLINE void __sincosf (float __x, float *__sinx, float *__cosx);
|
|
__MATH_INLINE void
|
|
__sincosf (float __x, float *__sinx, float *__cosx)
|
|
{
|
|
__sincos_code;
|
|
}
|
|
|
|
__MATH_INLINE void __sincosl (long double __x, long double *__sinx,
|
|
long double *__cosx);
|
|
__MATH_INLINE void
|
|
__sincosl (long double __x, long double *__sinx, long double *__cosx)
|
|
{
|
|
__sincos_code;
|
|
}
|
|
|
|
|
|
/* Optimized inline implementation, sometimes with reduced precision
|
|
and/or argument range. */
|
|
|
|
#define __expm1_code \
|
|
register long double __value; \
|
|
register long double __exponent; \
|
|
register long double __temp; \
|
|
__asm __volatile__ \
|
|
("fldl2e # e^x - 1 = 2^(x * log2(e)) - 1\n\t" \
|
|
"fmul %%st(1) # x * log2(e)\n\t" \
|
|
"fstl %%st(1)\n\t" \
|
|
"frndint # int(x * log2(e))\n\t" \
|
|
"fxch\n\t" \
|
|
"fsub %%st(1) # fract(x * log2(e))\n\t" \
|
|
"f2xm1 # 2^(fract(x * log2(e))) - 1\n\t" \
|
|
"fscale # 2^(x * log2(e)) - 2^(int(x * log2(e)))\n\t" \
|
|
: "=t" (__value), "=u" (__exponent) : "0" (__x)); \
|
|
__asm __volatile__ \
|
|
("fscale # 2^int(x * log2(e))\n\t" \
|
|
: "=t" (__temp) : "0" (1.0), "u" (__exponent)); \
|
|
__temp -= 1.0; \
|
|
return __temp + __value
|
|
__inline_mathcode_ (long double, __expm1l, __x, __expm1_code)
|
|
|
|
|
|
#define __exp_code \
|
|
register long double __value; \
|
|
register long double __exponent; \
|
|
__asm __volatile__ \
|
|
("fldl2e # e^x = 2^(x * log2(e))\n\t" \
|
|
"fmul %%st(1) # x * log2(e)\n\t" \
|
|
"fstl %%st(1)\n\t" \
|
|
"frndint # int(x * log2(e))\n\t" \
|
|
"fxch\n\t" \
|
|
"fsub %%st(1) # fract(x * log2(e))\n\t" \
|
|
"f2xm1 # 2^(fract(x * log2(e))) - 1\n\t" \
|
|
: "=t" (__value), "=u" (__exponent) : "0" (__x)); \
|
|
__value += 1.0; \
|
|
__asm __volatile__ \
|
|
("fscale" \
|
|
: "=t" (__value) : "0" (__value), "u" (__exponent)); \
|
|
return __value
|
|
__inline_mathcode (exp, __x, __exp_code)
|
|
__inline_mathcode_ (long double, __expl, __x, __exp_code)
|
|
|
|
|
|
__inline_mathcode (tan, __x, \
|
|
register long double __value; \
|
|
register long double __value2 __attribute__ ((unused)); \
|
|
__asm __volatile__ \
|
|
("fptan" \
|
|
: "=t" (__value2), "=u" (__value) : "0" (__x)); \
|
|
return __value)
|
|
|
|
|
|
#define __atan2_code \
|
|
register long double __value; \
|
|
__asm __volatile__ \
|
|
("fpatan\n\t" \
|
|
: "=t" (__value) : "0" (__x), "u" (__y) : "st(1)"); \
|
|
return __value
|
|
__inline_mathcode2 (atan2, __y, __x, __atan2_code)
|
|
__inline_mathcode2_ (long double, __atan2l, __y, __x, __atan2_code)
|
|
|
|
|
|
__inline_mathcode2 (fmod, __x, __y, \
|
|
register long double __value; \
|
|
__asm __volatile__ \
|
|
("1: fprem\n\t" \
|
|
"fnstsw %%ax\n\t" \
|
|
"sahf\n\t" \
|
|
"jp 1b" \
|
|
: "=t" (__value) : "0" (__x), "u" (__y) : "ax", "cc"); \
|
|
return __value)
|
|
|
|
|
|
__inline_mathcode2 (pow, __x, __y, \
|
|
register long double __value; \
|
|
register long double __exponent; \
|
|
long int __p = (long int) __y; \
|
|
if (__x == 0.0 && __y > 0.0) \
|
|
return 0.0; \
|
|
if (__y == (double) __p) \
|
|
{ \
|
|
long double __r = 1.0; \
|
|
if (__p == 0) \
|
|
return 1.0; \
|
|
if (__p < 0) \
|
|
{ \
|
|
__p = -__p; \
|
|
__x = 1.0 / __x; \
|
|
} \
|
|
while (1) \
|
|
{ \
|
|
if (__p & 1) \
|
|
__r *= __x; \
|
|
__p >>= 1; \
|
|
if (__p == 0) \
|
|
return __r; \
|
|
__x *= __x; \
|
|
} \
|
|
/* NOTREACHED */ \
|
|
} \
|
|
__asm __volatile__ \
|
|
("fyl2x" : "=t" (__value) : "0" (__x), "u" (1.0) : "st(1)"); \
|
|
__asm __volatile__ \
|
|
("fmul %%st(1) # y * log2(x)\n\t" \
|
|
"fst %%st(1)\n\t" \
|
|
"frndint # int(y * log2(x))\n\t" \
|
|
"fxch\n\t" \
|
|
"fsub %%st(1) # fract(y * log2(x))\n\t" \
|
|
"f2xm1 # 2^(fract(y * log2(x))) - 1\n\t" \
|
|
: "=t" (__value), "=u" (__exponent) : "0" (__y), "1" (__value)); \
|
|
__value += 1.0; \
|
|
__asm __volatile__ \
|
|
("fscale" \
|
|
: "=t" (__value) : "0" (__value), "u" (__exponent)); \
|
|
return __value)
|
|
|
|
|
|
__inline_mathop (sqrt, "fsqrt")
|
|
__inline_mathop_ (long double, __sqrtl, "fsqrt")
|
|
|
|
#if defined __GNUC__ && (__GNUC__ > 2 || __GNUC__ == 2 && __GNUC_MINOR__ >= 8)
|
|
__inline_mathcode_ (double, fabs, __x, return __builtin_fabs (__x))
|
|
__inline_mathcode_ (float, fabsf, __x, return __builtin_fabsf (__x))
|
|
__inline_mathcode_ (long double, fabsl, __x, return __builtin_fabsl (__x))
|
|
__inline_mathcode_ (long double, __fabsl, __x, return __builtin_fabsl (__x))
|
|
#else
|
|
__inline_mathop (fabs, "fabs")
|
|
__inline_mathop_ (long double, __fabsl, "fabs")
|
|
#endif
|
|
|
|
/* The argument range of this inline version is reduced. */
|
|
__inline_mathop (sin, "fsin")
|
|
/* The argument range of this inline version is reduced. */
|
|
__inline_mathop (cos, "fcos")
|
|
|
|
__inline_mathop (atan, "fld1; fpatan")
|
|
__inline_mathop (log, "fldln2; fxch; fyl2x")
|
|
__inline_mathop (log10, "fldlg2; fxch; fyl2x")
|
|
|
|
__inline_mathcode (asin, __x, return __atan2l (__x, __sqrtl (1.0 - __x * __x)))
|
|
__inline_mathcode (acos, __x, return __atan2l (__sqrtl (1.0 - __x * __x), __x))
|
|
|
|
__inline_mathcode_ (long double, __sgn1l, __x, return __x >= 0.0 ? 1.0 : -1.0)
|
|
|
|
|
|
/* The argument range of the inline version of sinhl is slightly reduced. */
|
|
__inline_mathcode (sinh, __x, \
|
|
register long double __exm1 = __expm1l (__fabsl (__x)); \
|
|
return 0.5 * (__exm1 / (__exm1 + 1.0) + __exm1) * __sgn1l (__x))
|
|
|
|
__inline_mathcode (cosh, __x, \
|
|
register long double __ex = __expl (__x); \
|
|
return 0.5 * (__ex + 1.0 / __ex))
|
|
|
|
__inline_mathcode (tanh, __x, \
|
|
register long double __exm1 = __expm1l (-__fabsl (__x + __x)); \
|
|
return __exm1 / (__exm1 + 2.0) * __sgn1l (-__x))
|
|
|
|
|
|
__inline_mathcode (floor, __x, \
|
|
register long double __value; \
|
|
__volatile unsigned short int __cw; \
|
|
__volatile unsigned short int __cwtmp; \
|
|
__asm __volatile ("fnstcw %0" : "=m" (__cw)); \
|
|
__cwtmp = (__cw & 0xf3ff) | 0x0400; /* rounding down */ \
|
|
__asm __volatile ("fldcw %0" : : "m" (__cwtmp)); \
|
|
__asm __volatile ("frndint" : "=t" (__value) : "0" (__x)); \
|
|
__asm __volatile ("fldcw %0" : : "m" (__cw)); \
|
|
return __value)
|
|
|
|
__inline_mathcode (ceil, __x, \
|
|
register long double __value; \
|
|
__volatile unsigned short int __cw; \
|
|
__volatile unsigned short int __cwtmp; \
|
|
__asm __volatile ("fnstcw %0" : "=m" (__cw)); \
|
|
__cwtmp = (__cw & 0xf3ff) | 0x0800; /* rounding up */ \
|
|
__asm __volatile ("fldcw %0" : : "m" (__cwtmp)); \
|
|
__asm __volatile ("frndint" : "=t" (__value) : "0" (__x)); \
|
|
__asm __volatile ("fldcw %0" : : "m" (__cw)); \
|
|
return __value)
|
|
|
|
#define __ldexp_code \
|
|
register long double __value; \
|
|
__asm __volatile__ \
|
|
("fscale" \
|
|
: "=t" (__value) : "0" (__x), "u" ((long double) __y)); \
|
|
return __value
|
|
|
|
__MATH_INLINE double ldexp (double __x, int __y);
|
|
__MATH_INLINE double
|
|
ldexp (double __x, int __y)
|
|
{
|
|
__ldexp_code;
|
|
}
|
|
|
|
|
|
/* Optimized versions for some non-standardized functions. */
|
|
#if defined __USE_ISOC9X || defined __USE_MISC
|
|
|
|
__inline_mathop(log2, "fld1; fxch; fyl2x")
|
|
|
|
__inline_mathcode (expm1, __x, __expm1_code)
|
|
|
|
/* We cannot rely on M_SQRT being defined. So we do it for ourself
|
|
here. */
|
|
# define __M_SQRT2 _Mldbl(1.41421356237309504880) /* sqrt(2) */
|
|
|
|
__inline_mathcode (log1p, __x, \
|
|
register long double __value; \
|
|
if (__fabsl (__x) >= 1.0 - 0.5 * __M_SQRT2) \
|
|
__value = logl (1.0 + __x); \
|
|
else \
|
|
__asm __volatile__ \
|
|
("fldln2\n\t" \
|
|
"fxch\n\t" \
|
|
"fyl2xp1" \
|
|
: "=t" (__value) : "0" (__x)); \
|
|
return __value)
|
|
|
|
|
|
/* The argument range of the inline version of asinhl is slightly reduced. */
|
|
__inline_mathcode (asinh, __x, \
|
|
register long double __y = __fabsl (__x); \
|
|
return (log1pl (__y * __y / (__sqrtl (__y * __y + 1.0) + 1.0) + __y) \
|
|
* __sgn1l (__x)))
|
|
|
|
__inline_mathcode (acosh, __x, \
|
|
return logl (__x + __sqrtl (__x - 1.0) * __sqrtl (__x + 1.0)))
|
|
|
|
__inline_mathcode (atanh, __x, \
|
|
register long double __y = __fabsl (__x); \
|
|
return (-0.5 * log1pl (-(__y + __y) / (1.0 + __y)) * \
|
|
__sgn1l (__x)))
|
|
|
|
|
|
/* The argument range of the inline version of hypotl is slightly reduced. */
|
|
__inline_mathcode2 (hypot, __x, __y, return __sqrtl (__x * __x + __y * __y))
|
|
|
|
__inline_mathcode(logb, __x, \
|
|
register long double __value; \
|
|
register long double __junk; \
|
|
__asm __volatile__ \
|
|
("fxtract\n\t" \
|
|
: "=t" (__junk), "=u" (__value) : "0" (__x)); \
|
|
return __value)
|
|
|
|
__MATH_INLINE float ldexpf (float __x, int __y);
|
|
__MATH_INLINE float
|
|
ldexpf (float __x, int __y)
|
|
{
|
|
__ldexp_code;
|
|
}
|
|
|
|
__MATH_INLINE long double ldexpl (long double __x, int __y);
|
|
__MATH_INLINE long double
|
|
ldexpl (long double __x, int __y)
|
|
{
|
|
__ldexp_code;
|
|
}
|
|
|
|
#endif
|
|
|
|
|
|
#ifdef __USE_MISC
|
|
|
|
__inline_mathcode2 (drem, __x, __y, \
|
|
register double __value; \
|
|
__asm __volatile__ \
|
|
("1: fprem1\n\t" \
|
|
"fstsw %%ax\n\t" \
|
|
"sahf\n\t" \
|
|
"jp 1b" \
|
|
: "=t" (__value) : "0" (__x), "u" (__y) : "ax", "cc"); \
|
|
return __value)
|
|
|
|
|
|
/* This function is used in the `isfinite' macro. */
|
|
__MATH_INLINE int __finite (double __x);
|
|
__MATH_INLINE int
|
|
__finite (double __x)
|
|
{
|
|
register int __result;
|
|
__asm__ __volatile__
|
|
("orl $0x800fffff, %0\n\t"
|
|
"incl %0\n\t"
|
|
"shrl $31, %0"
|
|
: "=q" (__result) : "0" (((int *) &__x)[1]) : "cc");
|
|
return __result;
|
|
}
|
|
|
|
/* Miscellaneous functions */
|
|
|
|
__inline_mathcode (__coshm1, __x, \
|
|
register long double __exm1 = __expm1l (__fabsl (__x)); \
|
|
return 0.5 * (__exm1 / (__exm1 + 1.0)) * __exm1)
|
|
|
|
__inline_mathcode (__acosh1p, __x, \
|
|
return log1pl (__x + __sqrtl (__x) * __sqrtl (__x + 2.0)))
|
|
|
|
#endif /* __USE_MISC */
|
|
|
|
/* Undefine some of the large macros which are not used anymore. */
|
|
#undef __expm1_code
|
|
#undef __exp_code
|
|
#undef __atan2_code
|
|
#undef __sincos_code
|
|
|
|
#endif /* Not gcc <= 2.7. */
|
|
#endif /* __NO_MATH_INLINES */
|
|
#endif /* __GNUC__ */
|
|
|
|
#endif /* _BITS_MATHINLINE_H */
|