On the Derivation of Quasi-Newton Formulas for Optimization in Function Spaces. (11th July 2020)

• Record Type: Journal Article
• Authors: Vuchkov, Radoslav G.; Petra, Cosmin G.; Petra, Noémi
• Is Part Of: Numerical functional analysis and optimization. Volume 41:Number 13(2020)
• Page Start: 1564
• Page End: 1587
• Publication Date: 2020-07-11
• Subjects: BFGS -- DFP -- optimization in infinite dimensions -- PDE-constrained optimization -- PSB -- Quasi-Newton -- SR1 -- variational problems; 90C53 -- 90C30 -- 65K10 -- 46N10 -- 35R30 -- 35Q93; Functional analysis -- Periodicals; Numerical analysis -- Periodicals; Mathematical optimization -- Periodicals; Numerical Analysis, Computer-Assisted; 515.705
• DOI: 10.1080/01630563.2020.1785496
• Languages: English
• ISSNs: 0163-0563

Abstract

Abstract: Newton's method is usually preferred when solving optimization problems due to its superior convergence properties compared to gradient-based or derivative-free optimization algorithms. However, deriving and computing second-order derivatives needed by Newton's method often is not trivial and, in some cases, not possible. In such cases quasi-Newton algorithms are a great alternative. In this paper, we provide a new derivation of well-known quasi-Newton formulas in an infinite-dimensional Hilbert space setting. It is known that quasi-Newton update formulas are solutions to certain variational problems over the space of symmetric matrices. In this paper, we formulate similar variational problems over the space of bounded symmetric operators in Hilbert spaces. By changing the constraints of the variational problem we obtain updates (for the Hessian and Hessian inverse) not only for the Broyden-Fletcher-Goldfarb-Shanno (BFGS) quasi-Newton method but also for Davidon–Fletcher–Powell (DFP), Symmetric Rank One (SR1), and Powell-Symmetric-Broyden (PSB). In addition, for an inverse problem governed by a partial differential equation (PDE), we derive DFP and BFGS "structured" secant formulas that explicitly use the derivative of the regularization and only approximates the second derivative of the misfit term. We show numerical results that demonstrate the desired mesh-independence property and superior performance of the resulting quasi-Newton methods.