Note this is written with the MPI API for use on a cluster. The external parameters are with init_params function in main.cpp. Each computing node will work on a block of size nx x ny, where ny = nx/(number of cores). ny must be even.
-
https://quantumtheory.physik.unibas.ch/people/bruder/Semesterprojekte2007/p6/index.html
-
Kosterlitz, J. M.; Thouless, D. J. (1973), "Ordering, metastability and phase transitions in two-dimensional systems", Journal of Physics C: Solid State Physics, 6 (7): 1181–1203