Efficient solution of symmetric systems of algebraic equations arising from Galerkin boundary element discretizations of 3D elastostatics problems. Part I: direct methods
Stefano Miccoli · Virtual Community of Pathological Anatomy (University of Castilla La Mancha) · 2002
In the present paper a natural extension of the Cholesky factorization, adapted to the peculiar matrix-structure of the solving system of Symmetric Galerkin Boundary Element Methods is presented. The implementation is described in some detail and preliminary timing figures are given, showing that the proposed approach is considerably faster than the Bunch-Kaufman factorization for non-definite matrices. Thread-based parallelism is exploited to show the scalability of the algorithm on moderately parallel shared memory multiprocessors.