@inproceedings{
kr{\"a}mer2024gradients,
title={Gradients of Functions of Large Matrices},
author={Nicholas Kr{\"a}mer and Pablo Moreno-Mu{\~n}oz and Hrittik Roy and S{\o}ren Hauberg},
booktitle={The Thirty-eighth Annual Conference on Neural Information Processing Systems},
year={2024},
url={https://openreview.net/forum?id=RL4FXrGcTw}
}
Tuning scientific and probabilistic machine learning models - for example, partial differential equations, Gaussian processes, or Bayesian neural networks - often relies on evaluating functions of matrices whose size grows with the data set or the number of parameters.
While the state-of-the-art for _evaluating_ these quantities is almost always based on Lanczos and Arnoldi iterations, the present work is the first to explain how to _differentiate_ these workhorses of numerical linear algebra efficiently.
To get there, we derive previously unknown adjoint systems for Lanczos and Arnoldi iterations, implement them in JAX, and show that the resulting code can compete with Diffrax when it comes to differentiating PDEs, GPyTorch for selecting Gaussian process models and beats standard factorisation methods for calibrating Bayesian neural networks.
All this is achieved without any problem-specific code optimisation.
Find the code at [link redacted] and install the library with *pip install [redacted]*.
Gradients of Functions of Large Matrices
Nicholas Krämer Pablo Moreno-Muñoz Hrittik Roy Søren Hauberg