On the Impact of Reordering in a Hierarchical Semi-Separable Compression Solver for Fractional Diffusion Problems

摘要

The performance of a hierarchical solver for systems of linear algebraic equations arising from finite elements (FEM) discretization of Fractional diffusion problems is the subject of this study. We consider the integral definition of Fractional Laplacian in a bounded domain, introduced through the Ritz potential. The problem is non-local and the related FEM system has a dense matrix. We utilize the Structured Matrix Package (STRUMPACK) and its implementation of a Hierarchical Semi-Separable compression in order to solve the system of linear equations. Our main aim is to improve the performance and accuracy of the method by proposing and analyzing 2 schemes for reordering of the unknowns. The numerical tests are run on the high performance cluster AVITOHOL at IICT-BAS.