Abstract:SUMMARY For software to fully exploit the computing power of emerging heterogeneous computers, not only must the required computational kernels be optimized for the specific hardware architectures but also an effective scheduling scheme is needed to utilize the available heterogeneous computational units and to hide the communication between them. As a case study, we develop a static scheduling scheme for the tridiagonalization of a symmetric dense matrix on multicore CPUs with multiple graphics processing units (GPUs) on a single compute node. We then parallelize and optimize the Basic Linear Algebra Subroutines (BLAS)‐2 symmetric matrix‐vector multiplication, and the BLAS‐3 low rank symmetric matrix updates on the GPUs. We demonstrate the good scalability of these multi‐GPU BLAS kernels and the effectiveness of our scheduling scheme on twelve Intel Xeon processors and three NVIDIA GPUs. We then integrate our hybrid CPU‐GPU kernel into computational kernels at higher‐levels of software stacks, that is, a shared‐memory dense eigensolver and a distributed‐memory sparse eigensolver. Our experimental results show that our kernels greatly improve the performance of these higher‐level kernels, not only reducing the solution time but also enabling the solution of larger‐scale problems. Because such symmetric eigenvalue problems arise in many scientific and engineering simulations, our kernels could potentially lead to new scientific discoveries. Furthermore, these dense linear algebra algorithms present algorithmic characteristics that can be found in other algorithms. Hence, they are not only important computational kernels on their own but also useful testbeds to study the performance of the emerging computers and the effects of the various optimization techniques. Copyright © 2013 John Wiley & Sons, Ltd.

An efficient hybrid tridiagonal divide-and-conquer algorithm on distributed memory architectures

Parallel Data Partitioning Strategy In Solving Large Scale Electromagnetic Scattering Problems

A Pre-splitting Green's Function Based Hybrid Fast Algorithm for Multiscale Problems

Tridiagonalization of a dense symmetric matrix on multiple GPUs and its application to symmetric eigenvalue problems

Parallel Tridiagonal Solver on Sunway Many-Core Processors*

A Hierarchical Tridiagonal System Solver for Heterogenous Supercomputers

Skew-Symmetric Matrix Decompositions on Shared-Memory Architectures

A Hybrid Parallel Tridiagonal Solver on Multi-core Architectures

swPTS: an efficient parallel Thomas split algorithm for tridiagonal systems on Sunway manycore processors

A Scalable Hybrid Algorithm for Solving Partial Differential Equations on a Cluster of CPU/GPU

An efficient multi-core implementation of a novel HSS-structured multifrontal solver using randomized sampling

A Novel Compute-Efficient Tridiagonal Solver for Many-Core Architectures

A Novel Multi-CPU/GPU Collaborative Computing Framework for SGD-based Matrix Factorization

Accelerating HPCG on Tianhe-2: A hybrid CPU-MIC algorithm

A Fast Tridiagonal Solver for Intel MIC Architecture

$O(N)$ distributed direct factorization of structured dense matrices using runtime systems

A sparsity-aware distributed-memory algorithm for sparse-sparse matrix multiplication

Extending DD-$α$AMG on heterogeneous machines

Hierarchical Jacobi Iteration for Structured Matrices on GPUs using Shared Memory

Scalable Hierarchical Parallel Algorithm for the Solution of Super Large-Scale Sparse Linear Equations

Parallel Solution of Triangular Systems in Direct Numerical Simulation