Fix bug in PermutedRadicalInverse()

Issue #52. Fix suggested by Andrew Kensler.
mmp · Jan 18, 2016 · e118a7a · dictoon · Jan 20, 2016 · e118a7a
1 parent b4491d4
commit e118a7a
Showing 1 changed file with 4 additions and 0 deletions.
diff --git a/src/core/montecarlo.h b/src/core/montecarlo.h
@@ -214,6 +214,10 @@ inline double PermutedRadicalInverse(uint32_t n, uint32_t base,
         n *= invBase;
         invBi *= invBase;
     }
+    // For the case where the permutation table permutes the digit 0 to
+    // another digit, account for the infinite sequence of that digit
+    // trailing at the end of the radical inverse value.
+    val += p[0] * base / (base-1.0f) * invBi;
     return val;
 }