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.

Caveats of three direct linear solvers for finite element analyses

High Performance Sparse Static Solver In Finite Element Analyses With Loop-Unrolling

Advances In Sparse Solver For Finite Element Analysis

A Comparison of Sparse Solvers for Severely Ill-Conditioned Linear Systems in Geophysical Marker-In-Cell Simulations

New High Performance Sparse Static Solver in Finite Element Analysis with Loop-Unrolling

Evaluating Accuracy and Efficiency of HPC Solvers for Sparse Linear Systems with Applications to PDEs

Scalability of high-performance PDE solvers

SOLVERS FOR SYSTEMS OF LARGE SPARSE LINEAR AND NONLINEAR EQUATIONS BASED ON MULTI-GPUS

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

Iterative Solvers for Linear Systems on Parallel Computers

Implementation of a Parallel Sparse Direct Solver on Vector Architecture

Natural Factor Based Solvers

Scalable parallel multilevel algorithms for solving partial differential equations

Discrete least-squares finite element methods

A Parallel DDM Preconditioner for Solving Finite Element Discretization of Elasticity Problem in Three Dimension

The Cell Index Concept in Fast Sparse Static and Eigenvalue Solvers for Finite Element Analysis

Performance of linear solvers in tensor-train format on current multicore architectures

Fast linear solvers for incompressible CFD simulations with compatible discrete operator schemes

On cost-efficient parallel iterative solvers for 3-D frequency-domain seismic multi-source viscoelastic anisotropic wave modeling

High-efficiency Improved Symmetric Successive Over-Relaxation Preconditioned Conjugate Gradient Method for Solving Large-Scale Finite Element Linear Equations

Computational complexity and memory usage for multi-frontal direct solvers in structured mesh finite elements