June 7 2010 2:12:48.651 PM NIEDERREITER_PRB FORTRAN90 version Test the NIEDERREITER routines. TEST01 NIEDERREITER computes the next element of a Niederreiter quasirandom sequence using base BASE. In this test, we call NIEDERREITER repeatedly. Using base BASE = 2 Using dimension DIM_NUM = 2 Seed Seed Niederreiter In Out 0 1 0.5000 0.5000 1 2 0.2500 0.7500 2 3 0.7500 0.2500 3 4 0.1250 0.6250 4 5 0.6250 0.1250 5 6 0.3750 0.3750 6 7 0.8750 0.8750 7 8 0.0625 0.9375 8 9 0.5625 0.4375 9 10 0.3125 0.1875 10 11 0.8125 0.6875 11 12 0.1875 0.3125 ...................... 95 96 0.0234 0.3984 96 97 0.5234 0.8984 97 98 0.2734 0.6484 98 99 0.7734 0.1484 99 100 0.1484 0.7734 100 101 0.6484 0.2734 101 102 0.3984 0.0234 102 103 0.8984 0.5234 103 104 0.0859 0.5859 104 105 0.5859 0.0859 105 106 0.3359 0.3359 106 107 0.8359 0.8359 107 108 0.2109 0.2109 108 109 0.7109 0.7109 109 110 0.4609 0.9609 110 111 0.9609 0.4609 Using dimension DIM_NUM = 3 Seed Seed Niederreiter In Out 0 1 0.5000 0.5000 0.7500 1 2 0.2500 0.7500 0.5625 2 3 0.7500 0.2500 0.3125 3 4 0.1250 0.6250 0.4375 4 5 0.6250 0.1250 0.6875 5 6 0.3750 0.3750 0.8750 6 7 0.8750 0.8750 0.1250 7 8 0.0625 0.9375 0.9531 8 9 0.5625 0.4375 0.2031 9 10 0.3125 0.1875 0.3906 10 11 0.8125 0.6875 0.6406 11 12 0.1875 0.3125 0.5156 ...................... 95 96 0.0234 0.3984 0.6484 96 97 0.5234 0.8984 0.3984 97 98 0.2734 0.6484 0.2109 98 99 0.7734 0.1484 0.9609 99 100 0.1484 0.7734 0.8359 100 101 0.6484 0.2734 0.0859 101 102 0.3984 0.0234 0.2734 102 103 0.8984 0.5234 0.5234 103 104 0.0859 0.5859 0.3203 104 105 0.5859 0.0859 0.5703 105 106 0.3359 0.3359 0.7578 106 107 0.8359 0.8359 0.0078 107 108 0.2109 0.2109 0.1328 108 109 0.7109 0.7109 0.8828 109 110 0.4609 0.9609 0.6953 110 111 0.9609 0.4609 0.4453 Using dimension DIM_NUM = 4 Seed Seed Niederreiter In Out 0 1 0.5000 0.5000 0.7500 0.8750 1 2 0.2500 0.7500 0.5625 0.7656 2 3 0.7500 0.2500 0.3125 0.1406 3 4 0.1250 0.6250 0.4375 0.5469 4 5 0.6250 0.1250 0.6875 0.4219 5 6 0.3750 0.3750 0.8750 0.2813 6 7 0.8750 0.8750 0.1250 0.6563 7 8 0.0625 0.9375 0.9531 0.2344 8 9 0.5625 0.4375 0.2031 0.8594 9 10 0.3125 0.1875 0.3906 0.9688 10 11 0.8125 0.6875 0.6406 0.0938 11 12 0.1875 0.3125 0.5156 0.6875 ...................... 95 96 0.0234 0.3984 0.6484 0.8516 96 97 0.5234 0.8984 0.3984 0.2266 97 98 0.2734 0.6484 0.2109 0.1172 98 99 0.7734 0.1484 0.9609 0.9922 99 100 0.1484 0.7734 0.8359 0.3359 100 101 0.6484 0.2734 0.0859 0.7109 101 102 0.3984 0.0234 0.2734 0.5703 102 103 0.8984 0.5234 0.5234 0.4453 103 104 0.0859 0.5859 0.3203 0.8984 104 105 0.5859 0.0859 0.5703 0.0234 105 106 0.3359 0.3359 0.7578 0.1328 106 107 0.8359 0.8359 0.0078 0.7578 107 108 0.2109 0.2109 0.1328 0.4141 108 109 0.7109 0.7109 0.8828 0.5391 109 110 0.4609 0.9609 0.6953 0.6797 110 111 0.9609 0.4609 0.4453 0.3047 TEST01 NIEDERREITER computes the next element of a Niederreiter quasirandom sequence using base BASE. In this test, we call NIEDERREITER repeatedly. Using base BASE = 3 Using dimension DIM_NUM = 2 Seed Seed Niederreiter In Out 0 1 0.3333 0.3333 1 2 0.6667 0.6667 2 3 0.1111 0.7778 3 4 0.4444 0.1111 4 5 0.7778 0.4444 5 6 0.2222 0.5556 6 7 0.5556 0.8889 7 8 0.8889 0.2222 8 9 0.0370 0.4815 9 10 0.3704 0.8148 10 11 0.7037 0.1481 11 12 0.1481 0.2593 ...................... 95 96 0.2634 0.2881 96 97 0.5967 0.6214 97 98 0.9300 0.9547 98 99 0.0782 0.2140 99 100 0.4115 0.5473 100 101 0.7449 0.8807 101 102 0.1893 0.9918 102 103 0.5226 0.3251 103 104 0.8560 0.6584 104 105 0.3004 0.4362 105 106 0.6337 0.7695 106 107 0.9671 0.1029 107 108 0.0165 0.2263 108 109 0.3498 0.5597 109 110 0.6831 0.8930 110 111 0.1276 0.6708 Using dimension DIM_NUM = 3 Seed Seed Niederreiter In Out 0 1 0.3333 0.3333 0.3333 1 2 0.6667 0.6667 0.6667 2 3 0.1111 0.7778 0.4444 3 4 0.4444 0.1111 0.7778 4 5 0.7778 0.4444 0.1111 5 6 0.2222 0.5556 0.8889 6 7 0.5556 0.8889 0.2222 7 8 0.8889 0.2222 0.5556 8 9 0.0370 0.4815 0.5926 9 10 0.3704 0.8148 0.9259 10 11 0.7037 0.1481 0.2593 11 12 0.1481 0.2593 0.7037 ...................... 95 96 0.2634 0.2881 0.6091 96 97 0.5967 0.6214 0.9424 97 98 0.9300 0.9547 0.2757 98 99 0.0782 0.2140 0.3128 99 100 0.4115 0.5473 0.6461 100 101 0.7449 0.8807 0.9794 101 102 0.1893 0.9918 0.4239 102 103 0.5226 0.3251 0.7572 103 104 0.8560 0.6584 0.0905 104 105 0.3004 0.4362 0.8683 105 106 0.6337 0.7695 0.2016 106 107 0.9671 0.1029 0.5350 107 108 0.0165 0.2263 0.8066 108 109 0.3498 0.5597 0.1399 109 110 0.6831 0.8930 0.4733 110 111 0.1276 0.6708 0.2510 Using dimension DIM_NUM = 4 Seed Seed Niederreiter In Out 0 1 0.3333 0.3333 0.3333 0.4444 1 2 0.6667 0.6667 0.6667 0.8889 2 3 0.1111 0.7778 0.4444 0.5679 3 4 0.4444 0.1111 0.7778 0.6790 4 5 0.7778 0.4444 0.1111 0.1235 5 6 0.2222 0.5556 0.8889 0.8025 6 7 0.5556 0.8889 0.2222 0.2469 7 8 0.8889 0.2222 0.5556 0.3580 8 9 0.0370 0.4815 0.5926 0.9383 9 10 0.3704 0.8148 0.9259 0.0494 10 11 0.7037 0.1481 0.2593 0.4938 11 12 0.1481 0.2593 0.7037 0.1728 ...................... 95 96 0.2634 0.2881 0.6091 0.8697 96 97 0.5967 0.6214 0.9424 0.3141 97 98 0.9300 0.9547 0.2757 0.4252 98 99 0.0782 0.2140 0.3128 0.8944 99 100 0.4115 0.5473 0.6461 0.0055 100 101 0.7449 0.8807 0.9794 0.4499 101 102 0.1893 0.9918 0.4239 0.1289 102 103 0.5226 0.3251 0.7572 0.5734 103 104 0.8560 0.6584 0.0905 0.6845 104 105 0.3004 0.4362 0.8683 0.3635 105 106 0.6337 0.7695 0.2016 0.8080 106 107 0.9671 0.1029 0.5350 0.2524 107 108 0.0165 0.2263 0.8066 0.3278 108 109 0.3498 0.5597 0.1399 0.4390 109 110 0.6831 0.8930 0.4733 0.8834 110 111 0.1276 0.6708 0.2510 0.5254 TEST01 NIEDERREITER computes the next element of a Niederreiter quasirandom sequence using base BASE. In this test, we call NIEDERREITER repeatedly. Using base BASE = 13 Using dimension DIM_NUM = 2 Seed Seed Niederreiter In Out 0 1 0.0769 0.0769 1 2 0.1538 0.1538 2 3 0.2308 0.2308 3 4 0.3077 0.3077 4 5 0.3846 0.3846 5 6 0.4615 0.4615 6 7 0.5385 0.5385 7 8 0.6154 0.6154 8 9 0.6923 0.6923 9 10 0.7692 0.7692 10 11 0.8462 0.8462 11 12 0.9231 0.9231 ...................... 95 96 0.4260 0.8876 96 97 0.5030 0.9645 97 98 0.5799 0.0414 98 99 0.6568 0.1183 99 100 0.7337 0.1953 100 101 0.8107 0.2722 101 102 0.8876 0.3491 102 103 0.9645 0.4260 103 104 0.0473 0.4320 104 105 0.1243 0.5089 105 106 0.2012 0.5858 106 107 0.2781 0.6627 107 108 0.3550 0.7396 108 109 0.4320 0.8166 109 110 0.5089 0.8935 110 111 0.5858 0.9704 Using dimension DIM_NUM = 3 Seed Seed Niederreiter In Out 0 1 0.0769 0.0769 0.0769 1 2 0.1538 0.1538 0.1538 2 3 0.2308 0.2308 0.2308 3 4 0.3077 0.3077 0.3077 4 5 0.3846 0.3846 0.3846 5 6 0.4615 0.4615 0.4615 6 7 0.5385 0.5385 0.5385 7 8 0.6154 0.6154 0.6154 8 9 0.6923 0.6923 0.6923 9 10 0.7692 0.7692 0.7692 10 11 0.8462 0.8462 0.8462 11 12 0.9231 0.9231 0.9231 ...................... 95 96 0.4260 0.8876 0.3491 96 97 0.5030 0.9645 0.4260 97 98 0.5799 0.0414 0.5030 98 99 0.6568 0.1183 0.5799 99 100 0.7337 0.1953 0.6568 100 101 0.8107 0.2722 0.7337 101 102 0.8876 0.3491 0.8107 102 103 0.9645 0.4260 0.8876 103 104 0.0473 0.4320 0.8166 104 105 0.1243 0.5089 0.8935 105 106 0.2012 0.5858 0.9704 106 107 0.2781 0.6627 0.0473 107 108 0.3550 0.7396 0.1243 108 109 0.4320 0.8166 0.2012 109 110 0.5089 0.8935 0.2781 110 111 0.5858 0.9704 0.3550 Using dimension DIM_NUM = 4 Seed Seed Niederreiter In Out 0 1 0.0769 0.0769 0.0769 0.0769 1 2 0.1538 0.1538 0.1538 0.1538 2 3 0.2308 0.2308 0.2308 0.2308 3 4 0.3077 0.3077 0.3077 0.3077 4 5 0.3846 0.3846 0.3846 0.3846 5 6 0.4615 0.4615 0.4615 0.4615 6 7 0.5385 0.5385 0.5385 0.5385 7 8 0.6154 0.6154 0.6154 0.6154 8 9 0.6923 0.6923 0.6923 0.6923 9 10 0.7692 0.7692 0.7692 0.7692 10 11 0.8462 0.8462 0.8462 0.8462 11 12 0.9231 0.9231 0.9231 0.9231 ...................... 95 96 0.4260 0.8876 0.3491 0.8107 96 97 0.5030 0.9645 0.4260 0.8876 97 98 0.5799 0.0414 0.5030 0.9645 98 99 0.6568 0.1183 0.5799 0.0414 99 100 0.7337 0.1953 0.6568 0.1183 100 101 0.8107 0.2722 0.7337 0.1953 101 102 0.8876 0.3491 0.8107 0.2722 102 103 0.9645 0.4260 0.8876 0.3491 103 104 0.0473 0.4320 0.8166 0.2012 104 105 0.1243 0.5089 0.8935 0.2781 105 106 0.2012 0.5858 0.9704 0.3550 106 107 0.2781 0.6627 0.0473 0.4320 107 108 0.3550 0.7396 0.1243 0.5089 108 109 0.4320 0.8166 0.2012 0.5858 109 110 0.5089 0.8935 0.2781 0.6627 110 111 0.5858 0.9704 0.3550 0.7396 TEST02 NIEDERREITER computes the next element of a Niederreiter quasirandom sequence using base BASE. In this test, we demonstrate how the SEED can be manipulated to skip ahead in the sequence, or to come back to any part of the sequence. Using base BASE = 2 Using dimension DIM_NUM = 3 Seed Seed Niederreiter In Out 0 1 0.5000 0.5000 0.7500 1 2 0.2500 0.7500 0.5625 2 3 0.7500 0.2500 0.3125 3 4 0.1250 0.6250 0.4375 4 5 0.6250 0.1250 0.6875 5 6 0.3750 0.3750 0.8750 6 7 0.8750 0.8750 0.1250 7 8 0.0625 0.9375 0.9531 8 9 0.5625 0.4375 0.2031 9 10 0.3125 0.1875 0.3906 10 11 0.8125 0.6875 0.6406 Jump ahead by increasing SEED: Seed Seed Niederreiter In Out 100 101 0.1875 0.3125 0.5156 101 102 0.6875 0.8125 0.2656 102 103 0.4375 0.5625 0.0781 103 104 0.9375 0.0625 0.8281 104 105 0.0313 0.5313 0.7344 Jump back by decreasing SEED: Seed Seed Niederreiter In Out 3 4 0.5313 0.0313 0.4844 4 5 0.2813 0.2813 0.1719 5 6 0.7813 0.7813 0.9219 6 7 0.1563 0.1563 0.7969 7 8 0.6563 0.6563 0.0469 8 9 0.4063 0.9063 0.3594 9 10 0.9063 0.4063 0.6094 10 11 0.0938 0.4688 0.2813 11 12 0.5938 0.9688 0.5313 12 13 0.3438 0.7188 0.8438 13 14 0.8438 0.2188 0.0938 Jump ahead by increasing SEED: Seed Seed Niederreiter In Out 98 99 0.2188 0.8438 0.2188 99 100 0.7188 0.3438 0.9688 100 101 0.4688 0.0938 0.6563 101 102 0.9688 0.5938 0.4063 102 103 0.0156 0.7969 0.4258 TEST02 NIEDERREITER computes the next element of a Niederreiter quasirandom sequence using base BASE. In this test, we demonstrate how the SEED can be manipulated to skip ahead in the sequence, or to come back to any part of the sequence. Using base BASE = 3 Using dimension DIM_NUM = 3 Seed Seed Niederreiter In Out 0 1 0.3333 0.3333 0.3333 1 2 0.6667 0.6667 0.6667 2 3 0.1111 0.7778 0.4444 3 4 0.4444 0.1111 0.7778 4 5 0.7778 0.4444 0.1111 5 6 0.2222 0.5556 0.8889 6 7 0.5556 0.8889 0.2222 7 8 0.8889 0.2222 0.5556 8 9 0.0370 0.4815 0.5926 9 10 0.3704 0.8148 0.9259 10 11 0.7037 0.1481 0.2593 Jump ahead by increasing SEED: Seed Seed Niederreiter In Out 100 101 0.1481 0.2593 0.7037 101 102 0.4815 0.5926 0.0370 102 103 0.8148 0.9259 0.3704 103 104 0.2593 0.7037 0.1481 104 105 0.5926 0.0370 0.4815 Jump back by decreasing SEED: Seed Seed Niederreiter In Out 3 4 0.9259 0.3704 0.8148 4 5 0.0741 0.9630 0.8519 5 6 0.4074 0.2963 0.1852 6 7 0.7407 0.6296 0.5185 7 8 0.1852 0.4074 0.2963 8 9 0.5185 0.7407 0.6296 9 10 0.8519 0.0741 0.9630 10 11 0.2963 0.1852 0.4074 11 12 0.6296 0.5185 0.7407 12 13 0.9630 0.8519 0.0741 13 14 0.0123 0.6790 0.3457 Jump ahead by increasing SEED: Seed Seed Niederreiter In Out 98 99 0.3457 0.0123 0.6790 99 100 0.6790 0.3457 0.0123 100 101 0.1235 0.4568 0.7901 101 102 0.4568 0.7901 0.1235 102 103 0.7901 0.1235 0.4568 TEST03 NIEDERREITER computes the next element of a Niederreiter quasirandom sequence using base BASE. Using base BASE = 2 Using dimension DIM_NUM = 20 Seed Seed Niederreiter In Out 0 1 0.5000 0.5000 0.7500 0.8750 0.8750 0.9375 0.9375 0.9375 0.9688 0.9688 0.9688 0.9688 0.9688 0.9688 0.9844 0.9844 0.9844 0.9844 0.9844 0.9844 1 2 0.2500 0.7500 0.5625 0.7656 0.7656 0.8789 0.8789 0.8789 0.9385 0.9385 0.9385 0.9385 0.9385 0.9385 0.9690 0.9690 0.9690 0.9690 0.9690 0.9690 2 3 0.7500 0.2500 0.3125 0.1406 0.1406 0.0664 0.0664 0.0664 0.0322 0.0322 0.0322 0.0322 0.0322 0.0322 0.0159 0.0159 0.0159 0.0159 0.0159 0.0159 3 4 0.1250 0.6250 0.4375 0.5469 0.6719 0.7617 0.8242 0.8242 0.8779 0.8779 0.8779 0.9092 0.9092 0.9092 0.9382 0.9382 0.9382 0.9382 0.9539 0.9539 4 5 0.6250 0.1250 0.6875 0.4219 0.2969 0.1992 0.1367 0.1367 0.0967 0.0967 0.0967 0.0654 0.0654 0.0654 0.0476 0.0476 0.0476 0.0476 0.0320 0.0320 TEST03 NIEDERREITER computes the next element of a Niederreiter quasirandom sequence using base BASE. Using base BASE = 2 Using dimension DIM_NUM = 29 Seed Seed Niederreiter In Out 0 1 0.5000 0.5000 0.7500 0.8750 0.8750 0.9375 0.9375 0.9375 0.9688 0.9688 0.9688 0.9688 0.9688 0.9688 0.9844 0.9844 0.9844 0.9844 0.9844 0.9844 0.9844 0.9844 0.9844 0.9922 0.9922 0.9922 0.9922 0.9922 0.9922 1 2 0.2500 0.7500 0.5625 0.7656 0.7656 0.8789 0.8789 0.8789 0.9385 0.9385 0.9385 0.9385 0.9385 0.9385 0.9690 0.9690 0.9690 0.9690 0.9690 0.9690 0.9690 0.9690 0.9690 0.9844 0.9844 0.9844 0.9844 0.9844 0.9844 2 3 0.7500 0.2500 0.3125 0.1406 0.1406 0.0664 0.0664 0.0664 0.0322 0.0322 0.0322 0.0322 0.0322 0.0322 0.0159 0.0159 0.0159 0.0159 0.0159 0.0159 0.0159 0.0159 0.0159 0.0079 0.0079 0.0079 0.0079 0.0079 0.0079 3 4 0.1250 0.6250 0.4375 0.5469 0.6719 0.7617 0.8242 0.8242 0.8779 0.8779 0.8779 0.9092 0.9092 0.9092 0.9382 0.9382 0.9382 0.9382 0.9539 0.9539 0.9539 0.9539 0.9539 0.9689 0.9689 0.9689 0.9689 0.9689 0.9689 4 5 0.6250 0.1250 0.6875 0.4219 0.2969 0.1992 0.1367 0.1367 0.0967 0.0967 0.0967 0.0654 0.0654 0.0654 0.0476 0.0476 0.0476 0.0476 0.0320 0.0320 0.0320 0.0320 0.0320 0.0236 0.0236 0.0236 0.0236 0.0236 0.0236 NIEDERREITER_PRB Normal end of execution. June 7 2010 2:12:48.717 PM