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.

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