Hierarchical Splitting for Weighted Graph-Laplacians


I. Georgiev, J. Kraus, S. Margenov
IPP, BAS, Sofia

Generalizing the approach of some previous works the authors present multilevel preconditioners for three-dimensional (3D) elliptic problems discretized by a family of Rannacher Turek non-conforming finite elements.
The first major contribution of the study is the derived estimates of the constant in the strengthened Cauchy-Bunyakowski-Schwarz (CBS) inequality which is shown to allow the efficient multilevel extension of the related two-level preconditioners. Representative numerical tests well illustrate the optimal complexity of the resulting iterative solver, also for the case of non-smooth coefficients. The second important achievement concerns the experimental study of AMLI solvers applied to the case of voxel FEM simulation. Here the coefficient jumps are resolved on the finest mesh only and therefore the classical CBS inequality based convergence theory is not directly applicable. The obtained results, however, demonstrate the efficiency of the proposed algorithms in this case also, as it is illustrated by an example of microstructure analysis of bones.