--- jsr166/src/main/java/util/SplittableRandom.java 2013/07/14 08:06:49 1.10
+++ jsr166/src/main/java/util/SplittableRandom.java 2013/07/16 12:32:05 1.11
@@ -51,7 +51,9 @@ import java.util.stream.DoubleStream;
* href="http://www.phy.duke.edu/~rgb/General/dieharder.php"> version
* 3.31.1.) These tests validate only the methods for certain
* types and ranges, but similar properties are expected to hold, at
- * least approximately, for others as well.
+ * least approximately, for others as well. The period
+ * (length of any series of generated values before it repeats) is at
+ * least 264.
*
*
Method {@link #split} constructs and returns a new
* SplittableRandom instance that shares no mutable state with the
@@ -156,13 +158,16 @@ public class SplittableRandom {
* SplittableRandom. Unlike other cases, this split must be
* performed in a thread-safe manner. We use
* AtomicLong.compareAndSet as the (typically) most efficient
- * mechanism. To bootstrap, we start off using System.nanotime(),
- * and update using another "genuinely random" constant
- * DEFAULT_SEED_GAMMA. The default constructor uses GAMMA_GAMMA,
- * not 0, for its splitSeed argument (addGammaModGeorge(0,
- * GAMMA_GAMMA) == GAMMA_GAMMA) to reflect that each is split from
- * this root generator, even though the root is not explicitly
- * represented as a SplittableRandom.
+ * mechanism. To bootstrap, we start off using a function of the
+ * current System time as seed, and update using another
+ * "genuinely random" constant DEFAULT_SEED_GAMMA. The default
+ * constructor uses GAMMA_GAMMA, not 0, for its splitSeed argument
+ * (addGammaModGeorge(0, GAMMA_GAMMA) == GAMMA_GAMMA) to reflect
+ * that each is split from this root generator, even though the
+ * root is not explicitly represented as a SplittableRandom. When
+ * establishing the initial seed, we use both
+ * System.currentTimeMillis and System.nanoTime(), to avoid
+ * regularities that may occur if using either alone.
*/
/**
@@ -180,6 +185,12 @@ public class SplittableRandom {
private static final long DEFAULT_SEED_GAMMA = 0xBD24B73A95FB84D9L;
/**
+ * The value 13 with 64bit sign bit set. Used in the signed
+ * comparison in addGammaModGeorge.
+ */
+ private static final long BOTTOM13 = 0x800000000000000DL;
+
+ /**
* The least non-zero value returned by nextDouble(). This value
* is scaled by a random value of 53 bits to produce a result.
*/
@@ -189,7 +200,8 @@ public class SplittableRandom {
* The next seed for default constructors.
*/
private static final AtomicLong defaultSeedGenerator =
- new AtomicLong(System.nanoTime());
+ new AtomicLong(mix64(System.currentTimeMillis()) ^
+ mix64(System.nanoTime()));
/**
* The seed, updated only via method nextSeed.
@@ -217,15 +229,16 @@ public class SplittableRandom {
* George < 2^64; thus we need only a conditional, not a loop,
* to be sure of getting a representable value.
*
- * @param s a seed value
+ * Because Java comparison operators are signed, we implement this
+ * by conceptually offsetting seed values downwards by 2^63, so
+ * 0..13 is represented as Long.MIN_VALUE..BOTTOM13.
+ *
+ * @param s a seed value, viewed as a signed long
* @param g a gamma value, 13 <= g (as unsigned)
*/
private static long addGammaModGeorge(long s, long g) {
long p = s + g;
- if (Long.compareUnsigned(p, g) >= 0)
- return p;
- long q = p - 13L;
- return (Long.compareUnsigned(p, 13L) >= 0) ? q : (q + g);
+ return (p >= s) ? p : ((p >= BOTTOM13) ? p : p + g) - 13L;
}
/**
@@ -264,7 +277,7 @@ public class SplittableRandom {
do { // ensure gamma >= 13, considered as an unsigned integer
s = addGammaModGeorge(s, GAMMA_GAMMA);
g = mix64(s);
- } while (Long.compareUnsigned(g, 13L) < 0);
+ } while (g >= 0L && g < 13L);
this.gamma = g;
this.nextSplit = s;
}
@@ -403,7 +416,7 @@ public class SplittableRandom {
/**
* Creates a new SplittableRandom instance using the specified
* initial seed. SplittableRandom instances created with the same
- * seed generate identical sequences of values.
+ * seed in the same program generate identical sequences of values.
*
* @param seed the initial seed
*/
@@ -553,7 +566,7 @@ public class SplittableRandom {
* (inclusive) and one (exclusive)
*/
public double nextDouble() {
- return (nextLong() >>> 11) * DOUBLE_UNIT;
+ return (mix64(nextSeed()) >>> 11) * DOUBLE_UNIT;
}
/**
@@ -569,7 +582,7 @@ public class SplittableRandom {
public double nextDouble(double bound) {
if (!(bound > 0.0))
throw new IllegalArgumentException("bound must be positive");
- double result = nextDouble() * bound;
+ double result = (mix64(nextSeed()) >>> 11) * DOUBLE_UNIT * bound;
return (result < bound) ? result : // correct for rounding
Double.longBitsToDouble(Double.doubleToLongBits(bound) - 1);
}
@@ -591,6 +604,15 @@ public class SplittableRandom {
return internalNextDouble(origin, bound);
}
+ /**
+ * Returns a pseudorandom {@code boolean} value.
+ *
+ * @return a pseudorandom {@code boolean} value
+ */
+ public boolean nextBoolean() {
+ return mix32(nextSeed()) < 0;
+ }
+
// stream methods, coded in a way intended to better isolate for
// maintenance purposes the small differences across forms.
@@ -854,7 +876,7 @@ public class SplittableRandom {
* approach. The long and double versions of this class are
* identical except for types.
*/
- static class RandomIntsSpliterator implements Spliterator.OfInt {
+ static final class RandomIntsSpliterator implements Spliterator.OfInt {
final SplittableRandom rng;
long index;
final long fence;
@@ -908,7 +930,7 @@ public class SplittableRandom {
/**
* Spliterator for long streams.
*/
- static class RandomLongsSpliterator implements Spliterator.OfLong {
+ static final class RandomLongsSpliterator implements Spliterator.OfLong {
final SplittableRandom rng;
long index;
final long fence;
@@ -963,7 +985,7 @@ public class SplittableRandom {
/**
* Spliterator for double streams.
*/
- static class RandomDoublesSpliterator implements Spliterator.OfDouble {
+ static final class RandomDoublesSpliterator implements Spliterator.OfDouble {
final SplittableRandom rng;
long index;
final long fence;