  We use an efficient method that eases the daunting task of simulating
dynamics in spin systems with long-range interaction. Our Monte Carlo
simulations of the long-range Ising model for the nonequilibrium phase ordering
dynamics in two spatial dimensions erform significantly faster than the
standard Metropolis approach and considerably more efficiently than the kinetic
Monte Carlo method. Importantly, this enables us to establish agreement with
the theoretical prediction for the time dependence of domain growth, in
contrast to previous numerical studies. This method can easily be generalized
to applications in other systems.
