示例#1
0
    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) + "");
        }
    }
示例#2
0
    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) + "");
        }
    }