From fa834172fd50fbd272cbc1c715ba6e48203a0f35 Mon Sep 17 00:00:00 2001 From: Taylor R Campbell Date: Thu, 25 Oct 2018 05:09:48 +0000 Subject: [PATCH] The bounds of Lemma 4 are <=, not <, as is this branch's condition. --- src/runtime/arith.scm | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/src/runtime/arith.scm b/src/runtime/arith.scm index b417b4843..a061722a5 100644 --- a/src/runtime/arith.scm +++ b/src/runtime/arith.scm @@ -2175,8 +2175,8 @@ USA. ;; ;; Since p = 2p/2 <= 1 <= 2*2p = 4p, the floating-point ;; evaluation of 1 - 2p is exact; the only error arises from - ;; division and log1p. First, note that if logistic(-1) < p < - ;; logistic(+1), (1 - 2p)/p lies in the bounds of Lemma 4. + ;; division and log1p. First, note that if logistic(-1) <= p + ;; <= logistic(+1), (1 - 2p)/p lies in the bounds of Lemma 4. ;; ;; If division has relative error d0 and log1p has relative ;; error d1, the outcome is -- 2.25.1