Fully Parallel Mesh I/O Using PETSc DMPlex with an Application to Waveform Modeling

Abstract

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) This page was generated automatically upon download from the ETH Zurich Research Collection. For more information, please consult the Terms of use.

Benchmark problems for transcranial ultrasound simulation: Intercomparison of compressional wave models

“…65 It utilizes a matrix-free implementation of the continuous-Galerkin spectral-element method 66 and an explicit second-order Newmark time-stepping scheme. The computational domain is discretized using unstructured conforming hexahedral meshes, 67 which enable the exact representation of interfaces and discontinuities in the tissue parameters. Absorbing boundaries are imposed using the first-order Sommerfeld radiation condition in addition to sponge layers.…”

Parallel Metric-Based Mesh Adaptation in PETSc using ParMmg

“…This enables algorithms to be written independently of mesh properties such as dimension and cell type [12]. DMPlex supports all of the operations required to handle meshing, such as generation, partitioning and distribution, creating of missing edges and faces, regular refinement, extrusion, traversal, selection, manipulation and I/O [13]. This work adds Riemannian metric utilities, to support metric-based adaptation.…”

Toward highly parallel loading of unstructured meshes

“…XDMF is used in many common visualization tools and e.g. PETSc [14] supports its parallel reading [15].…”

“…Elements SFC indices are used for computation of splitters for the second parallel sort (line 14). According to computed splitters, both nodes and elements are again redistributed to target processes (lines 15,16). Contrary to the first sort, SFC assures that spatially collocated elements and nodes are transferred to the same or a near process.…”