projects
/
musl
/ blobdiff
commit
grep
author
committer
pickaxe
?
search:
re
summary
|
shortlog
|
log
|
commit
|
commitdiff
|
tree
raw
|
inline
| side by side
math: sin cos cleanup
[musl]
/
src
/
math
/
scalbn.c
diff --git
a/src/math/scalbn.c
b/src/math/scalbn.c
index
42a2b29
..
003141e
100644
(file)
--- a/
src/math/scalbn.c
+++ b/
src/math/scalbn.c
@@
-2,8
+2,6
@@
double scalbn(double x, int n)
{
double scalbn(double x, int n)
{
- /* make sure result is stored as double on overflow or underflow */
- volatile double z;
double scale;
if (n > 1023) {
double scale;
if (n > 1023) {
@@
-13,8
+11,8
@@
double scalbn(double x, int n)
x *= 0x1p1023;
n -= 1023;
if (n > 1023) {
x *= 0x1p1023;
n -= 1023;
if (n > 1023) {
-
z = x * 0x1p1023
;
- return
z
;
+
STRICT_ASSIGN(double, x, x * 0x1p1023)
;
+ return
x
;
}
}
} else if (n < -1022) {
}
}
} else if (n < -1022) {
@@
-24,12
+22,12
@@
double scalbn(double x, int n)
x *= 0x1p-1022;
n += 1022;
if (n < -1022) {
x *= 0x1p-1022;
n += 1022;
if (n < -1022) {
-
z = x * 0x1p-1022
;
- return
z
;
+
STRICT_ASSIGN(double, x, x * 0x1p-1022)
;
+ return
x
;
}
}
}
INSERT_WORDS(scale, (uint32_t)(0x3ff+n)<<20, 0);
}
}
}
INSERT_WORDS(scale, (uint32_t)(0x3ff+n)<<20, 0);
-
z = x * scale
;
- return
z
;
+
STRICT_ASSIGN(double, x, x * scale)
;
+ return
x
;
}
}