Class Philox4x64
- All Implemented Interfaces:
ArbitrarilyJumpableUniformRandomProvider, RandomLongSource, JumpableUniformRandomProvider, LongJumpableUniformRandomProvider, RestorableUniformRandomProvider, UniformRandomProvider
This is a member of the Philox family of generators. Memory footprint is 384 bits and the period is 2258.
Jumping in the sequence is essentially instantaneous. This generator provides both subsequences and arbitrary jumps for easy parallelization.
References:
- Salmon, J.K. et al (2011) Parallel Random Numbers: As Easy as 1,2,3.
- Since:
- 1.7
-
Field Summary
FieldsModifier and TypeFieldDescriptionprivate final long[]Output buffer.private intOutput buffer index.private longCounter 0.private longCounter 1.private longCounter 2.private longCounter 3.private longKey low bits.private longKey high bits.private static final intThe base-2 logarithm of the period.private static final doubleThe period of 2^258 as a double.private static final intInternal buffer size.private static final longPhilox 64-bit mixing constant for counter 0.private static final longPhilox 64-bit mixing constant for counter 1.private static final longPhilox 64-bit constant for key 0.private static final longPhilox 64-bit constant for key 1.private static final intNumber of state variables.private static final double2^54. -
Constructor Summary
ConstructorsModifierConstructorDescriptionPhilox4x64(long[] seed) Creates a new instance given 6 long numbers containing, key (first two longs) and the counter (next 4 longs, low bits = first long).privatePhilox4x64(Philox4x64 source) Copy constructor. -
Method Summary
Modifier and TypeMethodDescriptioncopyAndJump(int skip, long[] increment) Copy the generator and advance the internal state.private voidFinish the jump of this generator.private static intgetBufferPositionIncrement(double distance) Gets the buffer position increment from the jump distance.private static long[]getCounterIncrement(double distance) Gets the counter increment from the jump distance.protected byte[]Creates a snapshot of the RNG state.private voidIncrement the counter by one.jump()Creates a copy of the UniformRandomProvider and then advances the state of the current instance.jump(double distance) Creates a copy of theArbitrarilyJumpableUniformRandomProviderand then advances the state cycle of the current instance by the specifieddistance.jumpPowerOfTwo(int logDistance) Creates a copy of theArbitrarilyJumpableUniformRandomProviderand then advances the state cycle of the current instance by a distance equal to 2logDistance.jumps(double distance) Returns an effectively unlimited stream of new random generators, each of which implements theArbitrarilyJumpableUniformRandomProviderinterface.longJump()Creates a copy of the JumpableUniformRandomProvider and then advances the state of the current instance.longnext()Return the next random value.private voidrand10()Perform 10 rounds, using counter0, counter1, counter2, counter3 as starting point.private voidsetState(long[] state) Copies the state from the array into the generator state.protected voidsetStateInternal(byte[] s) Resets the RNG to the givenstate.private static voidsingleRound(long[] counter, long key0, long key1) Performs a single round of philox.Methods inherited from class LongProvider
nextBoolean, nextInt, nextLong, resetCachedStateMethods inherited from class BaseProvider
checkIndex, checkStateSize, composeStateInternal, extendSeed, extendSeed, fillState, fillState, restoreState, saveState, splitStateInternal, toStringMethods inherited from class Object
clone, equals, finalize, getClass, hashCode, notify, notifyAll, wait, wait, waitMethods inherited from interface ArbitrarilyJumpableUniformRandomProvider
jumpsMethods inherited from interface JumpableUniformRandomProvider
jumps, jumpsMethods inherited from interface LongJumpableUniformRandomProvider
longJumps, longJumpsMethods inherited from interface UniformRandomProvider
doubles, doubles, doubles, doubles, ints, ints, ints, ints, longs, longs, longs, longs, nextBoolean, nextBytes, nextBytes, nextDouble, nextDouble, nextDouble, nextFloat, nextFloat, nextFloat, nextInt, nextInt, nextInt, nextLong, nextLong, nextLong
-
Field Details
-
PHILOX_M0
private static final long PHILOX_M0Philox 64-bit mixing constant for counter 0.- See Also:
-
PHILOX_M1
private static final long PHILOX_M1Philox 64-bit mixing constant for counter 1.- See Also:
-
PHILOX_W0
private static final long PHILOX_W0Philox 64-bit constant for key 0.- See Also:
-
PHILOX_W1
private static final long PHILOX_W1Philox 64-bit constant for key 1.- See Also:
-
PHILOX_BUFFER_SIZE
private static final int PHILOX_BUFFER_SIZEInternal buffer size.- See Also:
-
STATE_SIZE
private static final int STATE_SIZENumber of state variables.- See Also:
-
LOG_PERIOD
private static final int LOG_PERIODThe base-2 logarithm of the period.- See Also:
-
PERIOD
private static final double PERIODThe period of 2^258 as a double.- See Also:
-
TWO_POW_54
private static final double TWO_POW_542^54. Threshold for a double that cannot have the 2 least significant bits set when converted to a long.- See Also:
-
counter0
private long counter0Counter 0. -
counter1
private long counter1Counter 1. -
counter2
private long counter2Counter 2. -
counter3
private long counter3Counter 3. -
buffer
private final long[] bufferOutput buffer. -
key0
private long key0Key low bits. -
key1
private long key1Key high bits. -
bufferPosition
private int bufferPositionOutput buffer index. When at the end of the buffer the counter is incremented and the buffer regenerated.
-
-
Constructor Details
-
Philox4x64
public Philox4x64(long[] seed) Creates a new instance given 6 long numbers containing, key (first two longs) and the counter (next 4 longs, low bits = first long). The counter is not scrambled and may be used to create contiguous blocks with size a multiple of 4 longs. For example, setting seed[2] = 1 is equivalent to start with seed[2]=0 and callingnext()4 times.- Parameters:
seed- Array of size 6 defining key0,key1,counter0,counter1,counter2,counter3. If the size is smaller, zero values are assumed.
-
Philox4x64
-
-
Method Details
-
setState
private void setState(long[] state) Copies the state from the array into the generator state.- Parameters:
state- New state.
-
getStateInternal
protected byte[] getStateInternal()Creates a snapshot of the RNG state.- Overrides:
getStateInternalin classLongProvider- Returns:
- the internal state.
-
setStateInternal
protected void setStateInternal(byte[] s) Resets the RNG to the givenstate.- Overrides:
setStateInternalin classLongProvider- Parameters:
s- State (previously obtained by a call toBaseProvider.getStateInternal()).- See Also:
-
next
public long next()Return the next random value.- Specified by:
nextin interfaceRandomLongSource- Returns:
- the next random value.
-
incrementCounter
private void incrementCounter()Increment the counter by one. -
rand10
private void rand10()Perform 10 rounds, using counter0, counter1, counter2, counter3 as starting point. It updates the buffer member variable, but no others. -
singleRound
private static void singleRound(long[] counter, long key0, long key1) Performs a single round of philox.- Parameters:
counter- Counter, which will be updated after each call.key0- Key low bits.key1- Key high bits.
-
jump
Creates a copy of the UniformRandomProvider and then advances the state of the current instance. The copy is returned.The current state will be advanced in a single operation by the equivalent of a number of sequential calls to a method that updates the state of the provider. The size of the jump is implementation dependent.
Repeat invocations of this method will create a series of generators that are uniformly spaced at intervals of the output sequence. Each generator provides non-overlapping output for the length of the jump for use in parallel computations.
The jump size is the equivalent of 2130 calls to
nextLong(). It can provide up to 2128 non-overlapping subsequences.- Specified by:
jumpin interfaceJumpableUniformRandomProvider- Returns:
- A copy of the current state.
-
longJump
Creates a copy of the JumpableUniformRandomProvider and then advances the state of the current instance. The copy is returned.The current state will be advanced in a single operation by the equivalent of a number of sequential calls to a method that updates the state of the provider. The size of the long jump is implementation dependent.
Repeat invocations of this method will create a series of generators that are uniformly spaced at intervals of the output sequence. Each generator provides non-overlapping output for the length of the long jump for use in parallel computations.
The returned copy may be jumped
m / ntimes before overlap with the current instance wheremis the long jump length andnis the jump length of theJumpableUniformRandomProvider.jump()method.The jump size is the equivalent of 2194 calls to
nextLong(). It can provide up to 264 non-overlapping subsequences of length 2194; each subsequence can provide up to 264 non-overlapping subsequences of length 2130 using thejump()method.- Specified by:
longJumpin interfaceLongJumpableUniformRandomProvider- Returns:
- A copy of the current state.
-
jump
Description copied from interface:ArbitrarilyJumpableUniformRandomProviderCreates a copy of theArbitrarilyJumpableUniformRandomProviderand then advances the state cycle of the current instance by the specifieddistance. The copy is returned.The current state will be advanced in a single operation by the equivalent of a number of sequential calls to a method that updates the state cycle of the provider.
Repeat invocations of this method will create a series of generators that are uniformly spaced at intervals of the output sequence. Each generator provides non-overlapping output for the length specified by
distancefor use in parallel computations.- Specified by:
jumpin interfaceArbitrarilyJumpableUniformRandomProvider- Parameters:
distance- Distance to jump forward with the state cycle.- Returns:
- A copy of the current state.
-
jumpPowerOfTwo
Description copied from interface:ArbitrarilyJumpableUniformRandomProviderCreates a copy of theArbitrarilyJumpableUniformRandomProviderand then advances the state cycle of the current instance by a distance equal to 2logDistance. The copy is returned.The current state will be advanced in a single operation by the equivalent of a number of sequential calls to a method that updates the state cycle of the provider.
Repeat invocations of this method will create a series of generators that are uniformly spaced at intervals of the output sequence. Each generator provides non-overlapping output for the length specified by 2
logDistancefor use in parallel computations.- Specified by:
jumpPowerOfTwoin interfaceArbitrarilyJumpableUniformRandomProvider- Parameters:
logDistance- Base-2 logarithm of the distance to jump forward with the state cycle.- Returns:
- A copy of the current state.
-
jumps
Returns an effectively unlimited stream of new random generators, each of which implements theArbitrarilyJumpableUniformRandomProviderinterface. The generators are output at integer multiples of the specified jumpdistancein the generator's state cycle.- Specified by:
jumpsin interfaceArbitrarilyJumpableUniformRandomProvider- Parameters:
distance- Distance to jump forward with the state cycle.- Returns:
- a stream of random generators.
-
getBufferPositionIncrement
private static int getBufferPositionIncrement(double distance) Gets the buffer position increment from the jump distance.- Parameters:
distance- Jump distance.- Returns:
- the buffer position increment
-
getCounterIncrement
private static long[] getCounterIncrement(double distance) Gets the counter increment from the jump distance.- Parameters:
distance- Jump distance.- Returns:
- the counter increment
-
copyAndJump
Copy the generator and advance the internal state. The copy is returned.This method: (1) assumes that the arguments have been validated; and (2) regenerates the output buffer if required.
- Parameters:
skip- Amount to skip the buffer position in [0, 3].increment- Unsigned 256-bit increment, least significant bits first.- Returns:
- the copy
-
finishJump
private void finishJump()Finish the jump of this generator. Resets the cached state and regenerates the output buffer if required.
-