\usepackageamsmath

Block Multiplication of Matrices

(We needed this to derive the conditional distribution of a multivariate Gaussian).

Consider a matrix product AB. Partition the two outer dimensions (i.e. the rows of A and the columns of B) and the one inner dimension (i.e. the columns of A and the rows of B) arbitrarily. This defines a “block decomposition” of the product AB and of the factors A,B such that the blocks of AB are related to the blocks of A and B via the familiar formula for components of the product, i.e.

(AB)m,n=sAm,sBs,n.

Pictorially, we have the following:

blockmultnmatrices

Arithmetically, this is easy to prove by considering the formula above for the components of the product. The partitioning of the outer dimensions comes for free, while the partitioning of the inner dimension just corresponds to partitioning the summation:

(AB)m,n=sAm,sBs,n=isissi+1Am,sBs,n.

Zooming out to a categorical level, we can see that there is nothing peculiar about this situation. If, in an additive category, we have three objects X,Y,Z with biproduct decompositions, and a chain of morphisms:

XφBYφAZ

then this “block decomposition of matrices” finds expression as a formula in End(X,Z) using the injection and projection morphisms associated with each biproduct factor.

2 Replies to “Block Multiplication of Matrices”

Leave a Reply

Your email address will not be published. Required fields are marked *