public void SobolInitialValues(int max)
{
int d = 2;
var g = new Generator("SobolInitialValue");
var v = new bool[max + 1][];
v[0] = new bool[max];
v[0][0] = true;
foreach (var p in PolyGF2.GetAllPrimitive().Take(max - 1))
{
v[d - 1] = new bool[max];
int s = p.Degree;
var m = new int[s];
var V = new SquareMatrixGF2(d);
while (true)
{
for (int i = 0; i < s; i++)
{
m[i] = g.Int32(2 << i) | 1;
}
var dv = Subrandom.SobolDirectionVectors(max, p, m);
for (int i = 0; i < max; i++)
{
v[d - 1][i] = dv[i] > (ulong)long.MaxValue;
}
for (int row = 0; row < d; row++)
{
for (int col = 0; col < d; col++)
{
V[row, col] = v[row][col];
}
}
bool det = V.GetDeterminant();
if (det)
{
break;
}
}
Trace.Write($"{d}\t{s}\t{p.GetCode() / 2 ^ (1UL << (s-1))}\t");
for (int i = 0; i < s; i++)
{
Trace.Write($"{m[i]} ");
}
Trace.WriteLine("");
d++;
}
}
