Parallel and vectorized implementation of analytic evaluation of boundary integral operators. (November 2018)

• Record Type: Journal Article
• Authors: Zapletal, Jan; Of, Günther; Merta, Michal
• Is Part Of: Engineering analysis with boundary elements. Volume 96(2018)
• Page Start: 194
• Page End: 208
• Publication Date: 2018-11
• Subjects: Boundary element method -- Quadrature -- SIMD -- Vectorization -- Intel Xeon Phi -- Many-core; 65N38 -- 65Y05 -- 68W10; Boundary element methods -- Periodicals; Engineering mathematics -- Periodicals; Équations intégrales de frontière, Méthodes des -- Périodiques; Mathématiques de l'ingénieur -- Périodiques; Boundary element methods; Engineering mathematics; Periodicals; 620.00151
• DOI: 10.1016/j.enganabound.2018.08.015
• Languages: English
• ISSNs: 0955-7997

Abstract

Abstract: In this paper, we describe an efficient analytic evaluation of boundary integral operators. Firstly, we concentrate on a novel approach based on the simultaneous evaluation of all three linear shape functions defined on a boundary triangle. This results in a speedup of 2.35–3.15 times compared to the old approach of separate evaluations. In the second part we comment on the OpenMP parallelized and vectorized implementation of the suggested formulae. The employed code optimizations include techniques such as data alignment and padding, array-of-structures to structure-of-arrays data transformation, or unit-strided memory accesses. The presented scalability results, with respect both to the number of threads employed and the width of the SIMD register obtained on an Intel ® Xeon™ processor and two generations of the Intel ® Xeon Phi™ family (co)processors, validate the performed optimizations and show that vectorization needs to be an inherent part of modern scientific codes.