(We needed this to derive the conditional distribution of a multivariate Gaussian).
Consider a matrix product . Partition the two outer dimensions (i.e. the rows of
and the columns of
) and the one inner dimension (i.e. the columns of
and the rows of
) arbitrarily. This defines a “block decomposition” of the product
and of the factors
such that the blocks of
are related to the blocks of
and
via the familiar formula for components of the product, i.e.
.
Pictorially, we have the following:
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:
.
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 with biproduct decompositions, and a chain of morphisms:
then this “block decomposition of matrices” finds expression as a formula in using the injection and projection morphisms associated with each biproduct factor.
One Reply to “Block Multiplication of Matrices”