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) && __STDC_VERSION__ < 199901L && !defined(__cplusplus))
87 static inline long int
90 long int rounded = (long int)(d + 0.5);
92 if (d - floor(d) == 0.5) {
94 rounded += (d > 0) ? -1 : 1;
100 static inline long int
103 long int rounded = (long int)(f + 0.5f);
105 if (f - floorf(f) == 0.5f) {
106 if (rounded % 2 != 0)
107 rounded += (f > 0) ? -1 : 1;
113 static inline long long int
116 long long int rounded = (long long int)(d + 0.5);
118 if (d - floor(d) == 0.5) {
119 if (rounded % 2 != 0)
120 rounded += (d > 0) ? -1 : 1;
126 static inline long long int
129 long long int rounded = (long long int)(f + 0.5f);
131 if (f - floorf(f) == 0.5f) {
132 if (rounded % 2 != 0)
133 rounded += (f > 0) ? -1 : 1;
142 * signbit() is a macro on Linux. Not available on Windows.
145 #define signbit(x) ((x) < 0.0f)
150 #define M_PI (3.14159265358979323846)
154 #define M_E (2.7182818284590452354)
158 #define M_LOG2E (1.4426950408889634074)
162 #define FLT_MAX_EXP 128
166 #if defined(fpclassify)
167 /* ISO C99 says that fpclassify is a macro. Assume that any implementation
168 * of fpclassify, whether it's in a C99 compiler or not, will be a macro.
170 #elif defined(__cplusplus)
171 /* For C++, fpclassify() should be defined in <cmath> */
172 #elif defined(_MSC_VER)
173 /* Not required on VS2013 and above. Oddly, the fpclassify() function
174 * doesn't exist in such a form on MSVC. This is an implementation using
175 * slightly different lower-level Windows functions.
179 static inline enum {FP_NAN, FP_INFINITE, FP_ZERO, FP_SUBNORMAL, FP_NORMAL}
182 switch(_fpclass(x)) {
183 case _FPCLASS_SNAN: /* signaling NaN */
184 case _FPCLASS_QNAN: /* quiet NaN */
186 case _FPCLASS_NINF: /* negative infinity */
187 case _FPCLASS_PINF: /* positive infinity */
189 case _FPCLASS_NN: /* negative normal */
190 case _FPCLASS_PN: /* positive normal */
192 case _FPCLASS_ND: /* negative denormalized */
193 case _FPCLASS_PD: /* positive denormalized */
195 case _FPCLASS_NZ: /* negative zero */
196 case _FPCLASS_PZ: /* positive zero */
199 /* Should never get here; but if we do, this will guarantee
200 * that the pattern is not treated like a number.
206 #error "Need to include or define an fpclassify function"
210 #endif /* #define _C99_MATH_H_ */