log(a/b), not log(a)/log(b).
authorTaylor R Campbell <campbell@mumble.net>
Thu, 25 Oct 2018 05:20:53 +0000 (05:20 +0000)
committerTaylor R Campbell <campbell@mumble.net>
Thu, 25 Oct 2018 05:20:53 +0000 (05:20 +0000)
src/runtime/arith.scm

index 311fd583ebf2b90a56cfc4f84c7d4324c5cebec5..60459c49e23ebb2a818084f65e4f324b19fcd7c9 100644 (file)
@@ -2231,8 +2231,8 @@ USA.
 (define (logit-exp t)
   (guarantee-real t 'logit-exp)
   (cond ((<= t -37)
-         ;; e^t < eps/2, so since log(e^t)/log(1 - e^t) = t -
-         ;; log(1 - e^t), we have
+         ;; e^t < eps/2, so since log(e^t/(1 - e^t)) = t - log(1 -
+         ;; e^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)|.