More Web Proxy on the site http://driver.im/

research-article

Checkpointing Schemes for Adjoint Codes: : Application to the Meteorological Model Meso-NH

Author:

I. CharpentierAuthors Info & Claims

SIAM Journal on Scientific Computing, Volume 22, Issue 6

Pages 2135 - 2151

https://doi.org/10.1137/S1064827598343735

Published: 01 January 2001 Publication History

Abstract

The adjoint code of a nonlinear computer model calculates gradients along a trajectory that has to be known at integration time. When the storage of the whole trajectory requires too large an amount of memory, the calculation of the adjoint code is split and is done part by part from restart points called checkpoints. Griewank proposed a checkpointing method named Revolve, which provides an optimal logarithmic behavior with respect to time and memory requirement. In this work, some checkpointing schedules are proposed. Some of them correspond to special cases of Revolve.

The user's preference is essential to choose between time and memory requirements. This is a key point for adjoint codes of temporal models such as the meteorological model Meso-NH that may be used for weather forecasts. When the computational time is the top priority, a particular checkpointing scheme allows computation of the adjoint code with at most one extra integration of the model. The memory requirement behaves then as the square root of the number of iterations of the model. Checkpointing schemes are tested on adjoint simulations of Meso-NH.

References

[1]

I. Charpentier, The MesODiF package for gradients computations with the atmospheric model Meso‐NH, in Proceedings of the Air Pollution and Modeling Systems Conference (Paris, 1998), Environmental Modeling and Software, to appear.

[2]

Isabelle Charpentier, Mohammed Ghemires, Efficient adjoint derivatives: application to the meteorological model Meso‐NH, Optim. Methods Softw., 13 (2000), 35–63

[3]

J. Derber, Variational four‐dimensional analysis using quasi‐geostrophic constraints, Mon. Wea. Rev., 115 (1987), pp. 998–1008.

[4]

C. Faure and I. Charpentier, Comparing strategies for coding adjoints, Scientific Programming, to appear.

[5]

C. Faure and Y. Papegay, Odyssée Version 1.6, the User’s Reference Manual, Rapport technique INRIA RT‐0211, INRIA, Le Chesnay, France, 1997.

[6]

Andreas Griewank, On automatic differentiation, Math. Appl. (Japanese Ser.), Vol. 6, SCIPRESS, Tokyo, 1989, 83–107

[7]

A. Griewank, Achieving logarithmic growth of temporal and spatial complexity in reverse automatic differentiation, Optim. Methods Softw., 1 (1992), pp. 35–54.

[8]

A. Griewank, Evaluating Derivatives: Principles and Techniques of Algorithmic Differentiation, Frontiers Appl. Math. 19, SIAM, Philadelphia, 2000.

[9]

A. Griewank and A. Walther, Revolve: An Implementation of the Checkpointing for the Reverse or Adjoint Mode of Differentiation, preprint IOKOMO‐04‐1997, Technische Univerität Dresden, Dresden, Germany, 1997;

ACM Trans. Math. Software, to appear.

[10]

A. Griewank and A. Walther, Optimal Program Execution Reversal, in preparation.

[11]

J. Grimm, L. Pottier, and N. Rostaing‐Schmidt, Optimal time and minimum space‐time product for reversing a certain class of programs, in Computational Differentiation: Techniques, Applications, and Tools, M. Berz, C. Bischof, G. Corliss, and A. Griewank, eds., SIAM, Philadelphia, 1996, pp. 95–106.

[12]

J. P. Lafore, J. Stein, N. Ascencio, P. Bougeault, V. Ducrocq, J. Duron, C. Fischer, P. Héreil, P. Mascart, V. Masson, J. P. Pinty, J. L. Redelsperger, E. Richard, and J. Vilà‐Guerau de Arellano, The Meso‐NH atmospheric simulation system. Part I: Adiabatic formulation and control system, Ann. Geophysicae, 16 (1998), pp. 90–109.

[13]

F.‐X. Le Dimet and O. Talagrand, Variational algorithms for analysis and assimilation of meteorological observations: Theoretical aspects, Tellus A, 38 (1986), pp. 97–110.

[14]

J. Benamou, Domain decomposition, optimal control of systems governed by partial differential equations, and synthesis of feedback laws, J. Optim. Theory Appl., 102 (1999), 15–36

[15]

J. P. Lafore, J. Stein, N. Ascencio, P. Bougeault, V. Ducrocq, J. Duron, C. Fischer, P. Héreil, P. Mascart, V. Masson, J. P. Pinty, J. L. Redelsperger, E. Richard, and J. Vilà‐Guerau de Arellano, The Meso‐NH Atmospheric Simulation System: Scientific Documentation, Météo‐France and CNRS, Toulouse, France, 1995.

[16]

