Working f90 massively parallel source codes
for the CRAY T3E.
The codes heavily employs the CRAY-specific SHMEM routines.
All codes use Poisson asynchrony for the discrete-events (spin-flips)
and a conservative synchronization scheme [(Lubachevsky, J. Comput.
Phys.
75, 103 (1988)] for the parallel discrete-event simulations.
No warranties, no regrets. Use them at your own risk.
Appropriate acknowledgement of the source and the story behind it is
appreciated:
"Parallelization
of a Dynamic Monte Carlo Algorithm: A Partially Rejection-Free Conservative
Approach",
G. Korniss, M.A. Novotny, and P.A. Rikvold, J. Comput. Phys.153,
488 (1999), (paper on the implementation).
"From Massively
Parallel Algorithms and Fluctuating Time Horizons to Non-equilibrium Surface
Growth",
G. Korniss, Z. Toroczkai, M.A. Novotny, and P.A. Rikvold, Phys.
Rev. Lett. 84, 1351 (2000), (paper on scalability).
Continuous-time "Metropolis" algorithm for the two-dimensional kinetic
Ising model after an instateneous field reversal. The code can be
used to study hysteresis and metastability in large spatially-extended
bistable systems.
par_ctmp.f90
Partially rejection-free (n-fold-way-like) version of the same
algorithm as above.
par_ctnf.f90
Continuous-time "Metropolis" algorithm for the two-dimensional kinetic
Ising model in a periodic square-wave external field. The code can
be used to study the dynamic phase transition in large spatially-extended
bistable systems.
par_dpt_ctmp_J.f90
The "KISS" random number generator by G. Marsaglia and A. Zaman (1993).
kiss.f90