1996-03-06 05:41:30 +08:00
|
|
|
/* Declarations for math functions.
|
|
|
|
Copyright (C) 1991, 92, 93, 95, 96 Free Software Foundation, Inc.
|
1995-02-18 09:27:10 +08:00
|
|
|
This file is part of the GNU C Library.
|
|
|
|
|
|
|
|
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
|
Sat Apr 20 18:13:00 1996 Roland McGrath <roland@delasyd.gnu.ai.mit.edu>
* Makerules (no-whole-archive): Test $(have-no-whole-archive), not
$(libc_cv_ld_no_whole_archive).
Sat Apr 20 17:07:17 1996 Ulrich Drepper <drepper@cygnus.com>
* assert/assert.h, ctype/ctype.h, dirent/dirent.h, errno.h,
grp/grp.h, io/fcntl.h, io/sys/stat.h, io/utime.h, locale/locale.h,
math/math.h, misc/nlist.h, misc/sgtty.h, misc/sys/file.h,
misc/sys/ioctl.h, misc/sys/uio.h, posix/sys/times.h,
posix/sys/types.h, posix/sys/utsname.h, posix/sys/wait.h,
posix/tar.h, posix/wordexp.h, pwd/pwd.h, resource/sys/vlimit.h,
resource/sys/vtimes.h, setjmp/setjmp.h, signal/signal.h,
stdio-common/printf.h, stdlib/alloca.h, stdlib/stdlib.h,
string/string.h, sysdeps/generic/sigaction.h,
sysdeps/generic/sigset.h, sysdeps/generic/sys/ptrace.h,
sysdeps/generic/sys/ptrace.h, sysdeps/unix/bsd/osf/sigaction.h,
sysdeps/unix/sysv/linux/sys/ptrace.h,
sysdeps/unix/sysv/minix/sigaction.h,
sysdeps/unix/sysv/sco3.2.4/sigaction.h,
sysdeps/unix/sysv/sysv4/sigaction.h,
sysdeps/unix/sysv/sysv4/sigset.h, termios/termios.h,
time/sys/time.h, time/time.h: Fix copyright comment.
* sysdeps/generic/dl-sysdep.c (_dl_sysdep_start): If uids and gids are
not given in stack aux vector, fetch them with syscalls.
1996-04-21 08:27:19 +08:00
|
|
|
not, write to the Free Software Foundation, Inc., 675 Mass Ave,
|
1995-02-18 09:27:10 +08:00
|
|
|
Cambridge, MA 02139, USA. */
|
|
|
|
|
|
|
|
/*
|
|
|
|
* ANSI Standard: 4.5 MATHEMATICS <math.h>
|
|
|
|
*/
|
|
|
|
|
|
|
|
#ifndef _MATH_H
|
|
|
|
|
|
|
|
#define _MATH_H 1
|
|
|
|
#include <features.h>
|
|
|
|
|
|
|
|
__BEGIN_DECLS
|
|
|
|
|
|
|
|
#define __need_Emath
|
|
|
|
#include <errno.h>
|
|
|
|
|
1996-03-06 05:41:30 +08:00
|
|
|
/* Get machine-dependent HUGE_VAL value (returned on overflow).
|
|
|
|
On all IEEE754 machines, this is +Infinity. */
|
1995-02-18 09:27:10 +08:00
|
|
|
#include <huge_val.h>
|
|
|
|
|
|
|
|
/* Get machine-dependent NAN value (returned for some domain errors). */
|
|
|
|
#ifdef __USE_GNU
|
|
|
|
#include <nan.h>
|
|
|
|
#endif
|
|
|
|
|
|
|
|
|
1996-03-06 05:41:30 +08:00
|
|
|
/* The file <mathcalls.h> contains the prototypes for all the actual
|
|
|
|
math functions. These macros are used for those prototypes, so
|
|
|
|
we can easily declare each function as both `name' and `__name',
|
|
|
|
and can declare the float versions `namef' and `__namef'. */
|
|
|
|
|
|
|
|
#define __MATHCALL(function,suffix, args) \
|
|
|
|
__MATHDECL (_Mdouble_, function,suffix, args)
|
|
|
|
#define __MATHDECL(type, function,suffix, args) \
|
|
|
|
__MATHDECL_1(type, function,suffix, args); \
|
|
|
|
__MATHDECL_1(type, __##function,suffix, args)
|
|
|
|
#define __MATHDECL_1(type, function,suffix, args) \
|
|
|
|
extern type __MATH_PRECNAME(function,suffix) args
|
|
|
|
|
|
|
|
#define _Mdouble_ double
|
|
|
|
#define __MATH_PRECNAME(name,r) name##r
|
|
|
|
#include <mathcalls.h>
|
|
|
|
#undef _Mdouble_
|
|
|
|
#undef __MATH_PRECNAME
|
|
|
|
|
|
|
|
#ifdef __USE_MISC
|
Thu May 30 11:24:05 1996 Roland McGrath <roland@delasyd.gnu.ai.mit.edu>
* po/header.pot: Replace with exact boilerplate pinard dictates.
* sysdeps/i386/strtok.S (Lillegal_argument): Remove this code to set
errno and the check that jumped to it.
* sysdeps/mach/hurd/Makefile (errnos.d): Use $(sed-remove-objpfx).
Thu May 30 03:21:57 1996 Ulrich Drepper <drepper@cygnus.com>
* FAQ: Document need of gperf program for developers.
* elf/elf.h: Fix typos in comments.
* libio/stdio.h [!__STRICT_ANSI__ || _POSIX_SOURCE]: Add
prototypes for `ctermid' and `cuserid'.
* locale/programs/locale.c: Switch to user selected locale
before printing variables.
* math/Makefile [$(long-double-fcts)==yes]: Define long-m-routines
and long-c-routines. Only if the `long double' data type is
available we need to compile the functions.
(libm-routines): Add $(long-m-routines).
(routines): Remove isinfl, isnanl. Use new file s_isinfl and
s_isnanl instead if `long double' is available.
* math/math.h: Include <mathcalls.h> again to define `long double'
functions.
* math/math_private.h: Define data types, prototypes and access
macros for `long double'.
* stdlib/stdlib.h: Add prototypes for `strtoll' and `strtoull'.
[GCC2 && OPTIMIZE]: Define strto{,u}ll as inline function which
calls __strto{,u}q_internal.
* stdlib/strfmon.c: Replace PTR by `void *'.
* stdlib/strtoq.c: Define strtoll as weak alias.
* stdlib/strtouq.c: Define strtoull as weak alias.
* string/tester.c: Correct `strsep' test.
* sysdeps/generic/strsep.c: Make compatible with BSD version.
Trailing characters of skip set are not skipped. In this case
empty tokens are returned.
* sysdeps/i386/isinfl.c, sysdeps/i386/isnanl.c,
sysdeps/ieee754/isinf.c, sysdeps/ieee754/isinfl.c,
sysdeps/ieee754/isnan.c, sysdeps/ieee754/isnanl.c: Removed. We
now use the versions part of libm.
* sysdeps/i386/strsep.S: Removed. Generic C version is of
similar speed.
* sysdeps/i386/strtok.S: Remove support for `strsep'.
* sysdeps/libm-i387/e_acosl.S, sysdeps/libm-i387/s_ceill.S,
sysdeps/libm-i387/s_copysignl.S, sysdeps/libm-i387/s_finitel.S,
sysdeps/libm-i387/s_floorl.S, sysdeps/libm-i387/s_isinfl.c,
sysdeps/libm-i387/s_isnanl.c, sysdeps/libm-i387/s_nextafterl.c,
sysdeps/libm-i387/s_rintl.S, sysdeps/libm-i387/s_significandl.S:
New i387 specific math functions implementing `long double'
versions.
* sysdeps/libm-ieee754/s_ceill.c,
sysdeps/libm-ieee754/s_copysignl.c,
sysdeps/libm-ieee754/s_fabsl.c, sysdeps/libm-ieee754/s_finitel.c,
sysdeps/libm-ieee754/s_floorl.c, sysdeps/libm-ieee754/s_isinfl.c,
sysdeps/libm-ieee754/s_isnanl.c,
sysdeps/libm-ieee754/s_nextafterl.c,
sysdeps/libm-ieee754/s_rintl.c, sysdeps/libm-ieee754/s_scalbnl.c,
sysdeps/libm-ieee754/s_significandl.c: New generic `long double'
versions of libm functions.
* sysdeps/libm-i387/e_exp.S: Add a few comments to explain the
Intel FPU nonsense.
* sysdeps/libm-i387/s_ceil.S, sysdeps/libm-i387/s_ceilf.S,
sysdeps/libm-i387/s_floor.S, sysdeps/libm-i387/s_floorf.S: Correct
handling of local variables. The old version created a stack
frame but stored the values outside.
* sysdeps/libm-ieee754/s_isinf.c, sysdeps/libm-ieee754/s_isnan.c
[!NO_LONG_DOUBLE]: Define alias with `long double' versions name.
* login/pututline_r.c: Include sys/stat.h. Fix typos.
according to currently used locale for category LC_CTYPE by
inet_nsap_ntoa. Now in <arpa/inet.h>.
_IO_dup2 to contain complete parameter list.
1996-05-31 00:12:42 +08:00
|
|
|
|
1996-03-06 05:41:30 +08:00
|
|
|
/* Include the file of declarations again, this type using `float'
|
|
|
|
instead of `double' and appending f to each function name. */
|
|
|
|
|
|
|
|
#define _Mdouble_ float
|
|
|
|
#define __MATH_PRECNAME(name,r) name##f##r
|
|
|
|
#include <mathcalls.h>
|
|
|
|
#undef _Mdouble_
|
|
|
|
#undef __MATH_PRECNAME
|
Thu May 30 11:24:05 1996 Roland McGrath <roland@delasyd.gnu.ai.mit.edu>
* po/header.pot: Replace with exact boilerplate pinard dictates.
* sysdeps/i386/strtok.S (Lillegal_argument): Remove this code to set
errno and the check that jumped to it.
* sysdeps/mach/hurd/Makefile (errnos.d): Use $(sed-remove-objpfx).
Thu May 30 03:21:57 1996 Ulrich Drepper <drepper@cygnus.com>
* FAQ: Document need of gperf program for developers.
* elf/elf.h: Fix typos in comments.
* libio/stdio.h [!__STRICT_ANSI__ || _POSIX_SOURCE]: Add
prototypes for `ctermid' and `cuserid'.
* locale/programs/locale.c: Switch to user selected locale
before printing variables.
* math/Makefile [$(long-double-fcts)==yes]: Define long-m-routines
and long-c-routines. Only if the `long double' data type is
available we need to compile the functions.
(libm-routines): Add $(long-m-routines).
(routines): Remove isinfl, isnanl. Use new file s_isinfl and
s_isnanl instead if `long double' is available.
* math/math.h: Include <mathcalls.h> again to define `long double'
functions.
* math/math_private.h: Define data types, prototypes and access
macros for `long double'.
* stdlib/stdlib.h: Add prototypes for `strtoll' and `strtoull'.
[GCC2 && OPTIMIZE]: Define strto{,u}ll as inline function which
calls __strto{,u}q_internal.
* stdlib/strfmon.c: Replace PTR by `void *'.
* stdlib/strtoq.c: Define strtoll as weak alias.
* stdlib/strtouq.c: Define strtoull as weak alias.
* string/tester.c: Correct `strsep' test.
* sysdeps/generic/strsep.c: Make compatible with BSD version.
Trailing characters of skip set are not skipped. In this case
empty tokens are returned.
* sysdeps/i386/isinfl.c, sysdeps/i386/isnanl.c,
sysdeps/ieee754/isinf.c, sysdeps/ieee754/isinfl.c,
sysdeps/ieee754/isnan.c, sysdeps/ieee754/isnanl.c: Removed. We
now use the versions part of libm.
* sysdeps/i386/strsep.S: Removed. Generic C version is of
similar speed.
* sysdeps/i386/strtok.S: Remove support for `strsep'.
* sysdeps/libm-i387/e_acosl.S, sysdeps/libm-i387/s_ceill.S,
sysdeps/libm-i387/s_copysignl.S, sysdeps/libm-i387/s_finitel.S,
sysdeps/libm-i387/s_floorl.S, sysdeps/libm-i387/s_isinfl.c,
sysdeps/libm-i387/s_isnanl.c, sysdeps/libm-i387/s_nextafterl.c,
sysdeps/libm-i387/s_rintl.S, sysdeps/libm-i387/s_significandl.S:
New i387 specific math functions implementing `long double'
versions.
* sysdeps/libm-ieee754/s_ceill.c,
sysdeps/libm-ieee754/s_copysignl.c,
sysdeps/libm-ieee754/s_fabsl.c, sysdeps/libm-ieee754/s_finitel.c,
sysdeps/libm-ieee754/s_floorl.c, sysdeps/libm-ieee754/s_isinfl.c,
sysdeps/libm-ieee754/s_isnanl.c,
sysdeps/libm-ieee754/s_nextafterl.c,
sysdeps/libm-ieee754/s_rintl.c, sysdeps/libm-ieee754/s_scalbnl.c,
sysdeps/libm-ieee754/s_significandl.c: New generic `long double'
versions of libm functions.
* sysdeps/libm-i387/e_exp.S: Add a few comments to explain the
Intel FPU nonsense.
* sysdeps/libm-i387/s_ceil.S, sysdeps/libm-i387/s_ceilf.S,
sysdeps/libm-i387/s_floor.S, sysdeps/libm-i387/s_floorf.S: Correct
handling of local variables. The old version created a stack
frame but stored the values outside.
* sysdeps/libm-ieee754/s_isinf.c, sysdeps/libm-ieee754/s_isnan.c
[!NO_LONG_DOUBLE]: Define alias with `long double' versions name.
* login/pututline_r.c: Include sys/stat.h. Fix typos.
according to currently used locale for category LC_CTYPE by
inet_nsap_ntoa. Now in <arpa/inet.h>.
_IO_dup2 to contain complete parameter list.
1996-05-31 00:12:42 +08:00
|
|
|
|
|
|
|
/* Include the file of declarations again, this type using `long double'
|
|
|
|
instead of `double' and appending l to each function name. */
|
|
|
|
|
|
|
|
#define _Mdouble_ long double
|
|
|
|
#define __MATH_PRECNAME(name,r) name##l##r
|
|
|
|
#include <mathcalls.h>
|
|
|
|
#undef _Mdouble_
|
|
|
|
#undef __MATH_PRECNAME
|
|
|
|
|
|
|
|
#endif /* Use misc. */
|
1995-02-18 09:27:10 +08:00
|
|
|
|
|
|
|
|
|
|
|
#ifdef __USE_MISC
|
1996-03-06 05:41:30 +08:00
|
|
|
/* Support for various different standard error handling behaviors. */
|
1995-02-18 09:27:10 +08:00
|
|
|
|
1996-03-06 05:41:30 +08:00
|
|
|
typedef enum { _IEEE_ = -1, _SVID_, _XOPEN_, _POSIX_ } _LIB_VERSION_TYPE;
|
1995-02-18 09:27:10 +08:00
|
|
|
|
1996-03-06 05:41:30 +08:00
|
|
|
/* This variable can be changed at run-time to any of the values above to
|
|
|
|
affect floating point error handling behavior (it may also be necessary
|
|
|
|
to change the hardware FPU exception settings). */
|
|
|
|
extern _LIB_VERSION_TYPE _LIB_VERSION;
|
1995-02-18 09:27:10 +08:00
|
|
|
#endif
|
|
|
|
|
|
|
|
|
1996-03-06 05:41:30 +08:00
|
|
|
#ifdef __USE_SVID
|
|
|
|
/* In SVID error handling, `matherr' is called with this description
|
|
|
|
of the exceptional condition. */
|
|
|
|
struct exception
|
|
|
|
{
|
|
|
|
int type;
|
|
|
|
char *name;
|
|
|
|
double arg1;
|
|
|
|
double arg2;
|
|
|
|
double retval;
|
|
|
|
};
|
1995-02-18 09:27:10 +08:00
|
|
|
|
1996-03-06 05:41:30 +08:00
|
|
|
extern int matherr __P ((struct exception *));
|
1995-02-18 09:27:10 +08:00
|
|
|
|
1996-03-06 05:41:30 +08:00
|
|
|
#define X_TLOSS 1.41484755040568800000e+16
|
1995-02-18 09:27:10 +08:00
|
|
|
|
1996-03-06 05:41:30 +08:00
|
|
|
/* Types of exceptions in the `type' field. */
|
|
|
|
#define DOMAIN 1
|
|
|
|
#define SING 2
|
|
|
|
#define OVERFLOW 3
|
|
|
|
#define UNDERFLOW 4
|
|
|
|
#define TLOSS 5
|
|
|
|
#define PLOSS 6
|
1995-02-18 09:27:10 +08:00
|
|
|
|
1996-03-06 05:41:30 +08:00
|
|
|
/* SVID mode specifies returning this large value instead of infinity. */
|
|
|
|
#define HUGE FLT_MAX
|
|
|
|
#include <float.h> /* Defines FLT_MAX. */
|
1995-02-18 09:27:10 +08:00
|
|
|
|
|
|
|
#endif
|
|
|
|
|
|
|
|
|
|
|
|
#ifdef __USE_BSD
|
1996-06-03 12:46:40 +08:00
|
|
|
|
1995-02-18 09:27:10 +08:00
|
|
|
/* Some useful constants. */
|
1996-06-03 12:46:40 +08:00
|
|
|
#define M_E _Mldbl(2.7182818284590452354) /* e */
|
|
|
|
#define M_LOG2E _Mldbl(1.4426950408889634074) /* log 2e */
|
|
|
|
#define M_LOG10E _Mldbl(0.43429448190325182765) /* log 10e */
|
|
|
|
#define M_LN2 _Mldbl(0.69314718055994530942) /* log e2 */
|
|
|
|
#define M_LN10 _Mldbl(2.30258509299404568402) /* log e10 */
|
|
|
|
#define M_PI _Mldbl(3.14159265358979323846) /* pi */
|
|
|
|
#define M_PI_2 _Mldbl(1.57079632679489661923) /* pi/2 */
|
|
|
|
#define M_PI_4 _Mldbl(0.78539816339744830962) /* pi/4 */
|
|
|
|
#define M_1_PI _Mldbl(0.31830988618379067154) /* 1/pi */
|
|
|
|
#define M_2_PI _Mldbl(0.63661977236758134308) /* 2/pi */
|
|
|
|
#define M_2_SQRTPI _Mldbl(1.12837916709551257390) /* 2/sqrt(pi) */
|
|
|
|
#define M_SQRT2 _Mldbl(1.41421356237309504880) /* sqrt(2) */
|
|
|
|
#define M_SQRT1_2 _Mldbl(0.70710678118654752440) /* 1/sqrt(2) */
|
|
|
|
|
|
|
|
/* Our constants might specify more precision than `double' can represent.
|
|
|
|
Use `long double' constants in standard and GNU C, where they are
|
|
|
|
supported and the cast to `double'. */
|
|
|
|
#if __STDC__ - 0 || __GNUC__ - 0
|
|
|
|
#define _Mldbl(x) x##L
|
|
|
|
#else /* Traditional C. */
|
|
|
|
#define _Mldbl(x) x
|
|
|
|
#endif /* Standard or GNU C. */
|
|
|
|
|
1995-02-18 09:27:10 +08:00
|
|
|
#endif
|
|
|
|
|
|
|
|
|
1996-03-20 08:18:44 +08:00
|
|
|
/* Get machine-dependent inline versions (if there are any). */
|
|
|
|
#if defined (__NO_MATH_INLINES) || defined (__OPTIMIZE__)
|
|
|
|
#include <__math.h>
|
|
|
|
#endif
|
|
|
|
|
|
|
|
|
|
|
|
__END_DECLS
|
|
|
|
|
|
|
|
|
1995-02-18 09:27:10 +08:00
|
|
|
#endif /* math.h */
|