Axel Modave

GPU-accelerated discontinuous Galerkin methods on hybrid meshes

This work was also published as a Rice University thesis/dissertation: http://hdl.handle.net/1911/96151

A nodal discontinuous Galerkin method for reverse-time migration on GPU clusters

Improving both accuracy and computational performance of numerical tools is a major challenge for seismic imaging and generally requires specialized implementations to make full use of modern parallel architectures. We present a computational strategy for reverse-time migration (RTM) with accelerator-aided clusters. A new imaging condition computed from the pressure and velocity fields is introduced. The model solver is based on a high-order discontinuous Galerkin time-domain (DGTD) method for the pressure–velocity system with unstructured meshes and multirate local time stepping. We adopted the MPI+X approach for distributed programming where X is a threaded programming model. In this work we chose OCCA, a unified framework that makes use of major multithreading languages (e.g. CUDA and OpenCL) and offers the flexibility to run on several hardware architectures. DGTD schemes are suitable for efficient computations with accelerators thanks to localized element-to-element coupling and the dense algebraic operations required for each element. Moreover, compared to high-order finite-difference schemes, the thin halo inherent to DGTD method reduces the amount of data to be exchanged between MPI processes and storage requirements for RTM procedures. The amount of data to be recorded during simulation is reduced by storing only boundary values in memory rather than on disk and recreating the forward wavefields. Computational results are presented that indicate that these methods are strong scalable up to at least 32 GPUs for a three-dimensional RTM case.

GPU performance analysis of a nodal discontinuous Galerkin method for acoustic and elastic models

Finite element schemes based on discontinuous Galerkin methods possess features amenable to massively parallel computing accelerated with general purpose graphics processing units (GPUs). However, the computational performance of such schemes strongly depends on their implementation. In the past, several implementation strategies have been proposed. They are based exclusively on specialized compute kernels tuned for each operation, or they can leverage BLAS libraries that provide optimized routines for basic linear algebra operations. In this paper, we present and analyze up-to-date performance results for different implementations, tested in a unified framework on a single NVIDIA GTX980 GPU. We show that specialized kernels written with a one-node-per-thread strategy are competitive for polynomial bases up to the fifth and seventh degrees for acoustic and elastic models, respectively. For higher degrees, a strategy that makes use of the NVIDIA cuBLAS library provides better results, able to reach a net arithmetic throughput 35.7% of the theoretical peak value. HighlightsSeveral GPU implementations for time-domain wave simulations are compared.The numerical schemes are based on a high-order discontinuous finite element method.The implementations are profiled using the roofline model to highlight bottlenecks.The best implementation depends on the polynomial degree of the basis functions.

Perfectly matched layers for convex truncated domains with discontinuous Galerkin time domain simulations

This paper deals with the design of perfectly matched layers (PMLs) for transient acoustic wave propagation in generally-shaped convex truncated domains. After reviewing key elements to derive PML equations for such domains, we present two time-dependent formulations for the pressure-velocity system. These formulations are obtained by using a complex coordinate stretching of the time-harmonic version of the equations in a specific curvilinear coordinate system. The final PML equations are written in a general tensor form, which can easily be projected in Cartesian coordinates to facilitate implementation with classical discretization methods. Discontinuous Galerkin finite element schemes are proposed for both formulations. They are tested and compared using a three-dimensional benchmark with an ellipsoidal truncated domain. Our approach can be generalized to domains with corners.

3-D Modeling of Thin Sheets in the Discontinuous Galerkin Method for Transient Scattering Analysis

This paper presents a modeling of thin sheets. An interface condition based on analytical solution is used to avoid a fine mesh. This condition is integrated in a time-domain discontinuous Galerkin method to evaluate the shielding effectiveness. This approach is validated by a comparison with analytical solution. 2-D and 3-D cavities are simulated to illustrate the efficiency of the condition.