Diagonal scaling of stiffness matrices in the Galerkin boundary element method

Publisher: Cambridge University Press

E-ISSN: 1446-8735|42|1|141-150

ISSN: 1446-1811

Source: ANZIAM Journal, Vol.42, Iss.1, 2000-07, pp. : 141-150

Disclaimer: Any content in publications that violate the sovereignty, the constitution or regulations of the PRC is not accepted or approved by CNPIEC.

Previous Menu Next

Abstract

A boundary integral equation of the first kind is discretised using Galerkin's method with piecewise-constant trial functions. We show how the condition number of the stiffness matrix depends on the number of degrees of freedom and on the global mesh ratio. We also show that diagonal scaling eliminates the latter dependence. Numerical experiments confirm the theory, and demonstrate that in practical computations involving strong local mesh refinement, diagonal scaling dramatically improves the conditioning of the Galerkin equations.