<mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline" id="d1e10298" altimg="si178.svg"><mml:mi>p</mml:mi></mml:math>-multigrid methods and their comparison to <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline" id="d1e10303" altimg="si179.svg"><mml:mi>h</mml:mi></mml:math>-multigrid methods within Isogeometric Analysis
Roel Tielen, Matthias Möller, Dominik Göddeke, C. Vuik
Abstract
Over the years, Isogeometric Analysis has shown to be a successful alternative to the Finite Element Method (FEM). However, solving the resulting linear systems of equations efficiently remains a challenging task. In this paper, we consider a p-multigrid method, in which coarsening is applied in the spline degree p instead of the mesh width h, and compare it to h-multigrid methods. Since the use of classical smoothers (e.g. Gauss–Seidel) results in a p-multigrid/h-multigrid method with deteriorating performance for higher values of p, the use of an ILUT smoother is investigated as well. Numerical results and a spectral analysis indicate that the use of this smoother exhibits convergence rates essentially independent of h and p for both p-multigrid and h-multigrid methods. In particular, we compare both coarsening strategies (e.g. coarsening in h or p) adopting both smoothers for a variety of two and three dimensional benchmarks. Furthermore, the ILUT smoother is compared to a state-of-the-art smoother (Hofreither and Takacs 2017) using both coarsening strategies. Finally, the proposed p-multigrid method is used to solve linear systems resulting from THB-spline discretizations.