2016 Volume 23 Pages 1-11
The understanding and prediction of influences caused by seismic wave ranges from low to high frequency is important to mitigate the disaster of earthquakes by multilateral preparations. The structure of undergrounds causes the complex behavior of seismic wave propagation. The numerical simulation is an effective tool to predict the wave propagation. For metropolises on a large plain such as Kanto basin, however, the simulation requires enormous computing resources as well as techniques of high-performance computing. In this paper, we have developed the optimization method for the seismic simulation code on the Earth Simulator (ES2). We assume the finite difference method code ''k-fdm3d'' developed by KOZO KEKAKU ENGINEERING Inc.. Because the code adopt a stencil computation with the structured grid, we have optimized the memory-access and communication of the computation. For the reasons, the optimization method has applicability to simulation code on the other fields. The method enables a 570 million grids seismic simulation on ES2. The method also allows the simulation to achieve approximately 30％ to the peak performance of ES2. For the results, we have established the optimization method to calculate the high-accuracy prediction of seismic wave propagation applicable for a large plain.
