## Domain Decomposition Preconditioners for Hermite Collocation Problems

##### Abstract

Accelerating the convergence rate of Krylov subspace methods with parallelizable preconditioners is essential for obtaining effective iterative solvers for very large linear systems of equations. Substructuring provides a framework for constructing robust and parallel preconditioners for linear systems arising from the discretization of boundary value problems. Although collocation is a very general and effective discretization technique for many PDE problems, there has been relatively little work on preconditioners for collocation problems.

This thesis proposes two preconditioning methods for solving linear
systems of equations arising from Hermite bicubic collocation
discretization of elliptic partial differential equations on square
domains with mixed boundary conditions.
The first method, called *edge preconditioning*,
is based on a decomposition of the domain in parallel strips,
and the second, called *edge-vertex preconditioning*,
is based on a two-dimensional decomposition.
The preconditioners are derived in terms of two special rectangular
grids -- a coarse grid with diameter *H* and a hybrid coarse/fine
grid -- which together with the fine grid of diameter *h* provide the
framework for approximating the interface problem induced by
substructuring.

We show that the proposed methods are effective for nonsymmetric
indefinite problems, both from the point of view of the cost per
iteration and of the number of iterations. For an appropriate choice
of *H*, the edge preconditioner requires
*O(N)* arithmetic operations per
iteration, while the edge-vertex preconditioner requires
*O(N ^{ 4/3 })* operations,
where

*N*is the number of unknowns. For the edge-vertex preconditioner, the number of iterations is almost constant when

*h*and

*H*decrease such that

*H/h*is held constant and it increases very slowly with

*H*when

*h*is held constant. For both the edge- and edge-vertex preconditioners the number of iterations depends only weakly on

*h*when

*H*is constant. The edge-vertex preconditioner outperforms the edge-preconditioner for small enough

*H*. Numerical experiments illustrate the parallel efficiency of the preconditioners which is similar or even better than that provided by the well-known PETSc parallel software library for scientific computing.

##### Collections

- Doctoral Dissertations [11232]