Reduced Sum Implementation of the BURA Method for Spectral Fractional Diffusion Problems
Abstract
The numerical solution of spectral fractional diffusion problems in the form is studied, where is a selfadjoint elliptic operator in a bounded domain , and . The finite difference approximation of the problem leads to the system , where is a sparse, symmetric and positive definite (SPD) matrix, and is defined by its spectral decomposition. In the case of finite element approximation, is SPD with respect to the dot product associated with the mass matrix. The BURA method is introduced by the best uniform rational approximation of degree of in , denoted by . Then the approximation has the form , , thus requiring the solving of auxiliary linear systems with sparse SPD matrices. The BURA method has almost optimal computational complexity, assuming that an optimal PCG iterative solution method is applied to the involved auxiliary linear systems. The presented analysis shows that the absolute values of first % can be extremely large. In such a case the condition number of is practically equal to one. Obviously, such systems do not need preconditioning. The next question is if we can replace their solution by directly multiplying with . Comparative analysis of numerical results is presented as a proof-of-concept for the proposed RS-BURA method.
Cite
@article{arxiv.2105.09048,
title = {Reduced Sum Implementation of the BURA Method for Spectral Fractional Diffusion Problems},
author = {Stanislav Harizanov and Nikola Kosturski and Ivan Lirkov and Svetozar Margenov and Yavor Vutov},
journal= {arXiv preprint arXiv:2105.09048},
year = {2021}
}
Comments
8 pages, 4 figures