X-Git-Url: http://nsz.repo.hu/git/?a=blobdiff_plain;f=src%2Fmath%2Fpowl.c;h=690f2942a034ed02dad81dd92567bdde1d9ed446;hb=405a21d4cbce255eefd2b59626d8459d7d1108d4;hp=c829c18fcaeacbf0f963273f1932ada48f36bace;hpb=59c2bd747bce9b1e35c893b2ace9bea57523a901;p=libm diff --git a/src/math/powl.c b/src/math/powl.c index c829c18..690f294 100644 --- a/src/math/powl.c +++ b/src/math/powl.c @@ -69,7 +69,12 @@ #include "libm.h" -#if LDBL_MANT_DIG == 64 && LDBL_MAX_EXP == 16384 +#if LDBL_MANT_DIG == 53 && LDBL_MAX_EXP == 1024 +long double powl(long double x, long double y) +{ + return pow(x, y); +} +#elif LDBL_MANT_DIG == 64 && LDBL_MAX_EXP == 16384 /* Table size */ #define NXT 32 @@ -207,7 +212,7 @@ long double powl(long double x, long double y) if (y == 1.0L) return x; - // FIXME: this is wrong, see pow special cases in posix2008 + // FIXME: this is wrong, see pow special cases in c99 F.9.4.4 if (!isfinite(y) && (x == -1.0L || x == 1.0L) ) return y - y; /* +-1**inf is NaN */ if (x == 1.0L)