public void minval()
{
var samplesize = Pow2.T14;
var s1Range = Interval.closed(350.0, 1000.0);
var s1 = Random.Array <double>(samplesize, s1Range);
var s1Max = Observations.Load(s1).Max()[0];
NumericClaims.neq(s1Max, 0.0);
var zeroCount = s1.Count(x => x == 0);
Notify($"Found {zeroCount} zeroes");
}
public void CreateMt2203Generators()
{
var gencount = Pow2.T08;
var samplesize = Pow2.T16;
var seeds = Random.Array <uint>(gencount);
var streams = new MklRng[gencount];
for (var i = 0; i < gencount; i++)
{
streams[i] = rng.mt2203(seeds[i], i);
}
var bufferF64 = new double[samplesize];
var bufferU32 = new uint[samplesize];
var bufferI32 = new int[samplesize];
var ufRange = Interval.closed(1.0, 250.0);
for (var i = 0; i < gencount; i++)
{
var stream = streams[i];
sample.uniform(stream, ufRange, bufferF64);
Observations.Load(bufferF64, 1).Extrema();
var max = Observations.Load(bufferF64, 1).Max()[0];
NumericClaims.lteq(max, ufRange.Right);
NumericClaims.neq(max, 0);
sample.bits(stream, bufferU32);
sample.bernoulli(stream, .40, bufferI32);
for (var j = 0; j < samplesize; j++)
{
Claim.require(bufferI32[j] == 0 || bufferI32[j] == 1);
}
sample.gaussian(stream, .75, .75, bufferF64);
sample.laplace(stream, .5, .5, bufferF64);
}
for (var i = 0; i < gencount; i++)
{
streams[i].Dispose();
}
}
