Monte Carlo technique for very large ising models

Rebbi's multispin coding technique is improved and applied to the kinetic Ising model with size 600*600*600. We give the central part of our computer program (for a CDC Cyber 76), which will be helpful also in a simulation of smaller systems, and describe the other tricks necessary to go to large lattices. The magnetizationM atT=1.4*T_c is found to decay asymptotically as exp(-t/2.90) ift is measured in Monte Carlo steps per spin, and M(t = 0) = 1 initially.