Multilevel Interior Penalty Methods on GPUs
<|reference_start|>Multilevel Interior Penalty Methods on GPUs: We present a matrix-free multigrid method for high-order discontinuous Galerkin (DG) finite element methods with GPU acceleration. A performance analysis is conducted, comparing various data and compute layouts. Smoother implementations are optimized through localization and fast diagonalization techniques. Leveraging conflict-free access patterns in shared memory, arithmetic throughput of up to 39% of the peak performance on Nvidia A100 GPUs are achieved. Experimental results affirm the effectiveness of mixed-precision approaches and MPI parallelization in accelerating algorithms. Furthermore, an assessment of solver efficiency and robustness is provided across both two and three dimensions, with applications to Poisson problems.<|reference_end|>
@article{cui2024multilevel,
    title={Multilevel Interior Penalty Methods on GPUs},
    author={Cu Cui, Guido Kanschat},
    journal={arXiv preprint arXiv:2405.18982},
    year={2024},
    archivePrefix={arXiv},
    eprint={2405.18982},
    primaryClass={math.NA cs.NA}
}