Scalable discrete adjoint for the transient Savage–Hutter equation on the primal consistent adaptive grid

Computing the adjoint for a transient non-linear hyperbolic PDE, like the Savage–Hutter equation for modeling dry granular mass flows, is a challenge with often unaffordable computing and memory costs. The challenges are further compounded for an adaptive mesh refinement (AMR) parallel code. For legacy codes, addressing these difficulties are constrained with the original design and data structures which refactoring those to be well-suited for...