Abstract:Partial differential equations are discretized using the mortar finite element method, where the mortar space contains piecewise quadratic and -cubic functions. We prove the wellposedness of the resulting saddle point system under a theoretical framework by Braess et al. [SIAM J. Numer. Analy. (1999)]. For solving this saddle point system, we introduce two domain decomposition (DD) algorithms. Under the assumption that every local problem approximates the global problem well in the interior of its own partition, the required communication complexity is proportional to the size of the interface. With a globally conforming fine mesh, a DD algorithm by Bank et al. [SIAM Rev. (2003), SIAM J. Sci. Comput. (2004)] is used, in which both the interface solution and the interface residual of each local problem are sent to the neighboring processors. To accommodate globally nonconforming fine meshes, we propose a new algorithm where only the interface solution is communicated. Consequently, the communication complexity is reduced by 1/2. The convergence property of these DD algorithms are studied on one-dimensional two-grid meshes. If the interior residual of each local problem is also communicated, the DD algorithms are shown to converge in one iteration. Otherwise, these algorithms are shown to converge in two iterations. Numerical results are presented, on solving the Poisson equation and a convection diffusion equation with a highly adaptive mesh underlying the finite element discretization. The scalability and parallel efficiency of our DD algorithms are studied. As the number of processors increases, linear and logarithmic speed-up of the solving time is observed for the convection-diffusion equation and the Poisson equation respectively. Furthermore, several parallel algebraic multilevel algorithms are introduced. Three two-level and one three-level algorithms are proposed and implemented, together with an overlapping Schwartz method, which uses a three-level algorithm followed by a two-level one at each iteration. The scalability results indicate that these algorithms are effective in reducing the communication complexity during the solving stage. In summary, our algorithms improve the coordination among processors during the solving stage and may provide a powerful tool for scientific computing.

Automated MPI-X code generation for scalable finite-difference solvers

Large Scale Numerical Simulation Via Parallelization and Reconfigurable Computing Hardware

Automating GPU Scalability for Complex Scientific Models: Phonon Boltzman Transport Equation

Code generation for productive portable scalable finite element simulation in Firedrake

Scalability of high-performance PDE solvers

A multi-physics compiler for generating numerical solvers from differential equations

Mat2Stencil: A Modular Matrix-Based DSL for Explicit and Implicit Matrix-Free PDE Solvers on Structured Grid

Accelerating Dissipative Particle Dynamics Simulations on GPUs: Algorithms, Numerics and Applications

MPI+X:Massive Parallelization and Dynamic Load Balance of a Production-level Unstructured DSMC Solver

Scaling through abstractions -- high-performance vectorial wave simulations for seismic inversion with Devito

A shared compilation stack for distributed-memory parallelism in stencil DSLs

A Scalable Hybrid Algorithm for Solving Partial Differential Equations on a Cluster of CPU/GPU

Parallel Implicit Solvers for 2D Numerical Models on Structured Meshes

Automatic Code Generation and Optimization of Large-scale Stencil Computation on Many-core Processors.

Scalable parallel multilevel algorithms for solving partial differential equations

Automated translation and accelerated solving of differential equations on multiple GPU platforms

Devito: Towards a generic Finite Difference DSL using Symbolic Python

A Domain Specific Language for Performance Portable Molecular Dynamics Algorithms

An Application Driven Method for Assembling Numerical Schemes for the Solution of Complex Multiphysics Problems

Finite-element-wise Domain Decomposition Iterative Solvers with Polynomial Preconditioning.

Scaling up Probabilistic PDE Simulators with Structured Volumetric Information