2 * This file has no copyright assigned and is placed in the Public Domain.
3 * This file is part of the w64 mingw-runtime package.
4 * No warranty is given; refer to the file DISCLAIMER within this package.
10 #pragma GCC system_header
17 #pragma pack(push,_CRT_PACKING)
26 #ifndef __STRICT_ANSI__
28 #define DOMAIN _DOMAIN
30 #define OVERFLOW _OVERFLOW
31 #define UNDERFLOW _UNDERFLOW
37 #ifndef __STRICT_ANSI__
38 #define M_E 2.71828182845904523536
39 #define M_LOG2E 1.44269504088896340736
40 #define M_LOG10E 0.434294481903251827651
41 #define M_LN2 0.693147180559945309417
42 #define M_LN10 2.30258509299404568402
43 #define M_PI 3.14159265358979323846
44 #define M_PI_2 1.57079632679489661923
45 #define M_PI_4 0.785398163397448309616
46 #define M_1_PI 0.318309886183790671538
47 #define M_2_PI 0.636619772367581343076
48 #define M_2_SQRTPI 1.12837916709551257390
49 #define M_SQRT2 1.41421356237309504880
50 #define M_SQRT1_2 0.707106781186547524401
53 #ifndef __STRICT_ANSI__
54 /* See also float.h */
55 #ifndef __MINGW_FPCLASS_DEFINED
56 #define __MINGW_FPCLASS_DEFINED 1
57 #define _FPCLASS_SNAN 0x0001 /* Signaling "Not a Number" */
58 #define _FPCLASS_QNAN 0x0002 /* Quiet "Not a Number" */
59 #define _FPCLASS_NINF 0x0004 /* Negative Infinity */
60 #define _FPCLASS_NN 0x0008 /* Negative Normal */
61 #define _FPCLASS_ND 0x0010 /* Negative Denormal */
62 #define _FPCLASS_NZ 0x0020 /* Negative Zero */
63 #define _FPCLASS_PZ 0x0040 /* Positive Zero */
64 #define _FPCLASS_PD 0x0080 /* Positive Denormal */
65 #define _FPCLASS_PN 0x0100 /* Positive Normal */
66 #define _FPCLASS_PINF 0x0200 /* Positive Infinity */
74 #ifndef _EXCEPTION_DEFINED
75 #define _EXCEPTION_DEFINED
85 #ifndef _COMPLEX_DEFINED
86 #define _COMPLEX_DEFINED
99 extern double *_imp___HUGE
;
100 #define _HUGE (*_imp___HUGE)
104 #define HUGE_VAL _HUGE
106 #ifndef _CRT_ABS_DEFINED
107 #define _CRT_ABS_DEFINED
108 int __cdecl
abs(int _X
);
109 long __cdecl
labs(long _X
);
111 double __cdecl
acos(double _X
);
112 double __cdecl
asin(double _X
);
113 double __cdecl
atan(double _X
);
114 double __cdecl
atan2(double _Y
,double _X
);
115 #ifndef _SIGN_DEFINED
116 #define _SIGN_DEFINED
117 _CRTIMP
double __cdecl
_copysign (double _Number
,double _Sign
);
118 _CRTIMP
double __cdecl
_chgsign (double _X
);
120 double __cdecl
cos(double _X
);
121 double __cdecl
cosh(double _X
);
122 double __cdecl
exp(double _X
);
123 double __cdecl
expm1(double _X
);
124 double __cdecl
fabs(double _X
);
125 double __cdecl
fmod(double _X
,double _Y
);
126 double __cdecl
log(double _X
);
127 double __cdecl
log10(double _X
);
128 double __cdecl
pow(double _X
,double _Y
);
129 double __cdecl
sin(double _X
);
130 double __cdecl
sinh(double _X
);
131 double __cdecl
tan(double _X
);
132 double __cdecl
tanh(double _X
);
133 double __cdecl
sqrt(double _X
);
134 #ifndef _CRT_ATOF_DEFINED
135 #define _CRT_ATOF_DEFINED
136 double __cdecl
atof(const char *_String
);
137 double __cdecl
_atof_l(const char *_String
,_locale_t _Locale
);
140 _CRTIMP
double __cdecl
_cabs(struct _complex _ComplexA
);
141 double __cdecl
ceil(double _X
);
142 double __cdecl
floor(double _X
);
143 double __cdecl
frexp(double _X
,int *_Y
);
144 double __cdecl
_hypot(double _X
,double _Y
);
145 _CRTIMP
double __cdecl
_j0(double _X
);
146 _CRTIMP
double __cdecl
_j1(double _X
);
147 _CRTIMP
double __cdecl
_jn(int _X
,double _Y
);
148 double __cdecl
ldexp(double _X
,int _Y
);
149 #ifndef _CRT_MATHERR_DEFINED
150 #define _CRT_MATHERR_DEFINED
151 int __cdecl
_matherr(struct _exception
*_Except
);
153 double __cdecl
modf(double _X
,double *_Y
);
154 _CRTIMP
double __cdecl
_y0(double _X
);
155 _CRTIMP
double __cdecl
_y1(double _X
);
156 _CRTIMP
double __cdecl
_yn(int _X
,double _Y
);
158 #if(defined(_X86_) && !defined(__x86_64))
159 _CRTIMP
int __cdecl
_set_SSE2_enable(int _Flag
);
160 /* from libmingwex */
161 float __cdecl
_hypotf(float _X
,float _Y
);
164 float frexpf(float _X
,int *_Y
);
165 float __cdecl
ldexpf(float _X
,int _Y
);
166 long double __cdecl
ldexpl(long double _X
,int _Y
);
167 float __cdecl
acosf(float _X
);
168 float __cdecl
asinf(float _X
);
169 float __cdecl
atanf(float _X
);
170 float __cdecl
atan2f(float _X
,float _Y
);
171 float __cdecl
cosf(float _X
);
172 float __cdecl
sinf(float _X
);
173 float __cdecl
tanf(float _X
);
174 float __cdecl
coshf(float _X
);
175 float __cdecl
sinhf(float _X
);
176 float __cdecl
tanhf(float _X
);
177 float __cdecl
expf(float _X
);
178 float __cdecl
expm1f(float _X
);
179 float __cdecl
logf(float _X
);
180 float __cdecl
log10f(float _X
);
181 float __cdecl
modff(float _X
,float *_Y
);
182 float __cdecl
powf(float _X
,float _Y
);
183 float __cdecl
sqrtf(float _X
);
184 float __cdecl
ceilf(float _X
);
185 float __cdecl
floorf(float _X
);
186 float __cdecl
fmodf(float _X
,float _Y
);
187 float __cdecl
_hypotf(float _X
,float _Y
);
188 float __cdecl
fabsf(float _X
);
189 #if !defined(__ia64__)
190 /* from libmingwex */
191 float __cdecl
_copysignf (float _Number
,float _Sign
);
192 float __cdecl
_chgsignf (float _X
);
193 float __cdecl
_logbf(float _X
);
194 float __cdecl
_nextafterf(float _X
,float _Y
);
195 int __cdecl
_finitef(float _X
);
196 int __cdecl
_isnanf(float _X
);
197 int __cdecl
_fpclassf(float _X
);
201 __CRT_INLINE
long double __cdecl
fabsl (long double x
)
204 __asm__ ("fabs;" : "=t" (res
) : "0" (x
));
207 #define _hypotl(x,y) ((long double)_hypot((double)(x),(double)(y)))
208 #define _matherrl _matherr
209 __CRT_INLINE
long double _chgsignl(long double _Number
) { return _chgsign((double)(_Number
)); }
210 __CRT_INLINE
long double _copysignl(long double _Number
,long double _Sign
) { return _copysign((double)(_Number
),(double)(_Sign
)); }
211 __CRT_INLINE
float frexpf(float _X
,int *_Y
) { return ((float)frexp((double)_X
,_Y
)); }
213 #if !defined (__ia64__)
214 __CRT_INLINE
float __cdecl
fabsf (float x
)
217 __asm__ ("fabs;" : "=t" (res
) : "0" (x
));
221 __CRT_INLINE
float __cdecl
ldexpf (float x
, int expn
) { return (float) ldexp (x
, expn
); }
225 __CRT_INLINE
long double __cdecl
fabsl (long double x
)
228 __asm__ ("fabs;" : "=t" (res
) : "0" (x
));
231 __CRT_INLINE
long double modfl(long double _X
,long double *_Y
) {
232 double _Di
,_Df
= modf((double)_X
,&_Di
);
233 *_Y
= (long double)_Di
;
236 __CRT_INLINE
long double _chgsignl(long double _Number
) { return _chgsign(static_cast<double>(_Number
)); }
237 __CRT_INLINE
long double _copysignl(long double _Number
,long double _Sign
) { return _copysign(static_cast<double>(_Number
),static_cast<double>(_Sign
)); }
238 __CRT_INLINE
float frexpf(float _X
,int *_Y
) { return ((float)frexp((double)_X
,_Y
)); }
240 __CRT_INLINE
float __cdecl
fabsf (float x
)
243 __asm__ ("fabs;" : "=t" (res
) : "0" (x
));
246 __CRT_INLINE
float __cdecl
ldexpf (float x
, int expn
) { return (float) ldexp (x
, expn
); }
248 __CRT_INLINE
float acosf(float _X
) { return ((float)acos((double)_X
)); }
249 __CRT_INLINE
float asinf(float _X
) { return ((float)asin((double)_X
)); }
250 __CRT_INLINE
float atanf(float _X
) { return ((float)atan((double)_X
)); }
251 __CRT_INLINE
float atan2f(float _X
,float _Y
) { return ((float)atan2((double)_X
,(double)_Y
)); }
252 __CRT_INLINE
float ceilf(float _X
) { return ((float)ceil((double)_X
)); }
253 __CRT_INLINE
float cosf(float _X
) { return ((float)cos((double)_X
)); }
254 __CRT_INLINE
float coshf(float _X
) { return ((float)cosh((double)_X
)); }
255 __CRT_INLINE
float expf(float _X
) { return ((float)exp((double)_X
)); }
256 __CRT_INLINE
float floorf(float _X
) { return ((float)floor((double)_X
)); }
257 __CRT_INLINE
float fmodf(float _X
,float _Y
) { return ((float)fmod((double)_X
,(double)_Y
)); }
258 __CRT_INLINE
float logf(float _X
) { return ((float)log((double)_X
)); }
259 __CRT_INLINE
float log10f(float _X
) { return ((float)log10((double)_X
)); }
260 __CRT_INLINE
float modff(float _X
,float *_Y
) {
261 double _Di
,_Df
= modf((double)_X
,&_Di
);
265 __CRT_INLINE
float powf(float _X
,float _Y
) { return ((float)pow((double)_X
,(double)_Y
)); }
266 __CRT_INLINE
float sinf(float _X
) { return ((float)sin((double)_X
)); }
267 __CRT_INLINE
float sinhf(float _X
) { return ((float)sinh((double)_X
)); }
268 __CRT_INLINE
float sqrtf(float _X
) { return ((float)sqrt((double)_X
)); }
269 __CRT_INLINE
float tanf(float _X
) { return ((float)tan((double)_X
)); }
270 __CRT_INLINE
float tanhf(float _X
) { return ((float)tanh((double)_X
)); }
276 #define matherr _matherr
279 /* double __cdecl cabs(struct _complex _X); */
280 double __cdecl
hypot(double _X
,double _Y
);
281 _CRTIMP
double __cdecl
j0(double _X
);
282 _CRTIMP
double __cdecl
j1(double _X
);
283 _CRTIMP
double __cdecl
jn(int _X
,double _Y
);
284 _CRTIMP
double __cdecl
y0(double _X
);
285 _CRTIMP
double __cdecl
y1(double _X
);
286 _CRTIMP
double __cdecl
yn(int _X
,double _Y
);
290 #if (defined (__STDC_VERSION__) && __STDC_VERSION__ >= 199901L) \
291 || !defined __STRICT_ANSI__ || defined __GLIBCPP__
293 #define NAN (0.0F/0.0F)
294 #define HUGE_VALF (1.0F/0.0F)
295 #define HUGE_VALL (1.0L/0.0L)
296 #define INFINITY (1.0F/0.0F)
299 #define FP_NAN 0x0100
300 #define FP_NORMAL 0x0400
301 #define FP_INFINITE (FP_NAN | FP_NORMAL)
302 #define FP_ZERO 0x4000
303 #define FP_SUBNORMAL (FP_NORMAL | FP_ZERO)
304 /* 0x0200 is signbit mask */
308 We can't __CRT_INLINE float or double, because we want to ensure truncation
309 to semantic type before classification.
310 (A normal long double value might become subnormal when
311 converted to double, and zero when converted to float.)
314 extern int __cdecl
__fpclassifyf (float);
315 extern int __cdecl
__fpclassify (double);
317 __CRT_INLINE
int __cdecl
__fpclassifyl (long double x
){
319 __asm__ ("fxam; fstsw %%ax;" : "=a" (sw
): "t" (x
));
320 return sw
& (FP_NAN
| FP_NORMAL
| FP_ZERO
);
323 #define fpclassify(x) (sizeof (x) == sizeof (float) ? __fpclassifyf (x) \
324 : sizeof (x) == sizeof (double) ? __fpclassify (x) \
328 #define isfinite(x) ((fpclassify(x) & FP_NAN) == 0)
331 #define isinf(x) (fpclassify(x) == FP_INFINITE)
334 /* We don't need to worry about trucation here:
335 A NaN stays a NaN. */
337 __CRT_INLINE
int __cdecl
__isnan (double _x
)
341 "fstsw %%ax": "=a" (sw
) : "t" (_x
));
342 return (sw
& (FP_NAN
| FP_NORMAL
| FP_INFINITE
| FP_ZERO
| FP_SUBNORMAL
))
346 __CRT_INLINE
int __cdecl
__isnanf (float _x
)
350 "fstsw %%ax": "=a" (sw
) : "t" (_x
));
351 return (sw
& (FP_NAN
| FP_NORMAL
| FP_INFINITE
| FP_ZERO
| FP_SUBNORMAL
))
355 __CRT_INLINE
int __cdecl
__isnanl (long double _x
)
359 "fstsw %%ax": "=a" (sw
) : "t" (_x
));
360 return (sw
& (FP_NAN
| FP_NORMAL
| FP_INFINITE
| FP_ZERO
| FP_SUBNORMAL
))
365 #define isnan(x) (sizeof (x) == sizeof (float) ? __isnanf (x) \
366 : sizeof (x) == sizeof (double) ? __isnan (x) \
370 #define isnormal(x) (fpclassify(x) == FP_NORMAL)
372 /* 7.12.3.6 The signbit macro */
373 __CRT_INLINE
int __cdecl
__signbit (double x
) {
375 __asm__ ( "fxam; fstsw %%ax;": "=a" (stw
) : "t" (x
));
379 __CRT_INLINE
int __cdecl
__signbitf (float x
) {
381 __asm__ ("fxam; fstsw %%ax;": "=a" (stw
) : "t" (x
));
385 __CRT_INLINE
int __cdecl
__signbitl (long double x
) {
387 __asm__ ("fxam; fstsw %%ax;": "=a" (stw
) : "t" (x
));
391 #define signbit(x) (sizeof (x) == sizeof (float) ? __signbitf (x) \
392 : sizeof (x) == sizeof (double) ? __signbit (x) \
395 extern double __cdecl
exp2(double);
396 extern float __cdecl
exp2f(float);
397 extern long double __cdecl
exp2l(long double);
399 #define FP_ILOGB0 ((int)0x80000000)
400 #define FP_ILOGBNAN ((int)0x80000000)
401 extern int __cdecl
ilogb (double);
402 extern int __cdecl
ilogbf (float);
403 extern int __cdecl
ilogbl (long double);
405 extern double __cdecl
log1p(double);
406 extern float __cdecl
log1pf(float);
407 extern long double __cdecl
log1pl(long double);
409 extern double __cdecl
log2 (double);
410 extern float __cdecl
log2f (float);
411 extern long double __cdecl
log2l (long double);
413 extern double __cdecl
logb (double);
414 extern float __cdecl
logbf (float);
415 extern long double __cdecl
logbl (long double);
417 __CRT_INLINE
double __cdecl
logb (double x
)
420 __asm__ ("fxtract\n\t"
421 "fstp %%st" : "=t" (res
) : "0" (x
));
425 __CRT_INLINE
float __cdecl
logbf (float x
)
428 __asm__ ("fxtract\n\t"
429 "fstp %%st" : "=t" (res
) : "0" (x
));
433 __CRT_INLINE
long double __cdecl
logbl (long double x
)
436 __asm__ ("fxtract\n\t"
437 "fstp %%st" : "=t" (res
) : "0" (x
));
441 extern long double __cdecl
modfl (long double, long double*);
444 extern double __cdecl
scalbn (double, int);
445 extern float __cdecl
scalbnf (float, int);
446 extern long double __cdecl
scalbnl (long double, int);
448 extern double __cdecl
scalbln (double, long);
449 extern float __cdecl
scalblnf (float, long);
450 extern long double __cdecl
scalblnl (long double, long);
453 /* Implementations adapted from Cephes versions */
454 extern double __cdecl
cbrt (double);
455 extern float __cdecl
cbrtf (float);
456 extern long double __cdecl
cbrtl (long double);
458 __CRT_INLINE
float __cdecl
hypotf (float x
, float y
)
459 { return (float) hypot (x
, y
);}
460 extern long double __cdecl
hypotl (long double, long double);
462 extern long double __cdecl
powl (long double, long double);
463 extern long double __cdecl
expl(long double);
464 extern long double __cdecl
expm1l(long double);
465 extern long double __cdecl
coshl(long double);
466 extern long double __cdecl
fabsl (long double);
467 extern long double __cdecl
acosl(long double);
468 extern long double __cdecl
asinl(long double);
469 extern long double __cdecl
atanl(long double);
470 extern long double __cdecl
atan2l(long double,long double);
471 extern long double __cdecl
sinhl(long double);
472 extern long double __cdecl
tanhl(long double);
474 /* 7.12.8.1 The erf functions */
475 extern double __cdecl
erf (double);
476 extern float __cdecl
erff (float);
478 extern long double __cdecl erfl (long double);
481 /* 7.12.8.2 The erfc functions */
482 extern double __cdecl
erfc (double);
483 extern float __cdecl
erfcf (float);
485 extern long double __cdecl erfcl (long double);
488 /* 7.12.8.3 The lgamma functions */
489 extern double __cdecl
lgamma (double);
490 extern float __cdecl
lgammaf (float);
491 extern long double __cdecl
lgammal (long double);
493 /* 7.12.8.4 The tgamma functions */
494 extern double __cdecl
tgamma (double);
495 extern float __cdecl
tgammaf (float);
496 extern long double __cdecl
tgammal (long double);
498 extern long double __cdecl
ceill (long double);
499 extern long double __cdecl
floorl (long double);
500 extern long double __cdecl
frexpl(long double,int *);
501 extern long double __cdecl
log10l(long double);
502 extern long double __cdecl
logl(long double);
503 extern long double __cdecl
cosl(long double);
504 extern long double __cdecl
sinl(long double);
505 extern long double __cdecl
tanl(long double);
506 extern long double sqrtl(long double);
509 extern double __cdecl
nearbyint ( double);
510 extern float __cdecl
nearbyintf (float);
511 extern long double __cdecl
nearbyintl (long double);
514 /* round, using fpu control word settings */
515 __CRT_INLINE
double __cdecl
rint (double x
)
521 "fstl %0\n" : "=m" (retval
) : "m" (x
));
525 __CRT_INLINE
float __cdecl
rintf (float x
)
531 "fsts %0\n" : "=m" (retval
) : "m" (x
));
535 __CRT_INLINE
long double __cdecl
rintl (long double x
)
541 "fstt %0\n" : "=m" (retval
) : "m" (x
));
546 __CRT_INLINE
long __cdecl
lrint (double x
)
549 __asm__ __volatile__ \
551 "fistpl %0" : "=m" (retval
) : "m" (x
)); \
555 __CRT_INLINE
long __cdecl
lrintf (float x
)
558 __asm__ __volatile__ \
560 "fistpl %0" : "=m" (retval
) : "m" (x
)); \
564 __CRT_INLINE
long __cdecl
lrintl (long double x
)
567 __asm__ __volatile__ \
569 "fistpl %0" : "=m" (retval
) : "m" (x
)); \
573 __CRT_INLINE
long long __cdecl
llrint (double x
)
576 __asm__ __volatile__ \
578 "fistpll %0" : "=m" (retval
) : "m" (x
)); \
582 __CRT_INLINE
long long __cdecl
llrintf (float x
)
585 __asm__ __volatile__ \
587 "fistpll %0" : "=m" (retval
) : "m" (x
)); \
591 __CRT_INLINE
long long __cdecl
llrintl (long double x
)
594 __asm__ __volatile__ \
596 "fistpll %0" : "=m" (retval
) : "m" (x
)); \
601 /* round away from zero, regardless of fpu control word settings */
602 extern double __cdecl
round (double);
603 extern float __cdecl
roundf (float);
604 extern long double __cdecl
roundl (long double);
607 extern long __cdecl
lround (double);
608 extern long __cdecl
lroundf (float);
609 extern long __cdecl
lroundl (long double);
611 extern long long __cdecl
llround (double);
612 extern long long __cdecl
llroundf (float);
613 extern long long __cdecl
llroundl (long double);
616 /* round towards zero, regardless of fpu control word settings */
617 extern double __cdecl
trunc (double);
618 extern float __cdecl
truncf (float);
619 extern long double __cdecl
truncl (long double);
621 extern long double __cdecl
fmodl (long double, long double);
624 extern double __cdecl
remainder (double, double);
625 extern float __cdecl
remainderf (float, float);
626 extern long double __cdecl
remainderl (long double, long double);
629 extern double __cdecl
remquo(double, double, int *);
630 extern float __cdecl
remquof(float, float, int *);
631 extern long double __cdecl
remquol(long double, long double, int *);
634 extern double __cdecl
copysign (double, double); /* in libmoldname.a */
635 extern float __cdecl
copysignf (float, float);
636 extern long double __cdecl
copysignl (long double, long double);
638 /* 7.12.11.2 Return a NaN */
639 extern double __cdecl
nan(const char *tagp
);
640 extern float __cdecl
nanf(const char *tagp
);
641 extern long double __cdecl
nanl(const char *tagp
);
643 #ifndef __STRICT_ANSI__
644 #define _nan() nan("")
645 #define _nanf() nanf("")
646 #define _nanl() nanl("")
650 extern double __cdecl
nextafter (double, double); /* in libmoldname.a */
651 extern float __cdecl
nextafterf (float, float);
652 extern long double __cdecl
nextafterl (long double, long double);
654 /* 7.12.11.4 The nexttoward functions: TODO */
657 /* x > y ? (x - y) : 0.0 */
658 extern double __cdecl
fdim (double x
, double y
);
659 extern float __cdecl
fdimf (float x
, float y
);
660 extern long double __cdecl
fdiml (long double x
, long double y
);
663 NaN arguments are treated as missing data: if one argument is a NaN
664 and the other numeric, then these functions choose the numeric
668 extern double __cdecl
fmax (double, double);
669 extern float __cdecl
fmaxf (float, float);
670 extern long double __cdecl
fmaxl (long double, long double);
673 extern double __cdecl
fmin (double, double);
674 extern float __cdecl
fminf (float, float);
675 extern long double __cdecl
fminl (long double, long double);
678 /* return x * y + z as a ternary op */
679 extern double __cdecl
fma (double, double, double);
680 extern float __cdecl
fmaf (float, float, float);
681 extern long double __cdecl
fmal (long double, long double, long double);
684 #if 0 // gr: duplicate, see below
687 * With these functions, comparisons involving quiet NaNs set the FP
688 * condition code to "unordered". The IEEE floating-point spec
689 * dictates that the result of floating-point comparisons should be
690 * false whenever a NaN is involved, with the exception of the != op,
691 * which always returns true: yes, (NaN != NaN) is true).
696 #define isgreater(x, y) __builtin_isgreater(x, y)
697 #define isgreaterequal(x, y) __builtin_isgreaterequal(x, y)
698 #define isless(x, y) __builtin_isless(x, y)
699 #define islessequal(x, y) __builtin_islessequal(x, y)
700 #define islessgreater(x, y) __builtin_islessgreater(x, y)
701 #define isunordered(x, y) __builtin_isunordered(x, y)
705 __CRT_INLINE
int __cdecl
706 __fp_unordered_compare (long double x
, long double y
){
707 unsigned short retval
;
708 __asm__ ("fucom %%st(1);"
709 "fnstsw;": "=a" (retval
) : "t" (x
), "u" (y
));
713 #define isgreater(x, y) ((__fp_unordered_compare(x, y) \
715 #define isless(x, y) ((__fp_unordered_compare (y, x) \
717 #define isgreaterequal(x, y) ((__fp_unordered_compare (x, y) \
719 #define islessequal(x, y) ((__fp_unordered_compare(y, x) \
721 #define islessgreater(x, y) ((__fp_unordered_compare(x, y) \
722 & FP_SUBNORMAL) == 0)
723 #define isunordered(x, y) ((__fp_unordered_compare(x, y) \
730 #endif /* __STDC_VERSION__ >= 199901L */
731 #endif /* __NO_ISOCEXT */
736 template<class _Ty
> inline _Ty
_Pow_int(_Ty _X
,int _Y
) {
738 if(_Y
>= 0) _N
= (unsigned int)_Y
;
739 else _N
= (unsigned int)(-_Y
);
740 for(_Ty _Z
= _Ty(1);;_X
*= _X
) {
741 if((_N
& 1)!=0) _Z
*= _X
;
742 if((_N
>>= 1)==0) return (_Y
< 0 ? _Ty(1) / _Z
: _Z
);
752 * With these functions, comparisons involving quiet NaNs set the FP
753 * condition code to "unordered". The IEEE floating-point spec
754 * dictates that the result of floating-point comparisons should be
755 * false whenever a NaN is involved, with the exception of the != op,
756 * which always returns true: yes, (NaN != NaN) is true).
761 #define isgreater(x, y) __builtin_isgreater(x, y)
762 #define isgreaterequal(x, y) __builtin_isgreaterequal(x, y)
763 #define isless(x, y) __builtin_isless(x, y)
764 #define islessequal(x, y) __builtin_islessequal(x, y)
765 #define islessgreater(x, y) __builtin_islessgreater(x, y)
766 #define isunordered(x, y) __builtin_isunordered(x, y)
770 __CRT_INLINE
int __cdecl
771 __fp_unordered_compare (long double x
, long double y
){
772 unsigned short retval
;
773 __asm__ ("fucom %%st(1);"
774 "fnstsw;": "=a" (retval
) : "t" (x
), "u" (y
));
778 #define isgreater(x, y) ((__fp_unordered_compare(x, y) \
780 #define isless(x, y) ((__fp_unordered_compare (y, x) \
782 #define isgreaterequal(x, y) ((__fp_unordered_compare (x, y) \
784 #define islessequal(x, y) ((__fp_unordered_compare(y, x) \
786 #define islessgreater(x, y) ((__fp_unordered_compare(x, y) \
787 & FP_SUBNORMAL) == 0)
788 #define isunordered(x, y) ((__fp_unordered_compare(x, y) \
793 #endif /* End _MATH_H_ */