Matrix calculus
July 9, 2018 — July 24, 2023
We can generalize high school calculus, which is about scalar functions of a scalar argument, in various ways to handle matrix-valued functions or matrix-valued arguments and still look tidy. One could generalize this further, by going to full tensor calculus. But it happens that specifically matrix/vector operations are at a useful point of complexity for lots of algorithms. (I usually want this for higher-order gradient descent.)
I mention two convenient and popular formalisms for lazy matrix calculus. In practice, a mix of each is often useful.
1 Matrix differentials
🏗 I need to return to this and tidy it up with some examples.
A special case of tensor calculus; where the rank of the argument and value of the function is not too big. In this setting, we often get to cheat and use some handy shortcuts. Fun pain point: agreeing upon the layout of derivatives, numerator vs denominator.
If our problem is nice, this often gets us a low-fuss, compact, tidy solution even for some surprising cases where it seems that more general tensors would be more natural — for which, see below.
The Matrix Calculus You Need For Deep Learning (Parr and Howard 2018)
The rough-and-ready matrix differential notation is occasionally confusing, but it has a functional analysis interpretation.
Many quick recipes: the Matrix Cookbook (Petersen and Pedersen 2012)
Brookes’ Matrix Reference Manual
More expository but not as broad, Old and new matrix algebra useful for statistics (Minka 2000)
autodiff-focussed: Collected Matrix Derivative Results for Forward and Reverse Mode Algorithmic Differentiation (Giles 2008)
Alan Edelman’s lectures are v/ pedagogic on matrices and calculus (and go into useful random matrix theory)
2 Automating matrix calculus
- Use Mathematica+NCAlgebra to find matrix differentials.
- Soeren Laue, Matthias Mitterreiter, Joachim Giesen and Jens K. Mueller (Soeren Laue, Mitterreiter, and Giesen 2018) have a website MatrixCalculus.org which uses Ricci calculus to generate matrix calculus formulae. Bonus feature: generates both python and LaTeX code.
3 Indexed tensor calculus
Filed under multilinear algebra.