private Sobol(Sobol other)
{
v = other.v;
value = other.value;
skip = other.skip;
ruler = Sequences.Ruler(skip).GetEnumerator();
}
public void SobolSkip()
{
IUniformSource s1 = new Sobol(new[] { 0, 1, 3 }, 1, 3, 7);
var s2 = (Sobol)s1.Clone();
for (int i = 0; i < 1000; i++)
{
s1.Next();
}
s2.Skip(1000);
for (int i = 0; i < 1000; i++)
{
Assert.AreEqual(s1.Next(), s2.Next());
}
}
public void SobolRational()
{
IUniformSource source = new Sobol(new[] { 0, 1, 3 }, 1, 3, 7);
var g = new UniformGenerator(source);
Assert.AreEqual((Rational)0 / 1, g.Rational());
Assert.AreEqual((Rational)1 / 2, g.Rational());
Assert.AreEqual((Rational)1 / 4, g.Rational());
Assert.AreEqual((Rational)3 / 4, g.Rational());
Assert.AreEqual((Rational)1 / 8, g.Rational());
Assert.AreEqual((Rational)5 / 8, g.Rational());
Assert.AreEqual((Rational)3 / 8, g.Rational());
Assert.AreEqual((Rational)7 / 8, g.Rational());
Assert.AreEqual((Rational)11 / 16, g.Rational());
Assert.AreEqual((Rational)3 / 16, g.Rational());
Assert.AreEqual((Rational)15 / 16, g.Rational());
}
public void SobolDouble()
{
IUniformSource source = new Sobol(new[] { 0, 1, 3 }, 1, 3, 7);
var g = new UniformGenerator(source);
Assert.AreEqual(0.0, g.Double());
Assert.AreEqual(0.5, g.Double());
Assert.AreEqual(0.25, g.Double());
Assert.AreEqual(0.75, g.Double());
Assert.AreEqual(0.125, g.Double());
Assert.AreEqual(0.625, g.Double());
Assert.AreEqual(0.375, g.Double());
Assert.AreEqual(0.875, g.Double());
Assert.AreEqual(0.6875, g.Double());
Assert.AreEqual(0.1875, g.Double());
Assert.AreEqual(0.9375, g.Double());
}