During the initialization we pre-computed a maximum radius for each
object from its center of mass. We use this radius and the object's
initial and final position for this time step to compute a worst case
swept bounding volume. We also consider the apex of the object's
parabolic trajectory, should this occur during this time step. This
final volume is a discrete rectangular volume with faces parallel to
the global coordinate planes. An example of such a volume is show
above.
The computation of the swept volumes occur in parallel for the objects
on all processors. Since each processor maintains a distributed
subset of all of the objects, all computations are local for
determining the swept volumes.