Approximate tensor-product preconditioners for very high order discontinuous Galerkin methods

From MaRDI portal
Publication:1700740



Abstract: In this paper, we develop a new tensor-product based preconditioner for discontinuous Galerkin methods with polynomial degrees higher than those typically employed. This preconditioner uses an automatic, purely algebraic method to approximate the exact block Jacobi preconditioner by Kronecker products of several small, one-dimensional matrices. Traditional matrix-based preconditioners require mathcalO(p2d) storage and mathcalO(p3d) computational work, where p is the degree of basis polynomials used, and d is the spatial dimension. Our SVD-based tensor-product preconditioner requires mathcalO(pd+1) storage, mathcalO(pd+1) work in two spatial dimensions, and mathcalO(pd+2) work in three spatial dimensions. Combined with a matrix-free Newton-Krylov solver, these preconditioners allow for the solution of DG systems in linear time in p per degree of freedom in 2D, and reduce the computational complexity from mathcalO(p9) to mathcalO(p5) in 3D. Numerical results are shown in 2D and 3D for the advection and Euler equations, using polynomials of degree up to p=15. For many test cases, the preconditioner results in similar iteration counts when compared with the exact block Jacobi preconditioner, and performance is significantly improved for high polynomial degrees p.




Cites work


Cited in
(37)








This page was built for publication: Approximate tensor-product preconditioners for very high order discontinuous Galerkin methods

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q1700740)