English

Variational Inference Using a Differentiable Multigrid Linear Solver

Numerical Analysis 2026-08-01 v1 Medical Physics Computation

Abstract

Gradient-based Bayesian inference methods require efficient access to Jacobian and adjoint-Jacobian operators of high-dimensional forward models. While multigrid solvers provide near-optimal complexity for elliptic partial differential equations, they are rarely available in forms compatible with automatic differentiation (AD). We develop a differentiable multigrid solver for steady-state diffusion-absorption problems and derive its adjoint operations analytically through the full multigrid hierarchy. The resulting solver, DMGS, is implemented in C++ and interfaced with JAX to provide efficient Jacobian-vector and vector-Jacobian products for variational inference in the NIFTy framework. We validate the approach on a 3D inverse problem involving diffuse radiative transfer in tissue, reconstructing an effective radiative source from Monte Carlo-simulated data. The reconstruction reproduces the data at a reduced chi-squared of 1.1 and generalizes to 32 independent validation datasets. Benchmarks against a JAX-native multigrid implementation show comparable runtimes and consistently lower peak memory for the hand-derived adjoint, with modest reverse-mode overhead. These results establish differentiable multigrid solvers as practical building blocks for variational inference in PDE-constrained problems.

Cite

@article{arxiv.2608.00760,
  title  = {Variational Inference Using a Differentiable Multigrid Linear Solver},
  author = {Andrés Ramírez and Philipp Haim and David Föger and Torsten A. Enßlin and Philipp Frank and Philipp Gschwandtner and Dominik Jüstel and Philipp Mertsch and Vo Hong Minh Phan and Laurin Söding and Hanieh Zandinejad and Ralf Kissmann},
  journal= {arXiv preprint arXiv:2608.00760},
  year   = {2026}
}

Comments

18 pages, 13 figures, 1 table