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
Langue
Anglais
Affiliation de l'auteur
UNIVERSITY OF STRATHCLYDE GLASGOW GBR ; IRSTEA MONTPELLIER UMR G-EAU FRA ; UNIVERSITY OF STRATHCLYDE GLASGOW GBR
Année
2016
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.
Source
SIAM Journal on Scientific Computing, vol. 38, num. 5, p. A2934 - A2963

puce  Accés à la notice sur le site Irstea Publications / Display bibliographic record on Irstea Publications website

  Texte intégral / Full text

  Liste complète des notices de CemOA