### Abstract

Original language | English |
---|---|

Pages (from-to) | 53-62 |

Number of pages | 10 |

Journal | Computer Physics Communications |

Volume | 216 |

Early online date | 1 Mar 2017 |

DOIs | |

Publication status | Published - Jul 2017 |

Externally published | Yes |

### Fingerprint

### Cite this

}

*Computer Physics Communications*, vol. 216, pp. 53-62. https://doi.org/10.1016/j.cpc.2017.02.022

**Hybrid Multicore/vectorisation technique applied to the elastic wave equation on a staggered grid.** / Titarenko, Sofya; Hildyard, Mark.

Research output: Contribution to journal › Article

TY - JOUR

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

AU - Titarenko, Sofya

AU - Hildyard, Mark

PY - 2017/7

Y1 - 2017/7

N2 - 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.

AB - 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.

KW - OpenMP

KW - Vectorisation

KW - Multicore

KW - Elastic waves

KW - Staggered grid

UR - https://www.scopus.com/record/display.uri?eid=2-s2.0-85015364959&origin=resultslist&sort=plf-f&src=s&st1=Hybrid+Multicore%2fvectorisation+technique+applied+to+the+elastic+wave+equation+on+a+staggered+grid&st2=&sid=db53b3c1e0cdb317f659a38f1390a067&sot=b&sdt=b&sl=112&s=TITLE-ABS-KEY%28Hybrid+Multicore%2fvectorisation+technique+applied+to+the+elastic+wave+equation+on+a+staggered+grid%29&relpos=0&citeCnt=4&searchTerm=

U2 - 10.1016/j.cpc.2017.02.022

DO - 10.1016/j.cpc.2017.02.022

M3 - Article

VL - 216

SP - 53

EP - 62

JO - Computer Physics Communications

JF - Computer Physics Communications

SN - 0010-4655

ER -