MCSiMPI
While working as a post-graduate researcher in computational physics, I invented a new, optimally efficient parallisation scheme for Monte Carlo simulations of systems with short-range interactions. Unlike other schemes, it maintains the global Markov chain, making it possible to directly compare runs and results between serial and parallel computations and ensure the integrity of the simulation is not compromised when run at scale. It is written in C++ and uses MPI, and is designed so that the exact same algorithm is executed whether the simulation runs in two or three spatial dimensions.
The following animation shows a brief simulation of a ‘sticky hard disk fluid’, and shows how after the fluid has been mixed up, the disks rapidly aggregate to form a space filling structure. The simulation has a gravitational force applied, which leads to the structure sinking slightly before stablising.
The picture on the below shows a visualisation of the same simulation, but running in three dimensions.
