Preliminaries
The operators of a tensor can provide a wide lens into its structure. In particular, the operators can distinguish two tensors under the equivalence relation of change of bases. Throughout, we suppose is a -grid of elements of some field . For , we write for the -entry of . Then can be interpreted as a -multilinear map , where . Taking this interpretation a bit further, we can evaluate at , and we write for this evaluation.
The ring of adjoint operators of is the -algebra
where denotes the matrix transpose. There is a naive, super-quadratic-time algorithm to compute the adjoint operators of .
Lemma. 1 There exists a linear system of variables and equations whose kernel determines a basis for . This is solved using field operations.
This computation is the bottleneck of many computations aimed at the isomorphism problem for tensors. In some nice cases, we can construct a basis for in nearly linear time. Before we discuss special cases, we need to discuss special cases of slicing tensors. These have general descriptions for general tensors, but we keep the discussion to what is necessary for bilinear maps.
For , a -slice of is the tensor , and a -slice of is . A -slice is slightly different, resulting in a bilinear form . These -slices correspond to matrices in the -grid. For example, a -slice for corresponds to the matrix whose -entry is for .
For , denote the -slice of at by , so . The defining equations for can be divided into different equations based on the -slices.
If is the corresponding matrix from the -slice , then each of these equations can be represented as a Sylvester-like system:
The Sylvester equation determines given the matrices and .
Theorem 1. [Maglione–Wilson (2018)] There exists an algorithm that, given a -grid with an invertible - and -slice, constructs a basis for using field operations.
The key idea of Theorem 1 is to consider the structure of the linear system in Lemma 1 in terms of slices. Because of the invertibility conditions, we can clear most of the entries in the system. What remains is a homogeneous system whose solution is given by a formula in terms of certain slices.
Problems
We propose the following questions as a gradual approach to more general solutions. We think the approach taken in Theorem 1 might lead to an affirmative answer to the first question.
Problem 1. Does there exist a sub-quadratic-time or quadratic-time algorithm to construct a basis for , given a tensor ?
We suspect that an affirmative answer to Open Problem 1 would apply immediately to general tensors and general nuclei, improving the computational complexity in all cases. In addition, computations of centroids should also benefit from such an affirmative answer as they are, in some sense, the intersection of nuclei. It is not clear if such a solution would apply to derivations. The derivation algebra of a -grid is
The proof of Lemma 1 applies to , and so computing has the same complexity.
Problem 2 Does there exist a sub-quadratic-time or quadratic-time algorithm to construct a basis for , given a tensor ?
An affirmative answer to Open Problem 2 implies an affirmative answer to Open Problem 2. Assuming such an answer extends to general tensors, this implies that every operator set of a tensor whose annihilator is generated by linear polynomials can be constructed more efficiently than the naive approach.
Additional Questions?
If you have additional questions, feel free to reach out to a maintainer / contributor on the contact page.