J. P. Lafore, J. Stein, N. Ascencio, P. Bougeault, V. Ducrocq, J. Duron, C. Fischer, P. Héreil, P. Mascart, V. Masson, J. P. Pinty, J. L. Redelsperger, E. Richard, and J. Vilà‐Guerau de Arellano, The Meso‐NH Atmospheric Simulation System: Algorithmic Documentation, Météo‐France and CNRS, Toulouse, France, 1996.

[17]

J. Morgenstern, How to compute fast a function and all its derivatives, a variation on the theorem of Baur‐Strassen, SIGACT News, 16 (1985), pp. 60–62.

Digital Library

[18]

F. Rabier, H. Järvinen, E. Klinker, J‐F. Mahfouf, and A. Simmons, The ECMWF Operational Implementation of Four Dimensional Variational Assimilation. Part I: Experimental Results with Simplified Physics, ECMWF Research Department Tech. Memo. 271, ECMWF, Reading, UK, 1999.

[19]

A. Walther, Program Reversal Schedules for Single‐ and Multi‐Processor Machines, Thesis, Institute of Scientific Computing, Dresden University of Technology, Dresden, Germany, 1999.

Cited By

Zahr MPersson P(2016)An adjoint method for a high-order discretization of deforming domain conservation laws for optimization of flow problemsJournal of Computational Physics10.1016/j.jcp.2016.09.012326:C(516-543)Online publication date: 1-Dec-2016
https://dl.acm.org/doi/10.1016/j.jcp.2016.09.012
Liao W(2015)An adjoint-based Jacobi-type iterative method for elastic full waveform inversion problemApplied Mathematics and Computation10.1016/j.amc.2015.06.010267:C(56-70)Online publication date: 15-Sep-2015
https://dl.acm.org/doi/10.1016/j.amc.2015.06.010
Lemke MSchulze JSesterhenn J(2014)Adjoint-based reconstruction of an entropy source by discrete temperature measurementsInternational Journal of Computational Science and Engineering10.1504/IJCSE.2014.0645379:5/6(526-537)Online publication date: 1-Sep-2014
https://dl.acm.org/doi/10.1504/IJCSE.2014.064537

Index Terms

Checkpointing Schemes for Adjoint Codes: Application to the Meteorological Model Meso-NH

Index terms have been assigned to the content through auto-classification.

Recommendations

Symplectic Runge--Kutta Schemes for Adjoint Equations, Automatic Differentiation, Optimal Control, and More

The study of the sensitivity of the solution of a system of differential equations with respect to changes in the initial conditions leads to the introduction of an adjoint system, whose discretization is related to reverse accumulation in automatic ...
An improved PC scheme for nonlinear fractional differential equations

A novel computationally effective fractional predictorcorrector (PC) scheme is proposed to solve fractional differential equations involving Caputo derivative. The properties of the Caputo derivative are used to reduce the fractional differential ...
Matrix approach to discrete fractional calculus II: Partial fractional differential equations

A new method that enables easy and convenient discretization of partial differential equations with derivatives of arbitrary real order (so-called fractional derivatives) and delays is presented and illustrated on numerical solution of various types of ...

Comments

Please enable JavaScript to view thecomments powered by Disqus.

Information & Contributors

Information

Published In

cover image SIAM Journal on Scientific Computing

SIAM Journal on Scientific Computing Volume 22, Issue 6

2001

377 pages

ISSN:1064-8275

Issue’s Table of Contents

Copyright © 2001 Society for Industrial and Applied Mathematics.

Publisher

Society for Industrial and Applied Mathematics

United States

Publication History

Published: 01 January 2001

Author Tags

Author Tags

Qualifiers

Research-article

Contributors

Other Metrics

View Article Metrics

Bibliometrics & Citations

Bibliometrics

Article Metrics

3
Total Citations
View Citations
0
Total Downloads

Downloads (Last 12 months)0
Downloads (Last 6 weeks)0

Reflects downloads up to 09 Jan 2025

Other Metrics

View Author Metrics

Citations

Cited By

Zahr MPersson P(2016)An adjoint method for a high-order discretization of deforming domain conservation laws for optimization of flow problemsJournal of Computational Physics10.1016/j.jcp.2016.09.012326:C(516-543)Online publication date: 1-Dec-2016
https://dl.acm.org/doi/10.1016/j.jcp.2016.09.012
Liao W(2015)An adjoint-based Jacobi-type iterative method for elastic full waveform inversion problemApplied Mathematics and Computation10.1016/j.amc.2015.06.010267:C(56-70)Online publication date: 15-Sep-2015
https://dl.acm.org/doi/10.1016/j.amc.2015.06.010
Lemke MSchulze JSesterhenn J(2014)Adjoint-based reconstruction of an entropy source by discrete temperature measurementsInternational Journal of Computational Science and Engineering10.1504/IJCSE.2014.0645379:5/6(526-537)Online publication date: 1-Sep-2014
https://dl.acm.org/doi/10.1504/IJCSE.2014.064537

View Options

View options

Media

Figures

Other

Tables

View Issue’s Table of Contents