Approximating sparse Hessian matrices using large-scale linear least squares.

Saved in:
Bibliographic Details
Title: Approximating sparse Hessian matrices using large-scale linear least squares.
Authors: Fowkes, Jaroslav M., Gould, Nicholas I. M., Scott, Jennifer A.
Source: Numerical Algorithms; Aug2024, Vol. 96 Issue 4, p1675-1698, 24p
Subject Terms: HESSIAN matrices, SPARSE matrices, OPTIMIZATION algorithms, LEAST squares, PROBLEM solving
Abstract: Large-scale optimization algorithms frequently require sparse Hessian matrices that are not readily available. Existing methods for approximating large sparse Hessian matrices have limitations. To try and overcome these, we propose a novel approach that reformulates the problem as the solution of a large linear least squares problem. The least squares problem is sparse but can include a number of rows that contain significantly more entries than other rows and are regarded as dense. We exploit recent work on solving such problems using either the normal equations or an augmented system to derive a robust approach for computing approximate sparse Hessian matrices. Example sparse Hessians from the CUTEst test problem collection for optimization illustrate the effectiveness and robustness of the new method. [ABSTRACT FROM AUTHOR]
Copyright of Numerical Algorithms is the property of Springer Nature and its content may not be copied or emailed to multiple sites without the copyright holder's express written permission. Additionally, content may not be used with any artificial intelligence tools or machine learning technologies. However, users may print, download, or email articles for individual use. This abstract may be abridged. No warranty is given about the accuracy of the copy. Users should refer to the original published version of the material for the full abstract. (Copyright applies to all Abstracts.)
Database: Complementary Index
Description
Abstract:Large-scale optimization algorithms frequently require sparse Hessian matrices that are not readily available. Existing methods for approximating large sparse Hessian matrices have limitations. To try and overcome these, we propose a novel approach that reformulates the problem as the solution of a large linear least squares problem. The least squares problem is sparse but can include a number of rows that contain significantly more entries than other rows and are regarded as dense. We exploit recent work on solving such problems using either the normal equations or an augmented system to derive a robust approach for computing approximate sparse Hessian matrices. Example sparse Hessians from the CUTEst test problem collection for optimization illustrate the effectiveness and robustness of the new method. [ABSTRACT FROM AUTHOR]
ISSN:10171398
DOI:10.1007/s11075-023-01681-z