Sherman–Morrison formula


In linear algebra, the Sherman–Morrison formula, named after Jack Sherman and Winifred J. Morrison, computes the inverse of a "rank-1 update" to a matrix whose inverse has previously been computed. That is, given an invertible matrix and the outer product of vectors and the formula cheaply computes an updated matrix inverse
The Sherman–Morrison formula is a special case of the Woodbury formula. Though named after Sherman and Morrison, it appeared already in earlier publications.

Statement

Suppose is an invertible square matrix and are column vectors. Then is invertible if and only if. In this case,
Here, is the outer product of two vectors and. The general form shown here is the one published by Bartlett.

Proof

To prove that the backward direction is invertible is true, we verify the properties of the inverse. A matrix is the inverse of a matrix if and only if.
We first verify that the right hand side satisfies.
To end the proof of this direction, we need to show that in a similar way as above:
Reciprocally, if, then via the matrix determinant lemma,, so is not invertible.

Application

If the inverse of is already known, the formula provides a numerically cheap way to compute the inverse of corrected by the matrix . The computation is relatively cheap because the inverse of does not have to be computed from scratch, but can be computed by correcting .
Using unit columns for or, individual columns or rows of may be manipulated and a correspondingly updated inverse computed relatively cheaply in this way. In the general case, where is an -by- matrix and and are arbitrary vectors of dimension, the whole matrix is updated and the computation takes scalar multiplications. If is a unit column, the computation takes only scalar multiplications. The same goes if is a unit column. If both and are unit columns, the computation takes only scalar multiplications.
This formula also has application in theoretical physics. Namely, in quantum field theory, one uses this formula to calculate the propagator of a spin-1 field. The inverse propagator has the form. One uses the Sherman–Morrison formula to calculate the inverse of the inverse propagator—or simply the propagator—which is needed to perform any perturbative calculation involving the spin-1 field.
One of the issues with the formula is that little is known about its numerical stability. There are no published results concerning its error bounds. Anecdotal evidence suggests that the Woodbury matrix identity may diverge even for seemingly benign examples.

Alternative verification

Following is an alternate verification of the Sherman–Morrison formula using the easily verifiable identity
Let
then
Substituting gives

Generalization ([Woodbury matrix identity])

Given a square invertible matrix, an matrix, and a matrix, let be an matrix such that. Then, assuming is invertible, we have