Impact of mixed precision and storage layout on additive Schwarz smoothers

accepted for publication
  • The growing discrepancy between CPU computing power and memory bandwidth drives more and more numerical algorithms into a bandwidth-bound regime. One example is the overlapping Schwarz smoother, a highly effective building block for iterative multigrid solution of elliptic equations with higher order finite elements. Two options of reducing the required memory bandwidth are sparsity exploiting storage layouts and representing matrix entries with reduced precision in floating point or fixed point format. We investigate the impact of several options on storage demand and contraction rate, both analytically in the context of subspace correction methods and numerically at an example of solid mechanics. Both perspectives agree on the favourite scheme: fixed point representation of Cholesky factors in nested dissection storage.
Author:Jakob Schneck, Martin WeiserORCiD, Florian WendeORCiD
Document Type:Article
Parent Title (English):Numerical Linear Algebra with Applications
Year of first publication:2021