Affiliation:
1. Barcelona Supercomputing Center
2. Repsol USA
Abstract
Finite Difference (FD) is a widely used method to solve Partial Differential Equations (PDE). PDEs are the core of many simulations in different scientific fields, such as geophysics, astrophysics, etc. The typical FD solver performs stencil computations for the entire computational domain, thus solving the differential operators. In general terms, the stencil computation consists of a weighted accumulation of the contribution of neighbor points along the cartesian axis. Therefore, optimizing stencil computations is crucial in reducing the application execution time.
Stencil computation performance is bounded by two main factors: the memory access pattern and the inefficient reuse of the accessed data. We propose a novel algorithm, named
Semi-stencil
, that tackles these two problems. The main idea behind this algorithm is to change the way in which the stencil computation progresses within the computational domain. Instead of accessing all required neighbors and adding all their contributions at once, the Semi-stencil algorithm divides the computation into several updates. Then, each update gathers half of the axis neighbors, partially computing at the same time the stencil in a set of closely located points. As Semi-stencil progresses through the domain, the stencil computations are completed on precomputed points. This computation strategy improves the memory access pattern and efficiently reuses the accessed data.
Our initial target architecture was the Cell/B.E., where the Semi-stencil in a SPE was 44% faster than the naive stencil implementation. Since then, we have continued our research on emerging multicore architectures in order to assess and extend this work on homogeneous architectures. The experiments presented combine the Semi-stencil strategy with space- and time-blocking algorithms used in hierarchical memory architectures. Two x86 (Intel Nehalem and AMD Opteron) and two POWER (IBM POWER6 and IBM BG/P) platforms are used as testbeds, where the best improvements for a 25-point stencil range from 1.27 to 1.76× faster. The results show that this novel strategy is a feasible optimization method which may be integrated into auto-tuning frameworks. Also, since all current architectures are multicore based, we have introduced a brief section where scalability results on IBM POWER7-, Intel Xeon-, and MIC-based systems are presented. In a nutshell, the algorithm scales as well as or better than other stencil techniques. For instance, the scalability of Semi-stencil on MIC for a certain testcase reached 93.8 × over 244 threads.
Funder
Ministerio de Ciencia e Innovación
Seventh Framework Programme
Partnership for Advanced Computing in Europe AISBL
Publisher
Association for Computing Machinery (ACM)
Subject
Applied Mathematics,Software
Reference34 articles.
1. Software pipelining
2. Efficient Formalism for Large-ScaleAb InitioMolecular Dynamics based on Time-Dependent Density Functional Theory
3. ANAG. 2012. Chombo software package for amr applications. Applied Numerical Algorithms Group (ANAG) Lawrence Berkeley National Laboratory Berkeley CA. http://seesar.lbl.gov/anag/software.html. ANAG. 2012. Chombo software package for amr applications. Applied Numerical Algorithms Group (ANAG) Lawrence Berkeley National Laboratory Berkeley CA. http://seesar.lbl.gov/anag/software.html.
4. 3D seismic imaging through reverse-time migration on homogeneous and heterogeneous multi-core processors. Sci;Araya-Polo Mauricio;Program. Cell Process.,2008
Cited by
23 articles.
订阅此论文施引文献
订阅此论文施引文献,注册后可以免费订阅5篇论文的施引文献,订阅后可以查看论文全部施引文献
1. Scalable Distributed High-Order Stencil Computations;SC22: International Conference for High Performance Computing, Networking, Storage and Analysis;2022-11
2. Toward accelerated stencil computation by adapting tensor core unit on GPU;Proceedings of the 36th ACM International Conference on Supercomputing;2022-06-28
3. On the Transformation Optimization for Stencil Computation;Electronics;2021-12-23
4. DRStencil: Exploiting Data Reuse within Low-order Stencil on GPU;2021 IEEE 23rd Int Conf on High Performance Computing & Communications; 7th Int Conf on Data Science & Systems; 19th Int Conf on Smart City; 7th Int Conf on Dependability in Sensor, Cloud & Big Data Systems & Application (HPCC/DSS/SmartCity/DependSys);2021-12
5. Reducing redundancy in data organization and arithmetic calculation for stencil computations;Proceedings of the International Conference for High Performance Computing, Networking, Storage and Analysis;2021-11-13