From 9b0fcb441a44456c7b071c7cdaf90403f81ec05a Mon Sep 17 00:00:00 2001 From: Szabolcs Nagy Date: Fri, 6 Sep 2013 18:35:55 +0000 Subject: math: remove STRICT_ASSIGN macro gcc did not always drop excess precision according to c99 at assignments before version 4.5 even if -std=c99 was requested which caused badly broken mathematical functions on i386 when FLT_EVAL_METHOD!=0 but STRICT_ASSIGN was not used consistently and it is worked around for old compilers with -ffloat-store so it is no longer needed the new convention is to get the compiler respect c99 semantics and when excess precision is not harmful use float_t or double_t or to specialize code using FLT_EVAL_METHOD --- src/math/__rem_pio2.c | 2 +- src/math/__rem_pio2_large.c | 3 ++- src/math/__rem_pio2f.c | 2 +- src/math/exp.c | 4 ++-- src/math/exp2.c | 2 +- src/math/expf.c | 4 ++-- src/math/expm1.c | 2 +- src/math/expm1f.c | 2 +- src/math/log1p.c | 2 +- src/math/log1pf.c | 2 +- 10 files changed, 13 insertions(+), 12 deletions(-) (limited to 'src/math') diff --git a/src/math/__rem_pio2.c b/src/math/__rem_pio2.c index 0ef57fb5..9305be5f 100644 --- a/src/math/__rem_pio2.c +++ b/src/math/__rem_pio2.c @@ -112,7 +112,7 @@ int __rem_pio2(double x, double *y) uint32_t high; medium: /* Use a specialized rint() to get fn. Assume round-to-nearest. */ - STRICT_ASSIGN(double, fn, x*invpio2 + 0x1.8p52); + fn = x*invpio2 + 0x1.8p52; fn = fn - 0x1.8p52; // FIXME #ifdef HAVE_EFFICIENT_IRINT diff --git a/src/math/__rem_pio2_large.c b/src/math/__rem_pio2_large.c index 27b619cc..bb2dc43f 100644 --- a/src/math/__rem_pio2_large.c +++ b/src/math/__rem_pio2_large.c @@ -415,7 +415,8 @@ recompute: fw = 0.0; for (i=jz; i>=0; i--) fw += fq[i]; - STRICT_ASSIGN(double,fw,fw); + // TODO: drop excess precision here once double_t is used + fw = (double)fw; y[0] = ih==0 ? fw : -fw; fw = fq[0]-fw; for (i=1; i<=jz; i++) diff --git a/src/math/__rem_pio2f.c b/src/math/__rem_pio2f.c index 198ee91b..5bdeb529 100644 --- a/src/math/__rem_pio2f.c +++ b/src/math/__rem_pio2f.c @@ -44,7 +44,7 @@ int __rem_pio2f(float x, double *y) /* 25+53 bit pi is good enough for medium size */ if (ix < 0x4dc90fdb) { /* |x| ~< 2^28*(pi/2), medium size */ /* Use a specialized rint() to get fn. Assume round-to-nearest. */ - STRICT_ASSIGN(double, fn, x*invpio2 + 0x1.8p52); + fn = x*invpio2 + 0x1.8p52; fn = fn - 0x1.8p52; // FIXME #ifdef HAVE_EFFICIENT_IRINT diff --git a/src/math/exp.c b/src/math/exp.c index 2b182acd..9ea672fa 100644 --- a/src/math/exp.c +++ b/src/math/exp.c @@ -94,7 +94,7 @@ double exp(double x) return x; if (x > 709.782712893383973096) { /* overflow if x!=inf */ - STRICT_ASSIGN(double, x, 0x1p1023 * x); + x *= 0x1p1023; return x; } if (x < -708.39641853226410622) { @@ -113,7 +113,7 @@ double exp(double x) k = 1 - sign - sign; hi = x - k*ln2hi; /* k*ln2hi is exact here */ lo = k*ln2lo; - STRICT_ASSIGN(double, x, hi - lo); + x = hi - lo; } else if (hx > 0x3e300000) { /* if |x| > 2**-28 */ k = 0; hi = x; diff --git a/src/math/exp2.c b/src/math/exp2.c index 2e078fb0..e14adba5 100644 --- a/src/math/exp2.c +++ b/src/math/exp2.c @@ -340,7 +340,7 @@ double exp2(double x) if (ix >= 0x408ff000) { /* |x| >= 1022 or nan */ if (ix >= 0x40900000 && u.i>>63 == 0) { /* x >= 1024 or nan */ /* overflow */ - STRICT_ASSIGN(double, x, x * 0x1p1023); + x *= 0x1p1023; return x; } if (ix >= 0x7ff00000) /* -inf or -nan */ diff --git a/src/math/expf.c b/src/math/expf.c index 5572bbf6..16e9afe6 100644 --- a/src/math/expf.c +++ b/src/math/expf.c @@ -41,7 +41,7 @@ float expf(float x) if (hx >= 0x42aeac50) { /* if |x| >= -87.33655f or NaN */ if (hx >= 0x42b17218 && !sign) { /* x >= 88.722839f */ /* overflow */ - STRICT_ASSIGN(float, x, x * 0x1p127f); + x *= 0x1p127f; return x; } if (sign) { @@ -60,7 +60,7 @@ float expf(float x) k = 1 - sign - sign; hi = x - k*ln2hi; /* k*ln2hi is exact here */ lo = k*ln2lo; - STRICT_ASSIGN(float, x, hi - lo); + x = hi - lo; } else if (hx > 0x39000000) { /* |x| > 2**-14 */ k = 0; hi = x; diff --git a/src/math/expm1.c b/src/math/expm1.c index a7eb2c0b..ac1e61e4 100644 --- a/src/math/expm1.c +++ b/src/math/expm1.c @@ -155,7 +155,7 @@ double expm1(double x) hi = x - t*ln2_hi; /* t*ln2_hi is exact here */ lo = t*ln2_lo; } - STRICT_ASSIGN(double, x, hi - lo); + x = hi-lo; c = (hi-x)-lo; } else if (hx < 0x3c900000) { /* |x| < 2**-54, return x */ if (hx < 0x00100000) diff --git a/src/math/expm1f.c b/src/math/expm1f.c index 698ab45f..297e0b44 100644 --- a/src/math/expm1f.c +++ b/src/math/expm1f.c @@ -65,7 +65,7 @@ float expm1f(float x) hi = x - t*ln2_hi; /* t*ln2_hi is exact here */ lo = t*ln2_lo; } - STRICT_ASSIGN(float, x, hi - lo); + x = hi-lo; c = (hi-x)-lo; } else if (hx < 0x33000000) { /* when |x|<2**-25, return x */ if (hx < 0x00800000) diff --git a/src/math/log1p.c b/src/math/log1p.c index 9bed63c2..a71ac423 100644 --- a/src/math/log1p.c +++ b/src/math/log1p.c @@ -122,7 +122,7 @@ double log1p(double x) return x+x; if (k != 0) { if (hx < 0x43400000) { - STRICT_ASSIGN(double, u, 1.0 + x); + u = 1 + x; GET_HIGH_WORD(hu, u); k = (hu>>20) - 1023; c = k > 0 ? 1.0-(u-x) : x-(u-1.0); /* correction term */ diff --git a/src/math/log1pf.c b/src/math/log1pf.c index c38e0bcb..e6940d29 100644 --- a/src/math/log1pf.c +++ b/src/math/log1pf.c @@ -61,7 +61,7 @@ float log1pf(float x) return x+x; if (k != 0) { if (hx < 0x5a000000) { - STRICT_ASSIGN(float, u, 1.0f + x); + u = 1 + x; GET_FLOAT_WORD(hu, u); k = (hu>>23) - 127; /* correction term */ -- cgit v1.2.1