From 91cd80f1fffb7e10e4c1a94d89f4b4c67943b336 Mon Sep 17 00:00:00 2001
From: Taylor R Campbell <campbell@mumble.net>
Date: Sat, 27 Oct 2018 02:04:17 +0000
Subject: [PATCH] Use named constants rather than magic numbers.

---
 src/runtime/arith.scm | 22 +++++++++++-----------
 1 file changed, 11 insertions(+), 11 deletions(-)

diff --git a/src/runtime/arith.scm b/src/runtime/arith.scm
index 990af61f5..507005292 100644
--- a/src/runtime/arith.scm
+++ b/src/runtime/arith.scm
@@ -2025,8 +2025,8 @@ USA.
 ;;; log(1 + e^x)
 (define (log1pexp x)
   (guarantee-real x 'log1pexp)
-  (cond ((<= x -745) 0.)
-	((<= x -37) (exp x))
+  (cond ((<= x flo:subnormal-exponent-min-base-e) 0.)
+	((<= x flo:log-epsilon) (exp x))
 	((<= x 18) (log1p (exp x)))
 	((<= x 33.3) (+ x (exp (- x))))
 	(else x)))
@@ -2100,11 +2100,11 @@ USA.
 ;;; 1/(1 + e^{-x}) = e^x/(1 + e^x)
 (define (logistic x)
   (guarantee-real x 'logistic)
-  (cond ((<= x -745)
+  (cond ((<= x flo:subnormal-exponent-min-base-e)
          ;; e^x/(1 + e^x) < e^x < smallest positive float.  (XXX Should
          ;; raise inexact and underflow here.)
          0.)
-	((<= x -37)
+	((<= x flo:log-epsilon)
          ;; e^x < eps, so
          ;;
          ;;     |e^x - e^x/(1 + e^x)|/|e^x/(1 + e^x)|
@@ -2115,7 +2115,7 @@ USA.
          ;;      < eps.
          ;;
          (exp x))
-	((<= x 37)
+	((<= x (- flo:log-epsilon))
          ;; e^{-x} > 0, so 1 + e^{-x} > 1, and 0 < 1/(1 + e^{-x}) < 1;
          ;; further, since e^{-x} < 1 + e^{-x}, we also have 0 <
          ;; e^{-x}/(1 + e^{-x}) < 1.  Thus, if exp has relative error
@@ -2151,8 +2151,8 @@ USA.
          ;;
          (/ 1 (+ 1 (exp (- x)))))
 	(else
-         ;; If x > 37, then e^{-x} < eps, so the relative error of 1
-         ;; from 1/(1 + e^{-x}) is
+         ;; If x > -log eps, then e^{-x} < eps, so the relative error
+         ;; of 1 from 1/(1 + e^{-x}) is
          ;;
          ;;     |1/(1 + e^{-x}) - 1|/|1/(1 + e^{-x})|
          ;;      = |e^{-x}/(1 + e^{-x})|/|1/(1 + e^{-x})|
@@ -2240,16 +2240,16 @@ USA.
 ;;; log e^t/(1 - e^t) = logit(e^t)
 (define (logit-exp t)
   (guarantee-real t 'logit-exp)
-  (cond ((<= t -37)
-         ;; e^t < eps/2, so since log(e^t/(1 - e^t)) = t - log(1 -
-         ;; e^t), and |log(1 - e^t)| < 1 < |t|, we have
+  (cond ((<= t flo:log-epsilon)
+         ;; e^t < eps, so since log(e^t/(1 - e^t)) = t - log(1 - e^t),
+         ;; and |log(1 - e^t)| < 1 < |t|, we have
          ;;
          ;;     |t - log(e^t/(1 - e^t))|/|log(e^t/(1 - e^t))|
          ;;      = |log(1 - e^t)|/|t - log(1 - e^t)|
          ;;     <= |log(1 - e^t)|
          ;;     <= 1/(1 - e^t)
          ;;     <= 2|e^t|
-         ;;     <= eps.
+         ;;     <= 2 eps.
          ;;
          t)
         ((<= logit-exp-boundary-lo t logit-exp-boundary-hi)
-- 
2.25.1