1 /**************************************************************************
3 * Copyright 2007-2015 VMware, Inc.
6 * Permission is hereby granted, free of charge, to any person obtaining a
7 * copy of this software and associated documentation files (the
8 * "Software"), to deal in the Software without restriction, including
9 * without limitation the rights to use, copy, modify, merge, publish,
10 * distribute, sub license, and/or sell copies of the Software, and to
11 * permit persons to whom the Software is furnished to do so, subject to
12 * the following conditions:
14 * The above copyright notice and this permission notice (including the
15 * next paragraph) shall be included in all copies or substantial portions
18 * THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS
19 * OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF
20 * MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE AND NON-INFRINGEMENT.
21 * IN NO EVENT SHALL VMWARE AND/OR ITS SUPPLIERS BE LIABLE FOR
22 * ANY CLAIM, DAMAGES OR OTHER LIABILITY, WHETHER IN AN ACTION OF CONTRACT,
23 * TORT OR OTHERWISE, ARISING FROM, OUT OF OR IN CONNECTION WITH THE
24 * SOFTWARE OR THE USE OR OTHER DEALINGS IN THE SOFTWARE.
26 **************************************************************************/
29 * Wrapper for math.h which makes sure we have definitions of all the c99
38 #include "c99_compat.h"
43 /* This is to ensure that we get M_PI, etc. definitions */
44 #if !defined(_USE_MATH_DEFINES)
45 #error _USE_MATH_DEFINES define required when building with MSVC
49 #define isfinite(x) _finite((double)(x))
50 #define isnan(x) _isnan((double)(x))
51 #endif /* _MSC_VER < 1800 */
54 static inline double log2( double x
)
56 const double invln2
= 1.442695041;
57 return log( x
) * invln2
;
63 return x
>= 0.0 ? floor(x
+ 0.5) : ceil(x
- 0.5);
69 return x
>= 0.0f
? floorf(x
+ 0.5f
) : ceilf(x
- 0.5f
);
74 #include <float.h> // DBL_MAX
75 #define INFINITY (DBL_MAX + DBL_MAX)
79 #define NAN (INFINITY - INFINITY)
85 #if (defined(_MSC_VER) && _MSC_VER < 1800) || \
86 (!defined(_MSC_VER) && \
87 __STDC_VERSION__ < 199901L && \
88 (!defined(_XOPEN_SOURCE) || _XOPEN_SOURCE < 600) && \
89 !defined(__cplusplus))
91 static inline long int
94 long int rounded
= (long int)(d
+ 0.5);
96 if (d
- floor(d
) == 0.5) {
98 rounded
+= (d
> 0) ? -1 : 1;
104 static inline long int
107 long int rounded
= (long int)(f
+ 0.5f
);
109 if (f
- floorf(f
) == 0.5f
) {
110 if (rounded
% 2 != 0)
111 rounded
+= (f
> 0) ? -1 : 1;
117 static inline long long int
120 long long int rounded
= (long long int)(d
+ 0.5);
122 if (d
- floor(d
) == 0.5) {
123 if (rounded
% 2 != 0)
124 rounded
+= (d
> 0) ? -1 : 1;
130 static inline long long int
133 long long int rounded
= (long long int)(f
+ 0.5f
);
135 if (f
- floorf(f
) == 0.5f
) {
136 if (rounded
% 2 != 0)
137 rounded
+= (f
> 0) ? -1 : 1;
146 return powf(2.0f
, f
);
159 * signbit() is a macro on Linux. Not available on Windows.
162 #define signbit(x) ((x) < 0.0f)
167 #define M_PI (3.14159265358979323846)
171 #define M_E (2.7182818284590452354)
175 #define M_LOG2E (1.4426950408889634074)
179 #define FLT_MAX_EXP 128
183 #if defined(fpclassify)
184 /* ISO C99 says that fpclassify is a macro. Assume that any implementation
185 * of fpclassify, whether it's in a C99 compiler or not, will be a macro.
187 #elif defined(__cplusplus)
188 /* For C++, fpclassify() should be defined in <cmath> */
189 #elif defined(_MSC_VER)
190 /* Not required on VS2013 and above. Oddly, the fpclassify() function
191 * doesn't exist in such a form on MSVC. This is an implementation using
192 * slightly different lower-level Windows functions.
196 static inline enum {FP_NAN
, FP_INFINITE
, FP_ZERO
, FP_SUBNORMAL
, FP_NORMAL
}
199 switch(_fpclass(x
)) {
200 case _FPCLASS_SNAN
: /* signaling NaN */
201 case _FPCLASS_QNAN
: /* quiet NaN */
203 case _FPCLASS_NINF
: /* negative infinity */
204 case _FPCLASS_PINF
: /* positive infinity */
206 case _FPCLASS_NN
: /* negative normal */
207 case _FPCLASS_PN
: /* positive normal */
209 case _FPCLASS_ND
: /* negative denormalized */
210 case _FPCLASS_PD
: /* positive denormalized */
212 case _FPCLASS_NZ
: /* negative zero */
213 case _FPCLASS_PZ
: /* positive zero */
216 /* Should never get here; but if we do, this will guarantee
217 * that the pattern is not treated like a number.
223 #error "Need to include or define an fpclassify function"
227 #endif /* #define _C99_MATH_H_ */