projects
/
musl
/ blobdiff
commit
grep
author
committer
pickaxe
?
search:
re
summary
|
shortlog
|
log
|
commit
|
commitdiff
|
tree
raw
|
inline
| side by side
drop use of pthread_once for aio thread stack size init
[musl]
/
src
/
math
/
frexpl.c
diff --git
a/src/math/frexpl.c
b/src/math/frexpl.c
index
f9d90a6
..
3c1b553
100644
(file)
--- a/
src/math/frexpl.c
+++ b/
src/math/frexpl.c
@@
-1,20
+1,20
@@
-#include <math.h>
-#include <stdint.h>
-#include <float.h>
-
-#if LDBL_MANT_DIG == 64 && LDBL_MAX_EXP == 16384
-
-/* This version is for 80-bit little endian long double */
+#include "libm.h"
+#if LDBL_MANT_DIG == 53 && LDBL_MAX_EXP == 1024
long double frexpl(long double x, int *e)
{
long double frexpl(long double x, int *e)
{
- union { long double ld; uint16_t hw[5]; } y = { x };
- int ee = y.hw[4]&0x7fff;
+ return frexp(x, e);
+}
+#elif (LDBL_MANT_DIG == 64 || LDBL_MANT_DIG == 113) && LDBL_MAX_EXP == 16384
+long double frexpl(long double x, int *e)
+{
+ union ldshape u = {x};
+ int ee = u.i.se & 0x7fff;
if (!ee) {
if (x) {
if (!ee) {
if (x) {
- x = frexpl(x*0x1p
64
, e);
- *e -=
64
;
+ x = frexpl(x*0x1p
120
, e);
+ *e -=
120
;
} else *e = 0;
return x;
} else if (ee == 0x7fff) {
} else *e = 0;
return x;
} else if (ee == 0x7fff) {
@@
-22,16
+22,8
@@
long double frexpl(long double x, int *e)
}
*e = ee - 0x3ffe;
}
*e = ee - 0x3ffe;
-
y.hw[4]
&= 0x8000;
-
y.hw[4]
|= 0x3ffe;
- return
y.ld
;
+
u.i.se
&= 0x8000;
+
u.i.se
|= 0x3ffe;
+ return
u.f
;
}
}
-
-#else
-
-long double frexpl(long double x, int *e)
-{
- return frexp(x, e);
-}
-
#endif
#endif