Tue Oct 19 17:10:32 2021 sobol_test Test sobol(). sobol_test01 BITWISE_XOR returns the bitwise exclusive OR of two integers. I J BITXOR(I,J) 19 63 44 37 62 27 93 67 30 23 45 58 16 24 8 45 40 5 97 52 85 29 37 56 11 56 51 87 86 1 sobol_test02 i4_bit_hi1 returns the location of the high 1 bit. I i4_bit_hi1(I) 15 4 20 5 1 1 81 7 19 5 53 6 19 5 64 7 67 7 85 7 sobol_test03 i4_bit_lo0 returns the location of the low 0 bit. I i4_bit_lo0(I) 89 2 28 1 22 1 50 1 29 2 13 2 79 5 64 1 67 3 53 2 sobol_test04 i4_sobol returns the next element of a Sobol sequence. In this test, we call i4_sobol repeatedly. Using dimension DIM_NUM = 2 Seed Seed i4_sobol In Out 0 1 0.000000 0.000000 1 2 0.500000 0.500000 2 3 0.750000 0.250000 3 4 0.250000 0.750000 4 5 0.375000 0.375000 5 6 0.875000 0.875000 6 7 0.625000 0.125000 7 8 0.125000 0.625000 8 9 0.187500 0.312500 9 10 0.687500 0.812500 10 11 0.937500 0.062500 11 12 0.437500 0.562500 ...................... 95 96 0.054688 0.929688 96 97 0.039062 0.132812 97 98 0.539062 0.632812 98 99 0.789062 0.382812 99 100 0.289062 0.882812 100 101 0.414062 0.257812 101 102 0.914062 0.757812 102 103 0.664062 0.007812 103 104 0.164062 0.507812 104 105 0.226562 0.445312 105 106 0.726562 0.945312 106 107 0.976562 0.195312 107 108 0.476562 0.695312 108 109 0.351562 0.070312 109 110 0.851562 0.570312 110 111 0.601562 0.320312 Using dimension DIM_NUM = 3 Seed Seed i4_sobol In Out 0 1 0.000000 0.000000 0.000000 1 2 0.500000 0.500000 0.500000 2 3 0.750000 0.250000 0.750000 3 4 0.250000 0.750000 0.250000 4 5 0.375000 0.375000 0.625000 5 6 0.875000 0.875000 0.125000 6 7 0.625000 0.125000 0.375000 7 8 0.125000 0.625000 0.875000 8 9 0.187500 0.312500 0.312500 9 10 0.687500 0.812500 0.812500 10 11 0.937500 0.062500 0.562500 11 12 0.437500 0.562500 0.062500 ...................... 95 96 0.054688 0.929688 0.101562 96 97 0.039062 0.132812 0.929688 97 98 0.539062 0.632812 0.429688 98 99 0.789062 0.382812 0.179688 99 100 0.289062 0.882812 0.679688 100 101 0.414062 0.257812 0.304688 101 102 0.914062 0.757812 0.804688 102 103 0.664062 0.007812 0.554688 103 104 0.164062 0.507812 0.054688 104 105 0.226562 0.445312 0.742188 105 106 0.726562 0.945312 0.242188 106 107 0.976562 0.195312 0.492188 107 108 0.476562 0.695312 0.992188 108 109 0.351562 0.070312 0.117188 109 110 0.851562 0.570312 0.617188 110 111 0.601562 0.320312 0.867188 Using dimension DIM_NUM = 4 Seed Seed i4_sobol In Out 0 1 0.000000 0.000000 0.000000 0.000000 1 2 0.500000 0.500000 0.500000 0.500000 2 3 0.750000 0.250000 0.750000 0.250000 3 4 0.250000 0.750000 0.250000 0.750000 4 5 0.375000 0.375000 0.625000 0.125000 5 6 0.875000 0.875000 0.125000 0.625000 6 7 0.625000 0.125000 0.375000 0.375000 7 8 0.125000 0.625000 0.875000 0.875000 8 9 0.187500 0.312500 0.312500 0.687500 9 10 0.687500 0.812500 0.812500 0.187500 10 11 0.937500 0.062500 0.562500 0.937500 11 12 0.437500 0.562500 0.062500 0.437500 ...................... 95 96 0.054688 0.929688 0.101562 0.960938 96 97 0.039062 0.132812 0.929688 0.351562 97 98 0.539062 0.632812 0.429688 0.851562 98 99 0.789062 0.382812 0.179688 0.101562 99 100 0.289062 0.882812 0.679688 0.601562 100 101 0.414062 0.257812 0.304688 0.476562 101 102 0.914062 0.757812 0.804688 0.976562 102 103 0.664062 0.007812 0.554688 0.226562 103 104 0.164062 0.507812 0.054688 0.726562 104 105 0.226562 0.445312 0.742188 0.914062 105 106 0.726562 0.945312 0.242188 0.414062 106 107 0.976562 0.195312 0.492188 0.664062 107 108 0.476562 0.695312 0.992188 0.164062 108 109 0.351562 0.070312 0.117188 0.789062 109 110 0.851562 0.570312 0.617188 0.289062 110 111 0.601562 0.320312 0.867188 0.539062 sobol_test05 i4_sobol computes the next element of a Sobol sequence. 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 i4_sobol In Out 0 1 0.000000 0.000000 0.000000 1 2 0.500000 0.500000 0.500000 2 3 0.750000 0.250000 0.750000 3 4 0.250000 0.750000 0.250000 4 5 0.375000 0.375000 0.625000 5 6 0.875000 0.875000 0.125000 6 7 0.625000 0.125000 0.375000 7 8 0.125000 0.625000 0.875000 8 9 0.187500 0.312500 0.312500 9 10 0.687500 0.812500 0.812500 10 11 0.937500 0.062500 0.562500 Jump ahead by increasing SEED: Seed Seed i4_sobol In Out 100 101 0.414062 0.257812 0.304688 101 102 0.914062 0.757812 0.804688 102 103 0.664062 0.007812 0.554688 103 104 0.164062 0.507812 0.054688 104 105 0.226562 0.445312 0.742188 Jump back by decreasing SEED: Seed Seed i4_sobol In Out 3 4 0.250000 0.750000 0.250000 4 5 0.375000 0.375000 0.625000 5 6 0.875000 0.875000 0.125000 6 7 0.625000 0.125000 0.375000 7 8 0.125000 0.625000 0.875000 8 9 0.187500 0.312500 0.312500 9 10 0.687500 0.812500 0.812500 10 11 0.937500 0.062500 0.562500 11 12 0.437500 0.562500 0.062500 12 13 0.312500 0.187500 0.937500 13 14 0.812500 0.687500 0.437500 Jump back by decreasing SEED: Seed Seed i4_sobol In Out 98 99 0.789062 0.382812 0.179688 99 100 0.289062 0.882812 0.679688 100 101 0.414062 0.257812 0.304688 101 102 0.914062 0.757812 0.804688 102 103 0.664062 0.007812 0.554688 sobol_test Normal end of execution. Tue Oct 19 17:10:32 2021