Title | All-at-once Multigrid Methods for Optimality Systems Arising from Optimal Control Problems |
---|---|

Author | Stefan Takacs |

Funding | FWF Doctoral Program "Computational Mathematics" (W1214), Subproject DK 12 - Robust solvers for KKT systems (lead by W. Zulehner) |

Institution | Johannes Kepler University Linz |

Supervisor | Walter Zulehner |

External reviewer | Andrew Wathen (University of Oxford, UK) |

In this thesis we construct and analyze multigrid methods for solving the optimality system characterizing the solution of an optimal control problem. Originally multigrid methods were constructed for elliptic problems. However, the (discretized) optimality system is not elliptic. We make use of the fact that the matrix is a block-matrix with saddle point structure: on the one hand we have two blocks of variables representing state and control. These two blocks are already part of the optimal control problem itself. On the other hand, the conversion to the optimality system requires the introduction of Lagrange multipliers, which form the third block of variables.

There are several possibilities to use multigrid methods for constructing solvers for such saddle point problems. One approach to solve such problems is to apply multigrid methods in every step of an overall block-structured iterative method to equations in just one of these blocks of variables. Another approach, which we will follow here, is to apply the multigrid idea directly to the (reduced or not reduced) optimality system, which is called an all-at-once approach.

The choice of an appropriate smoother is the key issue in constructing such a multigrid method. The other part of the method—the coarse-grid correction—can be set up in a canonical way because we will use a framework of conforming geometric multigrid method. In this framework the smoother will be constructed such that the convergence rates are independent of the grid level. This leads to an overall computational complexity that is linear in the number of unknowns which is called an optimal convergence.

For a sub-class of the class of problems we will introduce in a first point, we will go one step further and construct methods where the convergence rates are also robust in a certain regularization or cost parameter which is part of the problem. Moreover, we will show this fact. Methods that do not take this into account typically show quite poor convergence rates if this parameter attains small values. For the analysis, we will introduce a general framework that is based on Hackbusch's splitting of the analysis into smoothing and approximation property. This allows to give general convergence theorems for the methods under investigation. A second point we consider is local Fourier analysis which allows to compute sharp bounds of the convergence rates.

- PhD thesis as PDF

- A local Fourier convergence analysis of a multigrid method using symbolic computation.
- Journal of Symbolic Computation, vol. 63, p. 1 – 20, 2014.
- DK-Report 2012-04 Supplemental material
- Convergence analysis of all-at-once multigrid methods for elliptic control problems under partial elliptic regularity.
- SIAM Journal on Numerical Analysis, vol. 51 (3), p. 1853 – 1874, 2013.
- DK-Report 2012-08
- Convergence analysis of multigrid methods with collective point smoothers for optimal control problems.
- Computing and Visualization in Science, vol. 14 (3), p. 131 – 141, 2011.
- DK-Report 2011-01
- Smoothing analysis of an all-at-once multigrid approach for optimal control problems using symbolic computation.
- In U. Langer and P. Paule (eds.): Numerical and Symbolic Scientific Computing: Progress and Prospects, p. 175 – 191, 2011.
- DK-Report 2010-09
- Multigrid methods for elliptic optimal control problems with Neumann boundary control.
- In G. Kreiss, P. Lötstedt, A. Malqvist, and M. Neytcheva (eds.): Numerical Mathematics and Advanced Applications 2009, p. 855 – 863, 2010.
- DK-Report 2009-01