diff options
Diffstat (limited to 'winsup/mingw/mingwex/math/ldexpl.c')
-rw-r--r-- | winsup/mingw/mingwex/math/ldexpl.c | 21 |
1 files changed, 13 insertions, 8 deletions
diff --git a/winsup/mingw/mingwex/math/ldexpl.c b/winsup/mingw/mingwex/math/ldexpl.c index e4477bf6c..19a3d56e3 100644 --- a/winsup/mingw/mingwex/math/ldexpl.c +++ b/winsup/mingw/mingwex/math/ldexpl.c @@ -1,14 +1,19 @@ #include <math.h> #include <errno.h> - - + long double ldexpl(long double x, int expn) { - if (isfinite (x) && x != 0.0L) - { - x = scalbnl (x , expn); - if (!isfinite (x) || x == 0.0) errno = ERANGE; - } - return x; + long double res; + if (!isfinite (x) || x == 0.0L) + return x; + + __asm__ ("fscale" + : "=t" (res) + : "0" (x), "u" ((long double) expn)); + + if (!isfinite (res) || res == 0.0L) + errno = ERANGE; + + return res; } |