Search

We study sweeping preconditioners for symmetric and positive definite block tridiagonal systems of linear equations. The algorithm provides an approximate inverse that can be used directly or in a preconditioned iterative scheme. These algorithms are based on replacing the Schur complements appearing in a block Gaussian elimination direct solve by hierarchical matrix approximations with reduced off-diagonal ranks. This involves developing low rank hierarchical approximations to inverses. We first provide a convergence analysis for the algorithm for reduced rank hierarchical inverse approximation. These results are then used to prove convergence and preconditioning estimates for the resulting sweeping preconditioner. Copyright © 2014 John Wiley & Sons, Ltd. [ABSTRACT FROM AUTHOR]

This paper analyses a two-level preconditioning scheme for H(curl) bilinear forms. The scheme utilizes an auxiliary problem on a related mesh that is more amenable for constructing optimal order multigrid methods. More specifically, we analyse the case when the auxiliary mesh only approximately covers the original domain. The latter assumption is important since it allows for easy construction of nested multilevel spaces on regular auxiliary meshes. Numerical experiments in both two and three space dimensions illustrate the optimal performance of the method. Published in 2007 by John Wiley & Sons, Ltd. [ABSTRACT FROM AUTHOR]

We present a negative-norm least-squares method for axisymmetric div–curl systems arising from Maxwell's equations for electrostatics and magnetostatics in three dimensions. The method approximates the solution in a two-dimensional meridian plane. To achieve this dimension reduction, we must work with weighted spaces in cylindrical co-ordinates. In this setting, a stable pair of approximation spaces is developed and analysed. We also report the results of some numerical experiments, which demonstrate a quasi-optimal convergence rate and a robustness with respect to the domain and coefficients. Copyright © 2006 John Wiley & Sons, Ltd. [ABSTRACT FROM AUTHOR]

We consider the Laplace equation under mixed boundary conditions on a polygonal domain Ω. Regularity estimates in terms of Sobolev norms of fractional order for this type of problem are proved. The analysis is based on new interpolation results and multilevel representation of norms on the Sobolev spaces Hα(Ω). The Fourier transform and the construction of extension operators to Sobolev spaces on ℝ2 are avoided in the proofs of the interpolation theorems. Copyright © 2002 John Wiley & Sons, Ltd. [ABSTRACT FROM AUTHOR]

In this paper, we consider approximation of a second-order elliptic problem defined on a domain in two-dimensional Euclidean space. Partitioning the domain into two subdomains, we consider a technique proposed by Wieners and Wohlmuth [9] for coupling mixed finite element approximation on one subdomain with a standard finite element approximation on the other. In this paper, we study the iterative solution of the resulting linear system of equations. This system is symmetric and indefinite (of saddle-point type). The stability estimates for the discretization imply that the algebraic system can be preconditioned by a block diagonal operator involving a preconditioner for H(div) (on the mixed side) and one for the discrete Laplacian (on the finite element side). Alternatively, we provide iterative techniques based on domain decomposition. Utilizing subdomain solvers, the composite problem is reduced to a problem defined only on the interface between the two subdomains. We prove that the interface problem is symmetric, positive definite and well conditioned and hence can be effectively solved by a conjugate gradient iteration. Copyright © 2001 John Wiley & Sons, Ltd. [ABSTRACT FROM AUTHOR]