Source-to-source adjoint Algorithmic Differentiation of an ice sheet model written in C
Algorithmic Differentiation (AD) has become a powerful tool to improve our understanding of the Earth System, because it can generate adjoint code which permits efficient calculation of gradients that are essential to sensitivity studies, inverse problems, parameter estimation and data assimilation....
Main Authors: | , |
---|---|
Format: | Article in Journal/Newspaper |
Language: | unknown |
Published: |
eScholarship, University of California
2018
|
Subjects: | |
Online Access: | https://escholarship.org/uc/item/9z43d6wv |
id |
ftcdlib:oai:escholarship.org/ark:/13030/qt9z43d6wv |
---|---|
record_format |
openpolar |
spelling |
ftcdlib:oai:escholarship.org/ark:/13030/qt9z43d6wv 2023-05-15T14:04:02+02:00 Source-to-source adjoint Algorithmic Differentiation of an ice sheet model written in C Hascoët, L Morlighem, M 829 - 843 2018-11-02 application/pdf https://escholarship.org/uc/item/9z43d6wv unknown eScholarship, University of California qt9z43d6wv https://escholarship.org/uc/item/9z43d6wv public Optimization Methods and Software, vol 33, iss 4-6 Ice sheet model ISSM Algorithmic Differentiation adjoint methods dynamic memory Tapenade Operations Research Applied Mathematics Numerical and Computational Mathematics Computation Theory and Mathematics article 2018 ftcdlib 2021-04-16T07:11:37Z Algorithmic Differentiation (AD) has become a powerful tool to improve our understanding of the Earth System, because it can generate adjoint code which permits efficient calculation of gradients that are essential to sensitivity studies, inverse problems, parameter estimation and data assimilation. Most source-to-source transformation tools, however, have been designed for FORTRAN and support for C remains limited. Here we use the Adjoinable Land Ice Flow model (ALIF), a C clone of the C++ Ice Sheet System Model (ISSM) and employ source-to-source AD to produce its adjoint code. We present the first running source-to-source adjoint of ALIF, and its application to basal drag inversion under Pine Island Glacier, West Antarctica. ALIF brought several challenges to AD tool development, such as the correct treatment of the context code, which does not compute the differentiable function, but controls this computation through the setup of data structures, including possible aliasing, as well as data-flow reversal in the presence of pointers and dynamic memory, which are ubiquitous in codes such as ISSM and ALIF. We present the strategies we have developed to overcome these challenges. Article in Journal/Newspaper Antarc* Antarctica Ice Sheet Pine Island Pine Island Glacier West Antarctica University of California: eScholarship West Antarctica Pine Island Glacier ENVELOPE(-101.000,-101.000,-75.000,-75.000) |
institution |
Open Polar |
collection |
University of California: eScholarship |
op_collection_id |
ftcdlib |
language |
unknown |
topic |
Ice sheet model ISSM Algorithmic Differentiation adjoint methods dynamic memory Tapenade Operations Research Applied Mathematics Numerical and Computational Mathematics Computation Theory and Mathematics |
spellingShingle |
Ice sheet model ISSM Algorithmic Differentiation adjoint methods dynamic memory Tapenade Operations Research Applied Mathematics Numerical and Computational Mathematics Computation Theory and Mathematics Hascoët, L Morlighem, M Source-to-source adjoint Algorithmic Differentiation of an ice sheet model written in C |
topic_facet |
Ice sheet model ISSM Algorithmic Differentiation adjoint methods dynamic memory Tapenade Operations Research Applied Mathematics Numerical and Computational Mathematics Computation Theory and Mathematics |
description |
Algorithmic Differentiation (AD) has become a powerful tool to improve our understanding of the Earth System, because it can generate adjoint code which permits efficient calculation of gradients that are essential to sensitivity studies, inverse problems, parameter estimation and data assimilation. Most source-to-source transformation tools, however, have been designed for FORTRAN and support for C remains limited. Here we use the Adjoinable Land Ice Flow model (ALIF), a C clone of the C++ Ice Sheet System Model (ISSM) and employ source-to-source AD to produce its adjoint code. We present the first running source-to-source adjoint of ALIF, and its application to basal drag inversion under Pine Island Glacier, West Antarctica. ALIF brought several challenges to AD tool development, such as the correct treatment of the context code, which does not compute the differentiable function, but controls this computation through the setup of data structures, including possible aliasing, as well as data-flow reversal in the presence of pointers and dynamic memory, which are ubiquitous in codes such as ISSM and ALIF. We present the strategies we have developed to overcome these challenges. |
format |
Article in Journal/Newspaper |
author |
Hascoët, L Morlighem, M |
author_facet |
Hascoët, L Morlighem, M |
author_sort |
Hascoët, L |
title |
Source-to-source adjoint Algorithmic Differentiation of an ice sheet model written in C |
title_short |
Source-to-source adjoint Algorithmic Differentiation of an ice sheet model written in C |
title_full |
Source-to-source adjoint Algorithmic Differentiation of an ice sheet model written in C |
title_fullStr |
Source-to-source adjoint Algorithmic Differentiation of an ice sheet model written in C |
title_full_unstemmed |
Source-to-source adjoint Algorithmic Differentiation of an ice sheet model written in C |
title_sort |
source-to-source adjoint algorithmic differentiation of an ice sheet model written in c |
publisher |
eScholarship, University of California |
publishDate |
2018 |
url |
https://escholarship.org/uc/item/9z43d6wv |
op_coverage |
829 - 843 |
long_lat |
ENVELOPE(-101.000,-101.000,-75.000,-75.000) |
geographic |
West Antarctica Pine Island Glacier |
geographic_facet |
West Antarctica Pine Island Glacier |
genre |
Antarc* Antarctica Ice Sheet Pine Island Pine Island Glacier West Antarctica |
genre_facet |
Antarc* Antarctica Ice Sheet Pine Island Pine Island Glacier West Antarctica |
op_source |
Optimization Methods and Software, vol 33, iss 4-6 |
op_relation |
qt9z43d6wv https://escholarship.org/uc/item/9z43d6wv |
op_rights |
public |
_version_ |
1766274971687976960 |