math: minor scalbn*.c simplification

diff --git a/src/math/scalbn.c b/src/math/scalbn.c
index 003141e3c14a9a8f59ec68407ca719c9f7d6de8b..1fec432f70caa913ef9e3817205c0facd2891c73 100644 (file)
--- a/src/math/scalbn.c
+++ b/src/math/scalbn.c
@@ -10,10 +10,8 @@ double scalbn(double x, int n)
                if (n > 1023) {
                        x *= 0x1p1023;
                        n -= 1023;
-                       if (n > 1023) {
-                               STRICT_ASSIGN(double, x, x * 0x1p1023);
-                               return x;
-                       }
+                       if (n > 1023)
+                               n = 1023;
                }
        } else if (n < -1022) {
                x *= 0x1p-1022;
@@ -21,10 +19,8 @@ double scalbn(double x, int n)
                if (n < -1022) {
                        x *= 0x1p-1022;
                        n += 1022;
-                       if (n < -1022) {
-                               STRICT_ASSIGN(double, x, x * 0x1p-1022);
-                               return x;
-                       }
+                       if (n < -1022)
+                               n = -1022;
                }
        }
        INSERT_WORDS(scale, (uint32_t)(0x3ff+n)<<20, 0);

diff --git a/src/math/scalbnf.c b/src/math/scalbnf.c
index f94b5d5992bb0ad8ed67e19323e68853201df045..c0eeaf840ae8c0459f4796d4914a674fea361bdd 100644 (file)
--- a/src/math/scalbnf.c
+++ b/src/math/scalbnf.c
@@ -10,10 +10,8 @@ float scalbnf(float x, int n)
                if (n > 127) {
                        x *= 0x1p127f;
                        n -= 127;
-                       if (n > 127) {
-                               STRICT_ASSIGN(float, x, x * 0x1p127f);
-                               return x;
-                       }
+                       if (n > 127)
+                               n = 127;
                }
        } else if (n < -126) {
                x *= 0x1p-126f;
@@ -21,10 +19,8 @@ float scalbnf(float x, int n)
                if (n < -126) {
                        x *= 0x1p-126f;
                        n += 126;
-                       if (n < -126) {
-                               STRICT_ASSIGN(float, x, x * 0x1p-126f);
-                               return x;
-                       }
+                       if (n < -126)
+                               n = -126;
                }
        }
        SET_FLOAT_WORD(scale, (uint32_t)(0x7f+n)<<23);

diff --git a/src/math/scalbnl.c b/src/math/scalbnl.c
index c605b8da3400a2e2a40d7e4f537974d2bb969f7f..7ad7688b6ccd2841cf54da09d3c59a3a2252006c 100644 (file)
--- a/src/math/scalbnl.c
+++ b/src/math/scalbnl.c
@@ -17,7 +17,7 @@ long double scalbnl(long double x, int n)
                        x *= 0x1p16383L;
                        n -= 16383;
                        if (n > 16383)
-                               return x * 0x1p16383L;
+                               n = 16383;
                }
        } else if (n < -16382) {
                x *= 0x1p-16382L;
@@ -26,7 +26,7 @@ long double scalbnl(long double x, int n)
                        x *= 0x1p-16382L;
                        n += 16382;
                        if (n < -16382)
-                               return x * 0x1p-16382L;
+                               n = -16382;
                }
        }
        scale.e = 1.0;