455.seismic
SPECaccel Benchmark Description

Benchmark Name

455.seismic

Benchmark Author

Dimitri Komatitsch and Roland Martin from University of Pau, France.

Initial port to OpenMP by Alexander Grund

Benchmark Program General Category

Seismic Wave Modeling

Benchmark Description

SEISMIC_CPML is a set of eleven open-source Fortran90 programs to solve the two-dimensional or three-dimensional isotropic or anisotropic elastic, viscoelastic or poroelastic wave equation using a finite-difference method with Convolutional or Auxiliary Perfectly Matched Layer (C-PML or ADE-PML) conditions. 455.seismic is an OpenMP enabled port of the "seismic_PML_Collino_3D_isotropic" solver, a 3D classical split PML program for an isotropic medium using a second-order finite-difference spatial operator, for comparison.

Input Description

The program takes in the X, Y, and Z size of the area to model, the number of receivers to use, the number of time steps to model, and finally how often (in timesteps) to output the image data.

Output Description

The "energy" data file contiain the accumlated total kinetic and potential energy after each timestep. The program also outputs simple image files showing the progress of the seismic wave at each time step.

Programming Language

Fortran

Known portability issues

None.

References

• More information about SEISMIC_CPML can be found at:https://github.com/geodynamics/seismic_cpml
• The steps involved in porting SEISMIC_CPML 3D Isotropic model to accelerators are show in the article 5x in 5 Hours: Porting a 3D Elastic Wave Simulator to GPUs Using the PGI Accelerator Model by Mathew Colgrove http://www.pgroup.com/lit/articles/insider/v4n1a3.htm