Matrix-free multigrid block-preconditioners for higher order discontinuous Galerkin discretisations

Efficient and suitably preconditioned iterative solvers for elliptic partial differential equations (PDEs) of the convection-diffusion type are used in all fields of science and engineering, including for example computational fluid dynamics, nuclear reactor simulations and combustion models. To ach...

Ausführliche Beschreibung

Gespeichert in:
Bibliographische Detailangaben
Hauptverfasser: Bastian, Peter (VerfasserIn) , Müller, Eike Hermann (VerfasserIn) , Müthing, Steffen (VerfasserIn) , Piatkowski, Stephan-Marian (VerfasserIn)
Dokumenttyp: Article (Journal)
Sprache:Englisch
Veröffentlicht: 4 June 2019
In: Journal of computational physics
Year: 2019, Jahrgang: 394, Pages: 417-439
ISSN:1090-2716
DOI:10.1016/j.jcp.2019.06.001
Online-Zugang:Verlag, Volltext: https://doi.org/10.1016/j.jcp.2019.06.001
Verlag, Volltext: http://www.sciencedirect.com/science/article/pii/S0021999119303973
Volltext
Verfasserangaben:Peter Bastian, Eike Hermann Müller, Steffen Müthing, Marian Piatkowski

MARC

LEADER 00000caa a2200000 c 4500
001 1672251664
003 DE-627
005 20220816214625.0
007 cr uuu---uuuuu
008 190827s2019 xx |||||o 00| ||eng c
024 7 |a 10.1016/j.jcp.2019.06.001  |2 doi 
035 |a (DE-627)1672251664 
035 |a (DE-599)KXP1672251664 
035 |a (OCoLC)1341238647 
040 |a DE-627  |b ger  |c DE-627  |e rda 
041 |a eng 
084 |a 27  |2 sdnb 
100 1 |a Bastian, Peter  |d 1964-  |e VerfasserIn  |0 (DE-588)1045599794  |0 (DE-627)774643374  |0 (DE-576)399098224  |4 aut 
245 1 0 |a Matrix-free multigrid block-preconditioners for higher order discontinuous Galerkin discretisations  |c Peter Bastian, Eike Hermann Müller, Steffen Müthing, Marian Piatkowski 
264 1 |c 4 June 2019 
300 |a 23 
336 |a Text  |b txt  |2 rdacontent 
337 |a Computermedien  |b c  |2 rdamedia 
338 |a Online-Ressource  |b cr  |2 rdacarrier 
500 |a Gesehen am 27.08.2019 
520 |a Efficient and suitably preconditioned iterative solvers for elliptic partial differential equations (PDEs) of the convection-diffusion type are used in all fields of science and engineering, including for example computational fluid dynamics, nuclear reactor simulations and combustion models. To achieve optimal performance, solvers have to exhibit high arithmetic intensity and need to exploit every form of parallelism available in modern manycore CPUs. This includes both distributed- or shared memory parallelisation between processors and vectorisation on individual cores. The computationally most expensive components of the solver are the repeatedapplications of the linear operator and the preconditioner. For discretisations based on higher-order Discontinuous Galerkin methods, sum-factorisation results in a dramatic reduction of the computational complexity of the operator application while, at the same time, the matrix-free implementation can run at a significant fraction of the theoretical peak floating point performance. Multigrid methods for high order methods often rely on block-smoothers to reduce high-frequency error components within one grid cell. Traditionally, this requires the assembly and expensive dense matrix solve in each grid cell, which counteracts any improvements achieved in the fast matrix-free operator application. To overcome this issue, we present a new matrix-free implementation of block-smoothers. Inverting the block matrices iteratively avoids storage and factorisation of the matrix and makes it is possible to harness the full power of the CPU. We implemented a hybrid multigrid algorithm with matrix-free block-smoothers in the high order Discontinuous Galerkin (DG) space combined with a low order coarse grid correction using algebraic multigrid where only low order components are explicitly assembled. The effectiveness of this approach is demonstrated by solving a set of representative elliptic PDEs of increasing complexity, including a convection dominated problem and the stationary SPE10 benchmark. 
650 4 |a Discontinuous Galerkin 
650 4 |a DUNE 
650 4 |a Elliptic PDE 
650 4 |a Matrix-free methods 
650 4 |a Multigrid 
650 4 |a Preconditioners 
700 1 |a Müller, Eike Hermann  |e VerfasserIn  |4 aut 
700 1 |a Müthing, Steffen  |e VerfasserIn  |0 (DE-588)1072062577  |0 (DE-627)826828302  |0 (DE-576)433524367  |4 aut 
700 1 |a Piatkowski, Stephan-Marian  |d 1987-  |e VerfasserIn  |0 (DE-588)1101485523  |0 (DE-627)859741532  |0 (DE-576)469918705  |4 aut 
773 0 8 |i Enthalten in  |t Journal of computational physics  |d Amsterdam : Elsevier, 1961  |g 394(2019), Seite 417-439  |h Online-Ressource  |w (DE-627)266892485  |w (DE-600)1469164-4  |w (DE-576)104193824  |x 1090-2716  |7 nnas  |a Matrix-free multigrid block-preconditioners for higher order discontinuous Galerkin discretisations 
773 1 8 |g volume:394  |g year:2019  |g pages:417-439  |g extent:23  |a Matrix-free multigrid block-preconditioners for higher order discontinuous Galerkin discretisations 
856 4 0 |u https://doi.org/10.1016/j.jcp.2019.06.001  |x Verlag  |x Resolving-System  |3 Volltext 
856 4 0 |u http://www.sciencedirect.com/science/article/pii/S0021999119303973  |x Verlag  |3 Volltext 
951 |a AR 
992 |a 20190827 
993 |a Article 
994 |a 2019 
998 |g 1101485523  |a Piatkowski, Stephan-Marian  |m 1101485523:Piatkowski, Stephan-Marian  |d 700000  |d 708000  |e 700000PP1101485523  |e 708000PP1101485523  |k 0/700000/  |k 1/700000/708000/  |p 4  |y j 
998 |g 1072062577  |a Müthing, Steffen  |m 1072062577:Müthing, Steffen  |d 700000  |d 708000  |e 700000PM1072062577  |e 708000PM1072062577  |k 0/700000/  |k 1/700000/708000/  |p 3 
998 |g 1045599794  |a Bastian, Peter  |m 1045599794:Bastian, Peter  |d 700000  |d 708000  |e 700000PB1045599794  |e 708000PB1045599794  |k 0/700000/  |k 1/700000/708000/  |p 1  |x j 
999 |a KXP-PPN1672251664  |e 3509124626 
BIB |a Y 
SER |a journal 
JSO |a {"id":{"doi":["10.1016/j.jcp.2019.06.001"],"eki":["1672251664"]},"relHost":[{"pubHistory":["1.1966 - 231.2012; Vol. 232.2013 -"],"language":["eng"],"type":{"media":"Online-Ressource","bibl":"periodical"},"part":{"pages":"417-439","text":"394(2019), Seite 417-439","volume":"394","extent":"23","year":"2019"},"title":[{"title":"Journal of computational physics","title_sort":"Journal of computational physics"}],"disp":"Matrix-free multigrid block-preconditioners for higher order discontinuous Galerkin discretisationsJournal of computational physics","recId":"266892485","id":{"issn":["1090-2716"],"eki":["266892485"],"zdb":["1469164-4"]},"origin":[{"publisherPlace":"Amsterdam ; Orlando, Fla.","dateIssuedDisp":"1961-","dateIssuedKey":"1961","publisher":"Elsevier ; Academic Press"}],"note":["Gesehen am 04.06.2020"],"physDesc":[{"extent":"Online-Ressource"}]}],"origin":[{"dateIssuedKey":"2019","dateIssuedDisp":"4 June 2019"}],"note":["Gesehen am 27.08.2019"],"person":[{"given":"Peter","role":"aut","display":"Bastian, Peter","roleDisplay":"VerfasserIn","family":"Bastian"},{"family":"Müller","roleDisplay":"VerfasserIn","role":"aut","display":"Müller, Eike Hermann","given":"Eike Hermann"},{"given":"Steffen","display":"Müthing, Steffen","role":"aut","family":"Müthing","roleDisplay":"VerfasserIn"},{"roleDisplay":"VerfasserIn","family":"Piatkowski","given":"Stephan-Marian","role":"aut","display":"Piatkowski, Stephan-Marian"}],"physDesc":[{"extent":"23 S."}],"language":["eng"],"type":{"bibl":"article-journal","media":"Online-Ressource"},"name":{"displayForm":["Peter Bastian, Eike Hermann Müller, Steffen Müthing, Marian Piatkowski"]},"title":[{"title":"Matrix-free multigrid block-preconditioners for higher order discontinuous Galerkin discretisations","title_sort":"Matrix-free multigrid block-preconditioners for higher order discontinuous Galerkin discretisations"}],"recId":"1672251664"} 
SRT |a BASTIANPETMATRIXFREE4201