June 29 2009 11:23:23.834 AM NIEDERREITER2_PRB FORTRAN90 version Test the NIEDERREITER2 routines. TEST01 NIEDERREITER2 computes the next element of a Niederreiter quasirandom sequence using base 2. In this test, we call NIEDERREITER2 repeatedly. Using dimension DIM_NUM = 2 Seed Seed Niederreiter2 In Out 0 1 0.0000 0.0000 1 2 0.5000 0.5000 2 3 0.7500 0.2500 3 4 0.2500 0.7500 4 5 0.3750 0.3750 5 6 0.8750 0.8750 6 7 0.6250 0.1250 7 8 0.1250 0.6250 8 9 0.1875 0.3125 9 10 0.6875 0.8125 10 11 0.9375 0.0625 11 12 0.4375 0.5625 ...................... 95 96 0.0547 0.9297 96 97 0.0391 0.1328 97 98 0.5391 0.6328 98 99 0.7891 0.3828 99 100 0.2891 0.8828 100 101 0.4141 0.2578 101 102 0.9141 0.7578 102 103 0.6641 0.0078 103 104 0.1641 0.5078 104 105 0.2266 0.4453 105 106 0.7266 0.9453 106 107 0.9766 0.1953 107 108 0.4766 0.6953 108 109 0.3516 0.0703 109 110 0.8516 0.5703 110 111 0.6016 0.3203 Using dimension DIM_NUM = 3 Seed Seed Niederreiter2 In Out 0 1 0.0000 0.0000 0.0000 1 2 0.5000 0.5000 0.7500 2 3 0.7500 0.2500 0.3125 3 4 0.2500 0.7500 0.5625 4 5 0.3750 0.3750 0.8750 5 6 0.8750 0.8750 0.1250 6 7 0.6250 0.1250 0.6875 7 8 0.1250 0.6250 0.4375 8 9 0.1875 0.3125 0.5156 9 10 0.6875 0.8125 0.2656 10 11 0.9375 0.0625 0.8281 11 12 0.4375 0.5625 0.0781 ...................... 95 96 0.0547 0.9297 0.1016 96 97 0.0391 0.1328 0.4648 97 98 0.5391 0.6328 0.7148 98 99 0.7891 0.3828 0.1523 99 100 0.2891 0.8828 0.9023 100 101 0.4141 0.2578 0.5898 101 102 0.9141 0.7578 0.3398 102 103 0.6641 0.0078 0.7773 103 104 0.1641 0.5078 0.0273 104 105 0.2266 0.4453 0.9492 105 106 0.7266 0.9453 0.1992 106 107 0.9766 0.1953 0.6367 107 108 0.4766 0.6953 0.3867 108 109 0.3516 0.0703 0.0742 109 110 0.8516 0.5703 0.8242 110 111 0.6016 0.3203 0.2617 Using dimension DIM_NUM = 4 Seed Seed Niederreiter2 In Out 0 1 0.0000 0.0000 0.0000 0.0000 1 2 0.5000 0.5000 0.7500 0.8750 2 3 0.7500 0.2500 0.3125 0.1406 3 4 0.2500 0.7500 0.5625 0.7656 4 5 0.3750 0.3750 0.8750 0.2813 5 6 0.8750 0.8750 0.1250 0.6563 6 7 0.6250 0.1250 0.6875 0.4219 7 8 0.1250 0.6250 0.4375 0.5469 8 9 0.1875 0.3125 0.5156 0.6875 9 10 0.6875 0.8125 0.2656 0.3125 10 11 0.9375 0.0625 0.8281 0.5781 11 12 0.4375 0.5625 0.0781 0.4531 ...................... 95 96 0.0547 0.9297 0.1016 0.5098 96 97 0.0391 0.1328 0.4648 0.2148 97 98 0.5391 0.6328 0.7148 0.8398 98 99 0.7891 0.3828 0.1523 0.0742 99 100 0.2891 0.8828 0.9023 0.9492 100 101 0.4141 0.2578 0.5898 0.4961 101 102 0.9141 0.7578 0.3398 0.6211 102 103 0.6641 0.0078 0.7773 0.3555 103 104 0.1641 0.5078 0.0273 0.7305 104 105 0.2266 0.4453 0.9492 0.5273 105 106 0.7266 0.9453 0.1992 0.4023 106 107 0.9766 0.1953 0.6367 0.6367 107 108 0.4766 0.6953 0.3867 0.2617 108 109 0.3516 0.0703 0.0742 0.8086 109 110 0.8516 0.5703 0.8242 0.1836 110 111 0.6016 0.3203 0.2617 0.9180 TEST02 NIEDERREITER2 computes the next element of a Niederreiter quasirandom sequence using base 2. 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 dimension DIM_NUM = 3 Seed Seed Niederreiter2 In Out 0 1 0.0000 0.0000 0.0000 1 2 0.5000 0.5000 0.7500 2 3 0.7500 0.2500 0.3125 3 4 0.2500 0.7500 0.5625 4 5 0.3750 0.3750 0.8750 5 6 0.8750 0.8750 0.1250 6 7 0.6250 0.1250 0.6875 7 8 0.1250 0.6250 0.4375 8 9 0.1875 0.3125 0.5156 9 10 0.6875 0.8125 0.2656 10 11 0.9375 0.0625 0.8281 Jump ahead by increasing SEED: Seed Seed Niederreiter2 In Out 100 101 0.4141 0.2578 0.5898 101 102 0.9141 0.7578 0.3398 102 103 0.6641 0.0078 0.7773 103 104 0.1641 0.5078 0.0273 104 105 0.2266 0.4453 0.9492 Jump back by decreasing SEED: Seed Seed Niederreiter2 In Out 3 4 0.2500 0.7500 0.5625 4 5 0.3750 0.3750 0.8750 5 6 0.8750 0.8750 0.1250 6 7 0.6250 0.1250 0.6875 7 8 0.1250 0.6250 0.4375 8 9 0.1875 0.3125 0.5156 9 10 0.6875 0.8125 0.2656 10 11 0.9375 0.0625 0.8281 11 12 0.4375 0.5625 0.0781 12 13 0.3125 0.1875 0.3906 13 14 0.8125 0.6875 0.6406 Jump ahead by increasing SEED: Seed Seed Niederreiter2 In Out 98 99 0.7891 0.3828 0.1523 99 100 0.2891 0.8828 0.9023 100 101 0.4141 0.2578 0.5898 101 102 0.9141 0.7578 0.3398 102 103 0.6641 0.0078 0.7773 NIEDERREITER2_PRB Normal end of execution. June 29 2009 11:23:23.841 AM