Abstract

Hierarchical ( $$\mathcal {H}$$ -) matrices provide a data-sparse way to approximate fully populated matrices. The two basic steps in the construction of an $$\mathcal {H}$$ -matrix are (a) the hierarchical construction of a matrix block partition, and (b) the blockwise approximation of matrix data by low rank matrices. In the context of finite element discretisations of elliptic boundary value problems, $$\mathcal {H}$$ -matrices can be used for the construction of preconditioners such as approximate $$\mathcal {H}$$ -LU factors. In this paper, we develop a new black box approach to construct the necessary partition. This new approach is based on the matrix graph of the sparse stiffness matrix and no longer requires geometric data associated with the indices like the standard clustering algorithms. The black box clustering and a subsequent $$\mathcal {H}$$ -LU factorisation have been implemented in parallel, and we provide numerical results in which the resulting black box $$\mathcal {H}$$ -LU factorisation is used as a preconditioner in the iterative solution of the discrete (three-dimensional) convection-diffusion equation.

Talk to us

Join us for a 30 min session where you can share your feedback and ask us any queries you have

Schedule a call

Disclaimer: All third-party content on this website/platform is and will remain the property of their respective owners and is provided on "as is" basis without any warranties, express or implied. Use of third-party content does not indicate any affiliation, sponsorship with or endorsement by them. Any references to third-party content is to identify the corresponding services and shall be considered fair use under The CopyrightLaw.