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.

Hierarchical subspace evolution method for super large parallel computing: A linear solver and an eigensolver as examples

Large Scale Numerical Simulation Via Parallelization and Reconfigurable Computing Hardware

Scalable Hierarchical Parallel Algorithm for the Solution of Super Large-Scale Sparse Linear Equations

The Hierarchical Subspace Iteration Method for Laplace--Beltrami Eigenproblems

A Parallel Function Evaluation Approach for Solution to Large-Scale Equation-Oriented Models.

Parallel Computation Method For Solving Large Scale Equation-Oriented Models

Hierarchical Node Method for Solving Large-Scale Sparse Linear Equations in Parallel

SEBSM-Based Iterative Method for Solving Large Systems of Linear Equations and Its Applications in Engineering Computation

Fast Calculation of Laplace-Beltrami Eigenproblems Via Subdivision Linear Subspace

A Parallel Iterative Method for Solving Periodical Block-Tridiagonal Linear Equations.

Parallel Computing of Large Eigenvalue Problems for Engineering Structures

Parallel Algorithms for Solving Semi Linear Elliptic Equations

Towards a Scalable Hierarchical High-order CFD Solver

The Improved Krylov Subspace Methods for Large and Sparse Linear Systems on Bulk Synchronous Parallel Architectures

Subspace recursive Fermi-operator expansion strategies for large-scale DFT eigenvalue problems on HPC architectures

A Hybrid Granularity Parallel Algorithm to Solve the Large Structure Eigenproblem

A Method for Solving Large Sparse Unsymmetric Linear Equations

Iterative Solvers for Linear Systems on Parallel Computers

Scalable parallel multilevel algorithms for solving partial differential equations

A new hybrid solver with two‐level parallel computing for large‐scale structural analysis

A Hierarchical Tridiagonal System Solver for Heterogenous Supercomputers