Merge pull request #1815 from ds26gte/log-of-small-nums

js-numbers.js: log() for nonintegral rational arguments = log(numr) - log(denr)

Author: ds26gte

src/js/base/js-numbers.js

@@ -802,6 +802,11 @@ define("pyret-base/js/js-numbers", function() {
     if (typeof(n) === 'number') {
       return Roughnum.makeInstance(Math.log(n), errbacks);
     }
+    if (isRational(n) && !isInteger(n)) {
+      return subtract(log(numerator(n, errbacks), errbacks),
+        log(denominator(n, errbacks), errbacks),
+        errbacks);
+    }
     var nFix = n.toFixnum();
     if (typeof(nFix) === 'number' && nFix !== Infinity) {
       return Roughnum.makeInstance(Math.log(nFix), errbacks);

tests/pyret/tests/test-numbers.arr

@@ -185,6 +185,9 @@ check:
   # Racket gives ~84709.80298615794, Wolfram ~84709.80298615795
   num-log(1e36789) is-roughly ~84709.80298615796
 
+  # Racket gives ~-171658.17010439213, Wolfram 171658.170104392139
+  num-log(1 / num-expt(9, num-expt(5, 7))) is-roughly ~-171658.17010439216
+
   # Racket, Wolfram match
   # commenting because arg calculation takes much time
   # num-log(num-expt(10, 1e5)) is-roughly ~230258.50929940457