Subscribe

Phoenix, Arizona

Nov. 15, 2003 to Nov. 21, 2003

ISBN: 1-58113-695-1

pp: 4

Dimitri Komatitsch , Caltech, Pasadena, California

Seiji Tsuboi , IFREE, JAMSTEC, Yokohama, Japan

Chen Ji , Caltech, Pasadena, California

Jeroen Tromp , Caltech, Pasadena, California

DOI Bookmark: http://doi.ieeecomputersociety.org/10.1109/SC.2003.10023

ABSTRACT

We use 1944 processors of the Earth Simulator to model seismic wave propagation resulting from large earthquakes. Simulations are conducted based upon the spectral-element method, a high-degree finite-element technique with an exactly diagonal mass matrix. We use a very large mesh with 5.5 billion grid points (14.6 billion degrees of freedom). We include the full complexity of the Earth, i.e., a three-dimensional wave-speed and density structure, a 3-D crustal model, ellipticity as well as topography and bathymetry. A total of 2.5 terabytes of memory is needed. Our implementation is purely based upon MPI, with loop vectorization on each processor. We obtain an excellent vectorization ratio of 99.3%, and we reach a performance of 5 teraflops (30% of the peak performance) on 38% of the machine. The very high resolution of the mesh allows us to perform fully three-dimensional calculations at seismic periods as low as 5 seconds.

CITATION

Dimitri Komatitsch,
Seiji Tsuboi,
Chen Ji,
Jeroen Tromp,
"A 14.6 billion degrees of freedom, 5 teraflops, 2.5 terabyte earthquake simulation on the Earth Simulator",

*SC*, 2003, SC Conference, SC Conference 2003, pp. 4, doi:10.1109/SC.2003.10023