Sara Faghih-Naini, Vadym Aizinger, Sebastian Kuckuk, Richard Angersbach, Harald Köstler
{"title":"p-adaptive discontinuous Galerkin method for the shallow water equations on heterogeneous computing architectures","authors":"Sara Faghih-Naini, Vadym Aizinger, Sebastian Kuckuk, Richard Angersbach, Harald Köstler","doi":"arxiv-2311.11348","DOIUrl":null,"url":null,"abstract":"Heterogeneous computing and exploiting integrated CPU-GPU architectures has\nbecome a clear current trend since the flattening of Moore's Law. In this work,\nwe propose a numerical and algorithmic re-design of a p-adaptive\nquadrature-free discontinuous Galerkin method (DG) for the shallow water\nequations (SWE). Our new approach separates the computations of the\nnon-adaptive (lower-order) and adaptive (higher-order) parts of the\ndiscretization form each other. Thereby, we can overlap computations of the\nlower-order and the higher-order DG solution components. Furthermore, we\ninvestigate execution times of main computational kernels and use automatic\ncode generation to optimize their distribution between the CPU and GPU. Several\nsetups, including a prototype of a tsunami simulation in a tide-driven flow\nscenario, are investigated, and the results show that significant performance\nimprovements can be achieved in suitable setups.","PeriodicalId":501256,"journal":{"name":"arXiv - CS - Mathematical Software","volume":"13 4","pages":""},"PeriodicalIF":0.0000,"publicationDate":"2023-11-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":null,"platform":"Semanticscholar","paperid":null,"PeriodicalName":"arXiv - CS - Mathematical Software","FirstCategoryId":"1085","ListUrlMain":"https://doi.org/arxiv-2311.11348","RegionNum":0,"RegionCategory":null,"ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"","JCRName":"","Score":null,"Total":0}
引用次数: 0
Abstract
Heterogeneous computing and exploiting integrated CPU-GPU architectures has
become a clear current trend since the flattening of Moore's Law. In this work,
we propose a numerical and algorithmic re-design of a p-adaptive
quadrature-free discontinuous Galerkin method (DG) for the shallow water
equations (SWE). Our new approach separates the computations of the
non-adaptive (lower-order) and adaptive (higher-order) parts of the
discretization form each other. Thereby, we can overlap computations of the
lower-order and the higher-order DG solution components. Furthermore, we
investigate execution times of main computational kernels and use automatic
code generation to optimize their distribution between the CPU and GPU. Several
setups, including a prototype of a tsunami simulation in a tide-driven flow
scenario, are investigated, and the results show that significant performance
improvements can be achieved in suitable setups.