1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
|
/* Declarations for math functions.
Copyright (C) 1991, 92, 93, 95, 96 Free Software Foundation, Inc.
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
not, write to the Free Software Foundation, Inc., 675 Mass Ave,
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>
/* Get machine-dependent HUGE_VAL value (returned on overflow).
On all IEEE754 machines, this is +Infinity. */
#include <huge_val.h>
/* Get machine-dependent NAN value (returned for some domain errors). */
#ifdef __USE_GNU
#include <nan.h>
#endif
/* 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
/* Include the file of declarations again, this time using `float'
instead of `double' and appending f to each function name. */
#ifndef _Mfloat_
#define _Mfloat_ float
#endif
#define _Mdouble_ _Mfloat_
#define __MATH_PRECNAME(name,r) name##f##r
#include <mathcalls.h>
#undef _Mdouble_
#undef __MATH_PRECNAME
/* Include the file of declarations again, this time using `long double'
instead of `double' and appending l to each function name. */
#ifndef _Mlong_double_
#define _Mlong_double_ long double
#endif
#define _Mdouble_ _Mlong_double_
#define __MATH_PRECNAME(name,r) name##l##r
#include <mathcalls.h>
#undef _Mdouble_
#undef __MATH_PRECNAME
#endif /* Use misc. */
#ifdef __USE_MISC
/* Support for various different standard error handling behaviors. */
typedef enum { _IEEE_ = -1, _SVID_, _XOPEN_, _POSIX_ } _LIB_VERSION_TYPE;
/* 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;
#endif
#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;
};
extern int matherr __P ((struct exception *));
#define X_TLOSS 1.41484755040568800000e+16
/* 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
/* SVID mode specifies returning this large value instead of infinity. */
#define HUGE FLT_MAX
#include <float.h> /* Defines FLT_MAX. */
#endif
#ifdef __USE_BSD
/* Some useful constants. */
#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. */
#endif
/* Get machine-dependent inline versions (if there are any). */
#if defined (__NO_MATH_INLINES) || defined (__OPTIMIZE__)
#include <__math.h>
#endif
__END_DECLS
#endif /* math.h */
|
