aboutsummaryrefslogtreecommitdiff
path: root/newlib/libm/math
diff options
context:
space:
mode:
authorJeff Johnston <jjohnstn@redhat.com>2021-06-04 14:36:38 -0400
committerJeff Johnston <jjohnstn@redhat.com>2021-06-04 14:42:58 -0400
commita9165ea07ce3a99bbfd294ada9c0fd54db738bb9 (patch)
tree8241f6063f952db54cb8a24dc8e2e83eac4e9b70 /newlib/libm/math
parentdfe5988f961ff97d283a9c460e75499db168163a (diff)
downloadnewlib-a9165ea07ce3a99bbfd294ada9c0fd54db738bb9.zip
newlib-a9165ea07ce3a99bbfd294ada9c0fd54db738bb9.tar.gz
newlib-a9165ea07ce3a99bbfd294ada9c0fd54db738bb9.tar.bz2
Fix rounding issues with sqrt/sqrtf
- compiler is sometimes optimizing out the rounding check in e_sqrt.c and ef_sqrt.c which uses two constants to create an inexact operation - there is a similar constant operation in s_tanh.c/sf_tanh.c - make the one and tiny constants volatile to stop this
Diffstat (limited to 'newlib/libm/math')
-rw-r--r--newlib/libm/math/e_sqrt.c4
-rw-r--r--newlib/libm/math/ef_sqrt.c4
-rw-r--r--newlib/libm/math/s_tanh.c4
-rw-r--r--newlib/libm/math/sf_tanh.c4
4 files changed, 8 insertions, 8 deletions
diff --git a/newlib/libm/math/e_sqrt.c b/newlib/libm/math/e_sqrt.c
index 78fc524..85dbdd1 100644
--- a/newlib/libm/math/e_sqrt.c
+++ b/newlib/libm/math/e_sqrt.c
@@ -86,9 +86,9 @@
#ifndef _DOUBLE_IS_32BITS
#ifdef __STDC__
-static const double one = 1.0, tiny=1.0e-300;
+static const volatile double one = 1.0, tiny=1.0e-300;
#else
-static double one = 1.0, tiny=1.0e-300;
+static double volatile one = 1.0, tiny=1.0e-300;
#endif
#ifdef __STDC__
diff --git a/newlib/libm/math/ef_sqrt.c b/newlib/libm/math/ef_sqrt.c
index 80e7f36..448e26e 100644
--- a/newlib/libm/math/ef_sqrt.c
+++ b/newlib/libm/math/ef_sqrt.c
@@ -16,9 +16,9 @@
#include "fdlibm.h"
#ifdef __STDC__
-static const float one = 1.0, tiny=1.0e-30;
+static const volatile float one = 1.0, tiny=1.0e-30;
#else
-static float one = 1.0, tiny=1.0e-30;
+static float volatile one = 1.0, tiny=1.0e-30;
#endif
#ifdef __STDC__
diff --git a/newlib/libm/math/s_tanh.c b/newlib/libm/math/s_tanh.c
index f786430..66018d7 100644
--- a/newlib/libm/math/s_tanh.c
+++ b/newlib/libm/math/s_tanh.c
@@ -73,9 +73,9 @@ PORTABILITY
#ifndef _DOUBLE_IS_32BITS
#ifdef __STDC__
-static const double one=1.0, two=2.0, tiny = 1.0e-300;
+static const volatile double one=1.0, two=2.0, tiny = 1.0e-300;
#else
-static double one=1.0, two=2.0, tiny = 1.0e-300;
+static double volatile one=1.0, two=2.0, tiny = 1.0e-300;
#endif
#ifdef __STDC__
diff --git a/newlib/libm/math/sf_tanh.c b/newlib/libm/math/sf_tanh.c
index 1eb44a2..214bdb4 100644
--- a/newlib/libm/math/sf_tanh.c
+++ b/newlib/libm/math/sf_tanh.c
@@ -16,9 +16,9 @@
#include "fdlibm.h"
#ifdef __STDC__
-static const float one=1.0, two=2.0, tiny = 1.0e-30;
+static const volatile float one=1.0, two=2.0, tiny = 1.0e-30;
#else
-static float one=1.0, two=2.0, tiny = 1.0e-30;
+static volatile float one=1.0, two=2.0, tiny = 1.0e-30;
#endif
#ifdef __STDC__