/* *
* Title:
* Description:
* Copyright:
* Notice:
* @Version
*/
package umontreal.iro.lecuyer.rng;
import umontreal.iro.lecuyer.util.BitVector;
import umontreal.iro.lecuyer.util.BitMatrix;
import java.io.Serializable;
import java.io.ObjectOutputStream;
import java.io.FileOutputStream;
import java.io.FileNotFoundException;
import java.io.IOException;
/**
* This class implements the {@link RandomStream} interface via inheritance
* from {@link RandomStreamBase}. The backbone generator is a Well
* Equidistributed Long period Linear Random Number Generator (WELL),
* proposed by F. Panneton.
* The implemented generator is the WELL607, which has a state
* size of 607 bits
* and a period length of approximatively
* 2607. The values of V, W and Z are 2150,
* 2250 and 2400 respectively (see {@link RandomStream} for their
* definition). The seed of the RNG, and the state of a stream at any given
* step, is a 19-dimensional vector of 32-bit integers.
* The output of nextValue has 32 bits of precision.
*
*
* This implementation requires the use of about
* 100K of memory to run. This memory is shared between all instances of the
* class, and is only loaded when the first instance is created.
*
*/
public class WELL607 extends WELL607base {
private static final long serialVersionUID = 70510L;
//La date de modification a l'envers, lire 10/05/2007
private static int[] curr_stream = {0xD6AFB71C, 0x82ADB18E, 0x326E714E,
0xB1EE42B6, 0xF1A834ED, 0x04AE5721,
0xC5EA2843, 0xFA04116B, 0x6ACE14EF,
0xCD5781A0, 0x6B1F731C, 0x7E3B8E3D,
0x8B34DE2A, 0x74EC15F5, 0x84EBC216,
0x83EA2C61, 0xE4A83B1E, 0xA5D82CB9,
0x9E1A6C89};
// P(z) = {0x987b2631, 0x2e33283d, 0x6a398474, 0xe9d24da1,
// 0x31235359, 0x6a2baf48, 0x7f97efd4, 0x468280f4,
// 0x7d9d9424, 0xa3238f8e, 0xe3edb4ef, 0x0e0a25f7,
// 0x92c4dff5, 0x55d0b8da, 0x7b982dec, 0xa06c078f,
// 0x38b65c31, 0xc8c3788d, 0x8000b200}
// Ce tableau represente les 512 coefficients du polynome suivant
// (z^(2^250) mod P(z)) mod 2
// P(z) est le polynome caracteristique du generateur.
static final int [] pw = new int[]
{0x83167621, 0x6b5515c8, 0x61a62bd2, 0xbceaa78f,
0xac04b304, 0x28a75ea4, 0xa9104058, 0x595ea53b,
0x35687e95, 0x7f8eca9b, 0x30beffb8, 0xc61e6111,
0x284ee30e, 0x4e9cd901, 0x659633ba, 0x344cc69e,
0xd6052ac1, 0x5d508b69, 0x62cf130};
// Ce tableau represente les 512 coefficients du polynome suivant
// (z^(2^400) mod P(z)) mod 2
// P(z) est le polynome caracteristique du generateur.
static final int [] pz = new int[]
{0x70b2bdee, 0x595828f1, 0x85a17885, 0x5100c7b2,
0xd3333da2, 0xb42857de, 0xf8a7a4a7, 0xabad2a33,
0xa2580cf, 0xf94c465e, 0x7df951d5, 0x35467053,
0xb3c9a4e, 0x6a33977, 0x443910e, 0xc25aec3d,
0xeb72e8c5, 0x8873b01, 0x7da57636};
/**
* Constructs a new stream.
*
*/
public WELL607() {
initialisation();
state = new int[BUFFER_SIZE];
stream = new int[R];
substream = new int[R];
for(int i = 0; i < R; i++)
stream[i] = curr_stream[i];
// advanceSeed(curr_stream, Apz);
advanceSeed(curr_stream, pz);
resetStartStream();
}
/**
* Constructs a new stream with the identifier name
* (used in the toString method).
*
* @param name name of the stream
*
*/
public WELL607 (String name) {
this();
this.name = name;
}
/**
* Sets the initial seed of the class WELL607 to the 19
* integers of the vector seed[0..18].
* This will be the initial seed of the next created stream.
* At least one of the integers must not be zero and if this integer is
* the last one, it must not be equal to 0x80000000.
*
* @param seed array of 19 elements representing the seed
*
*
*/
public static void setPackageSeed (int seed[]) {
verifySeed(seed);
for(int i = 0; i < R; i++)
curr_stream[i] = seed[i];
}
/**
* This method is discouraged for normal use.
* Initializes the stream at the beginning of a stream with the initial
* seed seed[0..18]. The seed must satisfy the same
* conditions as in setPackageSeed.
* This method only affects the specified stream; the others are not
* modified. Hence after calling this method, the beginning of the streams
* will no longer be spaced Z values apart.
* For this reason, this method should only be used in very exceptional cases;
* proper use of the reset... methods and of the stream constructor is
* preferable.
*
* @param seed array of 19 elements representing the seed
*
*
*/
public void setSeed (int seed[]) {
verifySeed(seed);
for(int i = 0; i < R; i++)
stream[i] = seed[i];
resetStartStream();
}
/**
* Returns the current state of the stream, represented as an
* array of 19 integers.
*
* @return the current state of the stream
*
*/
public int[] getState() {
return super.getState();
}
public void resetStartStream() {
for(int i = 0; i < R; i++)
substream[i] = stream[i];
resetStartSubstream();
}
public void resetStartSubstream() {
state_i= 0;
for(int i = 0; i < R; i++)
state[i] = substream[i];
}
public void resetNextSubstream() {
/// advanceSeed(substream, Apw);
advanceSeed(substream, pw);
resetStartSubstream();
}
public String toString() {
StringBuffer sb = new StringBuffer();
if(name == null)
sb.append("The state of this WELL607 is : {");
else
sb.append("The state of " + name + " is : {");
sb.append(super.stringState());
return sb.toString();
}
protected double nextValue() {
long result = nextInt();
if(result <= 0)
result += 0x100000000L;
return result * NORM;
}
/**
* This method is only meant to be used during the compilation process.
* It is used to create the resource file the class need in order to
* run.
*
*/
/**
* Clones the current generator and return its copy.
*
* @return A deep copy of the current generator
*
*/
public WELL607 clone() {
WELL607 retour = null;
retour = (WELL607)super.clone();
retour.state = new int[BUFFER_SIZE];
retour.substream = new int[R];
retour.stream = new int[R];
for (int i = 0; i>>19)) ^ (state[M1]^(state[M1]>>>11));
z2 = (state[M2]^(state[M2]<<(14))) ^ state[M3];
state[0] = z1 ^ z2;
state[BUFFER_SIZE - 1] = (z0^(z0>>>18)) ^
z1 ^ (state[0]^(state[0]<<5));
//put the state vector in vect (state_i = -1)
for(int j = 0; j < R; j++)
vect[j] = state[(j - 1) & MASK_STATE];
bv[i] = new BitVector(vect, NUM_BITS);
}
STp0 = (new BitMatrix(bv)).transpose();
STpw = STp0.power2e(w);
STpz = STpw.power2e(v);
try {
FileOutputStream fos = new FileOutputStream(args[0]);
ObjectOutputStream oos = new ObjectOutputStream(fos);
oos.writeObject(STpw);
oos.writeObject(STpz);
oos.close();
} catch(FileNotFoundException e) {
System.err.println("Couldn't create " + args[0]);
System.err.println(e);
} catch(IOException e) {
System.err.println(e);
}
}
}
Ads help maintain this website.