When one approximates elliptic equations by the spectral collocation method on the Chebyshev-Gauss-Lobatto (CGL) grid, the resulting coefficient matrix is dense and illconditioned. It is known that a good preconditioner, in the sense that the preconditioned system becomes well conditioned, can be constructed with finite difference on the CGL grid. However, there is a lack of an efficient solver...

This paper presents a preconditioner for non-overlapping Schwarz methods applied to the Helmholtz problem. Starting from a simple analytic example, we show how such a preconditioner can be designed by approximating the inverse of the iteration operator for a layered partitioning of the domain. The preconditioner works by propagating information globally by concurrently sweeping in both directio...

Aerodynamic flows involve features with a wide range of spatial and temporal scales which need to be resolved in order to accurately predict desired engineering quantities. While computational fluid dynamics (CFD) has advanced considerably in the past 30 years, the desire to perform more complex, higher-fidelity simulations remains. Present day CFD simulations are limited by the lack of an effi...

We develop a preconditioner for systems arising from space-time finite element discretizations of parabolic equations. The preconditioner is based on a transformation of the coupled system into block diagonal form and an efficient solution strategy for the arising 2 × 2 blocks. The suggested strategy makes use of an inexact factorization of the Schur complement of these blocks, for which unifor...

We present a domain decomposition method suitable for the drift-diiusion equations. The new scheme is applied to the linear systems resulting from Gummel's method and corresponds to a preconditioned conjugate gradient technique for the Schur complement of the linear systems. In designing the preconditioner two problems are addressed : anisotropic phenomena and large scale variations in the magn...

In this paper, the generalized shift-splitting preconditioner is implemented for saddle point problems with symmetric positive definite (1,1)-block and symmetric positive semidefinite (2,2)-block. The proposed preconditioner is extracted form a stationary iterative method which is unconditionally convergent. Moreover, a relaxed version of the proposed preconditioner is presented and some proper...

For solving large dense complex linear systems that arise in electromagnetic calculations, we perform experiments using a general purpose spectral low rank update preconditioner in the context of the GMRES method preconditioned by an approximate inverse preconditioner. The goal of the spectral preconditioner is to improve the convergence properties by shifting by one the smallest eigenvalues of...

A semi-Toeplitz preconditioner for the linearized Navier– Stokes equation for compressible flow is proposed and tested. The preconditioner is applied to the linear system of equations to be solved in each time step of an implicit method. The equations are solved with flat plate boundary conditions and are linearized around the Blasius solution. The grids are stretched in the normal direction to...

2 Hierarchical Preconditioner for Scalar Bases 3 2.1 General Hierarchical Preconditioner Formula . . . . . . . . . . . 3 2.2 Hierarchical Splitting in Case of Nodal Basis Functions . . . . . 4 2.3 Hierarchical Preconditioner Formula for Bases of Nodal Basis Functions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5 2.4 Implementation . . . . . . . . . . . . . . . . . . . . . . . . ...

In this paper, a relaxed Hermitian and skew-Hermitian splitting (RHSS) preconditioner is proposed for saddle point problems from the element-free Galerkin (EFG) discretization method. The EFG method is one of the most widely used meshfree methods for solving partial differential equations. The RHSS preconditioner is constructed much closer to the coefficient matrix than the well-known HSS preco...