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
){
322 "fstsw %0 \n" : "=a" (sw
): "m" (x
));
323 return sw
& (FP_NAN
| FP_NORMAL
| FP_ZERO
);
326 #define fpclassify(x) (sizeof (x) == sizeof (float) ? __fpclassifyf (x) \
327 : sizeof (x) == sizeof (double) ? __fpclassify (x) \
331 #define isfinite(x) ((fpclassify(x) & FP_NAN) == 0)
334 #define isinf(x) (fpclassify(x) == FP_INFINITE)
337 /* We don't need to worry about trucation here:
338 A NaN stays a NaN. */
340 __CRT_INLINE
int __cdecl
__isnan (double _x
)
346 "fstsw %0 \n" : "=a" (sw
) : "m" (_x
));
348 return (sw
& (FP_NAN
| FP_NORMAL
| FP_INFINITE
| FP_ZERO
| FP_SUBNORMAL
))
352 __CRT_INLINE
int __cdecl
__isnanf (float _x
)
358 "fstsw %0 \n" : "=a" (sw
) : "m" (_x
));
360 return (sw
& (FP_NAN
| FP_NORMAL
| FP_INFINITE
| FP_ZERO
| FP_SUBNORMAL
))
364 __CRT_INLINE
int __cdecl
__isnanl (long double _x
)
370 "fstsw %0 \n" : "=a" (sw
) : "m" (_x
));
372 return (sw
& (FP_NAN
| FP_NORMAL
| FP_INFINITE
| FP_ZERO
| FP_SUBNORMAL
))
377 #define isnan(x) (sizeof (x) == sizeof (float) ? __isnanf (x) \
378 : sizeof (x) == sizeof (double) ? __isnan (x) \
382 #define isnormal(x) (fpclassify(x) == FP_NORMAL)
384 /* 7.12.3.6 The signbit macro */
385 __CRT_INLINE
int __cdecl
__signbit (double x
) {
390 "fstsw %0 \n" : "=a" (stw
) : "m" (x
));
394 __CRT_INLINE
int __cdecl
__signbitf (float x
) {
399 "fstsw %0 \n" : "=a" (stw
) : "m" (x
));
403 __CRT_INLINE
int __cdecl
__signbitl (long double x
) {
408 "fstsw %0 \n" : "=a" (stw
) : "m" (x
));
412 #define signbit(x) (sizeof (x) == sizeof (float) ? __signbitf (x) \
413 : sizeof (x) == sizeof (double) ? __signbit (x) \
416 extern double __cdecl
exp2(double);
417 extern float __cdecl
exp2f(float);
418 extern long double __cdecl
exp2l(long double);
420 #define FP_ILOGB0 ((int)0x80000000)
421 #define FP_ILOGBNAN ((int)0x80000000)
422 extern int __cdecl
ilogb (double);
423 extern int __cdecl
ilogbf (float);
424 extern int __cdecl
ilogbl (long double);
426 extern double __cdecl
log1p(double);
427 extern float __cdecl
log1pf(float);
428 extern long double __cdecl
log1pl(long double);
430 extern double __cdecl
log2 (double);
431 extern float __cdecl
log2f (float);
432 extern long double __cdecl
log2l (long double);
434 extern double __cdecl
logb (double);
435 extern float __cdecl
logbf (float);
436 extern long double __cdecl
logbl (long double);
438 __CRT_INLINE
double __cdecl
logb (double x
)
441 __asm__ ("fxtract\n\t"
442 "fstp %%st" : "=t" (res
) : "0" (x
));
446 __CRT_INLINE
float __cdecl
logbf (float x
)
449 __asm__ ("fxtract\n\t"
450 "fstp %%st" : "=t" (res
) : "0" (x
));
454 __CRT_INLINE
long double __cdecl
logbl (long double x
)
457 __asm__ ("fxtract\n\t"
458 "fstp %%st" : "=t" (res
) : "0" (x
));
462 extern long double __cdecl
modfl (long double, long double*);
465 extern double __cdecl
scalbn (double, int);
466 extern float __cdecl
scalbnf (float, int);
467 extern long double __cdecl
scalbnl (long double, int);
469 extern double __cdecl
scalbln (double, long);
470 extern float __cdecl
scalblnf (float, long);
471 extern long double __cdecl
scalblnl (long double, long);
474 /* Implementations adapted from Cephes versions */
475 extern double __cdecl
cbrt (double);
476 extern float __cdecl
cbrtf (float);
477 extern long double __cdecl
cbrtl (long double);
479 __CRT_INLINE
float __cdecl
hypotf (float x
, float y
)
480 { return (float) hypot (x
, y
);}
481 extern long double __cdecl
hypotl (long double, long double);
483 extern long double __cdecl
powl (long double, long double);
484 extern long double __cdecl
expl(long double);
485 extern long double __cdecl
expm1l(long double);
486 extern long double __cdecl
coshl(long double);
487 extern long double __cdecl
fabsl (long double);
488 extern long double __cdecl
acosl(long double);
489 extern long double __cdecl
asinl(long double);
490 extern long double __cdecl
atanl(long double);
491 extern long double __cdecl
atan2l(long double,long double);
492 extern long double __cdecl
sinhl(long double);
493 extern long double __cdecl
tanhl(long double);
495 /* 7.12.8.1 The erf functions */
496 extern double __cdecl
erf (double);
497 extern float __cdecl
erff (float);
499 extern long double __cdecl erfl (long double);
502 /* 7.12.8.2 The erfc functions */
503 extern double __cdecl
erfc (double);
504 extern float __cdecl
erfcf (float);
506 extern long double __cdecl erfcl (long double);
509 /* 7.12.8.3 The lgamma functions */
510 extern double __cdecl
lgamma (double);
511 extern float __cdecl
lgammaf (float);
512 extern long double __cdecl
lgammal (long double);
514 /* 7.12.8.4 The tgamma functions */
515 extern double __cdecl
tgamma (double);
516 extern float __cdecl
tgammaf (float);
517 extern long double __cdecl
tgammal (long double);
519 extern long double __cdecl
ceill (long double);
520 extern long double __cdecl
floorl (long double);
521 extern long double __cdecl
frexpl(long double,int *);
522 extern long double __cdecl
log10l(long double);
523 extern long double __cdecl
logl(long double);
524 extern long double __cdecl
cosl(long double);
525 extern long double __cdecl
sinl(long double);
526 extern long double __cdecl
tanl(long double);
527 extern long double sqrtl(long double);
530 extern double __cdecl
nearbyint ( double);
531 extern float __cdecl
nearbyintf (float);
532 extern long double __cdecl
nearbyintl (long double);
535 /* round, using fpu control word settings */
536 __CRT_INLINE
double __cdecl
rint (double x
)
542 "fstl %0\n" : "=m" (retval
) : "m" (x
));
546 __CRT_INLINE
float __cdecl
rintf (float x
)
552 "fsts %0\n" : "=m" (retval
) : "m" (x
));
556 __CRT_INLINE
long double __cdecl
rintl (long double x
)
562 "fstt %0\n" : "=m" (retval
) : "m" (x
));
567 __CRT_INLINE
long __cdecl
lrint (double x
)
570 __asm__ __volatile__ \
572 "fistpl %0" : "=m" (retval
) : "m" (x
)); \
576 __CRT_INLINE
long __cdecl
lrintf (float x
)
579 __asm__ __volatile__ \
581 "fistpl %0" : "=m" (retval
) : "m" (x
)); \
585 __CRT_INLINE
long __cdecl
lrintl (long double x
)
588 __asm__ __volatile__ \
590 "fistpl %0" : "=m" (retval
) : "m" (x
)); \
594 __CRT_INLINE
long long __cdecl
llrint (double x
)
597 __asm__ __volatile__ \
599 "fistpll %0" : "=m" (retval
) : "m" (x
)); \
603 __CRT_INLINE
long long __cdecl
llrintf (float x
)
606 __asm__ __volatile__ \
608 "fistpll %0" : "=m" (retval
) : "m" (x
)); \
612 __CRT_INLINE
long long __cdecl
llrintl (long double x
)
615 __asm__ __volatile__ \
617 "fistpll %0" : "=m" (retval
) : "m" (x
)); \
621 #define FE_TONEAREST 0x0000
622 #define FE_DOWNWARD 0x0400
623 #define FE_UPWARD 0x0800
624 #define FE_TOWARDZERO 0x0c00
626 __CRT_INLINE
double trunc (double _x
)
629 unsigned short saved_cw
;
630 unsigned short tmp_cw
;
631 __asm__ ("fnstcw %0;" : "=m" (saved_cw
)); /* save FPU control word */
632 tmp_cw
= (saved_cw
& ~(FE_TONEAREST
| FE_DOWNWARD
| FE_UPWARD
| FE_TOWARDZERO
))
634 __asm__ ("fldcw %0;" : : "m" (tmp_cw
));
637 "fstl %0;" : "=m" (retval
) : "m" (_x
)); /* round towards zero */
638 __asm__ ("fldcw %0;" : : "m" (saved_cw
) ); /* restore saved control word */
643 /* round away from zero, regardless of fpu control word settings */
644 extern double __cdecl
round (double);
645 extern float __cdecl
roundf (float);
646 extern long double __cdecl
roundl (long double);
649 extern long __cdecl
lround (double);
650 extern long __cdecl
lroundf (float);
651 extern long __cdecl
lroundl (long double);
653 extern long long __cdecl
llround (double);
654 extern long long __cdecl
llroundf (float);
655 extern long long __cdecl
llroundl (long double);
658 /* round towards zero, regardless of fpu control word settings */
659 extern double __cdecl
trunc (double);
660 extern float __cdecl
truncf (float);
661 extern long double __cdecl
truncl (long double);
663 extern long double __cdecl
fmodl (long double, long double);
666 extern double __cdecl
remainder (double, double);
667 extern float __cdecl
remainderf (float, float);
668 extern long double __cdecl
remainderl (long double, long double);
671 extern double __cdecl
remquo(double, double, int *);
672 extern float __cdecl
remquof(float, float, int *);
673 extern long double __cdecl
remquol(long double, long double, int *);
676 extern double __cdecl
copysign (double, double); /* in libmoldname.a */
677 extern float __cdecl
copysignf (float, float);
678 extern long double __cdecl
copysignl (long double, long double);
680 /* 7.12.11.2 Return a NaN */
681 extern double __cdecl
nan(const char *tagp
);
682 extern float __cdecl
nanf(const char *tagp
);
683 extern long double __cdecl
nanl(const char *tagp
);
685 #ifndef __STRICT_ANSI__
686 #define _nan() nan("")
687 #define _nanf() nanf("")
688 #define _nanl() nanl("")
692 extern double __cdecl
nextafter (double, double); /* in libmoldname.a */
693 extern float __cdecl
nextafterf (float, float);
694 extern long double __cdecl
nextafterl (long double, long double);
696 /* 7.12.11.4 The nexttoward functions: TODO */
699 /* x > y ? (x - y) : 0.0 */
700 extern double __cdecl
fdim (double x
, double y
);
701 extern float __cdecl
fdimf (float x
, float y
);
702 extern long double __cdecl
fdiml (long double x
, long double y
);
705 NaN arguments are treated as missing data: if one argument is a NaN
706 and the other numeric, then these functions choose the numeric
710 extern double __cdecl
fmax (double, double);
711 extern float __cdecl
fmaxf (float, float);
712 extern long double __cdecl
fmaxl (long double, long double);
715 extern double __cdecl
fmin (double, double);
716 extern float __cdecl
fminf (float, float);
717 extern long double __cdecl
fminl (long double, long double);
720 /* return x * y + z as a ternary op */
721 extern double __cdecl
fma (double, double, double);
722 extern float __cdecl
fmaf (float, float, float);
723 extern long double __cdecl
fmal (long double, long double, long double);
726 #if 0 // gr: duplicate, see below
729 * With these functions, comparisons involving quiet NaNs set the FP
730 * condition code to "unordered". The IEEE floating-point spec
731 * dictates that the result of floating-point comparisons should be
732 * false whenever a NaN is involved, with the exception of the != op,
733 * which always returns true: yes, (NaN != NaN) is true).
738 #define isgreater(x, y) __builtin_isgreater(x, y)
739 #define isgreaterequal(x, y) __builtin_isgreaterequal(x, y)
740 #define isless(x, y) __builtin_isless(x, y)
741 #define islessequal(x, y) __builtin_islessequal(x, y)
742 #define islessgreater(x, y) __builtin_islessgreater(x, y)
743 #define isunordered(x, y) __builtin_isunordered(x, y)
747 __CRT_INLINE
int __cdecl
748 __fp_unordered_compare (long double x
, long double y
){
749 unsigned short retval
;
750 __asm__ ("fucom %%st(1);"
751 "fnstsw;": "=a" (retval
) : "t" (x
), "u" (y
));
755 #define isgreater(x, y) ((__fp_unordered_compare(x, y) \
757 #define isless(x, y) ((__fp_unordered_compare (y, x) \
759 #define isgreaterequal(x, y) ((__fp_unordered_compare (x, y) \
761 #define islessequal(x, y) ((__fp_unordered_compare(y, x) \
763 #define islessgreater(x, y) ((__fp_unordered_compare(x, y) \
764 & FP_SUBNORMAL) == 0)
765 #define isunordered(x, y) ((__fp_unordered_compare(x, y) \
772 #endif /* __STDC_VERSION__ >= 199901L */
773 #endif /* __NO_ISOCEXT */
778 template<class _Ty
> inline _Ty
_Pow_int(_Ty _X
,int _Y
) {
780 if(_Y
>= 0) _N
= (unsigned int)_Y
;
781 else _N
= (unsigned int)(-_Y
);
782 for(_Ty _Z
= _Ty(1);;_X
*= _X
) {
783 if((_N
& 1)!=0) _Z
*= _X
;
784 if((_N
>>= 1)==0) return (_Y
< 0 ? _Ty(1) / _Z
: _Z
);
794 * With these functions, comparisons involving quiet NaNs set the FP
795 * condition code to "unordered". The IEEE floating-point spec
796 * dictates that the result of floating-point comparisons should be
797 * false whenever a NaN is involved, with the exception of the != op,
798 * which always returns true: yes, (NaN != NaN) is true).
801 #endif /* End _MATH_H_ */