Abstract:Stencil computation arises from a large variety of scientific and engineering applications and often plays a critical role in the performance of extreme-scale simulations. Due to the memory bound nature, it is a challenging task to optimize stencil computation kernels on many leadership supercomputers, such as Sunway TaihuLight, which has relatively high computing throughput whilst relatively low data-moving capability. In this white paper, we show the efforts we have been making during the past two years in developing end-to-end implementation and optimization techniques for extreme-scale stencil computations on Sunway TaihuLight. We started with a work on optimizing the 3-D 2nd-order 13-point stencil for nonhydrostatic atmospheric dynamics simulation, which is an important part of the 2016 ACM Gordon Bell Prize winning work, and extended it in ways that can handle a broader range of realistic and challenging problems, such as the HPGMG benchmark that consists of memory-hungry stencils and the gaseous wave detonation simulation that relies on complex high-order stencils. The presented stencil computation paradigm on Sunway TaihuLight includes not only multilevel parallelization to exploit the parallelism on different hardware levels, but also systematic performance optimization techniques for communication, memory access, and computation. We show by extreme-scale tests that the proposed systematic stencil computation paradigm can successfully deliver remarkable performance on Sunway TaihuLight with ten million heterogeneous cores. In particular, we achieve an aggregate performance of 23.12 Pflops for the 3-D 5th-order WENO stencil computation in gaseous wave detonation simulation, which is the highest performance result for high-order stencil computations as far as we know, and an aggregate performance of solving over one trillion unknowns per second in the HPGMG benchmark, which ranks the first place in the HPGMG List of Nov 2017.

Dynamic Stencil: Effective Exploitation of Run-Time Resources in Reconfigurable Clusters

Computing Nodes In Reconfigurable Clusters Are Occupied And Released By Applications During Their

DRStencil: Exploiting Data Reuse Within Low-order Stencil on GPU

Mapping Stencils on Coarse-grained Reconfigurable Spatial Architecture

A Portable Framework for Accelerating Stencil Computations on Modern Node Architectures

Scaling and analyzing the stencil performance on multi-core and many-core architectures

Stencil Computations on AMD and Nvidia Graphics Processors: Performance and Tuning Strategies

HW/SW Co-Optimization for Stencil Computation: Beginning with a Customizable Core

Gamify Stencil Dwarf on Cloud for Democratizing Scientific Computing

Optimizing Three-Dimensional Stencil-Operations on Heterogeneous Computing Environments

Cache-Friendly Design for Complex Spatially-Variable Coefficient Stencils on Many-Core Architectures

Stencil-HMLS: A multi-layered approach to the automatic optimisation of stencil codes on FPGA

An approach of processor core customization for stencil computation

Graph-oriented Code Transformation Approach for Register-Limited Stencils on GPUs

Extreme-Scale Realistic Stencil Computations on Sunway TaihuLight with Ten Million Cores

Adapting Combined Tiling to Stencil Optimizations on Sunway Processor

Automatic Code Generation and Optimization of Large-scale Stencil Computation on Many-core Processors.

StencilMART: Predicting Optimization Selection for Stencil Computations Across GPUs

Program Optimization of Stencil Based Application on the GPU-Accelerated System

Employing polyhedral methods to optimize stencils on FPGAs with stencil-specific caches, data reuse, and wide data bursts

Revisiting Temporal Blocking Stencil Optimizations