trlib

• Verfasst von: Lenders, Felix [VerfasserIn]
• Verfasst von: Kirches, Christian [VerfasserIn]
• Verfasst von: Potschka, Andreas [VerfasserIn]
• Titel: trlib
• Titelzusatz: a vector-free implementation of the GLTR method for iterative solution of the trust region problem
• Verf.angabe: F. Lenders, C. Kirches, A. Potschka
• E-Jahr: 2018
• Jahr: 27 Mar 2018
• Umfang: 30 S.
• Fussnoten: Gesehen am 26.02.2020
• Titel Quelle: Enthalten in: Optimization methods & software
• Ort Quelle: London [u.a.] : Taylor & Francis, 1992
• Jahr Quelle: 2018
• Band/Heft Quelle: 33(2018), 3, Seite 420-449
• ISSN Quelle: 1029-4937
• DOI: doi:10.1080/10556788.2018.1449842
• Datenträger: Online-Ressource
• Sprache: eng
• Sach-SW: 35Q90
• Sach-SW: 65K05
• Sach-SW: 90C20
• Sach-SW: 90C30
• Sach-SW: 97N90
• Sach-SW: iterative method
• Sach-SW: Krylov subspace method
• Sach-SW: PDE-constrained optimization
• Sach-SW: trust region subproblem
• K10plus-PPN: 1691046248

Abstract

We describe trlib, a library that implements a variant of Generalized Lanczos method [Gould et al. in SIAM J. Opt. 9(2), 504-525, 1999] for solving the trust region problem. Our implementation has several distinct features that set it apart from preexisting ones. We implement both conjugate gradient (CG) and Lanczos iterations for assembly of Krylov subspaces. A vector- and matrix-free reverse communication interface allows the use of most general data structures, such as those arising after discretization of function space problems. The hard case of the trust region problem frequently arises in sequential methods for nonlinear optimization. In this implementation, we made an effort to fully address the hard case in an exact way by considering all invariant Krylov subspaces. We investigate the numerical performance of trlib on the full subset of unconstrained problems of the CUTEst benchmark set. In addition to this, interfacing the PDE discretization toolkit FEniCS with trlib using the vector-free reverse communication interface is demonstrated for a family of PDE-constrained control trust region problems adapted from the OPTPDE collection.