aboutsummaryrefslogtreecommitdiff
diff options
context:
space:
mode:
authorSteve Kargl <kargl@FreeBSD.org>2022-08-04 17:31:57 +0000
committerDimitry Andric <dim@FreeBSD.org>2022-08-07 09:03:20 +0000
commit162bdeaa7324c4c75b015530eac4bcbd6c58d639 (patch)
tree5f7de443375ec75623d549cfae3892015b9ce713
parentb38d534d58c7a7145e5779e811e85418cb9f2a89 (diff)
downloadsrc-162bdeaa7324c4c75b015530eac4bcbd6c58d639.tar.gz
src-162bdeaa7324c4c75b015530eac4bcbd6c58d639.zip
[libm] Correct comments in s_cbrt[l].c
Damian McGuckin <damianm at esi dot com dot au> noted that the accuracy claims in the code for cbrt(3) and cbrtl(3) were incorrect. Fix the comments to more accurately describe the accuracies. PR: 265603 MFC after: 3 days (cherry picked from commit 369ea0520a3061c07400d7cd32172efb6af39815)
-rw-r--r--lib/msun/src/s_cbrt.c2
-rw-r--r--lib/msun/src/s_cbrtl.c2
2 files changed, 2 insertions, 2 deletions
diff --git a/lib/msun/src/s_cbrt.c b/lib/msun/src/s_cbrt.c
index 0e609e110d9a..4353d340a1ac 100644
--- a/lib/msun/src/s_cbrt.c
+++ b/lib/msun/src/s_cbrt.c
@@ -108,7 +108,7 @@ cbrt(double x)
r=x/s; /* error <= 0.5 ulps; |r| < |t| */
w=t+t; /* t+t is exact */
r=(r-t)/(w+r); /* r-t is exact; w+r ~= 3*t */
- t=t+t*r; /* error <= 0.5 + 0.5/3 + epsilon */
+ t=t+t*r; /* error <= (0.5 + 0.5/3) * ulp */
return(t);
}
diff --git a/lib/msun/src/s_cbrtl.c b/lib/msun/src/s_cbrtl.c
index 2236c0fa0348..b15c96ecde03 100644
--- a/lib/msun/src/s_cbrtl.c
+++ b/lib/msun/src/s_cbrtl.c
@@ -136,7 +136,7 @@ cbrtl(long double x)
r=x/s; /* error <= 0.5 ulps; |r| < |t| */
w=t+t; /* t+t is exact */
r=(r-t)/(w+r); /* r-t is exact; w+r ~= 3*t */
- t=t+t*r; /* error <= 0.5 + 0.5/3 + epsilon */
+ t=t+t*r; /* error <= (0.5 + 0.5/3) * ulp */
t *= v.e;
RETURNI(t);