  We present a computational algorithm for computing short range forces between
particles. The algorithm has two distinguishing features. First, it is
optimized for multi-processor computers, and will use as many processors as are
available. Second, it is optimized for inhomogeneous, dynamic particle
distributions; for any distribution the computational load is distributed
evenly to all processors, and the communication time is less than 15\% of the
total run time.
  In this talk we present our new algorithm. We developed the program for a
grand-challenge problem in cosmology, simulation of the formation of
large-scale structure in the universe. This simulation, run on the Thinking
Machines Corporation CM-5, uses the particle-particle/particle-mesh (PPPM)
\cite{hock_east} algorithm. The particle-particle phase is computed using the
algorithm we describe in this paper. We discuss this and other applications.
