private static void test01()
//****************************************************************************80
//
// Purpose:
//
// TEST01 tests NIEDERREITER2.
//
// Licensing:
//
// This code is distributed under the GNU LGPL license.
//
// Modified:
//
// 02 March 2007
//
// Author:
//
// John Burkardt
//
{
int dim_num;
Niederreiter2.Niederreiter2Data data = new();
Console.WriteLine("");
Console.WriteLine("TEST01");
Console.WriteLine(" NIEDERREITER2 computes the next element of ");
Console.WriteLine(" a Niederreiter quasirandom sequence using base 2.");
Console.WriteLine("");
Console.WriteLine(" In this test, we call NIEDERREITER2 repeatedly.");
for (dim_num = 2; dim_num <= 4; dim_num++)
{
double[] r = new double[dim_num];
int seed = 0;
Console.WriteLine("");
Console.WriteLine(" Using dimension DIM_NUM = " + dim_num + "");
Console.WriteLine("");
Console.WriteLine(" Seed Seed Niederreiter2");
Console.WriteLine(" In Out");
Console.WriteLine("");
int i;
for (i = 0; i <= 110; i++)
{
int seed_in = seed;
Niederreiter2.niederreiter2(ref data, dim_num, ref seed, ref r);
int seed_out = seed;
switch (i)
{
case <= 11:
case >= 95:
{
string cout = seed_in.ToString(CultureInfo.InvariantCulture).PadLeft(6) + " ";
cout += seed_out.ToString(CultureInfo.InvariantCulture).PadLeft(6) + " ";
int j;
for (j = 0; j < dim_num; j++)
{
cout += r[j].ToString(CultureInfo.InvariantCulture).PadLeft(10) + " ";
}
Console.WriteLine(cout);
break;
}
case 12:
Console.WriteLine("......................");
break;
}
}
}
}
private static void test02()
//****************************************************************************80
//
// Purpose:
//
// TEST02 tests NIEDERREITER2.
//
// Licensing:
//
// This code is distributed under the GNU LGPL license.
//
// Modified:
//
// 02 March 2007
//
// Author:
//
// John Burkardt
//
{
const int DIM_NUM = 3;
int i;
int j;
double[] r = new double[DIM_NUM];
int seed_in;
int seed_out;
Niederreiter2.Niederreiter2Data data = new();
Console.WriteLine("");
Console.WriteLine("TEST02");
Console.WriteLine(" NIEDERREITER2 computes the next element of ");
Console.WriteLine(" a Niederreiter quasirandom sequence using base 2.");
Console.WriteLine("");
Console.WriteLine(" In this test, we demonstrate how the SEED can be");
Console.WriteLine(" manipulated to skip ahead in the sequence, or");
Console.WriteLine(" to come back to any part of the sequence.");
Console.WriteLine("");
Console.WriteLine(" Using dimension DIM_NUM = " + DIM_NUM + "");
int seed = 0;
Console.WriteLine("");
Console.WriteLine(" Seed Seed Niederreiter2");
Console.WriteLine(" In Out");
Console.WriteLine("");
for (i = 0; i <= 10; i++)
{
seed_in = seed;
Niederreiter2.niederreiter2(ref data, DIM_NUM, ref seed, ref r);
seed_out = seed;
string cout = seed_in.ToString(CultureInfo.InvariantCulture).PadLeft(6) + " ";
cout += seed_out.ToString(CultureInfo.InvariantCulture).PadLeft(6) + " ";
for (j = 0; j < DIM_NUM; j++)
{
cout += r[j].ToString(CultureInfo.InvariantCulture).PadLeft(10) + " ";
}
Console.WriteLine(cout);
}
Console.WriteLine("");
Console.WriteLine(" Jump ahead by increasing SEED:");
seed = 100;
Console.WriteLine("");
Console.WriteLine(" Seed Seed Niederreiter2");
Console.WriteLine(" In Out");
Console.WriteLine("");
for (i = 0; i <= 5; i++)
{
seed_in = seed;
Niederreiter2.niederreiter2(ref data, DIM_NUM, ref seed, ref r);
seed_out = seed;
seed_out = seed;
string cout = seed_in.ToString(CultureInfo.InvariantCulture).PadLeft(6) + " ";
cout += seed_out.ToString(CultureInfo.InvariantCulture).PadLeft(6) + " ";
for (j = 0; j < DIM_NUM; j++)
{
cout += r[j].ToString(CultureInfo.InvariantCulture).PadLeft(10) + " ";
}
Console.WriteLine(cout);
}
Console.WriteLine("");
Console.WriteLine(" Jump back by decreasing SEED:");
seed = 3;
Console.WriteLine("");
Console.WriteLine(" Seed Seed Niederreiter2");
Console.WriteLine(" In Out");
Console.WriteLine("");
for (i = 0; i <= 10; i++)
{
seed_in = seed;
Niederreiter2.niederreiter2(ref data, DIM_NUM, ref seed, ref r);
seed_out = seed;
seed_out = seed;
string cout = seed_in.ToString(CultureInfo.InvariantCulture).PadLeft(6) + " ";
cout += seed_out.ToString(CultureInfo.InvariantCulture).PadLeft(6) + " ";
for (j = 0; j < DIM_NUM; j++)
{
cout += r[j].ToString(CultureInfo.InvariantCulture).PadLeft(10) + " ";
}
Console.WriteLine(cout);
}
Console.WriteLine("");
Console.WriteLine(" Jump ahead by increasing SEED:");
seed = 98;
Console.WriteLine("");
Console.WriteLine(" Seed Seed Niederreiter2");
Console.WriteLine(" In Out");
Console.WriteLine("");
for (i = 1; i <= 5; i++)
{
seed_in = seed;
Niederreiter2.niederreiter2(ref data, DIM_NUM, ref seed, ref r);
seed_out = seed;
seed_out = seed;
string cout = seed_in.ToString(CultureInfo.InvariantCulture).PadLeft(6) + " ";
cout += seed_out.ToString(CultureInfo.InvariantCulture).PadLeft(6) + " ";
for (j = 0; j < DIM_NUM; j++)
{
cout += r[j].ToString(CultureInfo.InvariantCulture).PadLeft(10) + " ";
}
Console.WriteLine(cout);
}
}