Incomplete Iterative Solution of the Subdiffusion Problem
Abstract
In this work, we develop an efficient incomplete iterative scheme for the numerical solution of the subdiffusion model involving a Caputo derivative of order in time. It is based on piecewise linear Galerkin finite element method in space and backward Euler convolution quadrature in time and solves one linear algebraic system inexactly by an iterative algorithm at each time step. We present theoretical results for both smooth and nonsmooth solutions, using novel weighted estimates of the time-stepping scheme. The analysis indicates that with the number of iterations at each time level chosen properly, the error estimates are nearly identical with that for the exact linear solver, and the theoretical findings provide guidelines on the choice. Illustrative numerical results are presented to complement the theoretical analysis.
Cite
@article{arxiv.1906.06497,
title = {Incomplete Iterative Solution of the Subdiffusion Problem},
author = {Bangti Jin and Zhi Zhou},
journal= {arXiv preprint arXiv:1906.06497},
year = {2020}
}
Comments
25 pages