* mingwex/math/ldexpl.c (ldexpl): Call __asm__("fscale")

directly, rather than via scabnl.
This commit is contained in:
Danny Smith 2004-02-01 02:52:21 +00:00
parent 416bc45060
commit f6565cd1a6
2 changed files with 18 additions and 8 deletions

View File

@ -1,3 +1,8 @@
2004-02-01 Danny Smith <dannysmith@users.sourceforge.net>
* mingwex/math/ldexpl.c (ldexpl): Call __asm__("fscale")
directly, rather than via scabnl.
2004-02-01 Danny Smith <dannysmith@users.sourceforge.net>
* mingwex/math/powl.c (powl): Return infinity if

View File

@ -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;
}