public void CompareSequencesOfMersennes()
{
Random rnd = new Random();
for (int i = 0; i < 10; i++)
{
int seed = rnd.Next();
var m1 = new MT19937();
m1.init_genrand((ulong)seed);
var m2 = new RandomMT((ulong)seed);
var m3 = new SharpDevs.Randomization.MersenneTwister(seed);
var m4 = new Random(seed);
for (int j = 0; j < 10; j++)
{
var r1 = m1.genrand_int32();
var r2 = m2.RandomInt();
var r3 = m3.GetNextULong();
var r4 = m4.Next();
Console.WriteLine("{0}\t{1}\t{2}\t{3}", r1, r2, r3, r4);
//Assert.AreEqual(r1, r2);
//Assert.AreEqual(r2, r3);
}
Console.WriteLine();
}
}
public void CompareDistributionOfMersennes()
{
Random rnd = new Random();
int baseSeed = rnd.Next();
var m1 = new MT19937();
m1.init_genrand((ulong)baseSeed);
var m2 = new RandomMT((ulong)baseSeed);
var m3 = new SharpDevs.Randomization.MersenneTwister(baseSeed);
var m4 = new Random(baseSeed);
var maxRange = 1000;
var iterations = 10000000;
var dist1 = this.CalculateDistribution(() => m1.RandomRange(0, maxRange - 1), maxRange, iterations);
var dist2 = this.CalculateDistribution(() => m2.RandomRange(0, maxRange - 1), maxRange, iterations);
var dist3 = this.CalculateDistribution(() => m3.GetNext(0, maxRange), maxRange, iterations);
var dist4 = this.CalculateDistribution(() => m4.Next(0, maxRange), maxRange, iterations);
var distribution1 = this.AnalyzeDistribution(dist1, iterations);
var distribution2 = this.AnalyzeDistribution(dist2, iterations);
var distribution3 = this.AnalyzeDistribution(dist3, iterations);
var distribution4 = this.AnalyzeDistribution(dist4, iterations);
Console.WriteLine("{0}\t{1}\t{2}\t{3}",
distribution1, distribution2, distribution3, distribution4);
var expectedFill = (decimal)iterations / maxRange;
Console.WriteLine("Przy oczekiwanym napełnieniu: " + expectedFill);
Console.WriteLine("Co daje procentowo:");
Console.WriteLine("{0}\t{1}\t{2}\t{3}",
distribution1 / expectedFill * 100,
distribution2 / expectedFill * 100,
distribution3 / expectedFill * 100,
distribution4 / expectedFill * 100);
}