---
title: 'De Groote Equivalence: Tensors & Concurrency'
url: https://www.emergentmind.com/topics/de-groote-equivalence
type: topic
---

# De Groote Equivalence: Tensors & Concurrency

Searching arXiv for papers relevant to “de Groote equivalence,” covering both matrix multiplication tensor equivalence and Groote-style behavioural equivalences.
de Groote equivalence denotes, in modern tensor language, the equivalence relation on polyadic decompositions of matrix multiplication tensors generated by permutation transformations, scaling transformations, and trace transformations; two bilinear algorithms are regarded as the same when they differ only by these invariance operations. The available literature also uses closely related Groote-associated terminology in concurrency theory for divergence-blind stuttering equivalence on Kripke structures and branching bisimulation on labelled transition systems, together with the partition-refinement procedures used to decide them. The common theme is the identification of representations that are “essentially equivalent,” but the tensor-algebraic and behavioural settings are mathematically distinct [1902.03950][1603.05789].

## 1. Terminological scope

In the tensor-algebraic setting, de Groote equivalence is the equivalence relation attached to bilinear algorithms for matrix multiplication. Berger, Gesmundo, and de la Cruz make this explicit for matrix multiplication tensors and treat it as an orbit relation under natural symmetries of a polyadic decomposition [1902.03950].

In the concurrency-theoretic setting, the phrase can refer more loosely to Groote-style behavioural equivalences and their decision procedures. The relevant notions are divergence-blind stuttering equivalence on Kripke structures and branching bisimulation on labelled transition systems, as defined and algorithmically treated in work by Groote–Vaandrager and Groote–Wijs. In that context, the emphasis is not on a named algebraic orbit relation, but on standard weak behavioural equivalences and on how to compute them efficiently [1603.05789].

This suggests a useful disambiguation. In tensor theory, de Groote equivalence is an explicit equivalence relation on decompositions. In model checking and concurrency, Groote’s name is attached primarily to equivalence notions and refinement algorithms rather than to a single standardized term.

## 2. Matrix multiplication tensors and the equivalence relation

Let \(F\) be a field, and let
\[
\Phi_{m,p,n} \in \mathrm{Bil}(F^{m\times p}, F^{p\times n}; F^{m\times n})
\]
denote multiplication of an \(m\times p\) matrix by a \(p\times n\) matrix. Berger et al. identify \(\Phi_{m,p,n}\) with the matrix multiplication tensor of format \((m,p,n)\) [1902.03950].

A rank-1 bilinear map has the form
\[
\Phi(u,v) = f(u)\, g(v)\, w
\]
with \(f\in U^*\), \(g\in V^*\), and \(w\in W\). An \(F\)-term polyadic decomposition of a bilinear map is
\[
\Phi(u,v) = \sum_{r=1}^F f_r(u)\, g_r(v)\, w_r.
\]
For matrix multiplication, the dual functionals can be represented by matrices \(A_r\) and \(B_r\) through the trace pairing, and the output factors by matrices \(C_r\), yielding
\[
\Phi_{m,p,n}(A,B)
= \sum_{r=1}^F \operatorname{trace}(A_r A)\,\operatorname{trace}(B_r B)\,C_r.
\]
Each term corresponds to one active multiplication \(f_r(A)\cdot g_r(B)\); the bilinear rank is the smallest \(F\) admitting such a decomposition, hence the minimal number of active multiplications in any bilinear algorithm [1902.03950].

The defining invariance transformations are the following.

**Permutation transformations:** for \(\sigma\in S_F\),
\[
A'_r = A_{\sigma(r)},\qquad
B'_r = B_{\sigma(r)},\qquad
C'_r = C_{\sigma(r)}.
\]

**Scaling transformations:** for scalars \(\lambda_r,\mu_r,\nu_r\) with
\[
\lambda_r\mu_r\nu_r = 1,
\]
define
\[
A'_r = \lambda_r A_r,\qquad
B'_r = \mu_r B_r,\qquad
C'_r = \nu_r C_r.
\]

