Abstract:The International Journal of High Performance Computing Applications, Ahead of Print. Sparse linear iterative solvers are essential for many large-scale simulations. Much of the runtime of these solvers is often spent in the implicit evaluation of matrix polynomials via a sequence of sparse matrix-vector products. A variety of approaches has been proposed to make these polynomial evaluations explicit (i.e., fix the coefficients), e.g., polynomial preconditioners or s-step Krylov methods. Furthermore, it is nowadays a popular practice to approximate triangular solves by a matrix polynomial to increase parallelism. Such algorithms allow to evaluate the polynomial using a so-called matrix power kernel (MPK), which computes the product between a power of a sparse matrix [math] and a dense vector [math], i.e., [math], or a related operation. Recently we have shown that using the level-based formulation of sparse matrix-vector multiplications in the Recursive Algebraic Coloring Engine (RACE) framework we can perform temporal cache blocking of MPK to increase its performance. In this work, we demonstrate the application of this cache-blocking optimization in sparse iterative solvers. By integrating the RACE library into the Trilinos framework, we demonstrate the speedups achieved in (preconditioned) s-step GMRES, polynomial preconditioners, and algebraic multigrid (AMG). For MPK-dominated algorithms we achieve speedups of up to 3[math] on modern multi-core compute nodes. For algorithms with moderate contributions from subspace orthogonalization, the gain reduces significantly, which is often caused by the insufficient quality of the orthogonalization routines. Finally, we showcase the application of RACE-accelerated solvers in a real-world wind turbine simulation (Nalu-Wind) and highlight the new opportunities and perspectives opened up by RACE as a cache-blocking technique for MPK-enabled sparse solvers.

A High Performance Block Eigensolver for Nuclear Configuration Interaction Calculations.

Accelerating Nuclear Configuration Interaction Calculations through a Preconditioned Block Iterative Eigensolver

Optimizing Sparse Matrix-Multiple Vectors Multiplication for Nuclear Configuration Interaction Calculations

An Out-of-Core Eigensolver on SSD-equipped Clusters

Improving the Scalability of A Symmetric Iterative Eigensolver for Multi-Core Platforms

High Performance Optimizations For Nuclear Physics Code Mfdn On Knl

Hybrid Eigensolvers for Nuclear Configuration Interaction Calculations

Accelerating an Iterative Eigensolver for Nuclear Structure Configuration Interaction Calculations on GPUs Using OpenACC

Accelerating Eigenvalue Computation for Nuclear Structure Calculations via Perturbative Corrections

A High-Performance Implementation of Atomistic Spin Dynamics Simulations on x86 CPUs

A Novel Fully Hardware-Implemented SVD Solver Based on Ultra-Parallel BCV Jacobi Algorithm

Performance Optimization for Sparse A(T)Ax in Parallel on Multicore Cpu

Porting of the DBCSR library for Sparse Matrix-Matrix Multiplications to Intel Xeon Phi systems

Scaling of ab-initio nuclear physics calculations on multicore computer architectures

Eigenvalue Solvers for Modeling Nuclear Reactors on Leadership Class Machines

Large-scale parallel null space calculation for nuclear configuration interaction.

Parallel Efficient Sparse Matrix-Matrix Multiplication on Multicore Platforms

Accelerating Configuration Interaction Calculations for Nuclear Structure

Numerical eigen-spectrum slicing, accurate orthogonal eigen-basis, and mixed-precision eigenvalue refinement using OpenMP data-dependent tasks and accelerator offload

Multicore-Based Performance Optimization For Evaluating The Inverse Of Sparse Matrix

Algebraic temporal blocking for sparse iterative solvers on multi-core CPUs