private static void test02() //****************************************************************************80 // // Purpose: // // TEST02 tests VAN_DER_CORPUT_SEQUENCE. // // Licensing: // // This code is distributed under the GNU LGPL license. // // Modified: // // 09 August 2016 // // Author: // // John Burkardt // { const int N = 1; int i; double[] r = Array.Empty <double>(); Console.WriteLine(""); Console.WriteLine("TEST02"); Console.WriteLine(" VAN_DER_CORPUT_SEQUENCE computes several elements of"); Console.WriteLine(" a van der Corput sequence on a single call."); Console.WriteLine(""); Console.WriteLine(" In this test, we call VAN_DER_CORPUT_SEQUENCE once."); const int base_ = 2; VanDerCorput.van_der_corput_base_set(base_); const int seed = 0; VanDerCorput.van_der_corput_seed_set(seed); VanDerCorput.van_der_corput_sequence(N, ref r); Console.WriteLine(""); Console.WriteLine(" Element van der Corput"); Console.WriteLine(""); for (i = 0; i < N; i++) { Console.WriteLine(i.ToString(CultureInfo.InvariantCulture).PadLeft(6) + " " + r[i].ToString(CultureInfo.InvariantCulture).PadLeft(10) + ""); } }
private static void test045() //****************************************************************************80 // // Purpose: // // TEST045 tests VAN_DER_CORPUT, VAN_DER_CORPUT_SEED_GET, VAN_DER_CORPUT_SEED_SET. // // Licensing: // // This code is distributed under the GNU LGPL license. // // Modified: // // 09 August 2016 // // Author: // // John Burkardt // { int i; double r; Console.WriteLine(""); Console.WriteLine("TEST045"); Console.WriteLine(" VAN_DER_CORPUT_SEED_GET gets the current seed."); Console.WriteLine(" VAN_DER_CORPUT_SEED_SET sets the current seed."); Console.WriteLine(""); Console.WriteLine(" The van der Corput base is usually a prime, but this is"); Console.WriteLine(" not required."); Console.WriteLine(""); Console.WriteLine(" In this test, we compute a van der Corput sequence"); Console.WriteLine(" starting with the default seed, then check the seed,"); Console.WriteLine(" reset the seed, and recompute the sequence."); int base_ = 2; VanDerCorput.van_der_corput_base_set(base_); int seed = 0; VanDerCorput.van_der_corput_seed_set(seed); Console.WriteLine(""); Console.WriteLine(" All computations will use base " + base_ + "."); Console.WriteLine(""); seed = VanDerCorput.van_der_corput_seed_get(); Console.WriteLine(""); Console.WriteLine(" Set SEED = " + seed + ""); Console.WriteLine(""); Console.WriteLine(" Seed van der Corput"); Console.WriteLine(""); for (i = 1; i <= 10; i++) { seed = VanDerCorput.van_der_corput_seed_get(); r = VanDerCorput.van_der_corput(); Console.WriteLine(seed.ToString(CultureInfo.InvariantCulture).PadLeft(6) + " " + r.ToString(CultureInfo.InvariantCulture).PadLeft(10) + ""); } seed = VanDerCorput.van_der_corput_seed_get(); Console.WriteLine(""); Console.WriteLine(" Current seed is " + seed + ""); Console.WriteLine(""); seed = 100; VanDerCorput.van_der_corput_seed_set(seed); Console.WriteLine(""); Console.WriteLine(" Set SEED = " + seed + ""); Console.WriteLine(""); Console.WriteLine(" Seed van der Corput"); Console.WriteLine(""); for (i = 1; i <= 10; i++) { seed = VanDerCorput.van_der_corput_seed_get(); r = VanDerCorput.van_der_corput(); Console.WriteLine(seed.ToString(CultureInfo.InvariantCulture).PadLeft(6) + " " + r.ToString(CultureInfo.InvariantCulture).PadLeft(10) + ""); } seed = VanDerCorput.van_der_corput_seed_get(); Console.WriteLine(""); Console.WriteLine(" Current seed is " + seed + ""); Console.WriteLine(""); seed = 3; VanDerCorput.van_der_corput_seed_set(seed); Console.WriteLine(""); Console.WriteLine(" Reset seed to " + seed + ""); Console.WriteLine(""); Console.WriteLine(" Seed van der Corput"); Console.WriteLine(""); for (i = 1; i <= 10; i++) { seed = VanDerCorput.van_der_corput_seed_get(); r = VanDerCorput.van_der_corput(); Console.WriteLine(seed.ToString(CultureInfo.InvariantCulture).PadLeft(6) + " " + r.ToString(CultureInfo.InvariantCulture).PadLeft(10) + ""); } }