Abstract
This paper deals with a multilevel construction of hierarchical matrix approximations to the inverses of finite element stiffness matrices. Given a sequence of discretizations $A_{\ell} x_{\ell} = f_{\ell}$, $\ell=0,\ldots,n$, where $A_0 x_0 = f_0$ denotes the coarse grid problem, we will compute $A^{-1}_0$ exactly and then use interpolation to obtain an $\mathcal{H}$‐matrix approximation $A^{-\cal H}_{\ell+1}$ from the approximate $\mathcal{H}$‐matrix inverse $A^{-\cal H}_{\ell}$ on the next coarser grid. We develop an exact interpolation scheme for the inverse of tridiagonal matrices as they appear in the finite element discretization of one‐dimensional differential equations. We then generalize this approach to two spatial dimensions where these efficiently computed approximations to the inverse may serve as preconditioners in iterative solution methods. We illustrate this approach with some numerical tests for convection‐dominated convection‐diffusion problems.
Talk to us
Join us for a 30 min session where you can share your feedback and ask us any queries you have