Reversible random number generation for adjoint Monte Carlo simulation of the heat equation

arXiv2302.02778MaRDI QIDQ6425565FDOQ6425565

Giovanni Samaey, Author name not available (Why is that?), Frédéric Blondeel, Adam J. Lee, Andreas van Barel, Emil Løvbak

Publication date: 6 February 2023

Abstract: In PDE-constrained optimization, one aims to find design parameters that minimize some objective, subject to the satisfaction of a partial differential equation. A major challenges is computing gradients of the objective to the design parameters, as applying the chain rule requires computing the Jacobian of the design parameters to the PDE's state. The adjoint method avoids this Jacobian by computing partial derivatives of a Lagrangian. Evaluating these derivatives requires the solution of a second PDE with the adjoint differential operator to the constraint, resulting in a backwards-in-time simulation. Particle-based Monte Carlo solvers are often used to compute the solution to high-dimensional PDEs. However, such solvers have the drawback of introducing noise to the computed results, thus requiring stochastic optimization methods. To guarantee convergence in this setting, both the constraint and adjoint Monte Carlo simulations should simulate the same particle trajectories. For large simulations, storing full paths from the constraint equation for re-use in the adjoint equation becomes infeasible due to memory limitations. In this paper, we provide a reversible extension to the family of permuted congruential pseudorandom number generators (PCG). We then use such a generator to recompute these time-reversed paths for the heat equation, avoiding these memory issues.

Has companion code repository: https://gitlab.kuleuven.be/numa/public/mcqmc_2022_rrng

This page was built for publication: Reversible random number generation for adjoint Monte Carlo simulation of the heat equation

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q6425565)