private static void test02()
//****************************************************************************80
//
// Purpose:
//
// TEST02 tests FAURE.
//
// Licensing:
//
// This code is distributed under the GNU LGPL license.
//
// Modified:
//
// 23 January 2007
//
// Author:
//
// John Burkardt
//
{
const int DIM_NUM = 3;
int dim;
int i;
double[] r = new double[DIM_NUM];
int seed_in;
int seed_out;
Console.WriteLine("");
Console.WriteLine("TEST02");
Console.WriteLine(" FAURE computes the next element of a Faure sequence.");
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.");
int qs = Prime.prime_ge(DIM_NUM);
Console.WriteLine("");
Console.WriteLine(" Using dimension DIM_NUM = " + DIM_NUM + "");
Console.WriteLine(" The underlying base is QS = " + qs + "");
Console.WriteLine("");
Console.WriteLine(" Note that on the first call to FAURE, if");
Console.WriteLine(" SEED is negative, it is reset to a value that");
Console.WriteLine(" is the recommended starting point:");
int seed = -1;
Console.WriteLine("");
Console.WriteLine(" Seed Seed Faure");
Console.WriteLine(" In Out");
Console.WriteLine("");
FaureData data = new();
for (i = 1; i <= 5; i++)
{
seed_in = seed;
Faure.faure(ref data, DIM_NUM, ref seed, ref r);
seed_out = seed;
string cout = seed_in.ToString(CultureInfo.InvariantCulture).PadLeft(6) + " "
+ seed_out.ToString(CultureInfo.InvariantCulture).PadLeft(6) + " ";
for (dim = 0; dim < DIM_NUM; dim++)
{
cout += r[dim].ToString(CultureInfo.InvariantCulture).PadLeft(10) + " ";
}
Console.WriteLine(cout);
}
Console.WriteLine("");
Console.WriteLine(" However, if the input value of SEED is 0,");
Console.WriteLine(" then no initial skipping is done.");
seed = 0;
Console.WriteLine("");
Console.WriteLine(" Seed Seed Faure");
Console.WriteLine(" In Out");
Console.WriteLine("");
for (i = 1; i <= 10; i++)
{
seed_in = seed;
Faure.faure(ref data, DIM_NUM, ref seed, ref r);
seed_out = seed;
string cout = seed_in.ToString(CultureInfo.InvariantCulture).PadLeft(6) + " "
+ seed_out.ToString(CultureInfo.InvariantCulture).PadLeft(6) + " ";
for (dim = 0; dim < DIM_NUM; dim++)
{
cout += r[dim].ToString(CultureInfo.InvariantCulture).PadLeft(10) + " ";
}
Console.WriteLine(cout);
}
Console.WriteLine("");
Console.WriteLine(" Jump ahead by increasing SEED:");
Console.WriteLine("");
seed = 100;
Console.WriteLine("");
Console.WriteLine(" Seed Seed Faure");
Console.WriteLine(" In Out");
Console.WriteLine("");
for (i = 1; i <= 5; i++)
{
seed_in = seed;
Faure.faure(ref data, DIM_NUM, ref seed, ref r);
seed_out = seed;
string cout = seed_in.ToString(CultureInfo.InvariantCulture).PadLeft(6) + " "
+ seed_out.ToString(CultureInfo.InvariantCulture).PadLeft(6) + " ";
for (dim = 0; dim < DIM_NUM; dim++)
{
cout += r[dim].ToString(CultureInfo.InvariantCulture).PadLeft(10) + " ";
}
Console.WriteLine(cout);
}
Console.WriteLine("");
Console.WriteLine(" Jump back by decreasing SEED:");
Console.WriteLine("");
seed = 3;
Console.WriteLine("");
Console.WriteLine(" Seed Seed Faure");
Console.WriteLine(" In Out");
Console.WriteLine("");
for (i = 1; i <= 10; i++)
{
seed_in = seed;
Faure.faure(ref data, DIM_NUM, ref seed, ref r);
seed_out = seed;
string cout = seed_in.ToString(CultureInfo.InvariantCulture).PadLeft(6) + " "
+ seed_out.ToString(CultureInfo.InvariantCulture).PadLeft(6) + " ";
for (dim = 0; dim < DIM_NUM; dim++)
{
cout += r[dim].ToString(CultureInfo.InvariantCulture).PadLeft(10) + " ";
}
Console.WriteLine(cout);
}
}
private static void test03()
//****************************************************************************80
//
// Purpose:
//
// TEST03 tests FAURE.
//
// Licensing:
//
// This code is distributed under the GNU LGPL license.
//
// Modified:
//
// 04 June 2007
//
// Author:
//
// John Burkardt
//
{
const int dim_base = 10;
int dim_num;
Console.WriteLine("");
Console.WriteLine("TEST03");
Console.WriteLine(" FAURE computes the next element of a Faure sequence.");
Console.WriteLine("");
Console.WriteLine(" In this test, we try some large dimensions.");
FaureData data = new();
for (dim_num = dim_base; dim_num <= 6 * dim_base; dim_num += dim_base)
{
double[] r = new double[dim_num];
int seed = -1;
int qs = Prime.prime_ge(dim_num);
Console.WriteLine("");
Console.WriteLine(" Using dimension DIM_NUM = " + dim_num + "");
Console.WriteLine(" The underlying base is QS = " + qs + "");
Console.WriteLine("");
Console.WriteLine(" Seed Seed Faure");
Console.WriteLine(" In Out");
Console.WriteLine("");
int i;
for (i = 1; i <= 2; i++)
{
int seed_in = seed;
Faure.faure(ref data, dim_num, ref seed, ref r);
int seed_out = seed;
Console.WriteLine(" " + seed_in.ToString(CultureInfo.InvariantCulture).PadLeft(8)
+ " " + seed_out.ToString(CultureInfo.InvariantCulture).PadLeft(8) + "");
string cout = " ";
int dim;
for (dim = 0; dim < dim_num; dim++)
{
cout += " " + r[dim].ToString(CultureInfo.InvariantCulture).PadLeft(10);
if ((dim + 1) % 5 == 0 || dim + 1 == dim_num)
{
Console.WriteLine(cout);
}
switch ((dim + 1) % 5)
{
case 0 when dim + 1 < dim_num:
cout += " ";
break;
}
}
}
}
}
private static void test01()
//****************************************************************************80
//
// Purpose:
//
// TEST01 tests FAURE.
//
// Licensing:
//
// This code is distributed under the GNU LGPL license.
//
// Modified:
//
// 23 January 2007
//
// Author:
//
// John Burkardt
//
{
const int DIM_MAX = 4;
int dim_num;
double[] r = new double[DIM_MAX];
Console.WriteLine("");
Console.WriteLine("TEST01");
Console.WriteLine(" FAURE computes the next element of a Faure sequence.");
Console.WriteLine("");
Console.WriteLine(" In this test, we call FAURE repeatedly.");
FaureData data = new();
for (dim_num = 2; dim_num <= DIM_MAX; dim_num++)
{
int seed = -1;
int qs = Prime.prime_ge(dim_num);
Console.WriteLine("");
Console.WriteLine(" Using dimension DIM_NUM = " + dim_num + "");
Console.WriteLine(" The underlying base is QS = " + qs + "");
Console.WriteLine("");
Console.WriteLine(" Seed Seed Faure");
Console.WriteLine(" In Out");
Console.WriteLine("");
int i;
for (i = 1; i <= 10; i++)
{
int seed_in = seed;
Faure.faure(ref data, dim_num, ref seed, ref r);
int seed_out = seed;
string cout = seed_in.ToString(CultureInfo.InvariantCulture).PadLeft(6) + " "
+ seed_out.ToString(CultureInfo.InvariantCulture).PadLeft(6) + " ";
int dim;
for (dim = 0; dim < dim_num; dim++)
{
cout += r[dim].ToString(CultureInfo.InvariantCulture).PadLeft(10) + " ";
}
Console.WriteLine(cout);
}
}
}