Fully parallel mesh I/O using PETSC DMPLEX with an application to waveform modelling
Abstract
Large-scale PDE simulations using high-order finite-element methods on unstructured meshes are an indispensable tool in science and engineering. The widely used open-source PETSc library offers an efficient representation of generic unstructured meshes within its DMPlex module. This paper details our recent implementation of parallel mesh reading and topological interpolation (computation of intermediate codimensions from a cell-vertex mesh) into DMPlex. We apply these developments to seismic wave propagation scenarios on Mars as a sample application. The principal motivation is to overcome single-node memory limits and reach mesh sizes which were impossible before. Moreover, we demonstrate that scalability of I/O and topological interpolation goes beyond 12'000 cores, and memory imposed limits on maximum mesh size vanish. Show more
Publication status
publishedExternal links
Journal / series
SIAM Journal on Scientific ComputingVolume
Pages / Article No.
Publisher
Society for Industrial and Applied MathematicsSubject
Unstructured mesh; Directed acyclic graph; Partitioning; Topological interpolation; Parallel I/O; PETSc; DMPlex; Seismic waveform modeling; Spectral-element methodOrganisational unit
03971 - Fichtner, Andreas / Fichtner, Andreas03971 - Fichtner, Andreas / Fichtner, Andreas
Funding
175322 - Salvus for next-generation full-waveform ultrasonic imaging solutions (SNF)
823844 - Centre of Excellence for Exascale in Solid Earth (EC)
608881 - ETH Zurich Postdoctoral Fellowship Program II (EC)
More
Show all metadata