Abstract:We present a scalable approach to solve a class of elliptic partial differential equation (PDE)-constrained optimization problems with bound constraints. This approach utilizes a robust full-space interior-point (IP)-Gauss-Newton optimization method. To cope with the poorly-conditioned IP-Gauss-Newton saddle-point linear systems that need to be solved, once per optimization step, we propose two spectrally related preconditioners. These preconditioners leverage the limited informativeness of data in regularized PDE-constrained optimization problems. A block Gauss-Seidel preconditioner is proposed for the GMRES-based solution of the IP-Gauss-Newton linear systems. It is shown, for a large-class of PDE- and bound-constrained optimization problems, that the spectrum of the block Gauss-Seidel preconditioned IP-Gauss-Newton matrix is asymptotically independent of discretization and is not impacted by the ill-conditioning that notoriously plagues interior-point methods. We propose a regularization and log-barrier Hessian preconditioner for the preconditioned conjugate gradient (PCG)-based solution of the related IP-Gauss-Newton-Schur complement linear systems. The scalability of the approach is demonstrated on an example problem with bound and nonlinear elliptic PDE constraints. The numerical solution of the optimization problem is shown to require a discretization independent number of IP-Gauss-Newton linear solves. Furthermore, the linear systems are solved in a discretization and IP ill-conditioning independent number of preconditioned Krylov subspace iterations. The parallel scalability of preconditioner and linear system matrix applies, achieved with algebraic multigrid based solvers, and the aforementioned algorithmic scalability permits a parallel scalable means to compute solutions of a large class of PDE- and bound-constrained problems.

Computational Experience of an Interior-Point SQP Algorithm in a Parallel Branch-and-Bound Framework

A Parallel Method With Hybrid Algorithms For Mixed Integer Nonlinear Programming

An Efficient Interior-Point Algorithm with New Non-Monotone Line Search Filter Method for Nonlinear Constrained Programming

23 European Symposium on Computer Aided Process Engineering: A Parallel Method with Hybrid Algorithms for Mixed Integer Nonlinear Programming

A comparative study of SQP-type algorithms for nonlinear and nonconvex mixed-integer optimization

Integrating SQP and Branch-and-Bound for Mixed Integer Nonlinear Programming

On the Performance of SQP Methods for Nonlinear Optimization

A trust region SQP algorithm for mixed-integer nonlinear programming

A New SQP Algorithm for Large-Scale Nonlinear Programming

A QCQP-based splitting SQP algorithm for two-block nonconvex constrained optimization problems with application

A novel dual-decomposition method for non-convex mixed integer quadratically constrained quadratic problems

Quasi-Newton approaches to Interior Point Methods for quadratic problems

An Interior Point Method For Solving Optimization Problems

A branch-and-bound algorithm for parametric mixed-binary nonlinear programs

An Eigenvalue Decomposition Based Branch-and-bound Algorithm for Nonconvex Quadratic Programming Problems with Convex Quadratic Constraints

Monotone Splitting SQP Algorithms for Two-block Nonconvex Optimization Problems with General Linear Constraints and Applications

A New Branch and Bound Algorithm for Integer Quadratic Programming Problems

Simplex QP-based methods for minimizing a conic quadratic objective over polyhedra

A Scalable Interior-Point Gauss-Newton Method for PDE-Constrained Optimization with Bound Constraints

An SR1/BFGS SQP algorithm for nonconvex nonlinear programs with block-diagonal Hessian matrix