11 August 2016 8:44:44.678 PM HALTON_PRB: FORTRAN90 version Test the HALTON library HALTON_TEST HALTON returns the I-th element of an M-dimensional Halton sequence. I HALTON(I) Use M = 1 0 0.00000000 1 0.50000000 2 0.25000000 3 0.75000000 4 0.12500000 5 0.62500000 6 0.37500000 7 0.87500000 8 0.06250000 9 0.56250000 10 0.31250000 Use M = 2 0 0.00000000 0.00000000 1 0.50000000 0.33333333 2 0.25000000 0.66666667 3 0.75000000 0.11111111 4 0.12500000 0.44444444 5 0.62500000 0.77777778 6 0.37500000 0.22222222 7 0.87500000 0.55555556 8 0.06250000 0.88888889 9 0.56250000 0.03703704 10 0.31250000 0.37037037 Use M = 3 0 0.00000000 0.00000000 0.00000000 1 0.50000000 0.33333333 0.20000000 2 0.25000000 0.66666667 0.40000000 3 0.75000000 0.11111111 0.60000000 4 0.12500000 0.44444444 0.80000000 5 0.62500000 0.77777778 0.04000000 6 0.37500000 0.22222222 0.24000000 7 0.87500000 0.55555556 0.44000000 8 0.06250000 0.88888889 0.64000000 9 0.56250000 0.03703704 0.84000000 10 0.31250000 0.37037037 0.08000000 HALTON_SEQUENCE_TEST HALTON_SEQUENCE returns the elements I1 through I2 of an M-dimensional Halton sequence. HALTON_SEQUENCE(0,10,1,R): R: Col 1 2 3 4 5 Row 1: 0. 0.500000 0.250000 0.750000 0.125000 Col 6 7 8 9 10 Row 1: 0.625000 0.375000 0.875000 0.625000E-01 0.562500 Col 11 Row 1: 0.312500 HALTON_SEQUENCE(0,10,2,R): R: Col 1 2 3 4 5 Row 1: 0. 0.500000 0.250000 0.750000 0.125000 2: 0. 0.333333 0.666667 0.111111 0.444444 Col 6 7 8 9 10 Row 1: 0.625000 0.375000 0.875000 0.625000E-01 0.562500 2: 0.777778 0.222222 0.555556 0.888889 0.370370E-01 Col 11 Row 1: 0.312500 2: 0.370370 HALTON_SEQUENCE(0,10,3,R): R: Col 1 2 3 4 5 Row 1: 0. 0.500000 0.250000 0.750000 0.125000 2: 0. 0.333333 0.666667 0.111111 0.444444 3: 0. 0.200000 0.400000 0.600000 0.800000 Col 6 7 8 9 10 Row 1: 0.625000 0.375000 0.875000 0.625000E-01 0.562500 2: 0.777778 0.222222 0.555556 0.888889 0.370370E-01 3: 0.400000E-01 0.240000 0.440000 0.640000 0.840000 Col 11 Row 1: 0.312500 2: 0.370370 3: 0.800000E-01 HALTON_SEQUENCE(10,0,3,R): R: Col 1 2 3 4 5 Row 1: 0.312500 0.562500 0.625000E-01 0.875000 0.375000 2: 0.370370 0.370370E-01 0.888889 0.555556 0.222222 3: 0.800000E-01 0.840000 0.640000 0.440000 0.240000 Col 6 7 8 9 10 Row 1: 0.625000 0.125000 0.750000 0.250000 0.500000 2: 0.777778 0.444444 0.111111 0.666667 0.333333 3: 0.400000E-01 0.800000 0.600000 0.400000 0.200000 Col 11 Row 1: 0. 2: 0. 3: 0. HALTON_INVERSE_TEST HALTON_INVERSE inverts an element of a Halton sequence. I R=HALTON(I,3) HALTON_INVERSE(R,3) 0 0.00000000 0.00000000 0.00000000 0 1 0.50000000 0.33333333 0.20000000 1 2 0.25000000 0.66666667 0.40000000 2 3 0.75000000 0.11111111 0.60000000 3 4 0.12500000 0.44444444 0.80000000 4 5 0.62500000 0.77777778 0.04000000 5 6 0.37500000 0.22222222 0.24000000 6 7 0.87500000 0.55555556 0.44000000 7 8 0.06250000 0.88888889 0.64000000 8 9 0.56250000 0.03703704 0.84000000 9 10 0.31250000 0.37037037 0.08000000 10 HALTON_BASE_TEST HALTON_BASE returns the I-th element of an M-dimensional Halton sequence, using user-specified bases. M = 3 B: 2 3 5 0 0.00000000 0.00000000 0.00000000 1 0.50000000 0.33333333 0.20000000 2 0.25000000 0.66666667 0.40000000 3 0.75000000 0.11111111 0.60000000 4 0.12500000 0.44444444 0.80000000 5 0.62500000 0.77777778 0.04000000 6 0.37500000 0.22222222 0.24000000 7 0.87500000 0.55555556 0.44000000 8 0.06250000 0.88888889 0.64000000 9 0.56250000 0.03703704 0.84000000 10 0.31250000 0.37037037 0.08000000 M = 3 B: 3 10 2 0 0.00000000 0.00000000 0.00000000 1 0.33333333 0.10000000 0.50000000 2 0.66666667 0.20000000 0.25000000 3 0.11111111 0.30000000 0.75000000 4 0.44444444 0.40000000 0.12500000 5 0.77777778 0.50000000 0.62500000 6 0.22222222 0.60000000 0.37500000 7 0.55555556 0.70000000 0.87500000 8 0.88888889 0.80000000 0.06250000 9 0.03703704 0.90000000 0.56250000 10 0.37037037 0.01000000 0.31250000 HALTON_PRB Normal end of execution. 11 August 2016 8:44:44.679 PM