We report on a lattice based algorithm, completely vectorized for molecular
dynamics simulations. Its algorithmic complexity is of the order O(N), where
N is the number of particles. The algorithm works very effectively when the
particles have short range interaction, but it is applicable to each kind of
interaction. The code was tested on a Cray ymp el in a simulation of flowing
granular material.Comment: 9 pages, 6 figures, Late