Elsevier Science Home
Computer Physics Communications Program Library
Full text online from Science Direct
Programs in Physics & Physical Chemistry
CPC Home

[Licence| Download | New Version Template] aetf_v1_0.tar.gz(2879 Kbytes)
Manuscript Title: Optimized OpenCL implementation of the Elastodynamic Finite Integration Technique for viscoelastic media
Authors: M. Molero, Ursula Iturrarán-Viveros, S. Aparicio, M.G. Hernández
Program title: EFIT2D-PyOpenCL
Catalogue identifier: AETF_v1_0
Distribution format: tar.gz
Journal reference: Comput. Phys. Commun. 185(2014)2683
Programming language: Python.
Computer: Computers having GPU or Multicore CPU with OpenCL drivers.
Operating system: Multi-platform.
Has the code been vectorised or parallelized?: Yes.
RAM: 2 Gb
Keywords: EFIT, Kelvin-Voigt, GPUs, PyOpenCL, OpenCL.
Classification: 6.5.

External routines: numpy, scipy, matplotlib, glumpy, pyopencl

Nature of problem:
Development of parallel codes that are both scalable and portable for different processor architectures is a challenging task. To overcome this limitation we investigate the acceleration of the Elastodynamic Finite Integration Technique (EFIT) to model 2-D wave propagation in viscoelastic media by using modern parallel computing devices (PCDs), such as multi-core CPUs (central processing units) and GPUs (graphics processing units).

Solution method:
We choose the industry open standard Open Computing Language (OpenCL) and an open-source toolkit called PyOpenCL. The implementation is platform independent and can be used on AMD or NVIDIA GPUs as well as classical multi-core CPUs. The code is based on the Kelvin-Voigt mechanical model which has the gain of not requiring additional field variables. OpenCL performance can be in principle, improved once one can eliminate global memory access latency by using local memory. Our main contribution is the implementation of local memory and an analysis of performance of the local versus the global memory using eight different computing devices (including Kepler, one of the fastest and most efficient high performance computing technology) with various operating systems.

Restrictions:
Wave propagation simulation only in 2D Scenarios, OpenCL drivers needed.

Running time:
This code can process wave propagation simulations within a few minutes in a typical current computer with GPUs or multicore CPUs.