A multilevel approach for computing the limited-memory Hessian and its inverse in variational data assimilation
Brown, K.L. ; Gejadze, I. ; Ramage, A.
Type de document
Article de revue scientifique à comité de lecture
Affiliation de l'auteur
UNIVERSITY OF STRATHCLYDE GLASGOW GBR ; IRSTEA MONTPELLIER UMR G-EAU FRA ; UNIVERSITY OF STRATHCLYDE GLASGOW GBR
Résumé / Abstract
Use of data assimilation techniques is becoming increasingly common across many application areas. The inverse Hessian (and its square root) plays an important role in several different aspects of these processes. In geophysical and engineering applications, the Hessian-vector product is typically defined by sequential solution of a tangent linear and adjoint problem; for the inverse Hessian, however, no such definition is possible. Frequently, the requirement to work in a matrix-free environment means that compact representation schemes are employed. In this paper, we propose an enhanced approach based on a new algorithm for constructing a multilevel eigenvalue decomposition of a given operator, which results in a much more efficient compact representation of the inverse Hessian (and its square root). After introducing these multilevel approximations, we investigate their accuracy and demonstrate their efficiency (in terms of reducing memory requirements and/or computational time) using the example of preconditioning a Gauss-Newton minimization procedure.
SIAM Journal on Scientific Computing, vol. 38, num. 5, p. A2934 - A2963