-
Notifications
You must be signed in to change notification settings - Fork 129
/
Copy pathmath.h
1584 lines (1342 loc) · 52.5 KB
/
math.h
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
402
403
404
405
406
407
408
409
410
411
412
413
414
415
416
417
418
419
420
421
422
423
424
425
426
427
428
429
430
431
432
433
434
435
436
437
438
439
440
441
442
443
444
445
446
447
448
449
450
451
452
453
454
455
456
457
458
459
460
461
462
463
464
465
466
467
468
469
470
471
472
473
474
475
476
477
478
479
480
481
482
483
484
485
486
487
488
489
490
491
492
493
494
495
496
497
498
499
500
501
502
503
504
505
506
507
508
509
510
511
512
513
514
515
516
517
518
519
520
521
522
523
524
525
526
527
528
529
530
531
532
533
534
535
536
537
538
539
540
541
542
543
544
545
546
547
548
549
550
551
552
553
554
555
556
557
558
559
560
561
562
563
564
565
566
567
568
569
570
571
572
573
574
575
576
577
578
579
580
581
582
583
584
585
586
587
588
589
590
591
592
593
594
595
596
597
598
599
600
601
602
603
604
605
606
607
608
609
610
611
612
613
614
615
616
617
618
619
620
621
622
623
624
625
626
627
628
629
630
631
632
633
634
635
636
637
638
639
640
641
642
643
644
645
646
647
648
649
650
651
652
653
654
655
656
657
658
659
660
661
662
663
664
665
666
667
668
669
670
671
672
673
674
675
676
677
678
679
680
681
682
683
684
685
686
687
688
689
690
691
692
693
694
695
696
697
698
699
700
701
702
703
704
705
706
707
708
709
710
711
712
713
714
715
716
717
718
719
720
721
722
723
724
725
726
727
728
729
730
731
732
733
734
735
736
737
738
739
740
741
742
743
744
745
746
747
748
749
750
751
752
753
754
755
756
757
758
759
760
761
762
763
764
765
766
767
768
769
770
771
772
773
774
775
776
777
778
779
780
781
782
783
784
785
786
787
788
789
790
791
792
793
794
795
796
797
798
799
800
801
802
803
804
805
806
807
808
809
810
811
812
813
814
815
816
817
818
819
820
821
822
823
824
825
826
827
828
829
830
831
832
833
834
835
836
837
838
839
840
841
842
843
844
845
846
847
848
849
850
851
852
853
854
855
856
857
858
859
860
861
862
863
864
865
866
867
868
869
870
871
872
873
874
875
876
877
878
879
880
881
882
883
884
885
886
887
888
889
890
891
892
893
894
895
896
897
898
899
900
901
902
903
904
905
906
907
908
909
910
911
912
913
914
915
916
917
918
919
920
921
922
923
924
925
926
927
928
929
930
931
932
933
934
935
936
937
938
939
940
941
942
943
944
945
946
947
948
949
950
951
952
953
954
955
956
957
958
959
960
961
962
963
964
965
966
967
968
969
970
971
972
973
974
975
976
977
978
979
980
981
982
983
984
985
986
987
988
989
990
991
992
993
994
995
996
997
998
999
1000
/**
* This file has no copyright assigned and is placed in the Public Domain.
* This file is part of the mingw-w64 runtime package.
* No warranty is given; refer to the file DISCLAIMER.PD within this package.
*/
#ifndef _MATH_H_
#define _MATH_H_
#ifdef __GNUC__
#pragma GCC system_header
#endif /* __GNUC__ */
#include <crtdefs.h>
struct _exception;
#pragma pack(push,_CRT_PACKING)
#define _DOMAIN 1 /* domain error in argument */
#define _SING 2 /* singularity */
#define _OVERFLOW 3 /* range overflow */
#define _UNDERFLOW 4 /* range underflow */
#define _TLOSS 5 /* total loss of precision */
#define _PLOSS 6 /* partial loss of precision */
#ifndef __STRICT_ANSI__
#ifndef NO_OLDNAMES
#define DOMAIN _DOMAIN
#define SING _SING
#define OVERFLOW _OVERFLOW
#define UNDERFLOW _UNDERFLOW
#define TLOSS _TLOSS
#define PLOSS _PLOSS
#endif
#endif
#if !defined(__STRICT_ANSI__) || defined(_POSIX_C_SOURCE) || defined(_POSIX_SOURCE) || defined(_XOPEN_SOURCE) || defined(_GNU_SOURCE) || defined(_BSD_SOURCE) || defined(_USE_MATH_DEFINES)
#define M_E 2.7182818284590452354
#define M_LOG2E 1.4426950408889634074
#define M_LOG10E 0.43429448190325182765
#define M_LN2 0.69314718055994530942
#define M_LN10 2.30258509299404568402
#define M_PI 3.14159265358979323846
#define M_PI_2 1.57079632679489661923
#define M_PI_4 0.78539816339744830962
#define M_1_PI 0.31830988618379067154
#define M_2_PI 0.63661977236758134308
#define M_2_SQRTPI 1.12837916709551257390
#define M_SQRT2 1.41421356237309504880
#define M_SQRT1_2 0.70710678118654752440
#endif
#ifndef __STRICT_ANSI__
/* See also float.h */
#ifndef __MINGW_FPCLASS_DEFINED
#define __MINGW_FPCLASS_DEFINED 1
/* IEEE 754 classication */
#define _FPCLASS_SNAN 0x0001 /* Signaling "Not a Number" */
#define _FPCLASS_QNAN 0x0002 /* Quiet "Not a Number" */
#define _FPCLASS_NINF 0x0004 /* Negative Infinity */
#define _FPCLASS_NN 0x0008 /* Negative Normal */
#define _FPCLASS_ND 0x0010 /* Negative Denormal */
#define _FPCLASS_NZ 0x0020 /* Negative Zero */
#define _FPCLASS_PZ 0x0040 /* Positive Zero */
#define _FPCLASS_PD 0x0080 /* Positive Denormal */
#define _FPCLASS_PN 0x0100 /* Positive Normal */
#define _FPCLASS_PINF 0x0200 /* Positive Infinity */
#endif
#endif
#ifndef RC_INVOKED
#ifndef __mingw_types_compatible_p
#ifdef __cplusplus
extern "C++" {
template <typename type1, typename type2> struct __mingw_types_compatible_p {
static const bool result = false;
};
template <typename type1> struct __mingw_types_compatible_p<type1, type1> {
static const bool result = true;
};
template <typename type1> struct __mingw_types_compatible_p<const type1, type1> {
static const bool result = true;
};
template <typename type1> struct __mingw_types_compatible_p<type1, const type1> {
static const bool result = true;
};
}
#define __mingw_types_compatible_p(type1, type2) __mingw_types_compatible_p <type1, type2>::result
#else
#define __mingw_types_compatible_p(type1, type2) __builtin_types_compatible_p (type1, type2)
#endif
#endif
#ifndef __mingw_choose_expr
#ifdef __cplusplus
#define __mingw_choose_expr(C, E1, E2) ((C) ? E1 : E2)
#else
#define __mingw_choose_expr __builtin_choose_expr
#endif
#endif
#ifdef __cplusplus
extern "C" {
#endif
#ifndef __MINGW_SOFTMATH
#define __MINGW_SOFTMATH
/* IEEE float/double type shapes. */
typedef union __mingw_dbl_type_t {
double x;
unsigned long long val;
__C89_NAMELESS struct {
unsigned int low, high;
} lh;
} __mingw_dbl_type_t;
typedef union __mingw_flt_type_t {
float x;
unsigned int val;
} __mingw_flt_type_t;
typedef union __mingw_ldbl_type_t
{
long double x;
__C89_NAMELESS struct {
unsigned int low, high;
int sign_exponent : 16;
int res1 : 16;
int res0 : 32;
} lh;
} __mingw_ldbl_type_t;
typedef union __mingw_fp_types_t
{
long double *ld;
double *d;
float *f;
__mingw_ldbl_type_t *ldt;
__mingw_dbl_type_t *dt;
__mingw_flt_type_t *ft;
} __mingw_fp_types_t;
#endif
#ifndef _HUGE
extern double * __MINGW_IMP_SYMBOL(_HUGE);
#define _HUGE (* __MINGW_IMP_SYMBOL(_HUGE))
#endif
#ifdef __GNUC__
#define HUGE_VAL __builtin_huge_val()
#else
#define HUGE_VAL _HUGE
#endif /* __GNUC__ */
#ifndef _EXCEPTION_DEFINED
#define _EXCEPTION_DEFINED
struct _exception {
int type;
const char *name;
double arg1;
double arg2;
double retval;
};
void __mingw_raise_matherr (int typ, const char *name, double a1, double a2,
double rslt);
void __mingw_setusermatherr (int (__cdecl *)(struct _exception *));
_CRTIMP void __setusermatherr(int (__cdecl *)(struct _exception *));
#define __setusermatherr __mingw_setusermatherr
#endif
double __cdecl sin(double _X);
double __cdecl cos(double _X);
double __cdecl tan(double _X);
double __cdecl sinh(double _X);
double __cdecl cosh(double _X);
double __cdecl tanh(double _X);
double __cdecl asin(double _X);
double __cdecl acos(double _X);
double __cdecl atan(double _X);
double __cdecl atan2(double _Y,double _X);
double __cdecl exp(double _X);
double __cdecl log(double _X);
double __cdecl log10(double _X);
double __cdecl pow(double _X,double _Y);
double __cdecl sqrt(double _X);
double __cdecl ceil(double _X);
double __cdecl floor(double _X);
/* 7.12.7.2 The fabs functions: Double in C89 */
extern float __cdecl fabsf (float x);
extern long double __cdecl fabsl (long double);
extern double __cdecl fabs (double _X);
#ifndef __CRT__NO_INLINE
#if !defined (__ia64__)
__CRT_INLINE float __cdecl fabsf (float x)
{
#if defined(__x86_64__) || defined(__arm__)
return __builtin_fabsf (x);
#else
float res = 0.0F;
__asm__ __volatile__ ("fabs;" : "=t" (res) : "0" (x));
return res;
#endif
}
__CRT_INLINE long double __cdecl fabsl (long double x)
{
#ifdef __arm__
return __builtin_fabsl (x);
#else
long double res = 0.0l;
__asm__ __volatile__ ("fabs;" : "=t" (res) : "0" (x));
return res;
#endif
}
__CRT_INLINE double __cdecl fabs (double x)
{
#if defined(__x86_64__) || defined(__arm__)
return __builtin_fabs (x);
#else
double res = 0.0;
__asm__ __volatile__ ("fabs;" : "=t" (res) : "0" (x));
return res;
#endif
}
#endif
#endif
double __cdecl ldexp(double _X,int _Y);
double __cdecl frexp(double _X,int *_Y);
double __cdecl modf(double _X,double *_Y);
double __cdecl fmod(double _X,double _Y);
void __cdecl sincos (double __x, double *p_sin, double *p_cos);
void __cdecl sincosl (long double __x, long double *p_sin, long double *p_cos);
void __cdecl sincosf (float __x, float *p_sin, float *p_cos);
#ifndef _CRT_ABS_DEFINED
#define _CRT_ABS_DEFINED
int __cdecl abs(int _X);
long __cdecl labs(long _X);
#endif
#ifndef _CRT_ATOF_DEFINED
#define _CRT_ATOF_DEFINED
double __cdecl atof(const char *_String);
double __cdecl _atof_l(const char *_String,_locale_t _Locale);
#endif
#define EDOM 33
#define ERANGE 34
#if !defined(__STRICT_ANSI__) || defined(_POSIX_C_SOURCE) || defined(_POSIX_SOURCE) || defined(_XOPEN_SOURCE) || defined(_GNU_SOURCE) || defined(_BSD_SOURCE)
#ifndef _COMPLEX_DEFINED
#define _COMPLEX_DEFINED
struct _complex {
double x;
double y;
};
#endif
double __cdecl _cabs(struct _complex _ComplexA); /* Overriden to use our cabs. */
double __cdecl _hypot(double _X,double _Y);
_CRTIMP double __cdecl _j0(double _X);
_CRTIMP double __cdecl _j1(double _X);
_CRTIMP double __cdecl _jn(int _X,double _Y);
_CRTIMP double __cdecl _y0(double _X);
_CRTIMP double __cdecl _y1(double _X);
_CRTIMP double __cdecl _yn(int _X,double _Y);
#ifndef _CRT_MATHERR_DEFINED
#define _CRT_MATHERR_DEFINED
_CRTIMP int __cdecl _matherr (struct _exception *);
#endif
/* These are also declared in Mingw float.h; needed here as well to work
around GCC build issues. */
/* BEGIN FLOAT.H COPY */
/*
* IEEE recommended functions
*/
#ifndef _SIGN_DEFINED
#define _SIGN_DEFINED
_CRTIMP double __cdecl _chgsign (double _X);
_CRTIMP double __cdecl _copysign (double _Number,double _Sign);
_CRTIMP double __cdecl _logb (double);
_CRTIMP double __cdecl _nextafter (double, double);
_CRTIMP double __cdecl _scalb (double, long);
_CRTIMP int __cdecl _finite (double);
_CRTIMP int __cdecl _fpclass (double);
_CRTIMP int __cdecl _isnan (double);
#endif
/* END FLOAT.H COPY */
#if !defined(NO_OLDNAMES)
_CRTIMP double __cdecl j0 (double) __MINGW_ATTRIB_DEPRECATED_MSVC2005;
_CRTIMP double __cdecl j1 (double) __MINGW_ATTRIB_DEPRECATED_MSVC2005;
_CRTIMP double __cdecl jn (int, double) __MINGW_ATTRIB_DEPRECATED_MSVC2005;
_CRTIMP double __cdecl y0 (double) __MINGW_ATTRIB_DEPRECATED_MSVC2005;
_CRTIMP double __cdecl y1 (double) __MINGW_ATTRIB_DEPRECATED_MSVC2005;
_CRTIMP double __cdecl yn (int, double) __MINGW_ATTRIB_DEPRECATED_MSVC2005;
_CRTIMP double __cdecl chgsign (double);
/*
* scalb() is a GCC built-in.
* Exclude this _scalb() stub; the semantics are incompatible
* with the built-in implementation.
*
_CRTIMP double __cdecl scalb (double, long);
*
*/
_CRTIMP int __cdecl finite (double);
_CRTIMP int __cdecl fpclass (double);
#define FP_SNAN _FPCLASS_SNAN
#define FP_QNAN _FPCLASS_QNAN
#define FP_NINF _FPCLASS_NINF
#define FP_PINF _FPCLASS_PINF
#define FP_NDENORM _FPCLASS_ND
#define FP_PDENORM _FPCLASS_PD
#define FP_NZERO _FPCLASS_NZ
#define FP_PZERO _FPCLASS_PZ
#define FP_NNORM _FPCLASS_NN
#define FP_PNORM _FPCLASS_PN
#endif /* !defined (_NO_OLDNAMES) && !define (NO_OLDNAMES) */
#if(defined(_X86_) && !defined(__x86_64))
_CRTIMP int __cdecl _set_SSE2_enable(int _Flag);
#endif
#endif
#ifndef __NO_ISOCEXT
#if (defined (__STDC_VERSION__) && __STDC_VERSION__ >= 199901L) \
|| !defined __STRICT_ANSI__ || defined __cplusplus
#ifdef __GNUC__
#define HUGE_VALF __builtin_huge_valf()
#define HUGE_VALL __builtin_huge_vall()
#define INFINITY __builtin_inf()
#define NAN __builtin_nan("")
#else
extern const float __INFF;
#define HUGE_VALF __INFF
extern const long double __INFL;
#define HUGE_VALL __INFL
#define INFINITY HUGE_VALF
extern const double __QNAN;
#define NAN __QNAN
#endif /* __GNUC__ */
/* Use the compiler's builtin define for FLT_EVAL_METHOD to
set float_t and double_t. */
#if defined (__x86_64__) || defined(__FLT_EVAL_METHOD__)
# if defined (__x86_64__) || ( __FLT_EVAL_METHOD__== 0)
typedef float float_t;
typedef double double_t;
# elif (__FLT_EVAL_METHOD__ == 1)
typedef double float_t;
typedef double double_t;
# else /* (__FLT_EVAL_METHOD__ == 2) default ix87 FPU */
typedef long double float_t;
typedef long double double_t;
#endif
#else /* ix87 FPU default */
typedef long double float_t;
typedef long double double_t;
#endif
/* 7.12.3.1 */
/*
Return values for fpclassify.
These are based on Intel x87 fpu condition codes
in the high byte of status word and differ from
the return values for MS IEEE 754 extension _fpclass()
*/
#define FP_NAN 0x0100
#define FP_NORMAL 0x0400
#define FP_INFINITE (FP_NAN | FP_NORMAL)
#define FP_ZERO 0x4000
#define FP_SUBNORMAL (FP_NORMAL | FP_ZERO)
/* 0x0200 is signbit mask */
/*
We can't inline float or double, because we want to ensure truncation
to semantic type before classification.
(A normal long double value might become subnormal when
converted to double, and zero when converted to float.)
*/
extern int __cdecl __fpclassifyl (long double);
extern int __cdecl __fpclassifyf (float);
extern int __cdecl __fpclassify (double);
#ifndef __CRT__NO_INLINE
__CRT_INLINE int __cdecl __fpclassifyl (long double x) {
#if defined(__x86_64__) || defined(_AMD64_)
__mingw_fp_types_t hlp;
unsigned int e;
hlp.ld = &x;
e = hlp.ldt->lh.sign_exponent & 0x7fff;
if (!e)
{
unsigned int h = hlp.ldt->lh.high;
if (!(hlp.ldt->lh.low | h))
return FP_ZERO;
else if (!(h & 0x80000000))
return FP_SUBNORMAL;
}
else if (e == 0x7fff)
return (((hlp.ldt->lh.high & 0x7fffffff) | hlp.ldt->lh.low) == 0 ?
FP_INFINITE : FP_NAN);
return FP_NORMAL;
#elif defined(__arm__) || defined(_ARM_)
return __fpclassify(x);
#elif defined(__i386__) || defined(_X86_)
unsigned short sw;
__asm__ __volatile__ ("fxam; fstsw %%ax;" : "=a" (sw): "t" (x));
return sw & (FP_NAN | FP_NORMAL | FP_ZERO );
#endif
}
__CRT_INLINE int __cdecl __fpclassify (double x) {
#if defined(__x86_64__) || defined(_AMD64_) || defined(__arm__) || defined(_ARM_)
__mingw_fp_types_t hlp;
unsigned int l, h;
hlp.d = &x;
h = hlp.ldt->lh.high;
l = hlp.ldt->lh.low | (h & 0xfffff);
h &= 0x7ff00000;
if ((h | l) == 0)
return FP_ZERO;
if (!h)
return FP_SUBNORMAL;
if (h == 0x7ff00000)
return (l ? FP_NAN : FP_INFINITE);
return FP_NORMAL;
#elif defined(__i386__) || defined(_X86_)
unsigned short sw;
__asm__ __volatile__ ("fxam; fstsw %%ax;" : "=a" (sw): "t" (x));
return sw & (FP_NAN | FP_NORMAL | FP_ZERO );
#endif
}
__CRT_INLINE int __cdecl __fpclassifyf (float x) {
#if defined(__x86_64__) || defined(_AMD64_) || defined(__arm__) || defined(_ARM_)
__mingw_fp_types_t hlp;
hlp.f = &x;
hlp.ft->val &= 0x7fffffff;
if (hlp.ft->val == 0)
return FP_ZERO;
if (hlp.ft->val < 0x800000)
return FP_SUBNORMAL;
if (hlp.ft->val >= 0x7f800000)
return (hlp.ft->val > 0x7f800000 ? FP_NAN : FP_INFINITE);
return FP_NORMAL;
#elif defined(__i386__) || defined(_X86_)
unsigned short sw;
__asm__ __volatile__ ("fxam; fstsw %%ax;" : "=a" (sw): "t" (x));
return sw & (FP_NAN | FP_NORMAL | FP_ZERO );
#endif
}
#endif
#ifdef __STDC_WANT_DEC_FP__
#define __dfp_expansion(__call,__fin,x) \
__mingw_choose_expr ( \
__mingw_types_compatible_p (__typeof__ (x), _Decimal32), \
__call##d32(x), \
__mingw_choose_expr ( \
__mingw_types_compatible_p (__typeof__ (x), _Decimal64), \
__call##d64(x), \
__mingw_choose_expr ( \
__mingw_types_compatible_p (__typeof__ (x), _Decimal128), \
__call##d128(x), \
__fin)))
#else
#define __dfp_expansion(__call,__fin,x) __fin
#endif
#define fpclassify(x) \
__mingw_choose_expr ( \
__mingw_types_compatible_p (__typeof__ (x), double), \
__fpclassify(x), \
__mingw_choose_expr ( \
__mingw_types_compatible_p (__typeof__ (x), float), \
__fpclassifyf(x), \
__mingw_choose_expr ( \
__mingw_types_compatible_p (__typeof__ (x), long double), \
__fpclassifyl(x), \
__dfp_expansion(__fpclassify,(__builtin_trap(),0),x))))
/* 7.12.3.2 */
#define isfinite(x) ((fpclassify(x) & FP_NAN) == 0)
/* 7.12.3.3 */
#define isinf(x) (fpclassify(x) == FP_INFINITE)
/* 7.12.3.4 */
/* We don't need to worry about truncation here:
A NaN stays a NaN. */
extern int __cdecl __isnan (double);
extern int __cdecl __isnanf (float);
extern int __cdecl __isnanl (long double);
#ifndef __CRT__NO_INLINE
__CRT_INLINE int __cdecl __isnan (double _x)
{
#if defined(__x86_64__) || defined(_AMD64_) || defined(__arm__) || defined(_ARM_)
__mingw_fp_types_t hlp;
int l, h;
hlp.d = &_x;
l = hlp.dt->lh.low;
h = hlp.dt->lh.high & 0x7fffffff;
h |= (unsigned int) (l | -l) >> 31;
h = 0x7ff00000 - h;
return (int) ((unsigned int) h) >> 31;
#elif defined(__i386__) || defined(_X86_)
unsigned short sw;
__asm__ __volatile__ ("fxam;"
"fstsw %%ax": "=a" (sw) : "t" (_x));
return (sw & (FP_NAN | FP_NORMAL | FP_INFINITE | FP_ZERO | FP_SUBNORMAL))
== FP_NAN;
#endif
}
__CRT_INLINE int __cdecl __isnanf (float _x)
{
#if defined(__x86_64__) || defined(_AMD64_) || defined(__arm__) || defined(_ARM_)
__mingw_fp_types_t hlp;
int i;
hlp.f = &_x;
i = hlp.ft->val & 0x7fffffff;
i = 0x7f800000 - i;
return (int) (((unsigned int) i) >> 31);
#elif defined(__i386__) || defined(_X86_)
unsigned short sw;
__asm__ __volatile__ ("fxam;"
"fstsw %%ax": "=a" (sw) : "t" (_x));
return (sw & (FP_NAN | FP_NORMAL | FP_INFINITE | FP_ZERO | FP_SUBNORMAL))
== FP_NAN;
#endif
}
__CRT_INLINE int __cdecl __isnanl (long double _x)
{
#if defined(__x86_64__) || defined(_AMD64_)
__mingw_fp_types_t ld;
int xx, signexp;
ld.ld = &_x;
signexp = (ld.ldt->lh.sign_exponent & 0x7fff) << 1;
xx = (int) (ld.ldt->lh.low | (ld.ldt->lh.high & 0x7fffffffu)); /* explicit */
signexp |= (unsigned int) (xx | (-xx)) >> 31;
signexp = 0xfffe - signexp;
return (int) ((unsigned int) signexp) >> 16;
#elif defined(__arm__) || defined(_ARM_)
__isnan(_x);
#elif defined(__i386__) || defined(_X86_)
unsigned short sw;
__asm__ __volatile__ ("fxam;"
"fstsw %%ax": "=a" (sw) : "t" (_x));
return (sw & (FP_NAN | FP_NORMAL | FP_INFINITE | FP_ZERO | FP_SUBNORMAL))
== FP_NAN;
#endif
}
#endif
#define isnan(x) \
__mingw_choose_expr ( \
__mingw_types_compatible_p (__typeof__ (x), double), \
__isnan(x), \
__mingw_choose_expr ( \
__mingw_types_compatible_p (__typeof__ (x), float), \
__isnanf(x), \
__mingw_choose_expr ( \
__mingw_types_compatible_p (__typeof__ (x), long double), \
__isnanl(x), \
__dfp_expansion(__isnan,(__builtin_trap(),x),x))))
/* 7.12.3.5 */
#define isnormal(x) (fpclassify(x) == FP_NORMAL)
/* 7.12.3.6 The signbit macro */
extern int __cdecl __signbit (double);
extern int __cdecl __signbitf (float);
extern int __cdecl __signbitl (long double);
#ifndef __CRT__NO_INLINE
__CRT_INLINE int __cdecl __signbit (double x) {
#if defined(__x86_64__) || defined(_AMD64_) || defined(__arm__) || defined(_ARM_)
__mingw_fp_types_t hlp;
hlp.d = &x;
return ((hlp.dt->lh.high & 0x80000000) != 0);
#elif defined(__i386__) || defined(_X86_)
unsigned short stw;
__asm__ __volatile__ ( "fxam; fstsw %%ax;": "=a" (stw) : "t" (x));
return stw & 0x0200;
#endif
}
__CRT_INLINE int __cdecl __signbitf (float x) {
#if defined(__x86_64__) || defined(_AMD64_) || defined(__arm__) || defined(_ARM_)
__mingw_fp_types_t hlp;
hlp.f = &x;
return ((hlp.ft->val & 0x80000000) != 0);
#elif defined(__i386__) || defined(_X86_)
unsigned short stw;
__asm__ __volatile__ ("fxam; fstsw %%ax;": "=a" (stw) : "t" (x));
return stw & 0x0200;
#endif
}
__CRT_INLINE int __cdecl __signbitl (long double x) {
#if defined(__x86_64__) || defined(_AMD64_)
__mingw_fp_types_t ld;
ld.ld = &x;
return ((ld.ldt->lh.sign_exponent & 0x8000) != 0);
#elif defined(__arm__) || defined(_ARM_)
__signbit(x);
#elif defined(__i386__) || defined(_X86_)
unsigned short stw;
__asm__ __volatile__ ("fxam; fstsw %%ax;": "=a" (stw) : "t" (x));
return stw & 0x0200;
#endif
}
#endif
#define signbit(x) \
__mingw_choose_expr ( \
__mingw_types_compatible_p (__typeof__ (x), double), \
__signbit(x), \
__mingw_choose_expr ( \
__mingw_types_compatible_p (__typeof__ (x), float), \
__signbitf(x), \
__mingw_choose_expr ( \
__mingw_types_compatible_p (__typeof__ (x), long double), \
__signbitl(x), \
__dfp_expansion(__signbit,(__builtin_trap(),x),x))))
/* 7.12.4 Trigonometric functions: Double in C89 */
extern float __cdecl sinf(float _X);
extern long double __cdecl sinl(long double);
extern float __cdecl cosf(float _X);
extern long double __cdecl cosl(long double);
extern float __cdecl tanf(float _X);
extern long double __cdecl tanl(long double);
extern float __cdecl asinf(float _X);
extern long double __cdecl asinl(long double);
extern float __cdecl acosf (float);
extern long double __cdecl acosl (long double);
extern float __cdecl atanf (float);
extern long double __cdecl atanl (long double);
extern float __cdecl atan2f (float, float);
extern long double __cdecl atan2l (long double, long double);
/* 7.12.5 Hyperbolic functions: Double in C89 */
extern float __cdecl sinhf(float _X);
#ifndef __CRT__NO_INLINE
__CRT_INLINE float sinhf(float _X) { return ((float)sinh((double)_X)); }
#endif
extern long double __cdecl sinhl(long double);
extern float __cdecl coshf(float _X);
#ifndef __CRT__NO_INLINE
__CRT_INLINE float coshf(float _X) { return ((float)cosh((double)_X)); }
#endif
extern long double __cdecl coshl(long double);
extern float __cdecl tanhf(float _X);
#ifndef __CRT__NO_INLINE
__CRT_INLINE float tanhf(float _X) { return ((float)tanh((double)_X)); }
#endif
extern long double __cdecl tanhl(long double);
/* Inverse hyperbolic trig functions */
/* 7.12.5.1 */
extern double __cdecl acosh (double);
extern float __cdecl acoshf (float);
extern long double __cdecl acoshl (long double);
/* 7.12.5.2 */
extern double __cdecl asinh (double);
extern float __cdecl asinhf (float);
extern long double __cdecl asinhl (long double);
/* 7.12.5.3 */
extern double __cdecl atanh (double);
extern float __cdecl atanhf (float);
extern long double __cdecl atanhl (long double);
/* Exponentials and logarithms */
/* 7.12.6.1 Double in C89 */
extern float __cdecl expf(float _X);
#ifndef __CRT__NO_INLINE
__CRT_INLINE float expf(float _X) { return ((float)exp((double)_X)); }
#endif
extern long double __cdecl expl(long double);
/* 7.12.6.2 */
extern double __cdecl exp2(double);
extern float __cdecl exp2f(float);
extern long double __cdecl exp2l(long double);
/* 7.12.6.3 The expm1 functions */
/* TODO: These could be inlined */
extern double __cdecl expm1(double);
extern float __cdecl expm1f(float);
extern long double __cdecl expm1l(long double);
/* 7.12.6.4 Double in C89 */
extern float frexpf(float _X,int *_Y);
#ifndef __CRT__NO_INLINE
__CRT_INLINE float frexpf(float _X,int *_Y) { return ((float)frexp((double)_X,_Y)); }
#endif
extern long double __cdecl frexpl(long double,int *);
/* 7.12.6.5 */
#define FP_ILOGB0 ((int)0x80000000)
#define FP_ILOGBNAN ((int)0x80000000)
extern int __cdecl ilogb (double);
extern int __cdecl ilogbf (float);
extern int __cdecl ilogbl (long double);
/* 7.12.6.6 Double in C89 */
extern float __cdecl ldexpf(float _X,int _Y);
#ifndef __CRT__NO_INLINE
__CRT_INLINE float __cdecl ldexpf (float x, int expn) { return (float) ldexp ((double)x, expn); }
#endif
extern long double __cdecl ldexpl (long double, int);
/* 7.12.6.7 Double in C89 */
extern float __cdecl logf (float);
extern long double __cdecl logl(long double);
/* 7.12.6.8 Double in C89 */
extern float __cdecl log10f (float);
extern long double __cdecl log10l(long double);
/* 7.12.6.9 */
extern double __cdecl log1p(double);
extern float __cdecl log1pf(float);
extern long double __cdecl log1pl(long double);
/* 7.12.6.10 */
extern double __cdecl log2 (double);
extern float __cdecl log2f (float);
extern long double __cdecl log2l (long double);
/* 7.12.6.11 */
extern double __cdecl logb (double);
extern float __cdecl logbf (float);
extern long double __cdecl logbl (long double);
#ifndef __CRT__NO_INLINE
/* When compiling with gcc, always use gcc's builtins.
* The asm inlines below are kept here for future reference:
* they were written for gcc and do no error handling
* (exceptions/errno), therefore only valid if __FAST_MATH__
* is defined (-ffast-math) . */
#if 0 /*defined(__GNUC__) && defined(__FAST_MATH__)*/
__CRT_INLINE double __cdecl logb (double x)
{
#if defined(__x86_64__) || defined(_AMD64_) || defined(__arm__) || defined(_ARM_)
__mingw_fp_types_t hlp;
int lx, hx;
hlp.d = &x;
lx = hlp.dt->lh.low;
hx = hlp.dt->lh.high & 0x7fffffff; /* high |x| */
if ((hx | lx) == 0)
return -1.0 / fabs (x);
if (hx >= 0x7ff00000)
return x * x;
if ((hx >>= 20) == 0) /* IEEE 754 logb */
return -1022.0;
return (double) (hx - 1023);
#elif defined(__i386__) || defined(_X86_)
double res = 0.0;
__asm__ __volatile__ ("fxtract\n\t"
"fstp %%st" : "=t" (res) : "0" (x));
return res;
#endif
}
__CRT_INLINE float __cdecl logbf (float x)
{
#if defined(__x86_64__) || defined(_AMD64_) || defined(__arm__) || defined(_ARM_)
int v;
__mingw_fp_types_t hlp;
hlp.f = &x;
v = hlp.ft->val & 0x7fffffff; /* high |x| */
if (!v)
return (float)-1.0 / fabsf (x);
if (v >= 0x7f800000)
return x * x;
if ((v >>= 23) == 0) /* IEEE 754 logb */
return -126.0;
return (float) (v - 127);
#elif defined(__i386__) || defined(_X86_)
float res = 0.0F;
__asm__ __volatile__ ("fxtract\n\t"
"fstp %%st" : "=t" (res) : "0" (x));
return res;
#endif
}
__CRT_INLINE long double __cdecl logbl (long double x)
{
#if defined(__arm__) || defined(_ARM_)
__mingw_fp_types_t hlp;
int lx, hx;
hlp.d = &x;
lx = hlp.dt->lh.low;
hx = hlp.dt->lh.high & 0x7fffffff; /* high |x| */
if ((hx | lx) == 0)
return -1.0 / fabs (x);
if (hx >= 0x7ff00000)
return x * x;
if ((hx >>= 20) == 0) /* IEEE 754 logb */
return -1022.0;
return (double) (hx - 1023);
#elif defined(__x86_64__) || defined(_AMD64_) || defined(__i386__) || defined(_X86_)
long double res = 0.0l;
__asm__ __volatile__ ("fxtract\n\t"
"fstp %%st" : "=t" (res) : "0" (x));
return res;
#endif
}
#endif /* defined(__GNUC__) && defined(__FAST_MATH__) */
#endif /* __CRT__NO_INLINE */
/* 7.12.6.12 Double in C89 */
extern float __cdecl modff (float, float*);
extern long double __cdecl modfl (long double, long double*);
/* 7.12.6.13 */
extern double __cdecl scalbn (double, int);
extern float __cdecl scalbnf (float, int);
extern long double __cdecl scalbnl (long double, int);
extern double __cdecl scalbln (double, long);
extern float __cdecl scalblnf (float, long);
extern long double __cdecl scalblnl (long double, long);
/* 7.12.7.1 */
/* Implementations adapted from Cephes versions */
extern double __cdecl cbrt (double);
extern float __cdecl cbrtf (float);
extern long double __cdecl cbrtl (long double);
/* 7.12.7.3 */
extern double __cdecl hypot (double, double) __MINGW_ATTRIB_DEPRECATED_MSVC2005; /* in libmoldname.a */
extern float __cdecl hypotf (float x, float y);
#ifndef __CRT__NO_INLINE
__CRT_INLINE float __cdecl hypotf (float x, float y) { return (float) hypot ((double)x, (double)y);}
#endif
extern long double __cdecl hypotl (long double, long double);
/* 7.12.7.4 The pow functions. Double in C89 */
extern float __cdecl powf(float _X,float _Y);
#ifndef __CRT__NO_INLINE
__CRT_INLINE float powf(float _X,float _Y) { return ((float)pow((double)_X,(double)_Y)); }
#endif
extern long double __cdecl powl (long double, long double);
/* 7.12.7.5 The sqrt functions. Double in C89. */
extern float __cdecl sqrtf (float);
extern long double sqrtl(long double);
/* 7.12.8.1 The erf functions */
extern double __cdecl erf (double);
extern float __cdecl erff (float);
extern long double __cdecl erfl (long double);
/* 7.12.8.2 The erfc functions */
extern double __cdecl erfc (double);
extern float __cdecl erfcf (float);
extern long double __cdecl erfcl (long double);
/* 7.12.8.3 The lgamma functions */
extern double __cdecl lgamma (double);
extern float __cdecl lgammaf (float);
extern long double __cdecl lgammal (long double);
extern int signgam;
/* 7.12.8.4 The tgamma functions */
extern double __cdecl tgamma (double);
extern float __cdecl tgammaf (float);
extern long double __cdecl tgammal (long double);
/* 7.12.9.1 Double in C89 */
extern float __cdecl ceilf (float);
extern long double __cdecl ceill (long double);
/* 7.12.9.2 Double in C89 */
extern float __cdecl floorf (float);
extern long double __cdecl floorl (long double);
/* 7.12.9.3 */
extern double __cdecl nearbyint ( double);
extern float __cdecl nearbyintf (float);
extern long double __cdecl nearbyintl (long double);
/* 7.12.9.4 */
/* round, using fpu control word settings */
extern double __cdecl rint (double);
extern float __cdecl rintf (float);
extern long double __cdecl rintl (long double);
/* 7.12.9.5 */
extern long __cdecl lrint (double);
extern long __cdecl lrintf (float);
extern long __cdecl lrintl (long double);
__MINGW_EXTENSION long long __cdecl llrint (double);
__MINGW_EXTENSION long long __cdecl llrintf (float);
__MINGW_EXTENSION long long __cdecl llrintl (long double);
#ifndef __CRT__NO_INLINE
/* When compiling with gcc, always use gcc's builtins.
* The asm inlines below are kept here for future reference:
* they were written for gcc and do no error handling
* (exceptions/errno), therefore only valid if __FAST_MATH__
* is defined (-ffast-math) . */
#if 0 /*defined(__GNUC__) && defined(__FAST_MATH__)*/
__CRT_INLINE double __cdecl rint (double x)
{
double retval = 0.0;
__asm__ __volatile__ ("frndint;": "=t" (retval) : "0" (x));
return retval;
}
__CRT_INLINE float __cdecl rintf (float x)
{
float retval = 0.0;
__asm__ __volatile__ ("frndint;" : "=t" (retval) : "0" (x) );
return retval;
}
__CRT_INLINE long double __cdecl rintl (long double x)
{
long double retval = 0.0l;
__asm__ __volatile__ ("frndint;" : "=t" (retval) : "0" (x) );
return retval;
}
__CRT_INLINE long __cdecl lrint (double x)
{
long retval = 0;
__asm__ __volatile__ \
("fistpl %0" : "=m" (retval) : "t" (x) : "st"); \
return retval;
}
__CRT_INLINE long __cdecl lrintf (float x)
{
long retval = 0;
__asm__ __volatile__ \
("fistpl %0" : "=m" (retval) : "t" (x) : "st"); \
return retval;
}
__CRT_INLINE long __cdecl lrintl (long double x)
{
long retval = 0;
__asm__ __volatile__ \
("fistpl %0" : "=m" (retval) : "t" (x) : "st"); \
return retval;
}