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
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
|
#include <config.h> /* built by autoconf */
#ifndef NULL
#define NULL 0
#endif
/**********************************************************************/
/* SUPPORT COMPILATION WITH A C++ COMPILER (IF DESIRED) */
/**********************************************************************/
/* Support declarations of C linkage in C++, for functions not declared in
C headers the way they should be. */
#ifdef __cplusplus
# define __C_LINKAGE "C"
#else
# define __C_LINKAGE /* empty */
#endif
/**********************************************************************/
/* Include all the C headers we'll need. Because many platforms lack one
or more standard headers or function declarations, there are numerous
tests and substitutions here. */
/**********************************************************************/
/**********************************************************************/
/* If libxmi is being compiling as part of the libplot/libplotter package,
add support for multithreading, provided that libc includes pthread
functions and pthread.h is present. (This comes first, because defining
_REENTRANT may alter system header files.) */
/**********************************************************************/
#ifdef LIBPLOT
#ifdef PTHREAD_SUPPORT
#ifdef HAVE_PTHREAD_H
#define _REENTRANT
#include <pthread.h>
#endif
#endif
#endif
/**********************************************************************/
/* INCLUDE stdio.h, ctype.h, errno.h. (SUBSTITUTE AS NECESSARY.) */
/**********************************************************************/
#include <stdio.h>
#include <ctype.h> /* why is this needed? */
#include <errno.h>
/***************************************************************************/
/* INCLUDE math.h, limits.h. (SUBSTITUTE AS NECESSARY.) */
/***************************************************************************/
#ifdef __DJGPP__
/* for DJGPP math.h, must specify that -lm will be used;
thanks mdruiter@cs.vu.nl */
#define _USE_LIBM_MATH_H
#endif
/* Include math.h, and whichever other math-related header files we have */
#include <math.h>
#ifdef HAVE_LIMITS_H
#include <limits.h> /* for INT_MAX */
#endif
#ifdef HAVE_VALUES_H
#include <values.h> /* for MAXINT (backup) */
#endif
/* Bounds on integer datatypes (should be in limits.h, but may not be). */
#ifndef UINT_MAX
#ifdef __STDC__
#define UINT_MAX ((unsigned int)(~(0U)))
#else
#define UINT_MAX ((unsigned int)(~((unsigned int)0)))
#endif
#endif /* not UINT_MAX */
#ifndef INT_MAX
#ifdef MAXINT
#define INT_MAX MAXINT
#else
#define INT_MAX ((int)(~(1U << (8 * (int)sizeof(int) - 1))))
#endif
#endif /* not INT_MAX */
/* IBM's definition of INT_MAX is bizarre, in AIX 4.1 at least, and using
IROUND() below will yield a warning message unless we repair it */
#ifdef _AIX
#ifdef __GNUC__
#undef INT_MAX
#define INT_MAX ((int)(~(1U << (8 * (int)sizeof(int) - 1))))
#endif
#endif
/**********************************************************************/
/* INCLUDE stdlib.h, string.h. (SUBSTITUTE AS NECESSARY; if STDC_HEADERS
is defined then they're both present, and stdarg.h and float.h too.) */
/**********************************************************************/
#ifdef STDC_HEADERS
#include <stdlib.h> /* for getenv, atoi, atof, etc. */
#include <string.h> /* for memcpy, memmove, strchr, malloc, etc. */
#else /* not STDC_HEADERS, must do a LOT of declarations by hand */
#ifdef HAVE_SYS_STDTYPES_H
#include <sys/stdtypes.h> /* SunOS needs this for size_t */
#endif
/* supply declarations for functions declared in stdlib.h */
extern __C_LINKAGE char *getenv (const char *name);
extern __C_LINKAGE int atoi (const char *nptr);
extern __C_LINKAGE double atof (const char *nptr);
/* supply definitions in stdlib.h */
#define EXIT_FAILURE 1 /* Failing exit status. */
#define EXIT_SUCCESS 0 /* Successful exit status. */
/* determine how to declare (or define) functions declared in string.h */
#ifdef HAVE_STRCHR
#ifdef HAVE_STRING_H
#include <string.h>
#else
#ifdef HAVE_STRINGS_H
#include <strings.h>
#endif
#endif
#else /* don't have strchr, prefer strings.h */
#ifdef HAVE_STRINGS_H
#include <strings.h>
#else
#ifdef HAVE_STRING_H
#include <string.h>
#endif
#endif
#define strchr index
#define strrchr rindex
#endif /* not HAVE_STRCHR */
#ifndef HAVE_MEMCPY
#define memcpy(d, s, n) bcopy ((s), (d), (n))
#endif /* not HAVE_MEMCPY */
#ifndef HAVE_MEMMOVE
#define memmove(d, s, n) bcopy ((s), (d), (n))
#endif /* not HAVE_MEMMOVE */
#ifndef HAVE_STRCASECMP /* will use local version */
extern __C_LINKAGE int strcasecmp (const char *s1, const char *s2);
#endif /* not HAVE_STRCASECMP */
/* supply declarations for more functions declared in stdlib.h */
#ifdef HAVE_MALLOC_H
#include <malloc.h>
#else
extern __C_LINKAGE void * malloc (size_t size);
extern __C_LINKAGE void * realloc (void * ptr, size_t size);
extern __C_LINKAGE void * calloc (size_t nmemb, size_t size);
extern __C_LINKAGE void free (void * ptr);
#endif /* not HAVE_MALLOC_H */
#endif /* not STDC_HEADERS */
/**************************************************************************/
/* Support the `bool' datatype, which our code uses extensively. */
/**************************************************************************/
#ifndef __cplusplus
#ifndef HAVE_BOOL_IN_CC
#ifdef __STDC__
# include <stdbool.h>
#else /* not __STDC__, do things the old-fashioned way */
typedef int bool;
#define false 0
#define true 1
#endif
#endif /* not HAVE_BOOL_IN_CC */
#endif /* not __cplusplus */
/**************************************************************************/
/* Define numerical constants (unofficial, so may not be in math.h). */
/**************************************************************************/
#ifndef M_PI
#define M_PI 3.14159265358979323846264
#endif
#ifndef M_PI_2
#define M_PI_2 1.57079632679489661923
#endif
#ifndef M_SQRT2
#define M_SQRT2 1.41421356237309504880
#endif
#ifndef M_SQRT3
#define M_SQRT3 1.73205080756887719
#endif
/**************************************************************************/
/* Define misc. math macros (in GCC, can be evaluated more rapidly). */
/**************************************************************************/
#ifdef __GNUC__
#define DMAX(a,b) ({double _a = (a), _b = (b); _a > _b ? _a : _b; })
#define DMIN(a,b) ({double _a = (a), _b = (b); _a < _b ? _a : _b; })
#define IMAX(a,b) ({int _a = (a), _b = (b); _a > _b ? _a : _b; })
#define IMIN(a,b) ({int _a = (a), _b = (b); _a < _b ? _a : _b; })
#define UMAX(a,b) ({unsigned int _a = (a), _b = (b); _a > _b ? _a : _b; })
#define UMIN(a,b) ({unsigned int _a = (a), _b = (b); _a < _b ? _a : _b; })
#define IROUND(x) ({double _x = (x); int _i; \
if (_x >= INT_MAX) _i = INT_MAX; \
else if (_x <= -(INT_MAX)) _i = -(INT_MAX); \
else _i = (_x > 0.0 ? (int)(_x + 0.5) : (int)(_x - 0.5)); \
_i;})
#define FROUND(x) ({double _x = (x); float _f; \
if (_x >= FLT_MAX) _f = FLT_MAX; \
else if (_x <= -(FLT_MAX)) _f = -(FLT_MAX); \
else _f = _x; \
_f;})
#define FABS(x) ((x) >= 0.0 ? (x) : -(x))
#define ICEIL(x) ({double _x = (x); int _i = (int)_x; \
((_x == _i) || (_x < 0.0)) ? _i : _i + 1;})
#define IFLOOR(x) ({double _x = (x); int _i = (int)_x; \
((_x == _i) || (_x > 0.0)) ? _i : _i - 1;})
#else
#define DMAX(a,b) ((a) > (b) ? (a) : (b))
#define DMIN(a,b) ((a) < (b) ? (a) : (b))
#define IMAX(a,b) ((a) > (b) ? (a) : (b))
#define IMIN(a,b) ((a) < (b) ? (a) : (b))
#define UMAX(a,b) ((a) > (b) ? (a) : (b))
#define UMIN(a,b) ((a) < (b) ? (a) : (b))
#define IROUND(x) ((int) ((x) > 0 ? (x) + 0.5 : (x) - 0.5))
#define FROUND(x) ((float)(x))
#define FABS(x) ((x) >= 0.0 ? (x) : -(x))
#define ICEIL(x) ((int)ceil(x))
#define IFLOOR(x) ((int)floor(x))
#endif
|
