@brief Module containing subroutines for generating 1-D and 2-D Hammersley quasi-Monte Carlo sequences. @details The algorithm and code for generating the 1-D Hammersley sequence was developed by John Burkardt at the Florida State University. Visit "https://people.sc.fsu.edu/~jburkardt/f_src/hammersley/hammersley.html" for more information. The algorithm and code for generating the 2-D Hammersley sequence was developed by L. Carbajal at the Oak Ridge National Lab.
@brief For more info please visit "https://people.sc.fsu.edu/~jburkardt/f_src/hammersley/hammersley.html". ************80
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| integer(kind=4) | :: | n |
@brief Subroutine for generating a 2-D Hammersley sequence. @details This subroutine uses the algorithm for generating a 1-D Hammersley sequence. Each MPI process in KORC generates a (different) subset of pairs (X,Y) of a 2-D Hammersley sequence. The total number of pairs (X,Y) is NMPIS*N, where NMPIS is the number of MPI processes in the simulation and N is the number of particles followed by each MPI process. Each subset of pairs (X,Y) has N elements.
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| integer, | intent(in) | :: | ID | |||
| integer, | intent(in) | :: | NMPIS | |||
| real(kind=rp), | intent(inout), | DIMENSION(:) | :: | X | ||
| real(kind=rp), | intent(inout), | DIMENSION(:) | :: | Y |
@brief For more info please visit "https://people.sc.fsu.edu/~jburkardt/f_src/hammersley/hammersley.html". ************80
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| integer(kind=4) | :: | i | ||||
| integer(kind=4) | :: | m | ||||
| integer(kind=4) | :: | n | ||||
| real(kind=8) | :: | r(m) |
@brief For more info please visit "https://people.sc.fsu.edu/~jburkardt/f_src/hammersley/hammersley.html".
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| real(kind=8) | :: | r(m) | ||||
| integer(kind=4) | :: | m | ||||
| integer(kind=4) | :: | n | ||||
| integer(kind=4) | :: | i |
@brief For more info please visit "https://people.sc.fsu.edu/~jburkardt/f_src/hammersley/hammersley.html". ************80
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| integer(kind=4) | :: | i1 | ||||
| integer(kind=4) | :: | i2 | ||||
| integer(kind=4) | :: | m | ||||
| integer(kind=4) | :: | n | ||||
| real(kind=8) | :: | r(m,abs(i1-i2)+1) |
@brief For more info please visit "https://people.sc.fsu.edu/~jburkardt/f_src/hammersley/hammersley.html". ************80
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| integer(kind=4) | :: | m | ||||
| integer(kind=4) | :: | n | ||||
| real(kind=8) | :: | a(m,n) | ||||
| character(len=*) | :: | title |
@brief For more info please visit "https://people.sc.fsu.edu/~jburkardt/f_src/hammersley/hammersley.html". ************80
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| integer(kind=4) | :: | m | ||||
| integer(kind=4) | :: | n | ||||
| real(kind=8) | :: | a(m,n) | ||||
| integer(kind=4) | :: | ilo | ||||
| integer(kind=4) | :: | jlo | ||||
| integer(kind=4) | :: | ihi | ||||
| integer(kind=4) | :: | jhi | ||||
| character(len=*) | :: | title |
@brief For more info please visit "https://people.sc.fsu.edu/~jburkardt/f_src/hammersley/hammersley.html". ************80