A comparison of discrete schemes for numerical solution of parabolic problems with fractional power elliptic operators

The main aim of this article is to analyze the efficiency of general solvers for parabolic problems with fractional power elliptic operators. Such discrete schemes can be used in the cases of non-constant elliptic operators, non-uniform space meshes and general space domains. The stability results are proved for all algorithms and the accuracy of obtained approximations is estimated by solving well-known test problems. A modification of the second order splitting scheme is presented, it combines the splitting method to solve locally the nonlinear subproblem and the AAA algorithm to solve the nonlocal diffusion subproblem. Results of computational experiments are presented and analyzed