**Trace transformations:** for
\[
P \in \mathrm{GL}(m), \quad Q \in \mathrm{GL}(p), \quad R \in \mathrm{GL}(n),
\]
define
\[
A'_r = Q^{-1} A_r P,\qquad
B'_r = R^{-1} B_r Q,\qquad
C'_r = P^{-1} C_r R.
\]

Two \(F\)-term polyadic decompositions of \(\Phi_{m,p,n}\) are equivalent if one can be transformed into the other by a finite sequence of these operations. In tensor language, this is the natural orbit relation induced by the symmetric group on the summands, termwise rescaling, and compatible base changes in the three tensor factors [1902.03950].

## 3. Structural meaning and classical uniqueness results

The significance of de Groote equivalence is that it separates essential algorithmic content from presentation-dependent choices. Reordering multiplication steps does not change an algorithm. Termwise rescaling changes the normalization of intermediate linear forms but not the computed bilinear map. Trace transformations implement simultaneous coordinate changes compatible with matrix multiplication itself [1902.03950].

For two decompositions
\[
(A_{[F]},B_{[F]},C_{[F]})
\quad\text{and}\quad
(A'_{[F]},B'_{[F]},C'_{[F]}),
\]
equivalence is characterized by the existence of a permutation \(\sigma\in S_F\), invertible matrices \(P,Q,R\), and scalars \(\lambda_r,\mu_r,\nu_r\) satisfying \(\lambda_r\mu_r\nu_r=1\) such that, for all \(r\),
\[
\lambda_r A'_{\sigma(r)} = Q^{-1} A_r P,\qquad
\mu_r B'_{\sigma(r)} = R^{-1} B_r Q,\qquad
\nu_r C'_{\sigma(r)} = P^{-1} C_r R.
\]
This is the concrete form of de Groote equivalence used algorithmically in the tensor literature [1902.03950].

Its classical importance is illustrated by the \(2\times 2\) case. De Groote proved that for multiplication of \(2\times 2\) matrices with \(7\) active multiplications, all algorithms are essentially equivalent to Strassen’s algorithm. In the language above, optimal rank-\(7\) decompositions of \(\Phi_{2,2,2}\) form a single de Groote equivalence class [1902.03950].

The situation changes for larger formats. Berger et al. report that for multiplication of larger matrices, including \(2\times 3\) by \(3\times 2\) and \(3\times 3\) by \(3\times 3\), two decompositions are very likely to be essentially different. This aligns with earlier work showing inequivalent \(23\)-multiplication algorithms for \(3\times 3\) multiplication [1902.03950].

## 4. Decision procedures for tensor decompositions

Berger et al. give an algorithm for deciding whether two polyadic decompositions of a given matrix multiplication tensor are equivalent. The naive approach would enumerate all \(F!\) permutations and solve the nonlinear system defined by the equivalence conditions; the paper replaces this by a combination of linearization and permutation pruning [1902.03950].

For a fixed permutation \(\sigma\), the first step is to test whether the decompositions are \((\text{scaling}+\text{trace})\)-equivalent. Writing
\[
U = [\operatorname{vec}(A_1),\dots,\operatorname{vec}(A_F)] \in \mathbb{R}^{pm\times F},
\]
and similarly \(V\) and \(W\), the transformation
\[
A_r \mapsto Q^{-1}A_r P
\]
becomes, after vectorization,
\[
\operatorname{vec}(Q^{-1}A_r P) = (P^\top \otimes Q^{-1})\,\operatorname{vec}(A_r).
\]
Hence the \(A\)-factor equations can be written as
\[
(P^\top \otimes Q^{-1})\, U = U' \,\mathrm{diag}(\lambda_1,\dots,\lambda_F),
\]
which is linear in the entries of \(X:=P^\top\otimes Q^{-1}\) and the \(\lambda_r\). Analogous linear systems are solved for the \(B\)- and \(C\)-factors. The paper uses the clustering number \(cl(A)\) to control the solution space dimension; under the stated assumption that at least one of \(U,V,W\) and its primed counterpart has clustering number \(1\), the relevant solution space is one-dimensional up to scalar [1902.03950].

The permutation search is pruned through simultaneous similarity invariants. For each \(r\), define
\[
M_r := A_r B_r C_r,\qquad M'_r := A'_r B'_r C'_r.
\]
If the decompositions are equivalent, then after the correct permutation,
\[
M'_r = P^{-1}M_{\sigma(r)}P.
\]
Thus the two ordered families must be simultaneously similar. The recursive search builds a partial permutation and rejects branches by comparing eigenvalues of random linear combinations of the corresponding \(M_r\) and \(M'_r\); if the spectra differ, the branch cannot extend to a valid equivalence [1902.03950].

The worst-case complexity is bounded by
\[
O(F! \cdot [F\max\{mp,pn,nm\}]^3),
\]
but the reported empirical behavior is substantially better. In the experiments, the recursion depth never exceeded \(9\) even for \(F=23\), the expensive scaling-and-trace check was typically called only once per pair of decompositions, and the average time to test equivalence of two decompositions was at most about \(3\times 10^{-2}\) seconds in all cases tested [1902.03950].

## 5. Equivalence classes and discretizability

Berger et al. use the decision algorithm to study how many essentially different decompositions occur for several small matrix multiplication tensors. They generated \(N_s=10{,}000\) decompositions for each case and tested \(10{,}000\) random pairs for equivalence [1902.03950].

| Format \((m,p,n)\) | \(F\) | Equivalent random pairs |
|---|---:|---:|
| \((1,2,1)\) | 2 | 100% |
| \((2,1,2)\) | 4 | 0% |
| \((2,2,2)\) | 7 | 100% |
| \((2,3,2)\) | 11 | 0% |
| \((3,2,3)\) | 15 | 0% |
| \((3,3,3)\) | 23 | 0% |

These data confirm the exceptional status of the \(2\times 2\) rank-\(7\) case and suggest that for larger tensors, randomly obtained decompositions almost surely lie in different de Groote equivalence classes [1902.03950].

The same paper also studies whether a decomposition is equivalent to one with coefficients in a discrete set. A decomposition is called **discrete** if there exists \(q\in\mathbb{R}\) such that all entries of all factor matrices lie in \(q\mathbb{Z}\). It is **discretizable** if it is equivalent, via scaling and trace transformations, to some discrete decomposition. Strassen’s \(7\)-term algorithm and Laderman’s \(23\)-term algorithm are cited as discrete with entries in \(\{0,\pm1\}\) [1902.03950].

A necessary criterion for discretizability is obtained from the matrices
\[
M_r = A_r B_r C_r.
\]
If a decomposition is discretizable with parameter \(q\), then for every integer tuple \((\beta_1,\dots,\beta_F)\), the characteristic polynomial
\[
p(t;\beta_1,\dots,\beta_F)
=
\det\left(tI-\frac{1}{q^3}\sum_{r=1}^F \beta_r M_r\right)
\]
must have integer coefficients. Failure of this condition proves that the decomposition is not equivalent to any discrete decomposition with entries in \(q\mathbb{Z}\) [1902.03950].

For \(q=\tfrac12\), the reported experiments show that \(100\%\) of decompositions for \((1,2,1)\) and \((2,2,2)\) satisfy the criterion, whereas most decompositions for \((2,1,2)\), \((2,3,2)\), \((3,2,3)\), and \((3,3,3)\) do not. A plausible implication is that discrete or few-valued decompositions, although prominent in the constructive literature, are not representative of generic decompositions in larger formats [1902.03950].

## 6. Groote-style behavioural equivalences in concurrency theory

In concurrency theory, the Groote-associated equivalence notions discussed in the literature are divergence-blind stuttering equivalence on Kripke structures and branching bisimulation on labelled transition systems. A finite Kripke structure is written
\[
K = (S,\rightarrow,L),
\]
where \(S\) is finite, \(\rightarrow\subseteq S\times S\) is total, and \(L:S\to 2^{AP}\) assigns atomic propositions. The version of stuttering equivalence used in the algorithmic work is the divergence-blind variant associated with Browne–Clarke–Grünwald and with Groote–Vaandrager; it coincides with satisfaction of the same \(\mathrm{CTL}^*_{\neg X}\) formulas [1603.05789].

The partition-refinement framework distinguishes inert transitions, which remain inside the same current block, from non-inert transitions, which leave the block. Writing \(\mathit{out}_\tau(s)\) for the inert outgoing transitions of \(s\), a bottom state of a block \(B\) is characterized by
\[
s\in B \text{ is bottom} \iff |\mathit{out}_\tau(s)|=0.
\]
The global invariant maintained by the algorithm is stability with respect to constellations: all states in a block can reach the same constellations through a weak transition [1603.05789].

Groote–Wijs claimed that stuttering equivalence, and via standard encoding also branching bisimulation, could be computed in
\[
\mathcal{O}(m\log n)\text{ time and }\mathcal{O}(m)\text{ space},
\]
improving on the earlier \(\mathcal{O}(mn)\) algorithm of Groote–Vaandrager. The high-level method uses constellations, splitter blocks \(SpB\subset SpC\) with \(|SpB|\le \tfrac12|SpC|\), and a red/blue refinement routine `TrySplit` that attempts to process only the smaller half of each split [1603.05789].

The 2016 analysis shows, however, that the original Groote–Wijs description does not always meet the claimed complexity. Two counterexamples establish \(\Omega(md)\) behavior, where \(d\) is maximal out-degree:
\[
d := \max_{s\in S} |\mathit{out}(s)|.
\]
One counterexample targets `TrySplit'`, where full scans of outgoing transitions are triggered repeatedly; the other targets `PostprocessNewBottom`, where large lists \(S_C\) are repeatedly walked after splits. A third example shows that even a first-level fix to bottom-state postprocessing is not sufficient unless newly discovered bottom states are handled with additional care [1603.05789].

The corrective modifications preserve the underlying partition-refinement structure but change the cost accounting. In `TrySplit'`, the expensive test of whether a state still has a transition to the relevant constellation is deferred until the point at which its inert successors have all become blue, so that each such scan is charged only once within the amortized budget. In `PostprocessNewBottom`, work is localized to constellations reachable from the newly created block \(NewB\), rather than by globally scanning all \(S_C\) lists for the parent block. With these fixes, the paper concludes that stuttering equivalence on Kripke structures can indeed be computed in \(\mathcal{O}(m\log n)\) time and \(\mathcal{O}(m)\) space, and the same bound then applies to branching bisimulation through the standard encoding [1603.05789].

This concurrency-theoretic usage is not identical to tensor-theoretic de Groote equivalence, but it exhibits an analogous concern: identifying when two states or processes should count as behaviourally indistinguishable, and ensuring that the quotienting algorithm respects a tight complexity bound.

## 7. Conceptual synthesis

Across both settings, de Groote’s name is attached to mathematically natural quotient constructions. In matrix multiplication, the quotient is by explicit symmetries of a polyadic decomposition: permutation of rank-1 terms, termwise rescaling, and compatible base changes. The resulting equivalence classes capture when two bilinear algorithms are the same up to invariances that do not alter the underlying tensor [1902.03950].

In concurrency theory, Groote-associated equivalences identify states up to weak behavioural indistinguishability, either as divergence-blind stuttering equivalence on Kripke structures or as branching bisimulation on labelled transition systems. The associated algorithms are partition-refinement procedures based on inert transitions, bottom states, constellations, and smaller-half splitting [1603.05789].

The two literatures therefore embody different formal notions but a comparable methodological principle: replace syntactic or representational variation by an invariant notion of sameness, then study the algebraic or algorithmic structure of the resulting quotient. In the tensor setting, this leads to orbit classification, uniqueness results such as the rank-\(7\) \(2\times 2\) case, and practical equivalence testing. In the behavioural setting, it leads to refined minimization algorithms whose correctness and asymptotic complexity depend on careful control of refinement steps and data-structure updates [1902.03950][1603.05789].

Source: https://www.emergentmind.com/topics/de-groote-equivalence