An optimized treatment for algorithmic differentiation of an important glaciological fixed-point problem

Abstract : We apply an optimized method to the adjoint generation of a time-evolving land ice model through algorithmic differentiation (AD). The optimization involves a special treatment of the fixed-point iteration required to solve the nonlinear stress balance, which differs from a straightforward application of AD software, and leads to smaller memory requirements and in some cases shorter computation times of the adjoint. The optimization is done via implementation of the algorithm of 5 Christianson [1994] for reverse accumulation of fixed-point problems, with the AD tool OpenAD. For test problems, the optimized adjoint is shown to have far lower memory requirements, potentially enabling larger problem sizes on memory-limited machines. In the case of the land ice model, implementation of the algorithm allows further optimization by having the adjoint model solve a sequence of linear systems with identical (as opposed to varying) matrices, greatly improving per-10 formance. The methods introduced here will be of value to other efforts applying AD tools to ice models, particularly ones which solve a " hybrid " shallow ice / shallow shelf approximation to the Stokes equations.
Type de document :
Article dans une revue
Geoscientific Model Development, European Geosciences Union, 2016, 9 (5), pp.27. 〈http://www.geosci-model-dev.net/9/1891/2016/〉
Liste complète des métadonnées

https://hal.inria.fr/hal-01413295
Contributeur : Laurent Hascoet <>
Soumis le : vendredi 9 décembre 2016 - 15:54:20
Dernière modification le : vendredi 12 janvier 2018 - 11:02:49

Fichier

fp_adjoint_ice_revision1.pdf
Fichiers produits par l'(les) auteur(s)

Identifiants

  • HAL Id : hal-01413295, version 1

Collections

Citation

Daniel Goldberg, Sri Hari Krishna Narayanan, Laurent Hascoet, Jean Utke. An optimized treatment for algorithmic differentiation of an important glaciological fixed-point problem. Geoscientific Model Development, European Geosciences Union, 2016, 9 (5), pp.27. 〈http://www.geosci-model-dev.net/9/1891/2016/〉. 〈hal-01413295〉

Partager

Métriques

Consultations de la notice

109

Téléchargements de fichiers

44