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
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
|
/*
* bltMath.h --
*
* Copyright 1993-1998 Lucent Technologies, Inc.
*
* Permission to use, copy, modify, and distribute this software and
* its documentation for any purpose and without fee is hereby
* granted, provided that the above copyright notice appear in all
* copies and that both that the copyright notice and warranty
* disclaimer appear in supporting documentation, and that the names
* of Lucent Technologies any of their entities not be used in
* advertising or publicity pertaining to distribution of the software
* without specific, written prior permission.
*
* Lucent Technologies disclaims all warranties with regard to this
* software, including all implied warranties of merchantability and
* fitness. In no event shall Lucent Technologies be liable for any
* special, indirect or consequential damages or any damages
* whatsoever resulting from loss of use, data or profits, whether in
* an action of contract, negligence or other tortuous action, arising
* out of or in connection with the use or performance of this
* software.
*/
#ifndef _BLT_MATH_H
#define _BLT_MATH_H
#include <math.h>
#ifdef HAVE_FLOAT_H
#include <float.h>
#endif
#ifdef HAVE_IEEEFP_H
#include <ieeefp.h>
#endif /* HAVE_IEEEFP_H */
#ifndef M_PI
#define M_PI 3.14159265358979323846
#endif /* M_PI */
#ifndef M_PI_2
#define M_PI_2 1.57079632679489661923
#endif
#ifndef M_SQRT2
#define M_SQRT2 1.41421356237309504880
#endif /* M_SQRT2 */
#ifndef M_SQRT1_2
#define M_SQRT1_2 0.70710678118654752440
#endif /* M_SQRT1_2 */
#ifndef SHRT_MAX
#define SHRT_MAX 0x7FFF
#endif /* SHRT_MAX */
#ifndef SHRT_MIN
#define SHRT_MIN -(SHRT_MAX)
#endif /* SHRT_MAX */
#ifndef USHRT_MAX
#define USHRT_MAX 0xFFFF
#endif /* USHRT_MAX */
#ifndef INT_MAX
#define INT_MAX 2147483647
#endif /* INT_MAX */
#ifndef HAVE_FLOAT_H
/*
* ----------------------------------------------------------------------
*
* DBL_MIN, DBL_MAX --
*
* DBL_MAX and DBL_MIN are the largest and smaller double
* precision numbers that can be represented by the floating
* point hardware. If the compiler is ANSI, they can be found in
* float.h. Otherwise, we use HUGE_VAL or HUGE to determine
* them.
*
* ----------------------------------------------------------------------
*/
/*
* Don't want to include __infinity (definition of HUGE_VAL (SC1.x))
*/
#ifdef sun
#define DBL_MAX 1.7976931348623157E+308
#define DBL_MIN 2.2250738585072014E-308
#define DBL_EPSILON 2.2204460492503131e-16
#else
#ifndef DBL_EPSILON
#define DBL_EPSILON BLT_DBL_EPSILON
#endif
#ifdef HUGE_VAL
#define DBL_MAX HUGE_VAL
#define DBL_MIN (1/HUGE_VAL)
#else
#ifdef HUGE
#define DBL_MAX HUGE
#define DBL_MIN (1/HUGE)
#else
/*
* Punt: Assume relatively small values
*/
#define DBL_MAX 3.40282347E+38
#define DBL_MIN 1.17549435E-38
#endif /*HUGE*/
#endif /*HUGE_VAL*/
#endif /*sun*/
#endif /*!HAVE_FLOAT_H*/
/*
* ----------------------------------------------------------------------
*
* The following are macros replacing math library functions:
* "fabs", "fmod", "abs", "rint", and "exp10".
*
* Although many of these routines may exist in your math
* library, they aren't used in libtcl.a or libtk.a. This makes
* it difficult to dynamically load the BLT library as a shared
* object unless the math library is also shared (which isn't
* true on several systems). We can avoid the problem by
* replacing the "exotic" math routines with macros.
*
* ----------------------------------------------------------------------
*/
#undef ABS
#define ABS(x) (((x)<0)?(-(x)):(x))
#undef EXP10
#define EXP10(x) (pow(10.0,(x)))
#undef FABS
#define FABS(x) (((x)<0.0)?(-(x)):(x))
#undef SIGN
#define SIGN(x) (((x) < 0.0) ? -1 : 1)
/*
* Be careful when using the next two macros. They both assume the floating
* point number is less than the size of an int. That means, for example, you
* can't use these macros with numbers bigger than than 2^31-1.
*/
#undef FMOD
#define FMOD(x,y) ((x)-(((int)((x)/(y)))*y))
#undef ROUND
#define ROUND(x) ((int)((x) + (((x)<0.0) ? -0.5 : 0.5)))
#ifdef HAVE_FINITE
#define FINITE(x) finite(x)
#else
#ifdef HAVE_ISFINITE
#define FINITE(x) isfinite(x)
#else
#ifdef HAVE_ISNAN
#define FINITE(x) (!isnan(x))
#else
#define FINITE(x) (TRUE)
#endif /* HAVE_ISNAN */
#endif /* HAVE_ISFINITE */
#endif /* HAVE_FINITE */
extern double bltNaN;
#define DEFINED(x) (!isnan(x))
#define UNDEFINED(x) (isnan(x))
#define VALUE_UNDEFINED bltNaN
/*
* ----------------------------------------------------------------------
*
* On some systems "strdup" and "strcasecmp" are in the C library,
* but have no declarations in the C header files. Make sure we
* supply them here.
*
* ----------------------------------------------------------------------
*/
#ifdef NEED_DECL_STRDUP
extern char *strdup _ANSI_ARGS_((CONST char *s));
#endif /* NEED_DECL_STRDUP */
#ifndef HAVE_STRTOLOWER
extern void strtolower _ANSI_ARGS_((char *s));
#endif /* HAVE_STRTOLOWER */
#ifdef NEED_DECL_DRAND48
extern double drand48 _ANSI_ARGS_((void));
extern void srand48 _ANSI_ARGS_((long seed));
#endif /* NEED_DECL_DRAND48 */
#ifdef NEED_DECL_STRCASECMP
extern int strcasecmp _ANSI_ARGS_((CONST char *s1, CONST char *s2));
#endif /* NEED_DECL_STRCASECMP */
#endif /* BLT_MATH_H */
|
