I’m working on a project these days that uses four different kinds of matrix product, which made me wonder if there’s another kind of product out there that I could find some use for.

In the process of looking around for other matrix products, I ran across the **Kronecker sum**. I’ve seen Kronecker *products* many times, but I’d never heard of Kronecker *sums.*

The Kronecker sum is defined in terms of the Kronecker product, so if you’re not familiar with the latter, you can find a definition and examples here. Essentially, you multiply each scalar element of the first matrix by the second matrix *as a block matrix*.

The Kronecker product of an *m* × *n* matrix *A* and a *p* × *q* matrix *B* is a *mp* × *nq* matrix *K* = *A *⊗* B*. You could think of *K* as an *m* × *n* matrix whose entries are *p* × *q* blocks.

So, what is the Kronecker sum? It is defined for two square matrices, an *n* × *n* matrix *A* and an *m* × *m* matrix *B*. The sizes of the two matrices need not match, but the matrices do need to be square. The Kronecker sum of *A* and *B* is

*A* ⊕ *B* = *A* ⊗ *I*_{m} + *I*_{n} ⊗ *B*

where *I*_{m} and *I*_{n} are identity matrices of size *m* and *n* respectively.

Does this make sense dimensionally? The left side of the (ordinary) matrix addition is *nm* × *nm*, and so is the right side, so the addition makes sense.

However, the Kronecker sum is not commutative, and usually things called “sums” are commutative. Products are not always commutative, but it goes against convention to call a non-commutative operation a sum. Still, the Kronecker sum is kinda like a sum, so it’s not a bad name.

I don’t have any application in mind (*yet*) for the Kronecker sum, but presumably it was defined for a good reason, and maybe I’ll run an application, maybe even on the project alluded to at the beginning.

There are several identities involving Kronecker sums, and here’s one I found interesting:

exp( *A* ) ⊗ exp( *B* ) = exp( *A *⊕ *B* ).

If you haven’t seen the exponential of a matrix before, basically you stick your matrix into the power series for the exponential function.

## Examples

First, let’s define a couple matrices *A* and *B*.

We can compute the Kronecker sums

*S* = *A* ⊕ *B*

and

*T* = *B* ⊕ *A*

with Mathematica to show they are different.

A = {{1, 2}, {3, 4}} B = {{1, 0, 1}, {1, 2, 0}, {2, 0, 3}} S = KroneckerProduct[A, IdentityMatrix[3]] + KroneckerProduct[IdentityMatrix[2], B] T = KroneckerProduct[B, IdentityMatrix[2]] + KroneckerProduct[IdentityMatrix[3], A]

This shows

and so the two matrices are not equal.

We can compute the matrix exponentials of *A* and *B* with the Mathematica function `MatrixExp`

to see that

(I actually used `MatrixExp[N[A]]`

and similarly for *B* so Mathematica would compute the exponentials numerically rather than symbolically. The latter takes forever and it’s hard to read the result.)

Now we have

and so the two matrices are equal.

Even though the identity

exp( *A* ) ⊗ exp( *B* ) = exp( *A *⊕ *B* )

may look symmetrical, it’s not. The matrices on the left do not commute in general. And not only are *A* ⊕ *B* and *B* ⊕ *A* different in general, their exponentials are also different. For example

Kronecker sums can be useful for studying matrix equations (Riccati, Lyapunov, Sylvester) and matrix formulations of difference operators on tensor product grids. As for Kronecker products, there are nice theorems about the eigenvalues and eigenvectors of Kronecker sums. These applications and properties as well as the non-commutativity are discussed in a post by Nick Higham [1].

[1] https://nhigham.com/2020/08/25/what-is-the-kronecker-product/

Another interesting matrix product, related to the Kronecker product, is the Khatri-Rao product [1].

[1] https://en.wikipedia.org/wiki/Khatri%E2%80%93Rao_product

This is almost spooky: I hadn’t heard of the Khatri-Rao product before reading your comment, and this morning I found a use for it.

I was looking for a different Kronecker like sum, not sure what the correct name is, I would have also called it a (full?) Kronecker sum for any two matrices A (m x n) and B (o x p) yielding a matrix C (mo x np) via

A (+) B = kron(A, ones(o,p)) .+ kron(ones(m,n), B)

Do you know the name for this sum?