An MHSS-like iteration method for two-by-two linear systems with application to FDE optimization problems

Abstract

In this paper, we exploit the numerical solvers for a fractional differential equations (FDE) optimization problem with constraints given by fractional elliptic state equations. The discretized linear system can be rewritten as a two-by-two linear system. By making use of the strategy of modified Hermitian and skew-Hermitian splitting (MHSS) iteration method proposed by Bai, we propose an MHSS-like iteration method. Comparing with the MHSS iteration method, the MHSS-like method only requires one to solve the two-by-two linear systems by a sparse Cholesky factorization and a fast Fourier transform. Hence, the MHSS-like method has less workload than the MHSS method. The convergence properties and the quasi-optimal parameters are analyzed in detail. Numerical examples are used to testify the efficiency of the new method.