projects
/
musl
/ blobdiff
commit
grep
author
committer
pickaxe
?
search:
re
summary
|
shortlog
|
log
|
commit
|
commitdiff
|
tree
raw
|
inline
| side by side
floating point environment/exceptions support for mips
[musl]
/
src
/
math
/
expl.c
diff --git
a/src/math/expl.c
b/src/math/expl.c
index
9507fd2
..
b289ffe
100644
(file)
--- a/
src/math/expl.c
+++ b/
src/math/expl.c
@@
-102,13
+102,13
@@
long double expl(long double x)
if (x > MAXLOGL)
return INFINITY;
if (x < MINLOGL)
if (x > MAXLOGL)
return INFINITY;
if (x < MINLOGL)
- return 0.0
L
;
+ return 0.0;
/* Express e**x = e**g 2**n
* = e**g e**(n loge(2))
* = e**(g + n loge(2))
*/
/* Express e**x = e**g 2**n
* = e**g e**(n loge(2))
* = e**(g + n loge(2))
*/
- px = floorl(LOG2EL * x + 0.5
L
); /* floor() truncates toward -infinity. */
+ px = floorl(LOG2EL * x + 0.5); /* floor() truncates toward -infinity. */
n = px;
x -= px * C1;
x -= px * C2;
n = px;
x -= px * C1;
x -= px * C2;
@@
-120,8
+120,8
@@
long double expl(long double x)
xx = x * x;
px = x * __polevll(xx, P, 2);
x = px/(__polevll(xx, Q, 3) - px);
xx = x * x;
px = x * __polevll(xx, P, 2);
x = px/(__polevll(xx, Q, 3) - px);
- x = 1.0
L + ldexpl(x, 1)
;
- x =
ldexp
l(x, n);
+ x = 1.0
+ 2.0 * x
;
+ x =
scalbn
l(x, n);
return x;
}
#endif
return x;
}
#endif