By 苏剑林 | January 26, 2026
From "A Brief History of Linear Attention: From Imitation, Innovation to Feedback", we can see that the parallel form of DeltaNet involves an inverse matrix of the form $(I + KK^\top \odot M^-)^{-1}$. Recently, reader @Arch123 pointed out that experiments show the elements of this inverse matrix always lie within $[-1, 1]$, and asked if this could be mathematically proven or disproven.
In this article, we will prove that this conclusion holds strictly through two different methods.
Problem Description
First, let's restate the problem accurately. Let the matrix be $K=[k_1, k_2, \dots, k_n]^\top \in \mathbb{R}^{n \times d}$, where each $k_i \in \mathbb{R}^{d \times 1}$ is a column vector with a norm not exceeding 1. $M \in \mathbb{R}^{n \times n}$ is a lower triangular mask matrix, defined as:
\begin{equation}
M_{i,j} = \begin{cases} 1, & i \ge j \\ 0, & i < j \end{cases} \label{eq:1}
\end{equation}
$I$ is the identity matrix, and $M^- = M - I$. We want to prove that:
\begin{equation}
(I + KK^\top \odot M^-)^{-1} \in [-1, 1]^{n \times n} \label{eq:2}
\end{equation}
Let $X = I + KK^\top \odot M^-$, and denote its inverse matrix as $Y$. Writing $X$ explicitly gives:
\begin{equation}
X = I + KK^\top \odot M^- = \begin{pmatrix} 1 & 0 & 0 & \cdots & 0 & 0 \\ k_2^\top k_1 & 1 & 0 & \cdots & 0 & 0 \\ k_3^\top k_1 & k_3^\top k_2 & 1 & \cdots & 0 & 0 \\ \vdots & \vdots & \vdots & \ddots & \vdots & \vdots \\ k_{n-1}^\top k_1 & k_{n-1}^\top k_2 & k_{n-1}^\top k_3 & \cdots & 1 & 0 \\ k_n^\top k_1 & k_n^\top k_2 & k_n^\top k_3 & \cdots & k_n^\top k_{n-1} & 1 \end{pmatrix} \label{eq:3}
\end{equation}
Since the norm of each $k_i$ does not exceed 1, it is obvious that every element of $X$ lies in $[-1, 1]$. We need to prove that every element of $Y = X^{-1}$ also lies in $[-1, 1]$.
Mathematical Induction
Since this inverse matrix appears in DeltaNet, it must have a clear model background. Using this background might help us prove it faster, but we'll save that for the next section. In this section, we treat it purely as a mathematical problem.
Diagonal Properties
We are dealing with a lower triangular matrix. As discussed in "Efficient Inversion Methods for 'Diagonal + Low-Rank' Triangular Matrices", this class of matrices has excellent properties. First, it is closed under matrix multiplication. Second, it satisfies the property "the diagonal of the inverse is the inverse of the diagonal." This can be generalized to block lower triangular matrices: "the diagonal blocks of the inverse are the inverse of the diagonal blocks." From this, we can conclude that $Y[:n-1, :n-1] = (X[:n-1, :n-1])^{-1}$ and $Y[1:, 1:] = (X[1:, 1:])^{-1}$, where slicing is understood as in Numpy.
This property suggests using mathematical induction on $n$. Assuming the conclusion holds for any $(n-1) \times d$ matrix $K$, for matrix $X$, we only need to consider the following block partition:
\begin{equation}
X = \begin{pmatrix} X_{:n-1, :n-1} & 0 \\ \text{row}_n & 1 \end{pmatrix} \quad \text{or} \quad X = \begin{pmatrix} 1 & 0 \\ \text{col}_1 & X_{1:, 1:} \end{pmatrix} \label{eq:4}
\end{equation}
Using $Y[:n-1, :n-1] = (X[:n-1, :n-1])^{-1}$ and $Y[1:, 1:] = (X[1:, 1:])^{-1}$, we can prove that all elements of $Y$ except for the bottom-left $Y_{n,1}$ lie in $[-1, 1]$. Thus, we only need to prove $Y_{n,1} \in [-1, 1]$.
Adjugate Matrix
To prove $Y_{n,1} \in [-1, 1]$, we consider the explicit representation of $X^{-1}$ using the adjugate matrix. Since $X$ is lower triangular with 1s on the diagonal, its determinant is 1. According to the formula for the inverse via the adjugate matrix, we have:
\begin{equation}
Y_{n,1} = (-1)^{n+1} \det \begin{pmatrix} k_2^\top k_1 & 1 & 0 & \cdots & 0 \\ k_3^\top k_1 & k_3^\top k_2 & 1 & \cdots & 0 \\ k_4^\top k_1 & k_4^\top k_2 & k_4^\top k_3 & \cdots & 0 \\ \vdots & \vdots & \vdots & \ddots & \vdots \\ k_n^\top k_1 & k_n^\top k_2 & k_n^\top k_3 & \cdots & k_n^\top k_{n-1} \end{pmatrix} \label{eq:5}
\end{equation}
This matrix is derived from a lower triangular matrix where the super-diagonal elements have become 1. Calculating this determinant is more complex, but not overly so, as the last column has only two non-zero elements. Expanding the determinant along the last column reduces it to the subtraction of two determinants of order $n-1$, which still points toward induction.
Continued Product Form
We can try calculating the first few results manually:
\begin{align}
Y_{2,1} &= -k_2^\top k_1 \nonumber \\
Y_{3,1} &= \det \begin{pmatrix} k_2^\top k_1 & 1 \\ k_3^\top k_1 & k_3^\top k_2 \end{pmatrix} = k_2^\top k_1 k_3^\top k_2 - k_3^\top k_1 = -k_3^\top (I - k_2 k_2^\top) k_1 \nonumber \\
Y_{4,1} &= - \det \begin{pmatrix} k_2^\top k_1 & 1 & 0 \\ k_3^\top k_1 & k_3^\top k_2 & 1 \\ k_4^\top k_1 & k_4^\top k_2 & k_4^\top k_3 \end{pmatrix} = -k_4^\top (I - k_3 k_3^\top) (I - k_2 k_2^\top) k_1 \label{eq:6}
\end{align}
From this, we can conjecture:
\begin{equation}
Y_{n,1} = -k_n^\top (I - k_{n-1} k_{n-1}^\top) (I - k_{n-2} k_{n-2}^\top) \dots (I - k_2 k_2^\top) k_1 \label{eq:8}
\end{equation}
I invite readers to prove this conjecture using mathematical induction. Assuming the proof is complete, considering the condition $\|k_i\| \le 1$, we have $\|(I - k_i k_i^\top)x\|^2 = \|x\|^2 + (k_i^\top x)^2 (\|k_i\|^2 - 2) \le \|x\|^2$. Therefore:
\begin{equation}
|Y_{n,1}| \le \|k_n\| \times \|(I - k_{n-1} k_{n-1}^\top) \dots (I - k_2 k_2^\top) k_1\| \le \|k_n\| \times \|k_1\| \le 1 \label{eq:9}
\end{equation}
This completes the proof that $Y_{n,1} \in [-1, 1]$, and by induction, $Y \in [-1, 1]^{n \times n}$ holds for all $n$.
Dual Computation
The previous direct proof was somewhat "brute-force." Later, members of the FLA group reminded me that we can return to the original context of the inverse matrix—DeltaNet—and perform calculations in two different ways. By comparing the results, we can derive the explicit expression of the inverse matrix, leading to a more concise proof.
Inverse Matrix Form
We know the recursive form of DeltaNet is:
\begin{equation}
S_t = S_{t-1}(I - k_t k_t^\top) + v_t k_t^\top = S_{t-1} + (v_t - S_{t-1} k_t) k_t^\top \label{eq:10}
\end{equation}
where $S_0 = 0$. Looking at the second equality, let $u_t = v_t - S_{t-1} k_t$. Then:
\begin{equation}
S_t = S_{t-1} + u_t k_t^\top = \sum_{i=1}^t u_i k_i^\top \label{eq:11}
\end{equation}
So,
\begin{equation}
u_t = v_t - S_{t-1} k_t = v_t - \left(\sum_{i=1}^{t-1} u_i k_i^\top \right) k_t = v_t - \sum_{i=1}^{t-1} u_i (k_i^\top k_t) \label{eq:12}
\end{equation}
Defining $U = [u_1, u_2, \dots, u_n]^\top$ and $V = [v_1, v_2, \dots, v_n]^\top$, the above equation is equivalent to $U = V - (KK^\top \odot M^-) U$, or $U = (I + KK^\top \odot M^-)^{-1} V$. This is the origin of the inverse matrix $(I + KK^\top \odot M^-)^{-1}$, which is key to parallel computing in DeltaNet and subsequent works like GDN and KDA.
Recursive Expansion
Now looking at the first equality, expanding it recursively gives:
\begin{equation}
S_t = v_1 k_1^\top H_{2 \to t} + \dots + v_{t-1} k_{t-1}^\top H_{t \to t} + v_t k_t^\top = \sum_{i=1}^t v_i k_i^\top H_{i+1 \to t} \label{eq:13}
\end{equation}
Here $H_{r \to t} \triangleq (I - k_r k_r^\top) (I - k_{r+1} k_{r+1}^\top) \dots (I - k_t k_t^\top)$, with the convention $H_{t+1 \to t} = I$. Substituting this into the definition of $u_t$:
\begin{equation}
u_t = v_t - S_{t-1} k_t = v_t - \sum_{i=1}^{t-1} v_i k_i^\top H_{i+1 \to t-1} k_t \label{eq:14}
\end{equation}
Comparison of Results
Again letting $Y = (I + KK^\top \odot M^-)^{-1}$, then $U = YV$ implies $u_t = \sum_{i=1}^n Y_{t,i} v_i$. By comparing with the equation above, we can identify:
\begin{equation}
Y_{t,i} = \begin{cases} 0, & i > t \\ 1, & i = t \\ -k_i^\top H_{i+1 \to t-1} k_t, & i < t \end{cases} \label{eq:15}
\end{equation}
Therefore, we only need to prove $|k_i^\top H_{i+1 \to t-1} k_t| \le 1$. Given $\|k_i\| \le 1$ and the property $\|(I - k_j k_j^\top)x\| \le \|x\|$ established in the first proof, it follows immediately that:
\begin{equation}
|k_i^\top H_{i+1 \to t-1} k_t| \le \|k_i^\top\| \times \|H_{i+1 \to t-1} k_t\| \le \|k_i^\top\| \times \|k_t\| \le 1 \label{eq:16}
\end{equation}
Thus, the theorem is proven.
Summary
This article provides two proofs for the boundedness of the core inverse matrix elements in DeltaNet.