English

Smoothing analysis of two robust multigrid methods for elliptic optimal control problems

Numerical Analysis 2022-03-25 v1 Numerical Analysis Optimization and Control

Abstract

In this paper we study and compare two multigrid relaxation schemes with coarsening by two, three, and four for solving elliptic sparse optimal control problems with control constraints. First, we perform a detailed local Fourier analysis (LFA) of a well-known collective Jacobi relaxation (CJR) scheme, where the optimal smoothing factors are derived. This insightful analysis reveals that the optimal relaxation parameters depend on mesh size and regularization parameters, which was not investigated in literature. Second, we propose and analyze a new mass-based Braess-Sarazin relaxation (BSR) scheme, which is proven to provide smaller smoothing factors than the CJR scheme when αch4\alpha\ge ch^4 for a small constant cc. Here α\alpha is the regularization parameter and hh is the spatial mesh step size. These schemes are successfully extended to control-constrained cases through the semi-smooth Newton method. Coarsening by three or four with BSR is competitive with coarsening by two. Numerical examples are presented to validate our theoretical outcomes. The proposed inexact BSR (IBSR) scheme, where two preconditioned conjugate gradients iterations are applied to the Schur complement system, yields better computational efficiency than the CJR scheme.

Keywords

Cite

@article{arxiv.2203.13066,
  title  = {Smoothing analysis of two robust multigrid methods for elliptic optimal control problems},
  author = {Yunhui He and Jun Liu},
  journal= {arXiv preprint arXiv:2203.13066},
  year   = {2022}
}

Comments

23 pages; 2 tables and 5 figures

R2 v1 2026-06-24T10:24:40.984Z