Hybrid Multicore/vectorisation technique applied to the elastic wave equation on a staggered grid

Sofya Titarenko, Mark Hildyard

Research output: Contribution to journalArticlepeer-review

10 Citations (Scopus)


In modern physics it has become common to find the solution of a problem by solving numerically a set of PDEs. Whether solving them on a finite difference grid or by a finite element approach, the main calculations are often applied to a stencil structure. In the last decade it has become usual to work with so called big data problems where calculations are very heavy and accelerators and modern architectures are widely used. Although CPU and GPU clusters are often used to solve such problems, parallelisation of any calculation ideally starts from a single processor optimisation. Unfortunately, it is impossible to vectorise a stencil structured loop with high level instructions. In this paper we suggest a new approach to rearranging the data structure which makes it possible to apply high level vectorisation instructions to a stencil loop and which results in significant acceleration. The suggested method allows further acceleration if shared memory APIs are used. We show the effectiveness of the method by applying it to an elastic wave propagation problem on a finite difference grid. We have chosen Intel architecture for the test problem and OpenMP (Open Multi-Processing) since they are extensively used in many applications.
Original languageEnglish
Pages (from-to)53-62
Number of pages10
JournalComputer Physics Communications
Early online date1 Mar 2017
Publication statusPublished - Jul 2017
Externally publishedYes


Dive into the research topics of 'Hybrid Multicore/vectorisation technique applied to the elastic wave equation on a staggered grid'. Together they form a unique fingerprint.

Cite this