355.seismic
SPEC ACCEL Benchmark Description

Benchmark Name

355.seismic

Benchmark Author

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

Ported to OpenACC by Mathew Colgrove

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. 355.seismic is an OpenACC 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

Threading Model

OpenACC

Known portability issues

None.

References

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

Last updated: $Date: 2015-03-02 15:15:22 -0500 (Mon, 02 Mar 2015) $