A new MOX report entitled “A multigrid method for PDE-constrained optimization with uncertain inputs” by Ciaramella, G.; Nobile, F.; Vanzan, T. has appeared in the MOX Report Collection.
The report can be donwloaded at the following link:
https://www.mate.polimi.it/biblioteca/add/qmox/20/2023.pdf
Abstract: We present a multigrid algorithm to solve efficiently the large saddle-point systems of equations that typically arise in PDE-constrained optimization under uncertainty. The algorithm is based on a collective smoother that at each iteration sweeps over the nodes of the computational mesh, and solves a reduced saddle-point system whose size depends on the number N of samples used to discretized the probability space. We show that this reduced system can be solved with optimal O(N) complexity. We test the multigrid method on three problems: a linear-quadratic problem for which the multigrid method is used to solve directly the linear optimality system; a nonsmooth problem with box constraints and L 1 -norm penalization on the control, in which the multigrid scheme is used within a semismooth Newton iteration; a risk-adverse problem with the smoothed CVaR risk measure where the multigrid method is called within a preconditioned Newton iteration. I! n all cas es, the multigrid algorithm exhibits very good performances and robustness with respect to all parameters of interest.
