/// <summary>
/// Sobol sequence.
/// </summary>
public static IEnumerable <Rational> SobolR(int polynomialCode, params int[] initialValues)
{
var primitivePolynomial = PolyGF2.FromCode(2UL * (ulong)polynomialCode ^ 1UL ^ (1UL << initialValues.Length));
return(SobolR(primitivePolynomial, initialValues));
}
/// <summary>
/// Creates a <see cref="Sobol"/> with specified initial
/// values and primitive polynomial coefficients packed into
/// bits, with leading and trailing coefficient omitted.
/// </summary>
/// <param name="polynomialCode">Bit 0 corresponds to x^1.</param>
/// <param name="initialValues">Initial values.</param>
public Sobol(int polynomialCode, params int[] initialValues)
: this(PolyGF2.FromCode(2UL * (ulong)polynomialCode ^ 1UL ^ (1UL << initialValues.Length))
.Exponents.ToArray(), initialValues)
{
}