A system, method and computer program product for seismic imaging implements a seismic imaging algorithm utilizing Reverse Time Migration technique requiring large communication bandwidth and low latency to convert a parallel problem into one solved using massive domain partitioning. Several aspects of the imaging problem, including very regular and local communication patterns, balanced compute and communication requirements, scratch data handling and multiple-pass approaches. The partitioning of the velocity model into processing blocks allows each sub-problem to fit in a local cache, increasing locality and bandwidth and reducing latency.