Stabilised rounded addition of hierarchical matrices
Mario Bebendorf and Wolfgang Hackbusch
Contact the author: Please use for correspondence this email.
Submission date: 12. Jan. 2007
published in: Numerical linear algebra with applications, 14 (2007) 5, p. 407-423
DOI number (of the published article): 10.1002/nla.525
Download full preprint: PDF (318 kB), PS ziped (353 kB)
The efficiency of hierarchical matrices is based on the approximate evaluation of usual matrix operations. The introduced approximation error may, however, lead to a loss of important matrix properties. In this article we present a technique which preserves the positive definitness of a matrix independently of the approximation quality. The importance of this technique is illustrated by an elliptic mixed boundary value problem with tiny Dirichlet part.