From f6565cd1a68b4a4b7c43b588513cd802597bc5d3 Mon Sep 17 00:00:00 2001 From: Danny Smith Date: Sun, 1 Feb 2004 02:52:21 +0000 Subject: [PATCH] * mingwex/math/ldexpl.c (ldexpl): Call __asm__("fscale") directly, rather than via scabnl. --- winsup/mingw/ChangeLog | 5 +++++ winsup/mingw/mingwex/math/ldexpl.c | 21 +++++++++++++-------- 2 files changed, 18 insertions(+), 8 deletions(-) diff --git a/winsup/mingw/ChangeLog b/winsup/mingw/ChangeLog index f32087fac..bd57ce6b4 100644 --- a/winsup/mingw/ChangeLog +++ b/winsup/mingw/ChangeLog @@ -1,3 +1,8 @@ +2004-02-01 Danny Smith + + * mingwex/math/ldexpl.c (ldexpl): Call __asm__("fscale") + directly, rather than via scabnl. + 2004-02-01 Danny Smith * mingwex/math/powl.c (powl): Return infinity if 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 #include - - + 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; } -- 2.43.5