On an efficient line smoother for the p-multigrid {\gamma}-cycle
Abstract
As part of the development of a Poisson solver for the spectral element discretization used in the GeoFluid Object Workbench (GeoFLOW) code, we propose a solver for the linear system arising from a Gauss-Legendre-Lobatto global spectral method. We precondition using a p-multigrid {\gamma}-cycle with highly-vectorizable smoothers, that we refer to as line smoothers. Our smoothers are restrictions of spectral and finite element discretizations to low-order one-dimensional problems along lines, that are solved by a reformulation of cyclic reduction as a direct multigrid method. We illustrate our method with numerical experiments showing the apparent boundedness of the iteration count for a fixed residual reduction over a range of moderately deformed domains, right hand sides and Dirichlet boundary conditions.
Cite
@article{arxiv.2504.10710,
title = {On an efficient line smoother for the p-multigrid {\gamma}-cycle},
author = {José Pablo Lucero Lorca and Duane Rosenberg and Isidora Jankov and Conor McCoid and Martin Jakob Gander},
journal= {arXiv preprint arXiv:2504.10710},
year = {2025}
}