20 August 2016 9:43:48.149 AM HAMMERSLEY_PRB: FORTRAN90 version Test the HAMMERSLEY library HAMMERSLEY_TEST HAMMERSLEY returns the I-th element of an M-dimensional Hammersley sequence. I HAMMERSLEY(I) Use M = 1 N = 16 0 0.00000000 1 0.06250000 2 0.12500000 3 0.18750000 4 0.25000000 5 0.31250000 6 0.37500000 7 0.43750000 8 0.50000000 9 0.56250000 10 0.62500000 Use M = 2 N = 16 0 0.00000000 0.00000000 1 0.06250000 0.50000000 2 0.12500000 0.25000000 3 0.18750000 0.75000000 4 0.25000000 0.12500000 5 0.31250000 0.62500000 6 0.37500000 0.37500000 7 0.43750000 0.87500000 8 0.50000000 0.06250000 9 0.56250000 0.56250000 10 0.62500000 0.31250000 Use M = 3 N = 16 0 0.00000000 0.00000000 0.00000000 1 0.06250000 0.50000000 0.33333333 2 0.12500000 0.25000000 0.66666667 3 0.18750000 0.75000000 0.11111111 4 0.25000000 0.12500000 0.44444444 5 0.31250000 0.62500000 0.77777778 6 0.37500000 0.37500000 0.22222222 7 0.43750000 0.87500000 0.55555556 8 0.50000000 0.06250000 0.88888889 9 0.56250000 0.56250000 0.03703704 10 0.62500000 0.31250000 0.37037037 HAMMERSLEY_SEQUENCE_TEST HAMMERSLEY_SEQUENCE returns the elements I1 through I2 of an M-dimensional Hammersley sequence. HAMMERSLEY_SEQUENCE(0,10,1N,R): R: Col 1 2 3 4 5 Row 1: 0. 0.625000E-01 0.125000 0.187500 0.250000 Col 6 7 8 9 10 Row 1: 0.312500 0.375000 0.437500 0.500000 0.562500 Col 11 Row 1: 0.625000 HAMMERSLEY_SEQUENCE(0,10,2N,R): R: Col 1 2 3 4 5 Row 1: 0. 0.625000E-01 0.125000 0.187500 0.250000 2: 0. 0.500000 0.250000 0.750000 0.125000 Col 6 7 8 9 10 Row 1: 0.312500 0.375000 0.437500 0.500000 0.562500 2: 0.625000 0.375000 0.875000 0.625000E-01 0.562500 Col 11 Row 1: 0.625000 2: 0.312500 HAMMERSLEY_SEQUENCE(0,10,3N,R): R: Col 1 2 3 4 5 Row 1: 0. 0.625000E-01 0.125000 0.187500 0.250000 2: 0. 0.500000 0.250000 0.750000 0.125000 3: 0. 0.333333 0.666667 0.111111 0.444444 Col 6 7 8 9 10 Row 1: 0.312500 0.375000 0.437500 0.500000 0.562500 2: 0.625000 0.375000 0.875000 0.625000E-01 0.562500 3: 0.777778 0.222222 0.555556 0.888889 0.370370E-01 Col 11 Row 1: 0.625000 2: 0.312500 3: 0.370370 HAMMERSLEY_SEQUENCE(10,0,3,N,R): R: Col 1 2 3 4 5 Row 1: 0.625000 0.562500 0.500000 0.437500 0.375000 2: 0.312500 0.562500 0.625000E-01 0.875000 0.375000 3: 0.370370 0.370370E-01 0.888889 0.555556 0.222222 Col 6 7 8 9 10 Row 1: 0.312500 0.250000 0.187500 0.125000 0.625000E-01 2: 0.625000 0.125000 0.750000 0.250000 0.500000 3: 0.777778 0.444444 0.111111 0.666667 0.333333 Col 11 Row 1: 0. 2: 0. 3: 0. HAMMERSLEY_INVERSE_TEST HAMMERSLEY_INVERSE inverts an element of a Hammersley sequence. I R=HAMMERSLEY(I,3) HAMMERSLEY_INVERSE(R,3) 0 0.00000000 0.00000000 0.00000000 0 1 0.06250000 0.50000000 0.33333333 1 2 0.12500000 0.25000000 0.66666667 2 3 0.18750000 0.75000000 0.11111111 3 4 0.25000000 0.12500000 0.44444444 4 5 0.31250000 0.62500000 0.77777778 5 6 0.37500000 0.37500000 0.22222222 6 7 0.43750000 0.87500000 0.55555556 7 8 0.50000000 0.06250000 0.88888889 8 9 0.56250000 0.56250000 0.03703704 9 10 0.62500000 0.31250000 0.37037037 10 HAMMERSLEY_PRB Normal end of execution. 20 August 2016 9:43:48.150 AM