5b01d53a83af269094b601d30a9155ffb27c432f
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 __STDC_VERSION__ < 199901L && (!defined(__cplusplus) || defined(_MSC_VER))
86 static inline long int
89 long int rounded
= (long int)(d
+ 0.5);
91 if (d
- floor(d
) == 0.5) {
93 rounded
+= (d
> 0) ? -1 : 1;
99 static inline long int
102 long int rounded
= (long int)(f
+ 0.5f
);
104 if (f
- floorf(f
) == 0.5f
) {
105 if (rounded
% 2 != 0)
106 rounded
+= (f
> 0) ? -1 : 1;
112 static inline long long int
115 long long int rounded
= (long long int)(d
+ 0.5);
117 if (d
- floor(d
) == 0.5) {
118 if (rounded
% 2 != 0)
119 rounded
+= (d
> 0) ? -1 : 1;
125 static inline long long int
128 long long int rounded
= (long long int)(f
+ 0.5f
);
130 if (f
- floorf(f
) == 0.5f
) {
131 if (rounded
% 2 != 0)
132 rounded
+= (f
> 0) ? -1 : 1;
141 * signbit() is a macro on Linux. Not available on Windows.
144 #define signbit(x) ((x) < 0.0f)
149 #define M_PI (3.14159265358979323846)
153 #define M_E (2.7182818284590452354)
157 #define M_LOG2E (1.4426950408889634074)
161 #define FLT_MAX_EXP 128
165 #if defined(fpclassify)
166 /* ISO C99 says that fpclassify is a macro. Assume that any implementation
167 * of fpclassify, whether it's in a C99 compiler or not, will be a macro.
169 #elif defined(__cplusplus)
170 /* For C++, fpclassify() should be defined in <cmath> */
171 #elif defined(_MSC_VER)
172 /* Not required on VS2013 and above. Oddly, the fpclassify() function
173 * doesn't exist in such a form on MSVC. This is an implementation using
174 * slightly different lower-level Windows functions.
178 static inline enum {FP_NAN
, FP_INFINITE
, FP_ZERO
, FP_SUBNORMAL
, FP_NORMAL
}
181 switch(_fpclass(x
)) {
182 case _FPCLASS_SNAN
: /* signaling NaN */
183 case _FPCLASS_QNAN
: /* quiet NaN */
185 case _FPCLASS_NINF
: /* negative infinity */
186 case _FPCLASS_PINF
: /* positive infinity */
188 case _FPCLASS_NN
: /* negative normal */
189 case _FPCLASS_PN
: /* positive normal */
191 case _FPCLASS_ND
: /* negative denormalized */
192 case _FPCLASS_PD
: /* positive denormalized */
194 case _FPCLASS_NZ
: /* negative zero */
195 case _FPCLASS_PZ
: /* positive zero */
198 /* Should never get here; but if we do, this will guarantee
199 * that the pattern is not treated like a number.
205 #error "Need to include or define an fpclassify function"
209 #endif /* #define _C99_MATH_H_ */