Fast algorithm for the simulation of Ising models

Using a new microcanonical algorithm efficiently vectorized on a Cray XMP, we reach a simulation speed of 1.5 nsec per update of one spin, three times faster than the best previous method known to us. Data for the nonlinear relaxation with conserved energy are presented for the two-dimensional Ising model.