Guarantees for existence of a best canonical polyadic approximation of a noisy low-rank tensor

Abstract: The canonical polyadic decomposition (CPD) of a low rank tensor plays a major role in data analysis and signal processing by allowing for unique recovery of underlying factors. However, it is well known that the low rank CPD approximation problem is ill-posed. That is, a tensor may fail to have a best rank $R$ CPD approximation when $R>1$.This article gives deterministic bounds for the existence of best low rank tensor approximations over $\mathbb{K}=\mathbb{R}$ or $\mathbb{K}=\mathbb{C}$. More precisely, given a tensor $\mathcal{T} \in \mathbb{K}^{I \times I \times I}$ of rank $R \leq I$, we compute the radius of a Frobenius norm ball centered at $\mathcal{T}$ in which best $\mathbb{K}$-rank $R$ approximations are guaranteed to exist. In addition we show that every $\mathbb{K}$-rank $R$ tensor inside of this ball has a unique canonical polyadic decomposition. This neighborhood may be interpreted as a neighborhood of "mathematical truth" in with CPD approximation and computation is well-posed.In pursuit of these bounds, we describe low rank tensor decomposition as a ``joint generalized eigenvalue" problem. Using this framework, we show that, under mild assumptions, a low rank tensor which has rank strictly greater than border rank is defective in the sense of algebraic and geometric multiplicities for joint generalized eigenvalues. Bounds for existence of best low rank approximations are then obtained by establishing perturbation theoretic results for the joint generalized eigenvalue problem. In this way we establish a connection between existence of best low rank approximations and the tensor spectral norm. In addition we solve a "tensor Procrustes problem" which examines orthogonal compressions for pairs of tensors.The main results of the article are illustrated by a variety of numerical experiments.

\section{Introduction}

Tensors, i.e. multiindexed arrays, are natural generalizations of matrices and are commonplace  in fields such as machine learning and signal processing where data sets often have inherent higher-order structure \cite{Setal17,Cetal15}. Tensor decompositions are powerful tools which can be used either to recover source signals responsible for measured data or to compress data \cite{CJ10}.

One of the most popular decompositions for tensors is the canonical polyadic decomposition (CPD) which expresses a given tensor as a sum of {\it rank one tensors}. A key property of this decomposition is that, under mild conditions, the CPD of a low rank tensor is essentially unique \cite{DD13,K77,SB00}. This essential uniqueness makes CPD an indispensable tool in applications such as indpendent compenent analysis where one wants to recover source signals from measured data \cite{KB09,DDV00a}. In the case that the data has a natural tensor structure, CPD can allow for unique recovery of the source signals.

In applications one typically works with a measurement of a {\it low rank tensor}, that is, a tensor which is representable by a small number of source signals. This measurement is  always corrupted by measurement error, thus the measurement itself does not truly have low rank. As a consequence, in order to obtain the signals which represent the tensor, a practitioner must compute a low rank approximation to the measured tensor.

Although computation of a best low rank approximation is a natural problem, it is not necessarily a well-posed problem. In particular, an arbitrary tensor may not have a best rank (less than or equal to) $R$ approximation when $R>1$ \cite{KHL89,SL08}. Furthermore, even if a best rank $ {\mathbb{R}} $ approximation exists, it can fail to be unique. In fact, when considering CPD over $ {\mathbb{R}} $, nonexistence of a best low rank approximation occurs with positive probability \cite{S06,TKD88}. In contrast, it has recently been shown that  a unique best rank $R$ approximation exists with probability one when working over $ {\mathbb{C}} $ \cite{QML19}. See \cite{FO14,FS15} for further discussion of uniqueness of best approximations.\footnote{Note that uniqueness of a best rank $R$ approximation is a separate issue from uniqueness of the CPD of that approximation. }

The issue of potential nonexistence of best low rank approximations is not one which should be taken lightly. Should a tensor fail to have a best low rank approximation, one could of course choose to simply take a near optimal low rank approximation. However, when a best low rank approximation fails to exist, near optimal low rank approximations exhibit a phenomena known as {\it diverging components} \cite{KDS08,S08}. Roughly, the optimal way to represent a tensor which does not have a best low rank approximation is to express the tensor as destructive interference between signals which, in the limit, have infinite magnitude. Stated in an alternative fashion, computation of low rank approximations is unstable in a neighborhood around a tensor which does not have a best low rank approximation \cite{V17}.

Although computation of a best rank $R$ approximation can suffer from issues of nonexistence, tensor approximations typically work surprisingly well in practice. The reason for this is closely linked to the low rank structure of tensors which naturally arise in applications. In this article we explain this connection and quantify how much error can be present in the measurement of a low rank tensor before one encounters issues of nonexistence. Additionally, we further develop the perspective that low rank tensor approximation is in fact well-posed in a neighborhood of ``mathematical truth" for tensors occurring in practical applications. In particular, we provide deterministic bounds for the existence of best low rank approximations for tensors of order three.

More precisely, given a low rank tensor of order three, we compute a deterministic lower bound for the radius of a Frobenius norm ball centered at the tensor in which best low rank approximations are guaranteed to exist and in which every rank $R$ tensor has a unique CPD. While \cite[Theorem 5.1]{AGHKT14} gives a probabilistic lower bound for the radius of such a neighborhood which may be applied to a restricted class of rank $R$ symmetric tensors, the computation of a deterministic lower bound which may be applied to a large collection of rank $R$ tensors is novel.

 To obtain these bounds, we develop perturbation theoretic bounds for a ``joint generalized eigenvalue decomposition" which, with mild assumptions, is equivalent to the CPD of a low rank tensor. In particular, we show that low border rank tensors which have rank strictly greater than border rank can be identified with matrix tuples which are defective in the sense of algebraic and geometric multiplicities for the joint generalized eigenvalue problem. Given a low rank tensor, we then compute lower bounds on the distance to a tensor which is defective in this sense.

\subsection{Brief statement of results}

Given a tensor ${\mathcal T} \in  {\mathbb{K}} ^{R \times R \times K}$, we call a nonzero vector $ {\vec{x}}  \in  {\mathbb{K}} ^R$ a \df{joint generalized eigenvector} (JGE vector) of ${\mathcal T}$ if there exists a $ {\vec{\lambda}}  \in  {\mathbb{K}} ^K$ and a nonzero vector $ {\vec{y}}  \in  {\mathbb{K}} ^R$ such that
\begin{equation}
\label{eq:JGEVdefIntro}
 {\mathbf{T}} (:,:,\ell)  {\vec{x}}  =  {\lambda} _\ell  {\vec{y}}
\end{equation}
for all $\ell = 1, \dots, K$, and we call $\mathrm{span}( {\vec{\lambda}} )$ the \df{joint generalized eigenvalue} of ${\mathcal T}$ corresponding to $ {\vec{x}} $. In Section \ref{sec:JGEmultdefs} we formally define appropriate notions of characteristic polynomials and algebraic multiplicity for JGE vectors. For now we roughly define the \df{algebraic multiplicity} of $\mathrm{span}( {\vec{\lambda}} )$ to be the degree to which $ {\vec{\lambda}} $ ``divides" the characteristic polynomial of ${\mathcal T}$, and we define the \df{geometric multiplicity} of $\mathrm{span}( {\vec{\lambda}} )$ to be the number of linearly independent JGE vectors corresponding to $\mathrm{span}( {\vec{\lambda}} )$.

  \begin{theorem}
\label{theorem:CPDiffJGEVDSimple}
Let ${\mathcal T} \in  {\mathbb{K}} ^{R \times R \times K}$ and assume there is some vector $ {\vec{v}}  \in  {\mathbb{K}} ^K$ such that ${\mathcal T} \cdot_3  {\vec{v}} $ is invertible. Then ${\mathcal T}$ has $ {\mathbb{K}} $-rank $R$ if and only if ${\mathcal T}$ has a basis of JGE eigenvectors.

Furthermore, if ${\mathcal S} \in  {\mathbb{K}} ^{R \times R \times K}$ is a border $ {\mathbb{K}} $-rank $R$ tensor which has $ {\mathbb{K}} $-rank strictly greater than $R$ and there is a vector $ {\vec{u}}  \in  {\mathbb{K}} ^K$ such that ${\mathcal S} \cdot_3  {\vec{u}} $ is invertible, then ${\mathcal S}$ has a joint generalized eigenvalue with algebraic multiplicity strictly greater than geometric multiplicity.
\end{theorem}

 Theorem \ref{theorem:CPDiffJGEVDSimple} is a simplified statement of the upcoming Theorem \ref{theorem:rankvsmult}. We remark that Theorem \ref{theorem:CPDiffJGEVDSimple} can be viewed as a reformulation of algebraic-geometric tensor results in the language of the joint generalized eigenvalue problem, e.g. see \cite{Stra83,LM17}. While this result is known in a different language, the reformulation in terms of joint generalized eigenvalues illustrates the connection between perturbation theory for joint generalized eigenvalues and existence of best rank $R$ tensor approximations. Using this connection we give computable deterministic guarantees for the existence of best low rank tensor approximations.

\begin{theorem}
\label{theorem:MultiplePencilBoundSimple}
Let ${\mathcal T} \in  {\mathbb{K}} ^{R \times R \times K}$ be a tensor of $ {\mathbb{K}} $-rank $R$ and assume there is some vector $ {\vec{v}}  \in  {\mathbb{K}} ^K$ such that ${\mathcal T} \cdot_3  {\vec{v}} $ is invertible. For each $k=1,\dots, \lfloor K/2 \rfloor$, let $[\![ {\mathbf{A}} _k, {\mathbf{B}} _k, {\mathbf{C}} _k]\!]$ be a CPD of the matrix subpencil $({\mathcal T}(:,:,2k-1),{\mathcal T}(:,:,2k))$ of ${\mathcal T}$ where the columns of $ {\mathbf{C}} _k$ have unit Euclidean norm. Define
\[
\epsilon_k := \frac{\sigma_{\min} ( {\mathbf{A}} _k) \sigma_{\min} ( {\mathbf{B}} _k) \min_{i \neq j} \chi ( {\mathbf{C}} _k(:,i),  {\mathbf{C}} _k(:,j))}{2} \qquad \mathrm{for \ } k=1,\dots,\lfloor K/2 \rfloor
\]
and set $\epsilon = \|(\epsilon_1,\dots,\epsilon_{\lfloor K/2 \rfloor})\|_2$. If ${\mathcal W} \in  {\mathbb{K}} ^{R \times R \times K}$ satisfies $\|{\mathcal W}-{\mathcal T}\|_{\mathrm{F}} < \epsilon/2$ then ${\mathcal W}$ has a best $ {\mathbb{K}} $-rank $R$ approximation. Furthermore, any best $ {\mathbb{K}} $-rank $R$ approximation of ${\mathcal W}$ has a unique CPD, and, with probability one, the best $ {\mathbb{K}} $-rank $R$ approximation of ${\mathcal W}$ is unique. Here $\chi(\cdot,\cdot)$ denotes the chordal distance between unit vectors and $\sigma_{\min} (\cdot)$ denotes the smallest (possibly zero) singular value of a matrix.
\end{theorem}

The full statement of this result, Theorem \ref{theorem:MultiplePencilBound}, allows for more general collections of subpencils to be considered than those considered above. Additionally, a slight modification of this result may be applied to a measured tensor ${{\mathcal M}}'$ in an attempt to guarantee that ${{\mathcal M}}'$ has a best rank $R$ approximation, see Theorem \ref{theorem:MLSVDExistenceBound}.

\subsection{An overview of our approach}

The primary goal in this article is to give deterministic bounds for the existence of best low rank tensor approximations. A simple yet important idea for our approach is the following: If ${\hat{\cT}}'$ is a best border rank $R$ approximation of ${{\mathcal M}}' = {\mathcal T}' +  {\mathcal N} ' \in  {\mathbb{K}} ^{I_1 \times I_2 \times I_3}$ where ${\mathcal T}'$ has rank $R$, then the Frobenius distance from ${\mathcal T}'$ to ${\hat{\cT}}'$ cannot exceed $2 \| {\mathcal N} '\|_\mathrm{F}$. Therefore, to show that every tensor in a ball of radius $\epsilon/2$ around ${\mathcal T}'$ has a best rank $R$ approximation, it is sufficient to show that every border rank $R$ tensor in a $\epsilon$ radius ball around ${\mathcal T}'$ has rank $R$.

We temporarily restrict our discussion to the case of $R \times R \times K$ tensors which are ``slice mix invertible". In this setting, using the language of joint generalized eigenvalues, Theorem \ref{theorem:rankvsmult} shows that slice mix invertible tensors of size $R \times R \times K$ with border rank $R$ and rank strictly greater than $R$ have a joint generalized eigenvalue with algebraic multiplicity strictly greater than geometric multiplicity. We call such tensors defective. These defective tensors can be seen as a natural extension to the tensor setting of matrices which have a nontrivial Jordan block.

From this point our goal is to guarantee that a best border rank $R$ approximation is nondefective, hence has rank equal to border rank. Mirroring the matrix case, a natural sufficient condition for a tensor to be nondefective is that the tensor has distinct joint generalized eigenvalues. This motivates the study and development of perturbation theoretic bounds for the joint generalized eigenvalue problem, see Theorems \ref{theorem:TensMDBound} and \ref{theorem:TensBauerFike}.

Theorem \ref{theorem:TensMDBound} in particular may be used to guarantee that a perturbation of a matrix pencil is nondefective. By applying this result to a collection of subpencils of the signal tensor, we obtain a full tensor based bound which may be used to guarantee the existence of a best rank $R$ approximation of the measured tensor, see Theorems \ref{theorem:MultiplePencilBoundSimple}, \ref{theorem:MultiplePencilBound}, and \ref{theorem:MLSVDExistenceBound}.

The more general setting where ${{\mathcal M}}' = {\mathcal T}' + {\mathcal N} '$ has size $I_1 \times I_2 \times I_3$ and ${\mathcal T}'$ has rank $R \leq \min \{I_1,I_2\}$ is treated by a reduction to the $R \times R \times K$ setting. As above let ${\hat{\cT}}'$ be a best border rank $R$ approximation of ${{\mathcal M}}'={\mathcal T}'+ {\mathcal N} '$. Numerically, we wish to work with orthogonal compressions of ${\mathcal T}'$ and ${\hat{\cT}}'$, as this puts us in the $R \times R \times K$ setting discussed above. However, there can be adverse effects if one independently computes the compressions of ${\mathcal T}'$ and ${\hat{\cT}}'$.

Theorem \ref{theorem:TensorProcrustes} discusses how to orthogonally compress a \textit{pair}\footnote{Computing an orthogonal compression or a low multilinear rank approximation of a \textit{single} tensor is a standard preprocessing in many applications, e.g., see \cite{Setal17}. } of tensors and shows that there exist orthogonal compressions ${\mathcal T}$ and ${\hat{\cT}}$ of ${\mathcal T}'$ and ${\hat{\cT}}'$, respectively, so that
\begin{equation}
\label{eq:ProcrustesEqIntro}
\|{\mathcal T}-{\hat{\cT}}\|_\mathrm{F} \leq \|{\mathcal T}'-{\hat{\cT}}'\|_\mathrm{F}.
\end{equation}
From this point we may use our results in the $R \times R \times K$ setting to guarantee that ${\hat{\cT}}$ has rank $R$, hence ${\hat{\cT}}'$ has rank $R$.

\subsection{Connections to existing literature}

We now briefly discuss the relationship between our results and existing literature.

\subsubsection{Constrained existence results}

While our results  are novel in providing computable deterministic bounds for the existence of best unconstrained low rank tensor approximations, there are several known existence results in the setting of \textit{constrained} tensor decompositions. For example, \cite[Theorem 25]{LC14} shows that if one places constraints on the coherence of the factor matrices for the approximation of ${{\mathcal M}}'$, then the constrained best low rank approximation problem always has a minimizer. As a special case, this result implies that best low rank approximations exist if one restricts to considering factor matrices which have orthogonal columns. Also see \cite{Hars1984d,Setal12} and references therein.

 Another commonly considered constraint is nonnegativity. \cite[Theorem 6.1]{LC09}, see also \cite{QCL16}, shows that for nonnegative tensors, best low nonnegative rank approximations always exist. That is, a best low rank approximation exists if one constrains their factors to be nonnegative.

As previously mentioned, it is in general known that if ${{\mathcal M}}'$ does not have a best rank $R$ approximation then, near optimal approximations exhibit diverging components. Thus if one places constraints on the factors which prevent diverging components from occurring, then a best approximation with respect to those conditions will always exist. While in some applications it may be natural to consider such constraints, it is important to separately consider the \textit{unconstrained} problem, as constraints will not always be applicable. Additionally, as discussed in \cite{LC14}, imposing constraints to ``fix" ill-posedness may be misleading, as the problem at hand may in fact be ill-posed. For example, CPD may not be a proper data model for the problem, or the rank in consideration may be too low.

\subsubsection{Existence in tubular neighborhoods}

In the language of algebraic geometry, it is well known that given a smooth point on an algebraic variety, there is an open neighborhood around that point in which unique best approximations to the variety exist, e.g. see \cite{FS13,Aba80,Sot16,Lan12}. Translating to this terminology, the ``simple" (see Section \ref{sec:JGEdefs}) rank $R$ tensors we consider are smooth points on the $R$th secant variety of the Segre variety, thus there is necessarily some neighborhood around a simple rank $R$ tensor in which best rank $R$ approximations exist. This leads to a so called \textit{tubular neighborhood} around the set of simple rank $R$ tensors in which best rank $R$ approximations are guaranteed to exist.

 Though it was already known that there is  \textit{some neighborhood} around a given simple rank $R$ tensor in which best rank $R$ approximations are guaranteed to exist, until now there has been no known method for computing a lower bound for the \textit{radius} of such a neighborhood. That is, until now there has been no known method for providing a \textit{deterministic bound on noise levels} which may be used to guarantee that a measured tensor has a best rank $R$ approximation.

\subsubsection{The joint eigenvalue problem}

A problem which is closely related to our joint generalized eigenvalue is the \textit{joint eigenvalue problem} where one studies joint eigenvalues of a tuple of commuting matrices. In particular, by performing a modal multiplication in the first or second mode, a rank $R$ slice mix invertible tensor ${\mathcal T} \in  {\mathbb{K}} ^{R \times R \times K}$ can be transformed to a tensor whose frontal slices form a collection of commuting matrices.

There are indeed results in literature for the joint eigenvalue problem which may be used to show that if $( {\mathbf{T}} _1,\dots, {\mathbf{T}} _k) \in ( {\mathbb{K}} ^{R \times R})^K \cong  {\mathbb{K}} ^{R \times R \times K}$ is a simultaneously diagonalizable collection of matrices whose span contains a rank $R$ matrix, then there exists \textit{some neighborhood} around this tuple in which \textit{other tuples of commuting matrices} are simultaneously diagonalizable, e.g., see \cite[Theorem 3.5]{KP02}; however, the results in this direction that we are aware of do not provide a method for computing a lower bound for the \textit{radius} of such a neighborhood. Furthermore, though there is much literature, e.g. see \cite{KSS12}, for the joint eigenvalue problem, there are two significant shortcomings in attempting to directly apply this literature to our setting.

First, we essentially wish to determine if some tensor ${\hat{\cT}} \in  {\mathbb{K}} ^{R \times R \times K}$ is diagonalizable in the sense of the joint generalized eigenvalue problem based on its distance to a slice mix invertible tensor ${\mathcal T}$ which itself is diagonalizable in the sense of the joint generalized eigenvalue problem.  Here ${\hat{\cT}}$ is thought of as a best low rank approximation to some noisy observed tensor ${\mathcal T}+ {\mathcal N} $ for which ${\mathcal T}$ is the signal portion. However, \textit{a transformation which makes the frontal slices of ${\mathcal T}$ commute may not make the frontal slices of ${\hat{\cT}}$ commute and there may exist no invertible transformation which simultaneously converts both the frontal slices of ${\mathcal T}$ and the frontal slices of ${\hat{\cT}}$ to tuples of commuting matrices}. Thus, many tensors which can be treated with our results are missed by converting to the commuting case.

Second, treating a generalized eigenvalue problem\footnote{If $K=2$, then our notions of joint generalized eigenvectors and joint generalized eigenvalues coincide exactly with the (classical) notions of generalized eigenvalues and generalized eigenvectors as discussed in \cite{SS90}.} by transforming it to a classical eigenvalue problem can  have adverse effects on stability. This is illustrated by the discussion in \cite[Section 7.7.1]{GV96}. These difficulties motivate studying perturbation theory for the joint generalized eigenvalue problem in its own right.

\subsubsection{Perturbation theory for the generalized eigenvalue problem}

One of our main results, Theorem \ref{theorem:TensMDBound}, gives a lower bound on the radius of a neighborhood around a given diagonalizable matrix pencil in which all matrix pencils are diagonalizable. This result extends the ``spectral variation" bound of \cite[Chapter VI, Theorem 2.7]{SS90} to the much stronger notion of a ``matching distance" bound. Additionally, our extension is given in the tensor spectral norm. The proof of this extension may be unsurprising to experts in perturbation theory for classical generalized eigenvalues; however, the authors are unaware of an equivalent  extension in the literature.

In general the literature on perturbation theory for generalized eigenvalues and generalized eigenvectors is rich, for example see \cite{GV96,Naka11,Cang03}. We note that other bounds presented in the literature often use information from both the original matrix pencil and its perturbation. While such bounds can provide helpful insight into how generalized eigenvalues are affected by a specific perturbation, they are not well-suited to our needs.

\subsection{Organization}

Section \ref{sec:notationdefinitions} formally introduces our notation and definitions and gives formal statements of our main results. In Section \ref{sec:JGEbasics} we discuss basic theory for the joint generalized eigenvalue problem. Section \ref{sec:perturbation} discusses perturbation theoretic results for the joint generalized eigenvalue problem, while Section \ref{sec:Bounds} gives deterministic bounds for existence of a best low rank tensor approximation using results from Sections \ref{sec:JGEbasics} and \ref{sec:perturbation}. Numerical experiments illustrating our results are given in Section \ref{sec:figures}.

\section{Notation, Definitions, and Main Results}
\label{sec:notationdefinitions}

We now give our basic definitions and notations as well as formal statements of the main results of the article. Let $ {\mathbb{K}} $ denote $ {\mathbb{R}} $ or $ {\mathbb{C}} $. We denote scalars, vectors, matrices, and tensors with entries in $ {\mathbb{K}} $ by lower case $(m)$, bold lower case $(\vec{m})$, bold upper case $( {\mathbf{M}} )$, and calligraphic script $({{\mathcal M}})$, respectively. Additionally, script font $(\mathscr{M})$ is used to denote a subspace of $ {\mathbb{K}} ^I$. Given a matrix $ {\mathbf{M}}  \in  {\mathbb{K}} ^{I_1 \times I_2}$, let $ {\mathbf{M}} ^\mathrm{T}, {\mathbf{M}} ^\mathrm{H}, {\mathbf{M}} ^\dagger$ denote the transpose, Hermitian, and Moore-Penrose pseudo inverse of $ {\mathbf{M}} $, respectively.  If $ {\mathbf{M}} $ is invertible, let $ {\mathbf{M}} ^{-1}$ and $ {\mathbf{M}} ^{-\mathrm{T}}$ denote the inverse and inverse transpose of $ {\mathbf{M}} $. Additionally let $\| {\mathbf{M}} \|_2$ and $\| {\mathbf{M}} \|_\mathrm{F}$ denote the spectral and Frobenius norms of $ {\mathbf{M}} $.

Let $\langle \cdot, \cdot \rangle$ denote the usual inner product on $ {\mathbb{K}} ^R$ which is conjugate linear in the second term. We make frequent use of the linear form
\[
\langle  {\vec{\lambda}} ,  {\bar{\boldsymbol{\gamma}}}  \rangle = \lambda_1 \gamma_1 + \lambda_2 \gamma_2 +\cdots + \lambda_K \gamma_K.
\]
Here $ {\vec{\lambda}} , {\boldsymbol{\gamma}}  \in  {\mathbb{K}} ^K$ and $ {\bar{\boldsymbol{\gamma}}} $ denotes the complex conjugate of $ {\boldsymbol{\gamma}} $. Note that $\langle \cdot,\cdot \rangle$ denotes the usual inner product on $ {\mathbb{K}} ^K$ which is conjugate linear in the second term.\footnote{The linear form we are interested in is linear in $ {\boldsymbol{\gamma}} $ while the inner product $\langle \cdot,\cdot \rangle$ is conjugate linear in the second argument. Hence, the complex conjugate of $ {\boldsymbol{\gamma}} $ appears in the expression $\langle  {\vec{\lambda}} ,  {\bar{\boldsymbol{\gamma}}}  \rangle$. We use the notation $\langle \cdot,\cdot \rangle$ as we often consider some geometric relationship between $ {\vec{\lambda}} $ and $ {\boldsymbol{\gamma}} $.}

A \df{tensor} is a multiindexed array ${{\mathcal M}}'$ with entries in $ {\mathbb{K}} $. In this article our primary interest lies with tensors of order $3$. That is, we study tensors ${{\mathcal M}}' \in  {\mathbb{K}} ^{I_1 \times I_2 \times I_3}$.  Often we will restrict our study to low rank (pencil-like) tensors having size $R \times R \times K$ and rank $R$. To the reader keep track of context, we label tensors of size $I_1 \times I_2 \times I_3$ with the symbol $'$, while tensors of size $R \times R \times K$ not labelled in this way. So for example, ${{\mathcal M}}' \in  {\mathbb{K}} ^{I_1 \times I_2 \times I_3}$ while ${{\mathcal M}} \in  {\mathbb{K}} ^{R \times R \times K}$.

Define the \df{outer product} of vectors $ {\vec{a}}  \in  {\mathbb{K}} ^{I_1},  {\vec{b}}  \in  {\mathbb{K}} ^{I_2}$ and  $ {\vec{c}}  \in  {\mathbb{K}} ^{I_3}$, denoted $ {\vec{a}}  \otimes  {\vec{b}}  \otimes  {\vec{c}} $, to be the $I_1 \times I_2 \times I_3$ tensor whose $i,j,k$th entry is given by
\[
( {\vec{a}}  \otimes  {\vec{b}}  \otimes  {\vec{c}} )_{ijk}= a_i b_j c_k.
\]
A tensor which can be expressed as the outer product of nonzero vectors is called a \df{rank one tensor}.

Every tensor may be expressed as a sum of rank one tensors. That is, given a tensor ${{\mathcal M}}' \in  {\mathbb{K}} ^{I_1 \times I_2 \times I_3}$ there exists a collection of vector tuples  $\{( {\vec{a}} _r,  {\vec{b}} _r,  {\vec{c}} _r) \}_{r=1}^R \subset  {\mathbb{K}} ^{I_1} \times  {\mathbb{K}} ^{I_2} \times  {\mathbb{K}} ^{I_3}$ such that
\begin{equation}
{{\mathcal M}}' = \sum_{r=1}^R  {\vec{a}} _r  \otimes  {\vec{b}} _r \otimes  {\vec{c}} _r.
\label{eq:CPDdef}
\end{equation}
In the case that $R$ is as small as possible, the tensor ${{\mathcal M}}'$ is said to have $ {\mathbb{K}} $-\df{rank} $R$, and the decomposition in equation \eqref{eq:CPDdef} is called a \df{canonical polyadic decomposition} (CPD) of ${{\mathcal M}}'$. We let $\mbox{rank}_ {\mathbb{K}}  ({{\mathcal M}}')$ denote the rank of ${{\mathcal M}}'$ over $ {\mathbb{K}} $. A fundamental question in the study of tensors and in many applications is to determine the CPD of a given tensor over either $ {\mathbb{R}} $ or $ {\mathbb{C}} $. We emphasize that the field in which the entries of the tensor lie may differ from the field over which the decomposition is taken. In particular, if $ {{\mathcal M}}'$ is a tensor with real entries, then both real and complex rank are interesting to consider, and these ranks are not necessarily equal, e.g. see \cite{BK99,B13,K89}.

Let $ {\mathbf{A}}  \in  {\mathbb{K}} ^{I_1 \times R}$ be the matrix whose $r$th column is the vector $ {\vec{a}} _r$, and similarly define $ {\mathbf{B}}  \in  {\mathbb{K}} ^{I_2 \times R}$ and $ {\mathbf{C}}  \in  {\mathbb{K}} ^{I_3 \times R}$. Write
\[
{{\mathcal M}}'= {\cp{ {\mathbf{A}} , {\mathbf{B}} , {\mathbf{C}} }}
\]
as a compact notation for the CPD of ${{\mathcal M}}'$. The matrices $ {\mathbf{A}} , {\mathbf{B}} ,$ and $ {\mathbf{C}} $ are called the \df{factor matrices} of ${{\mathcal M}}'$. As discussed in \cite{Setal17} one has
\[
{{\mathcal M}}'(:,:,k) =  {\mathbf{A}}  D_k ( {\mathbf{C}} )  {\mathbf{B}} ^\mathrm{T} \qquad \mathrm{for \ all \ } k=1,\dots, I_3.
\]
Here $D_k ( {\mathbf{C}} ) \in  {\mathbb{K}} ^{R \times R}$ is the diagonal matrix whose diagonal entries are given by the $k$th row of $ {\mathbf{C}} $.

Say a CPD ${{\mathcal M}}'={\cp{ {\mathbf{A}} , {\mathbf{B}} , {\mathbf{C}} }}$ is \df{unique} if for any other CPD ${{\mathcal M}}'=[\![\tilde {\mathbf{A}} ,\tilde {\mathbf{B}} ,\tilde {\mathbf{C}} ]\!]$, there exist a permutation matrix $\Pi$ and invertible diagonal matrices $ {\mathbf{D}} _ {\mathbf{A}} ,  {\mathbf{D}} _ {\mathbf{B}} , {\mathbf{D}} _ {\mathbf{C}} $ such that
\[
 {\mathbf{A}}   {\mathbf{D}} _ {\mathbf{A}}  \Pi = \tilde {\mathbf{A}}  {\quad \quad \mathrm{and} \quad \quad}   {\mathbf{B}}   {\mathbf{D}} _ {\mathbf{B}}  \Pi = \tilde {\mathbf{B}}  {\quad \quad \mathrm{and} \quad \quad}   {\mathbf{C}}   {\mathbf{D}} _ {\mathbf{C}}  \Pi = \tilde {\mathbf{C}} .
\]
That is, a CPD is unique if the only indeterminacies in the CPD are permutation and scaling/counter scaling of the columns of the factor matrices.

Central to the potential nonexistence of a best low rank tensor approximation is the fact that the set of tensors of rank less than or equal to $R$ is, in general, not a closed set \cite{KHL89,SL08}. For example, there are tensors having rank $3$ which are a limit of rank $2$ tensors. Of course, such a tensor cannot have a best rank $2$ approximation. Say a ${{\mathcal M}}'$ tensor has \df{border $ {\mathbb{K}} $-rank $R$} if ${{\mathcal M}}'$ is a limit of $ {\mathbb{K}} $-rank $R$ tensors and there is no integer $L<R$ such that ${{\mathcal M}}'$ is a limit of $ {\mathbb{K}} $-rank $L$ tensors. That is, $R$ should be the smallest integer such that ${{\mathcal M}}'$ is a limit of $ {\mathbb{K}} $-rank $R$ tensors. We let $ {\mathbb{K}} $-$\underline{\mbox{rank}} ({{\mathcal M}}')$ denote the border $ {\mathbb{K}} $-rank of ${{\mathcal M}}'$.

We will make frequent use of the fact that the set of border $ {\mathbb{K}} $-rank $R$ tensors is a closed set. This in particular implies that a tensor always has a (not necessarily unique) best border $ {\mathbb{K}} $-rank $R$ approximation.

\subsection{Approximately low rank tensors}
\label{sec:Compress}
 In this article we are mainly interested in tensors which have (approximately) low rank in the following sense. Stated briefly, for tensors ${{\mathcal M}}' \in  {\mathbb{K}} ^{I_1 \times I_2 \times I_3}$ occuring in applications, one often has ${{\mathcal M}}' = {\mathcal T}'+ {\mathcal N} '$ where ${\mathcal T}'$ is a rank $R$ tensor with rank equal to multilinear rank and $ {\mathcal N} '$ is noise. That is, ${{\mathcal M}}'$ is a noisy measurement of low rank signal tensor ${\mathcal T}'$. Our main result, Theorem \ref{theorem:MultiplePencilBound}, uses information from ${\mathcal T}'$ to compute an upper bound $\epsilon$ on the Frobenius norm of the noise $\| {\mathcal N} '\|_\mathrm{F}$ such that ${{\mathcal M}}'$ is guaranteed to have a best rank $R$ approximation if $\| {\mathcal N} '\|_\mathrm{F} < \epsilon$. A detailed explanation of our perspective follows.

Given a tensor ${{\mathcal M}}' \in  {\mathbb{K}} ^{I_1 \times I_2 \times I_3}$, define
\begin{equation}
\label{eq:MLrankDef}
\begin{array}{ccc}
R_1 ({{\mathcal M}}'):= \dim \mathrm{span} \{{{\mathcal M}}'(:,k,\ell)\}_{\forall k,\ell}, \\
R_2 ({{\mathcal M}}'):= \dim \mathrm{span} \{{{\mathcal M}}'(j,:,\ell)\}_{\forall j,\ell}, \\
R_3 ({{\mathcal M}}'):= \dim \mathrm{span} \{{{\mathcal M}}'(j,k,:)\}_{\forall j,k}. \\
\end{array}
\end{equation}
The integer triple $(R_1 ({{\mathcal M}}'), R_2 ({{\mathcal M}}'), R_3 ({{\mathcal M}}'))$ is called the \df{multilinear rank} of ${{\mathcal M}}'$. In the case where $R_i ({{\mathcal M}}') = I_i$ for each $i=1,2,3$ we say ${{\mathcal M}}'$ has \df{full multilinear rank}. The vectors in equation \eqref{eq:MLrankDef} used to define each $R_i ({{\mathcal M}}')$ are called the \df{mode-$i$ fibers} of ${{\mathcal M}}'$, and we let ${{\mathcal M}}'_{(i)}$ denote a matrix whose columns are equal to the mode-$i$ fibers of ${{\mathcal M}}'$. The matrix ${{\mathcal M}}'_{(i)}$ is called a \df{mode-$i$ unfolding} of ${{\mathcal M}}'$.\footnote{There is freedom in the ordering of the columns of ${{\mathcal M}}'_{(i)}$, so in this sense we have not well-defined ${{\mathcal M}}'_{(i)}$. However, the properties of mode-$i$ unfoldings that we are most interested in are invariant under permutations of columns, so this ambiguity makes little difference for our purpose so long as one is consistent in the ordering of columns.}

It is not difficult to show that one always has
\begin{equation}
\label{eq:RankVsMLRank}
\mbox{rank}_ {\mathbb{K}}  ({{\mathcal M}}') \geq \max \{R_1 ({{\mathcal M}}'), R_2 ({{\mathcal M}}'), R_3 ({{\mathcal M}}')\}.
\end{equation}
For generic tensors this inequality is strict; however, for the signal portion of many tensors occurring in applications we have equality. In particular, in applications the tensor ${{\mathcal M}}'$ is often a noisy measurement of some low rank signal tensor ${\mathcal T}'$ which satisfies  $\mbox{rank}_{ {\mathbb{K}} } ({\mathcal T}') \leq \min\{I_1,I_2,I_3\}$ and, up to a permutation of indices,
\begin{equation}
\label{eq:LowRankNotion}
R=\mbox{rank}_ {\mathbb{K}}  ({\mathcal T}') = R_{1} ({\mathcal T}') = R_{2} ({\mathcal T}') \geq R_{3} ({\mathcal T}')=K.
\end{equation}
That is, we think of ${{\mathcal M}}'$ as ${{\mathcal M}}'={\mathcal T}'+ {\mathcal N} '$ where ${\mathcal T}'$ has $ {\mathbb{K}} $-rank $R$ and multilinear rank $(R,R,K)$ and $ {\mathcal N} '$ is noise. Note that our results still hold if the multilinear rank of ${\mathcal T}'$ in one mode is less than $R$, and allowing this increases flexibility, hence $R_3 ({\mathcal T}')=K$ is allowed to be less than or equal to $R$.

For later use we note that the multilinear ranks of a tensor not only lower bound the rank of the tensor, but also the border rank of the tensor. That is, one has
\begin{equation}
\label{eq:BRankVsMLRank}
 {\mathbb{K}} \text{-}\underline{\mbox{rank}} ({{\mathcal M}}') \geq \max \{R_1 ({{\mathcal M}}'), R_2 ({{\mathcal M}}'), R_3 ({{\mathcal M}}')\}.
\end{equation}
This fact follows from lower semi-continuity of matrix rank, as tensors in a sufficiently small neighborhood around ${{\mathcal M}}'$ must have multilinear rank greater than or equal to that of ${{\mathcal M}}'$. It follows that the multilinear rank of ${{\mathcal M}}'$ provides a lower bound for the rank of tensors in a neighborhood of ${{\mathcal M}}'$.

\subsection{Orthogonal compressions and modal products}
In the case where $R_i({\mathcal T}') < I_i$ for some $i$ the tensor ${\mathcal T}'$ maybe be orthogonally compressed to a tensor ${\mathcal T}\in  {\mathbb{K}} ^{R_1 ({\mathcal T}') \times R_2 ({\mathcal T}') \times R_3 ({\mathcal T}')}$. In particular, for $i=1,2,3$ let $ {\mathbf{V}} _i \in  {\mathbb{K}} ^{I_i \times R_i ({\mathcal T}')}$ be a matrix whose columns form an orthonormal basis for the span of the mode-$i$ fibers of ${\mathcal T}'$. Then the tensor
\begin{equation}
\label{eq:orthcompr}
{\mathcal T}' = {\mathcal T}' \cdot_1  {\mathbf{V}} _1^\mathrm{H} \cdot_2  {\mathbf{V}} _2^\mathrm{H} \cdot_3  {\mathbf{V}} _3^\mathrm{H} \in  {\mathbb{K}} ^{R_1 ({\mathcal T}') \times R_2 ({\mathcal T}') \times R_3 ({\mathcal T}')}
\end{equation}
contains the same algebraic and geometric information as ${\mathcal T}$, and we call this tensor an \df{orthogonal compression} of ${\mathcal T}'$. One may recover ${\mathcal T}'$ from ${\mathcal T}$ using
\begin{equation}
\label{eq:orthrecover}
{\mathcal T}'=({\mathcal T}' \cdot_1  {\mathbf{V}} _1^\mathrm{H} \cdot_2  {\mathbf{V}} _2^\mathrm{H} \cdot_3  {\mathbf{V}} _3^\mathrm{H}) \cdot_1  {\mathbf{V}} _1 \cdot_2  {\mathbf{V}} _2  \cdot_3  {\mathbf{V}} _3.
\end{equation}

In equations \eqref{eq:orthcompr} and \eqref{eq:orthrecover} above, the product ${\mathcal T}' \cdot_1  {\mathbf{V}} _1^\mathrm{H}$ denotes the \df{mode-1 product} between the tensor ${\mathcal T}'$ and the matrix $ {\mathbf{V}} _1^\mathrm{H}$. In general, the mode-1 product between a tensor ${{\mathcal M}}' \in  {\mathbb{K}} ^{I_1 \times I_2 \times I_3}$ and a matrix $ {\mathbf{U}}  \in  {\mathbb{K}} ^{J_1 \times I_1}$ is the tensor ${{\mathcal M}}' \cdot_1  {\mathbf{U}}  \in  {\mathbb{K}} ^{J_1 \times I_2 \times I_3}$ with mode-$1$ fibers given by
\[
({{\mathcal M}}' \cdot_1  {\mathbf{U}} )(:,k,\ell) =  {\mathbf{U}} ({{\mathcal M}}' (:,k,\ell)) \qquad \mathrm{for \ all \ } k,\ell.
\]
The mode-$2$ and mode-$3$ products $\cdot_2$ and $\cdot_3$ are defined analogously.

A common first step in CPD computation in applications is to compute\footnote{To compute an orthogonal compression or an (approximately) best low multilinear rank approximation of a tensor it is common to take the truncated core of a multilinear singular value decomposition of the tensor \cite{DDV00b}.}  the core of a multilinear $(R_1 ({\mathcal T}'), R_2 ({\mathcal T}'), R_3 ({\mathcal T}'))$ rank approximation to ${{\mathcal M}}'={\mathcal T}'+ {\mathcal N} ' \in  {\mathbb{K}} ^{I_1 \times I_2 \times I_3}$, and it will be convenient for us to have notation for such tensors. Recalling that we typically assume ${\mathcal T}'$ has multilinear rank $R \times R \times K$, we let ${\mathcal W}'$ denote a tensor of size $I_1 \times I_2 \times I_3$ which has multilinear rank $(R,R,K)$ and we let ${\mathcal W} \in  {\mathbb{K}} ^{R \times R \times K}$ be an orthogonal compression of ${\mathcal W}'$. In particular, ${\mathcal W}$ is thought of as the core of a low multilinear rank approximation to ${{\mathcal M}}'$, while ${\mathcal W}'$ itself is thought of as a low multilinear rank approximation to ${{\mathcal M}}'$.

\subsubsection{Implicit subpencils of tensors}

We frequently consider subpencils of a given tensor. That is, we will consider pairs of matrices which are formed from the frontal slices of a given tensor. The subpencils considered in this article differ slightly from the standard terminology in that we allow subpencils which are ``revealed" by an orthogonal change of basis in the third mode of the tensor. We call such a subpencil an \df{implicit subpencil}.

Defined precisely, an implicit subpencil of $ {\mathcal S}$ is an $R \times R \times 2$ tensor of the form ${\mathcal S} = {\mathcal W} \cdot_3  {\mathbf{V}}  \in  {\mathbb{K}} ^{R \times R \times 2}$ where $ {\mathbf{V}}  \in  {\mathbb{K}} ^{2 \times K}$ has rank $2$ and has orthonormal rows. Since we always work with implicit subpencils, we drop the distinction ``implicit" in the remainder of the article.

\subsection{The joint generalized eigenvalue problem}
\label{sec:JGEdefs}

 We now discuss the joint generalized eigenvalue problem in greater detail. Note that we will interchangeably view ${\mathcal T} \in  {\mathbb{K}} ^{R \times R \times K}$ as either a tensor or a $K$-tuple of $R \times R$ matrices tuple via the map
\[
{\mathcal T} \to ({\mathcal T}(:,:,1), {\mathcal T}(:,:,2), \cdots ,{\mathcal T}(:,:,K)).
\]
For notational convenience we denote $ {\mathbf{T}} _k={\mathcal T}(:,:,k)$.

Recall that given a tensor (or a matrix tuple) ${\mathcal T} \in  {\mathbb{K}} ^{R \times R \times K}$, we call a nonzero vector $ {\vec{x}}  \in  {\mathbb{K}} ^R$ a joint generalized eigenvector of ${\mathcal T}$ if there exists a $ {\vec{\lambda}}  \in  {\mathbb{K}} ^K$ and a nonzero vector $ {\vec{y}}  \in  {\mathbb{K}} ^R$ such that
\begin{equation}
\label{eq:JGEVdef}
 {\mathbf{T}} _\ell  {\vec{x}}  =  {\lambda} _\ell  {\vec{y}}
\end{equation}
for all $\ell = 1, \dots, K$. Note that assuming $ {\vec{y}}  \neq  {\vec{0}} $ implies that $ {\lambda} _\ell=0$ if and only if $ {\mathbf{T}} _\ell  {\vec{x}}  =  {\vec{0}} $.

There are natural scaling indeterminates in equation \eqref{eq:JGEVdef}. In particular, if the triple $( {\vec{x}} , {\vec{\lambda}} , {\vec{y}} )$ is a solution to equation \eqref{eq:JGEVdef}, then so is $(\alpha  {\vec{x}} ,\beta  {\vec{\lambda}} , \gamma  {\vec{y}} )$ for any nonzero scalars $\alpha,\beta,\gamma \in  {\mathbb{K}} $ satisfying $\alpha  = \beta \gamma$. For this reason, if $( {\vec{\lambda}} , {\vec{x}} )$ is a solution to \eqref{eq:JGEVdef}, then we say the one-dimensional subspace ${{\mathscr L}}:= \mathrm{span} ( {\vec{\lambda}} ) \subset  {\mathbb{K}} ^K$ is a \df{joint generalized eigenvalue} (JGE value) of ${\mathcal T}$ and that $({{\mathscr L}}, {\vec{x}} )$ is a \df{joint generalized eigenpair} of ${\mathcal T}$. We typically use ${{\mathscr L}} \subset  {\mathbb{K}} ^K$ to denote a JGE value, while $ {\vec{\lambda}}  \in  {\mathbb{K}} ^K$ typically denotes a vector of unit norm such that $\mathrm{span}( {\vec{\lambda}} )={{\mathscr L}}$. Say a tensor has a \df{joint generalized eigenbasis} if its joint generalized eigenvectors form a basis for $ {\mathbb{K}} ^R$. We call the set of JGE values of a tensor the \df{(joint generalized) spectrum} of the tensor. An equivalent formulation of joint generalized eigenvectors is used in \cite{ZLP18} for blind source extraction.

We note that some authors have considered other notions of eigenvalues and eigenvectors for a tensor. These are not equivalent to and should not be confused with joint generalized eigenvalues. The eigenvalues and eigenvectors considered by other authors are intended to allow a spectral theory for general tensors. See \cite{S16,L05,QCC18}.

\subsection{Multiplicities of JGE values}
\label{sec:JGEmultdefs}
One of the main themes of this article is that low border rank tensors which have rank strictly greater than border rank can be identified with matrix tuples which are defective in the sense of algebraic and geometric multiplicities for joint generalized eigenvalues. These multiplicities we now define. To simplify the exposition, we only consider tensors of size $R \times R \times K$ when discussing multiplicities. While it is possible to define multiplicities of JGE values for general third order tensors, certain pathologies arise in the $I_1 \times I_2 \times I_3$ case which must be treated with care. For our purposes, it will not be necessary to consider multiplicities in the $I_1 \times I_2 \times I_3$ setting.

Given a tensor ${\mathcal T} \in  {\mathbb{K}} ^{R \times R \times K}$, define the polynomial $p_{\mathcal T} (\gamma)$ in the indeterminate $ {\boldsymbol{\gamma}} =(\gamma_1,\dots, \gamma_K)$ by
\[
p_{\mathcal T} ( {\boldsymbol{\gamma}} ):=\det \left( \sum_{k=1}^K \gamma_k  {\mathbf{T}} _k \right).
\]
We call $p_{\mathcal T} ( {\boldsymbol{\gamma}} )$ the \df{characteristic polynomial} of ${\mathcal T}$.
Given a nonzero $ {\vec{\lambda}}  \in  {\mathbb{K}} ^K$, the \df{algebraic multiplicity} of ${{\mathscr L}}=\mathrm{span}( {\vec{\lambda}} )$ as a JGE value of the tensor ${\mathcal T} \in  {\mathbb{K}} ^{R \times R \times K}$, denoted $\mathrm{am}({{\mathscr L}})$, is the largest integer $m$ such that $\langle  {\vec{\lambda}} ,  {\bar{\boldsymbol{\gamma}}}  \rangle ^m$ divides the characteristic polynomial of ${\mathcal T}$. If $p_{\mathcal T}$ factors as a product of $R$  (not necessarily distinct) linear terms, then we say that ${\mathcal T}$ has $R$ JGE values counting algebraic multiplicity.

Define the \df{geometric multiplicity} of ${{\mathscr L}}$ as a JGE value of ${\mathcal T} \in  {\mathbb{K}} ^{R \times R \times K}$, denoted $\mathrm{gm}({{\mathscr L}})$, to be the dimension of the span of the set of JGE vectors of ${\mathcal T}$ which correspond to $\mathrm{span}( {\vec{\lambda}} )$. That is,
\[
\mathrm{gm}({{\mathscr L}})= \dim \mathrm{span} (\{ {\vec{x}}  |\  ({{\mathscr L}}, {\vec{x}} ) \mathrm{\ is \ a \ JGE \ pair \ of \ } {\mathcal T}\}).
\]
Similar to the cases of classical and generalized eigenvalues, the algebraic multiplicity of a JGE value is always greater than or equal to  the geometric multiplicity of said JGE value, see Proposition \ref{proposition:alggreatergeo}. However, we warn that it is possible for a JGE value to have nonzero algebraic multiplicity while having geometric multiplicity equal to zero, see Example \ref{exa:geomultzero}.

\begin{definition}
Say a tensor ${\mathcal T} \in {R \times R \times K}$ is \df{defective} if it has a JGE value ${{\mathscr L}}$ with $\mathrm{am}({{\mathscr L}})>\mathrm{gm}({{\mathscr L}})$.
\end{definition}

A degenerate case can occur in which $p_{{\mathcal T}} ( {\boldsymbol{\gamma}} )$ is identically equal to zero. In this case, $\langle   {\vec{\lambda}} , {\bar{\boldsymbol{\gamma}}}  \rangle^m$ divides $p_{{\mathcal T}} ( {\boldsymbol{\gamma}} )$ for all $ {\vec{\lambda}}  \in  {\mathbb{K}} ^K$ and for all integers $m$, hence we may reasonably say that $\mathrm{span} ( {\vec{\lambda}} )$ is a JGE value of ${\mathcal T}$ with infinite algebraic multiplicity for all $ {\vec{\lambda}}  \in  {\mathbb{K}} ^K$. In this sense, if $p_{{\mathcal T}} ( {\boldsymbol{\gamma}} ) \equiv 0$ then the tensor ${\mathcal T}$ is necessarily defective.

It is straightforward to show that $p_{{\mathcal T}} (\gamma)$ is not identically zero if and only if the span of the frontal slices of ${\mathcal T}$ contains an invertible matrix. We call such a tensor a \df{slice mix invertible} tensor. See Example \ref{exa:pequiv0} and Lemma \ref{lemma:invertslice} for further discussion.

\begin{definition}
Say ${\mathcal T} \in {R \times R \times K}$ is \df{simple} if ${\mathcal T}$ is slice mix invertible and has $R$ distinct JGE values each having algebraic multiplicity equal to $1$. In other words, ${\mathcal T}$ is simple if $p_{{\mathcal T}}( {\boldsymbol{\gamma}} )$ is not identically zero and factors as a product of $R$ distinct linear terms.
\end{definition}

We are now in position to formally state the connection between the joint generalized eigenvalue decomposition and the canonical polyadic decomposition.

\begin{theorem}
\label{theorem:rankvsmult}
Let ${\mathcal T}' \in  {\mathbb{K}} ^{I_1 \times I_2 \times I_3}$ and assume ${\mathcal T}'$ has multilinear rank $(R,R,K)$. Let ${\mathcal T} \in  {\mathbb{K}} ^{R \times R \times K}$ be an orthogonal compression of ${\mathcal T}'$. Then we have the following.

\begin{enumerate}
\item \label{it:RankIFFBasis} ${\mathcal T}'$ has $ {\mathbb{K}} $-rank $R$ if and only if the JGE vectors of ${\mathcal T}$ span $ {\mathbb{K}} ^R$.  Furthermore, if ${\mathcal T}={\cp{ {\mathbf{A}} , {\mathbf{B}} , {\mathbf{C}} }}$ has $ {\mathbb{K}} $-rank $R$, then the JGE values of ${\mathcal T}$ are equal to the spans of columns of the factor matrix $ {\mathbf{C}} $, and the columns of $ {\mathbf{B}} ^{-\mathrm{T}}$ are JGE vectors of ${\mathcal T}$.

\item \label{it:NonderogatoryImpliesRankR} If ${\mathcal T}'$ has border $ {\mathbb{K}} $-rank $R$ and ${\mathcal T}$ is simple then ${\mathcal T}'$ has $ {\mathbb{K}} $-rank $R$.

\item \label{it:BRankLessRankImpliesDefective} If ${\mathcal T}'$ has border $ {\mathbb{K}} $-rank $R$ but $ {\mathbb{K}} $-rank strictly greater than $R$ and if ${\mathcal T}$ is slice mix invertible, then ${\mathcal T}$ must have a JGE value ${{\mathscr L}} \subset  {\mathbb{K}} ^K$ which satisfies
\[
\mathrm{am}({{\mathscr L}}) > \mathrm{gm}({{\mathscr L}}) \geq 1.
\]
In this case, ${\mathcal T}$ is defective in the sense of joint generalized eigenvalues.
\end{enumerate}
\end{theorem}

\begin{proof}
The proof of Theorem \ref{theorem:rankvsmult} is given in Section \ref{sec:rankvsmultproof}.
\end{proof}

Theorem \ref{theorem:rankvsmult} leads to a strategy for determining the radius of an open ball centered at a given rank $R$ tensor in which best $ {\mathbb{K}} $-rank $R$ approximations are guaranteed to exist. In particular, it is sufficient to guarantee that nearby border $ {\mathbb{K}} $-rank $R$ tensors are simple.

\begin{proposition}
\label{proposition:existenceCor}
Let ${\mathcal T}' \in  {\mathbb{K}} ^{I_1 \times I_2 \times I_3}$ and assume ${\mathcal T}'$ has $ {\mathbb{K}} $-rank $R$ and multilinear rank $(R,R,K)$, and let ${\mathcal T} \in  {\mathbb{K}} ^{R \times R \times K}$ be an orthogonal compression of ${\mathcal T}'$. If every border $ {\mathbb{K}} $-rank $R$ tensor in an open ball of Frobenius radius $\epsilon$ around ${\mathcal T}$ is simple, then every tensor in an open ball of radius $\epsilon/2$ around ${\mathcal T}'$ has a best $ {\mathbb{K}} $-rank $R$ approximation. In addition, every $ {\mathbb{K}} $-rank $R$ tensor inside the open ball of Frobenius radius $\epsilon$ centered at ${\mathcal T}'$ has a unique CPD. Furthermore, the set of tensors inside this open ball for which the best rank $R$ approximation is not unique has Lebesgue measure zero.

In particular, with these assumptions the ball of radius $\epsilon$ around ${\mathcal T}'$ does not contain a tensor with border $ {\mathbb{K}} $-rank $R$ and $ {\mathbb{K}} $-rank strictly greater than $R$.
\end{proposition}
\begin{proof} See Section \ref{sec:existenceCorProof}.
\end{proof}

 As previously mentioned, the existence of some neighborhood in which best rank $R$ approximations exist is known. The main motivation for Proposition \ref{proposition:existenceCor} is that a lower bound for the radius of a ball  centered at a simple rank $R$ tensor in which border rank $R$-tensors are simple can be computed using perturbation theoretic results for joint generalized eigenvalues, e.g., see Theorems \ref{theorem:MultiplePencilBound} and \ref{theorem:TensMDBound}.

We also note that \cite{QML19, FO14} together show that the set of tensors which do not have a unique best border $ {\mathbb{K}} $-rank $R$ approximation has measure $0$. The uniqueness of best $ {\mathbb{K}} $-rank $R$ approximations discussed in Proposition \ref{proposition:existenceCor} follows quickly from these results after having shown existence of best $ {\mathbb{K}} $-rank $R$ approximations.\footnote{In fact it is known that there is a neighborhood in which best rank $R$ approximations exist \textit{and are unique}. In regards to uniqueness, the conclusion of Proposition \ref{proposition:existenceCor} is slightly weaker; however, the fact that a lower bound for $\epsilon$ is computable is a significant advantage.}

\subsection{Deterministic bounds for existence of a best rank $R$ approximation}

We now explain how to compute the $\epsilon$ appearing in Proposition \ref{proposition:existenceCor}. For general $K$, Proposition \ref{proposition:DefectiveSubpencils} shows that if a tensor is defective, then any subpencil of the tensor must be defective. Hence by examining any given subpencil, one may use Theorem \ref{theorem:TensMDBound} to compute the radius of a ball around a given tensor in which best $ {\mathbb{K}} $-rank $R$ approximations exist. A bound obtained from a single pencil may then be improved by combining information from a collection of non-overlapping subpencils of the tensor. In this way one arrives at the following multiple pencil based bound for the existence of a best $ {\mathbb{K}} $-rank $R$ approximation.

\begin{theorem}
\label{theorem:MultiplePencilBound}
Let ${\mathcal T}' \in  {\mathbb{K}} ^{I_1 \times I_2 \times I_3}$ be a tensor of $ {\mathbb{K}} $-rank $R$ and multilinear rank $(R,R,K)$, and let ${\mathcal T} \in  {\mathbb{K}} ^{R \times R \times K}$ be an orthogonal compression of ${\mathcal T}'$. Let $ {\mathbf{U}}  \in  {\mathbb{K}} ^{K \times K}$ be a unitary matrix and set ${\mathcal S} = {\mathcal T} \cdot_3  {\mathbf{U}} $. For each $i=1,\dots, \lfloor K/2 \rfloor$, let $\epsilon_i \geq 0$ be a nonnegative real number such that every matrix pencil in an open Frobenius norm ball of radius $\epsilon_i$ around the pencil
\[
( {\mathbf{S}} _{2i-1}, {\mathbf{S}} _{2i})
\]
is simple and set $\epsilon = ||(\epsilon_1, \dots, \epsilon_{ \lfloor K/2 \rfloor})||_2$. Then every tensor in an open ball of Frobenius radius $\epsilon/2$ centered at ${\mathcal T}$ has a best $ {\mathbb{K}} $-rank $R$ approximation and every rank $R$ tensor inside the open ball of Frobenius radius $\epsilon$ centered at ${\mathcal T}$ has a unique CPD. Furthermore, the set of tensors inside this open ball for which the best rank $R$ approximation is not unique has Lebesgue measure zero.

In particular, each $\epsilon_i$ may be computed using Theorem \ref{theorem:TensMDBound}.
\end{theorem}
\begin{proof}
See Section \ref{sec:MultiplePencilBoundProof}.
\end{proof}

\begin{remark}
As an alternative to Theorem \ref{theorem:TensMDBound}, one could use \cite[Corollary 3]{DK87} to compute the  individual $\epsilon_i$ appearing in Theorem \ref{theorem:MultiplePencilBound}. We observe in experiments that the  $\epsilon_i$ computed by Theorem \ref{theorem:TensMDBound} are typically about the same or slightly larger than those computed using \cite[Corollary 3]{DK87}. Additionally, it is much more computationally intensive to obtain the $\epsilon_i$ using \cite[Corollary 3]{DK87}. However, for particularly ill-conditioned tensors, \cite[Corollary 3]{DK87} may outperform Theorem \ref{theorem:TensMDBound}.
\end{remark}

Though Theorem \ref{theorem:MultiplePencilBound} explains how to compute a lower bound for the radius of a neighborhood around a given rank $R$ tensor in which best rank $R$ approximations are guaranteed to exist, this result may be difficult to directly use in applications as one will not have access to the low rank signal tensor one is interested in.

To remedy this difficulty, we present a variation on Theorem \ref{theorem:TensMDBound} which can be directly applied to (a low multilinear rank approximation of) a measured tensor. This variation, Theorem \ref{theorem:MLSVDExistenceBound}, only requires access to the measured tensor to use.

\begin{theorem}
\label{theorem:MLSVDExistenceBound}

Given ${{\mathcal M}}' \in  {\mathbb{K}} ^{I_1 \times I_2 \times I_3}$ and an integer $R \leq \min\{I_1,I_2\}$, set $K=\min\{R_3({{\mathcal M}}'),R\}$. Let ${\mathcal W}'$ be a multilinear rank $(R,R,K)$ approximation of ${{\mathcal M}}'$ and let ${\mathcal W} \in  {\mathbb{K}} ^{R \times R \times K}$ be an orthogonal compression of ${\mathcal W}'$. Also let $ {\mathbf{U}}  \in  {\mathbb{K}} ^{K \times K}$ be a unitary matrix and set ${\mathcal S} = {\mathcal W} \cdot_3  {\mathbf{U}} $. For each $i=1,\dots, \lfloor K/2 \rfloor$, let $\epsilon_i \geq 0$ be a nonnegative real number such that every matrix pencil in an open Frobenius norm ball of radius $\epsilon_i$ around the pencil
\[
( {\mathbf{S}} _{2i-1}, {\mathbf{S}} _{2i})
\]
is simple and set $\epsilon = ||(\epsilon_1, \dots, \epsilon_{ \lfloor K/2 \rfloor})||_2$. If there exists some $ {\mathbb{K}} $-rank $R$ tensor ${\tilde{\cT}}' \in  {\mathbb{K}} ^{I_1 \times I_2 \times I_3} $ such that
\begin{equation}
\label{eq:MLSVDexistIneq}
\|{{\mathcal M}}'-{\tilde{\cT}}'\|_\mathrm{F} < \epsilon - \|{{\mathcal M}}'-{\mathcal W}'\|_\mathrm{F}
\end{equation}
then ${{\mathcal M}}'$ has a best $ {\mathbb{K}} $-rank $R$ approximation and  any best $ {\mathbb{K}} $-rank $R$ approximation of ${{\mathcal M}}'$ has a unique CPD. Furthermore, the set of tensors ${{\mathcal M}}'$ which satisfy these assumptions and do not have a unique best $ {\mathbb{K}} $-rank $R$ approximation has Lebesgue measure zero.

In particular, each $\epsilon_i$ may be computed using Theorem \ref{theorem:TensMDBound}.
\end{theorem}
\begin{proof}
The proof of Theorem \ref{theorem:MLSVDExistenceBound} is nearly identical to the proof of Theorem \ref{theorem:MultiplePencilBound}. For completeness, details of the proof are given in the supplementary materials.
\end{proof}

Theorem \ref{theorem:MLSVDExistenceBound} is most naturally applicable to tensors which have low or approximately low multilinear rank. Indeed if this theorem is applied to a tensor which does not approximately have low multilinear rank, then the right hand side of equation \eqref{eq:MLSVDexistIneq} can be negative in which case no conclusion can be made.

However, as previously discussed, a common approach to computing a CPD in application is to first compute a core ${\mathcal W}$ of a low multilinear rank approximation of ${{\mathcal M}}'$, then to compute a low rank CPD approximation of ${\mathcal W}$. In this approach, a critical question to answer is whether or not the core tensor ${\mathcal W}$ has a best low rank approximation. To address this question, Theorem \ref{theorem:MLSVDExistenceBound} can be directly applied to ${\mathcal W}$. Since ${\mathcal W}$ already has low multilinear rank, if one uses Theorem \ref{theorem:MLSVDExistenceBound} in this way, then inequality \eqref{eq:MLSVDexistIneq} simplifies to checking if there is some $ {\mathbb{K}} $-rank $R$ tensor ${\tilde{\cT}} \in  {\mathbb{K}} ^{R \times R \times K}$ such that
\[
\|{\mathcal W}-{\tilde{\cT}}\|_\mathrm{F} < \epsilon
\]
where $\epsilon$ is computed as described in the theorem statement.
The tensor ${\tilde{\cT}}$ used in this theorem can for instance be obtained by applying any standard CPD algorithm to ${\tilde{\cT}}$. The resulting approximation error being bounded by $\epsilon$ is sufficient to make sure the approximation problem for ${\mathcal W}$ is well-posed. The fact that ${\tilde{\cT}}$ itself may not be optimal is not an issue.

Numerical experiments which illustrate Theorem \ref{theorem:MLSVDExistenceBound} are given in Section \ref{sec:existenceImagesBigDim}.

\subsection{Deterministic bounds when $K=2$}
The following theorem explains how to compute the $\epsilon_i$ appearing in Theorems \ref{theorem:MultiplePencilBound} and \ref{theorem:MLSVDExistenceBound}. Before stating the result we give three brief definitions. Given two one-dimensional subspaces ${{\mathscr L}}_1,{{\mathscr L}}_2 \subseteq  {\mathbb{K}} ^K$, the \df{chordal metric}\footnote{The chordal metric is a popular choice of metric for studying generalized eigenvalues as it accounts for the natural scaling invariance in the generalized eigenvalue problem. See e.g. \cite[Chapter VI]{SS90}.} between ${{\mathscr L}}_1$ and ${{\mathscr L}}_2$, denoted $\chi({{\mathscr L}}_1,{{\mathscr L}}_2)$, is equal to the sine of the angle between ${{\mathscr L}}_1$ and ${{\mathscr L}}_2$.

Next, for $i=1,2$ let ${\mathcal T}_i \in  {\mathbb{K}} ^{R \times R \times K}$ be a tensor which has $R$ JGE values counting algebraic multiplicity and let $\{{{\mathscr L}}_{i,r}\}_{r=1}^R$ be the spectrum of ${\mathcal T}_i$. Define the \df{matching distance} between ${\mathcal T}_1$ and ${\mathcal T}_2$, denoted ${\mathrm{md}}[{\mathcal T}_1,{\mathcal T}_2]$, by
\begin{equation}
\label{eq:MatchingDistanceDef}
{\mathrm{md}}[{\mathcal T}_1,{\mathcal T}_2]= \min_{\pi \in
S_R} \max_{k=1, \dots, K} \chi({{\mathscr L}}_{1,k}, {{\mathscr L}}_{2,\pi (k)}).
\end{equation}
Here $S_R$ is the group of permutations of $\{1, \dots, R\}$. We note that the matching distance is a metric between the spectrum ${\mathcal T}_1$ and the spectrum of ${\mathcal T}_2$.

Finally, for a tensor ${{\mathcal M}}' \in K^{I_1 \times I_2 \times I_3}$, define the \df{spectral norm} of ${{\mathcal M}}'$, denoted $\spnorm{{{\mathcal M}}'}$  by
\[
\|{{\mathcal M}}'\|_{\mathrm{sp}}:= \max_{\| {\vec{x}} \|_2=\| {\vec{y}} \|_2=\| {\vec{z}} \|_2=1} |{{\mathcal M}}' \cdot_1  {\vec{x}}  \cdot_2  {\vec{y}}  \cdot_3  {\vec{z}} |
\]
where the vectors $ {\vec{x}} , {\vec{y}} , {\vec{z}} $ are of compatible dimension. Equivalently, the spectral norm of ${{\mathcal M}}'$ is the largest singular value of ${{\mathcal M}}'$ which is in turn given by the Frobenius norm of a best rank $1$ approximation of ${{\mathcal M}}'$. See \cite{FL18,FMPS13,FO14,QHX21} for more information about best rank $1$ approximations of tensors, tensor singular values, and the spectral norm of tensors.

\begin{theorem}
\label{theorem:TensMDBound}
Let ${\mathcal T}$ and ${\mathcal W}$ be $R \times R \times 2$ tensors with entries in $ {\mathbb{K}} $. Assume that ${\mathcal T}$ has $ {\mathbb{K}} $-rank $R$ with CPD ${\cp{ {\mathbf{A}} , {\mathbf{B}} , {\mathbf{C}} }}$ where the columns of $ {\mathbf{C}} $ are normalized to have unit Euclidean norm. Let $\{{{\mathscr L}}_r\}_{r=1}^R$ be the spectrum of ${\mathcal T}$. If
\[
\spnorm{{\mathcal T}-{\mathcal W}} < \frac{\sigma_{\min} ( {\mathbf{A}} ) \sigma_{\min} ( {\mathbf{B}} ) \min_{i \neq j} \chi ({{\mathscr L}}_i, {{\mathscr L}}_j)}{2},
\]
then ${\mathcal W}$ is a  simple  tensor of $ {\mathbb{K}} $-rank $R$ and satisfies
\begin{equation}
\label{eq:MatchingDistanceBound}
{\mathrm{md}}[{\mathcal T},{\mathcal W}]  \leq \frac{\spnorm{{\mathcal T}-{\mathcal W}}}{\sigma_{\min} ( {\mathbf{A}} ) \sigma_{\min} ( {\mathbf{B}} )}.
\end{equation}
\end{theorem}

Since the Frobenius norm is an upper bound for the spectral norm of a tensor, Theorem \ref{theorem:TensMDBound} illustrates that if the pencil $( {\mathbf{S}} _{2i-1}, {\mathbf{S}} _{2i})\cong {\mathcal S}_i$ appearing in Theorem \ref{theorem:MultiplePencilBound} has CPD ${\mathcal S}_i= [\![ {\mathbf{A}} _i, {\mathbf{B}} _i, {\mathbf{C}} _i]\!]$ where the columns of $ {\mathbf{C}} _i$ have unit Euclidean norm, then one may take
\[
\epsilon_i = \frac{\sigma_{\min} ( {\mathbf{A}} _i ) \sigma_{\min} ( {\mathbf{B}} _i ) \min_{k \neq j} \chi ( {\mathbf{C}} _i(:,k),  {\mathbf{C}} _i(:,j))}{2}.
\]
Here $ {\mathbf{C}} _i(:,k)$ is the $k$th column of $ {\mathbf{C}} _i$. Note that one can use a generalized eigenvalue decomposition to determine if an $R \times R \times 2$ tensor is slice mix invertible and has rank $R$ and, if so, compute a CPD of the tensor. This approach is often called Jennrich's algorithm \cite{Harsh70}.

In Theorem \ref{theorem:TensBauerFike} we give a perturbation theoretic bound for the \textit{joint} generalized eigenvalues of tensors ${\mathcal T}$ and ${\mathcal W}$ of size $R \times R \times K$. However, the measurement of variation between the joint generalized spectra of ${\mathcal T}$ and ${\mathcal W}$ used in Theorem \ref{theorem:TensBauerFike} is notably weaker than the matching distance which is used in Theorem \ref{theorem:TensMDBound}.

\begin{remark}
In the language of the classical generalized eigenvalue problem and matrix pencils, Theorem \ref{theorem:TensMDBound} gives a bound which guarantees that the matrix pencil $( {\mathbf{W}} _1, {\mathbf{W}} _2)$ is diagonalizable over $ {\mathbb{K}} $.
\end{remark}

\section{Basic results for the joint generalized eigenvalue problem}
\label{sec:JGEbasics}

We now develop several basic results for joint generalized eigenvalues. These results help reformulate CPD in terms of the joint generalized eigenvalue problem. This reformulation allows us to apply our upcoming perturbation theoretic results to the CPD setting.

\subsection{A simplifying assumption}

As a simplifying assumption throughout the section we restrict to the case where the tensor ${\mathcal T}$ of interest has size $R \times R \times K$ for some $K \leq R$. In addition, we typically assume that ${\mathcal T}$ is slice mix invertible. The assumption that ${\mathcal T}$ is slice mix invertible is closely related to the tensor having multilinear rank $(R,R,K)$ and rank $R$, see Lemma \ref{lemma:invertslice}.

As described in Section \ref{sec:Compress}, orthogonal compressions preserve geometry, so the assumption that ${\mathcal T}$ has multilinear rank $R \times R \times K$ is not restrictive. The assumption that ${\mathcal T}$ has a linear combination of frontal slices which is invertible holds for tensors with generic entries.

We remark that, without assuming ${\mathcal T} \in  {\mathbb{K}} ^{R \times R \times K}$ has rank $R$, the assumption that ${\mathcal T}$ has multilinear rank $(R,R,K)$ is not sufficient for ${\mathcal T}$ to be slice mix invertible. In fact, it is possible for a tensor to have multilinear rank $(R,R,K)$ and border $ {\mathbb{K}} $-rank $R$ and not be slice mix invertible. This is illustrated by the following example.

\begin{exa}
\label{exa:pequiv0}
Let ${\mathcal T} \in  {\mathbb{R}} ^{3 \times 3 \times 3}$ be the tensor with frontal slices
\[
 {\mathbf{T}} _1= \begin{pmatrix}
1 & 0 & 0 \\
0 & 0 & 0 \\
0 & 0 & 1
\end{pmatrix}
\quad
 {\mathbf{T}} _2= \begin{pmatrix}
0 & 2 & 0 \\
0 & 0 & -2 \\
0 & 0 & 0
\end{pmatrix}
\quad
 {\mathbf{T}} _3= \begin{pmatrix}
0 & 0 & -1 \\
0 & 0 & 0 \\
0 & 0 & 0
\end{pmatrix}.
\]
Then ${\mathcal T}$ has multilinear rank $(3,3,3)$ and border $ {\mathbb{R}} $-rank $3$, but ${\mathcal T}$ is not slice mix invertible. That is, the span of the frontal slices of ${\mathcal T}$ does not contain an invertible matrix. To see that ${\mathcal T}$ has border $ {\mathbb{R}} $-rank $3$, for each $n$ define the factor matrices
\[
 {\mathbf{A}} _n = \begin{pmatrix}
1 & 1 & 1 \\ 0 & \frac{1}{n} & \frac{1}{n^4} \\
0 & 0 & \frac{1}{n^2}
\end{pmatrix}
\quad
 {\mathbf{B}} _n = \begin{pmatrix}
1 & 0 & 0 \\ -\frac{1}{n^2} & n & 0 \\ -n^2 & -n^2 & n^2
\end{pmatrix}
\quad
 {\mathbf{C}} _n = \begin{pmatrix}
1 & \frac{1}{n^3} & 1 \\ \frac{4}{n} & \frac{2}{n} & \frac{6}{n} \\
\frac{1}{n^2} & 0 & 0
\end{pmatrix}
\]
and let ${\mathcal T}^n$ be the $ {\mathbb{K}} $-rank $3$ tensor ${\mathcal T}^n = [\![ {\mathbf{A}} _n,  {\mathbf{B}} _n,  {\mathbf{C}} _n]\!]$. Then $\lim {\mathcal T}^n = {\mathcal T}$, hence the border $ {\mathbb{K}} $--rank of ${\mathcal T}$ is at most $3$. That ${\mathcal T}$ has border $ {\mathbb{K}} $--rank $3$ then follows from the fact that ${\mathcal T}$ has multilinear rank $(3,3,3)$  together with equation \eqref{eq:BRankVsMLRank}. As a consequence of the forthcoming Lemma \ref{lemma:invertslice}, ${\mathcal T}$ has $ {\mathbb{K}} $-rank strictly greater than $3$.
\end{exa}

\begin{lemma}
\label{lemma:invertslice}
Let ${\mathcal T} \in  {\mathbb{K}} ^{R \times R \times K}$ be a tensor with $ {\mathbb{K}} $-rank $R$ and full multilinear rank, and let ${\mathcal T}={\cp{ {\mathbf{A}} , {\mathbf{B}} , {\mathbf{C}} }}$ be a CPD of ${\mathcal T}$. Then ${\mathcal T}$ is slice mix invertible and the factor matrices $ {\mathbf{A}} , {\mathbf{B}}  \in  {\mathbb{K}} ^{R \times R}$ are invertible.
\end{lemma}
\begin{proof}
Assuming that ${\mathcal T}$ has multilinear rank $(R,R,K)$ implies that each mode-$i$ unfolding of ${\mathcal T}$ has full row rank. By choosing an appropriate ordering of columns, the mode-$1$ unfolding of ${\mathcal T}$ is given by
\[
 {\mathbf{A}}  ( {\mathbf{C}}  \odot  {\mathbf{B}} )^\mathrm{T} \in  {\mathbb{K}} ^{R \times RK},
\]
where $\odot$ denotes the Khatri-Rao (column-wise Kronecker) product, e.g. see \cite{Cetal15}. ${\mathcal T}$ having multilinear rank $(R,R,K)$ implies that the mode-$1$ unfolding has rank $R$ from which it follows that $ {\mathbf{A}}  \in  {\mathbb{K}} ^{R \times R}$ has rank $R$ and is invertible. A similar argument shows that $ {\mathbf{B}} $ is invertible.

We now show that ${\mathcal T}$ is slice mix invertible. Using $ {\mathbf{T}} _i =  {\mathbf{A}}  D_i ( {\mathbf{C}} )  {\mathbf{B}} ^\mathrm{T}$ where $D_i ( {\mathbf{C}} )$ is the diagonal matrix whose entries are given by the $i$th row of $ {\mathbf{C}} $, to show that ${\mathcal T}$ has a linear combination of frontal slices which is invertible, it is sufficient to show that there is a linear combination of rows of $ {\mathbf{C}} $ which has no entry equal to zero. To this end, note that since ${\mathcal T}$ has $ {\mathbb{K}} $-rank $R$, the factor matrix $ {\mathbf{C}} $ cannot have a column with all entries equal to zero. Therefore, for each $r=1,\dots,R$, the set of $\alpha \in  {\mathbb{K}} ^K$ such that
\[
\sum_{k=1}^K  \alpha_k  {\mathbf{C}} (r,k)=0
\]
is a nontrivial Zariski closed subset of $ {\mathbb{K}} ^K$. It follows that for $\alpha \in  {\mathbb{K}} ^K$ outside of a set of measure zero, one has
\[
\sum_{k=1}^K  \alpha_k  {\mathbf{C}} (r,k) \neq 0 \mathrm{\ for \ all \ } r=1,\dots,R
\]
from which the result follows.
\end{proof}

\subsection{Distinct JGE values correspond to linearly independent JGE vectors}

A perhaps unsurprising yet important fact is that a set of joint generalized eigenvectors corresponding to a set of distinct joint generalized eigenvalues is linearly independent.

\begin{lemma}
\label{lemma:DistinctENsImpliesEigenbasis}
Let ${\mathcal T} \in  {\mathbb{K}} ^{R \times R \times K}$ and assume ${\mathcal T}$ is slice mix invertible. Let $\{({{\mathscr L}}_j, {\vec{x}} _j)\}_{j=1}^\ell$ be a collection of JGE pairs of ${\mathcal T}$ such that the ${{\mathscr L}}_j$ are all distinct. Then $\{ {\vec{x}} _j\}_{j=1}^\ell$ is a linearly independent set. As a consequence, if ${\mathcal T}$ is simple and each of its JGE values has geometric multiplicity at least one, then ${\mathcal T}$ has a JGE basis.
\end{lemma}
\begin{proof}
 We treat the problem by reducing to the $K=2$ setting. To this end, for $j=1,\dots,\ell$ let $ {\vec{\lambda}} _j \in  {\mathbb{K}} ^K$ be a vector such that $\mathrm{span}( {\vec{\lambda}} _j)={{\mathscr L}}_j$, and let $ {\mathbf{L}}  \in  {\mathbb{K}} ^{K \times \ell}$ be the matrix whose $j$th column is equal to $ {\vec{\lambda}} _j$. Since the ${{\mathscr L}}_j$ are distinct, the matrix $ {\mathbf{L}} $ has Kruskal rank at least $2$. That is, no pair of columns of $ {\mathbf{L}} $ is linearly dependent. It follows that there exists a matrix $ {\mathbf{V}}  \in  {\mathbb{K}} ^{2 \times K}$ such that $ {\mathbf{V}}   {\mathbf{L}} $ has Kruskal rank $2$.

Set ${\mathcal S} = {\mathcal T}  \cdot_3  {\mathbf{V}}  \in  {\mathbb{K}} ^{R \times R \times 2}$ and let $ {\mathbf{S}} _1$ and $ {\mathbf{S}} _2$ be the frontal slices of ${\mathcal S}$. It is straightforward to show that each $ {\vec{x}} _j$ is a generalized eigenvector of the matrix pencil $( {\mathbf{S}} _1,  {\mathbf{S}} _2)$ which corresponds to the generalized eigenvalue $\mathrm{span}( {\mathbf{V}}   {\vec{\lambda}} _j)$. Furthermore, since the matrix $ {\mathbf{V}}   {\mathbf{L}} $ has Kruskal rank $2$, the generalized eigenvalues $\{\mathrm{span}( {\mathbf{V}}   {\vec{\lambda}} _j)\}_{j=1}^\ell$ are all distinct. As a consequence of \cite[Chapter VI, Theorem 1.11]{SS90}, distinct generalized eigenvalues correspond to linearly independent eigenvectors, so it follows that $\{ {\vec{x}} _j\}_{j=1}^\ell$ is a linearly independent set.
\end{proof}

It is possible for a joint generalized eigenvalue of a slice mix invertible tensor to have nonzero algebraic multiplicity while having geometric multiplicity equal to zero, hence, the assumption that the joint generalized eigenvalues of ${\mathcal T}$ have nonzero geometric multiplicity in the second part of of Lemma \ref{lemma:DistinctENsImpliesEigenbasis} is necessary.

\begin{exa}
\label{exa:geomultzero}
Let ${\mathcal T} \in  {\mathbb{R}} ^{3 \times 3 \times 3}$ be the tensor with frontal slices
\[
 {\mathbf{T}} _1= \begin{pmatrix}
1 & 0 & 0 \\
1 & 1 & 0 \\
3 & 4 & 1
\end{pmatrix}
\quad
 {\mathbf{T}} _2= \begin{pmatrix}
1 & 0 & 0 \\
1 & 2 & 0 \\
4 & 5 & -1
\end{pmatrix}
\quad
 {\mathbf{T}} _3= \begin{pmatrix}
1 & 0 & 0 \\
5 & 3 & 0 \\
7 & 8 & 5
\end{pmatrix}.
\]
Then
\[
p_{\mathcal T}(\gamma)=\det \left(\sum_{k=1}^3 \gamma_k  {\mathbf{T}} _k \right)= (\gamma_1 + \gamma_2 + \gamma_3) (\gamma_1 + 2 \gamma_2 + 3 \gamma_3) (\gamma_1 - \gamma_2 + 5 \gamma_3)
\]
hence ${\mathcal T}$ is slice mix invertible and has three distinct JGE values counting algebraic multiplicity, namely $\mathrm{span} ((1,1,1))$, $\mathrm{span} ((1,2,3)),$ and $\mathrm{span} ((1,-1,5))$. However, the JGE values $\mathrm{span} ((1,1,1))$ and $\mathrm{span} ((1,2,3))$ do not have a corresponding JGE vector.
\end{exa}

In Example \ref{exa:geomultzero}, one could reasonably consider left joint generalized eigenvalues of ${\mathcal T}$ in which case the JGE value $\mathrm{span} ((1,1,1))$ in Example \ref{exa:geomultzero} corresponds to the left JGE vector $(1,0,0) \in  {\mathbb{R}} ^3$. However, the JGE values $\mathrm{span} ((1,-1,5))$ and $\mathrm{span} ((1,2,3))$ do not have left JGE vectors. The JGE value $\mathrm{span} ((1,2,3))$ is of particular note since it neither has left nor right JGE vectors.

Assuming that a tensor ${\mathcal T}$ is slice mix invertible, it is not difficult to show that ${\mathcal T}$ has a basis of right JGE vectors if and only if it has a basis of left JGE vectors. However, this fact and left JGE vectors in general have little impact on our results, so we do not discuss them further.

\subsection{Algebraic and geometric multiplicities}

We now more carefully examine algebraic and geometric multiplicities for joint generalized eigenvalues. In particular we show that algebraic multiplicity is always greater than or equal to geometric multiplicity of a joint generalized eigenvalue. As an immediate consequence, if the tensor ${\mathcal T} \in  {\mathbb{K}} ^{R \times R \times K}$ has a JGE pair $({{\mathscr L}}, {\vec{x}} )$, then $\mathrm{am}({{\mathscr L}}) \geq 1$.

Before stating our result we recall that if ${\mathcal T}$ is not slice mix invertible, then $p_{{\mathcal T}}$ is identically zero in which case $\langle  {\vec{\lambda}} , {\boldsymbol{\gamma}}  \rangle^m$ divides $p_{{\mathcal T}} ( {\boldsymbol{\gamma}} )$ for all $ {\vec{\lambda}}  \in  {\mathbb{K}} ^K$ and all $m \in \mathbb{N}$. In this case one can say that $\mathrm{span}( {\vec{\lambda}} )$ has infinite algebraic multiplicity as a JGE value of ${\mathcal T}$ for all $ {\vec{\lambda}}  \in  {\mathbb{K}} ^K$. For this reason, in the upcoming proposition it is sufficient to restrict to the case where ${\mathcal T}$ is slice mix invertible.

\begin{proposition}
\label{proposition:alggreatergeo}
Let ${\mathcal T} \in  {\mathbb{K}} ^{R \times R \times K}$ and assume ${\mathcal T}$ is slice mix invertible. Then for any nonzero $ {\vec{\lambda}}  \in  {\mathbb{K}} ^K$ one has $\mathrm{am}(\mathrm{span}( {\vec{\lambda}} )) \geq \mathrm{gm}(\mathrm{span}( {\vec{\lambda}} ))$. In particular, if $(\mathrm{span}( {\vec{\lambda}} ), {\vec{x}} )$ is a JGE pair of ${\mathcal T}$, then there is an integer $m>0$ such that the characteristic polynomial of ${\mathcal T}$ factors as
\begin{equation}
\label{eq:gmamFactor}
p_{\mathcal T} ( {\boldsymbol{\gamma}} ) = \langle  {\vec{\lambda}} , \overline{ {\boldsymbol{\gamma}} } \rangle^m g( {\boldsymbol{\gamma}} )
\end{equation}
where $g$ is a nonzero polynomial in the indeterminate $ {\boldsymbol{\gamma}} $.

Furthermore, if ${\mathcal T}$ has  border $ {\mathbb{K}} $-rank $R$, then $p_{\mathcal T}( {\boldsymbol{\gamma}} )$ factors into a product of $R$ linear terms. That is, one has
\[
p_{\mathcal T} ( {\boldsymbol{\gamma}} ) = \Pi_{r=1}^R \langle  {\vec{\lambda}} _r, \overline{ {\boldsymbol{\gamma}} } \rangle
\]
hence ${\mathcal T}$ has $R$ JGE values counting algebraic multiplicity.
\end{proposition}
\begin{proof}

Suppose $\mathrm{gm}(  {\vec{\lambda}} ) = m$ and let $ {\vec{x}} _1, \dots,  {\vec{x}} _m$ be the corresponding JGE vectors. For $i=1, \dots, m$ let $ {\vec{y}} _i$ be the nonzero vector corresponding to the eigenpair $(\mathrm{span}( {\vec{\lambda}} ), {\vec{x}} _i)$, as in equation \eqref{eq:JGEVdef}. Note that the assumption that ${\mathcal T}$ is slice mix invertible guarantees that $\{ {\vec{y}} _1, \dots,  {\vec{y}} _m\}$ is a linearly independent set. To show that $p_{{\mathcal T}}$ factors as desired, we will use the $ {\vec{x}} _i$ and $ {\vec{y}} _i$ to construct matrices $ {\mathbf{X}} $ and $ {\mathbf{Q}} $ which simultaneously block upper triangularize the frontal slices $ {\mathbf{T}} _k$ with $\lambda_k  {\mathbf{I}} _m$ on the first diagonal block.

To simultaneously block upper diagonalize the $ {\mathbf{T}} _k$, let $ {\mathbf{X}} _1 \in  {\mathbb{R}} ^{R \times m}$ be a matrix whose $i$th column is $ {\vec{x}} _i$, and let $ {\mathbf{X}} _2 \in  {\mathbb{R}} ^{R \times (R-m)}$ be any matrix such that the block matrix $ {\mathbf{X}} = ( {\mathbf{X}} _1 \  {\mathbf{X}} _2)$ is invertible. Similarly let $ {\mathbf{Y}} _1  \in  {\mathbb{R}} ^{R \times m}$ be the matrix whose $i$th column is $ {\vec{y}} _i$, and let $ {\mathbf{Q}} _2 \in  {\mathbb{R}} ^{(R-m) \times R}$ be a matrix whose rows form a basis for the orthogonal complement of the range of $ {\mathbf{Y}} _1$. Set
\[
 {\mathbf{Q}}  =\begin{pmatrix}
 {\mathbf{Y}} _1^\dagger  \\
 {\mathbf{Q}} _2
\end{pmatrix} .
\]

With this set up, for each $k=1, \dots, K$ one has
\[
 {\mathbf{Q}}   {\mathbf{T}} _k  {\mathbf{X}}  =  {\mathbf{Q}}  ( {\lambda} _k  {\mathbf{Y}} _1 \  {\mathbf{T}} _k  {\mathbf{X}} _2) = \begin{pmatrix}
 {\lambda} _k   {\mathbf{Y}} _1^\dagger  {\mathbf{Y}} _1 &  {\mathbf{Y}} _1^\dagger  {\mathbf{T}} _k  {\mathbf{X}} _2 \\
\lambda_k  {\mathbf{Q}} _2  {\mathbf{Y}} _1 &   {\mathbf{Q}} _2  {\mathbf{T}} _k  {\mathbf{X}} _2
\end{pmatrix} = \begin{pmatrix}
 {\lambda} _k   {\mathbf{I}} _m &  {\mathbf{Y}} _1^\dagger  {\mathbf{T}} _k  {\mathbf{X}} _2 \\
\mathbf{0} &  {\mathbf{Q}} _2  {\mathbf{T}} _k  {\mathbf{X}} _2
\end{pmatrix} .
\]
The matrices $ {\mathbf{X}} $ and $ {\mathbf{Q}} $ are both invertible by construction so using
\[
\det\left( {\mathbf{Q}} \right) \det \left( \sum_{k=1}^K \gamma_k  {\mathbf{T}} _k \right) \det( {\mathbf{X}} )  =\det \left( \sum_{k=1}^K \gamma_k  {\mathbf{Q}}   {\mathbf{T}} _k  {\mathbf{X}}  \right) =
\det \left( \sum_{k=1}^K \gamma_k \begin{pmatrix}
 {\lambda} _k   {\mathbf{I}} _m &  {\mathbf{Y}} _1^\dagger  {\mathbf{T}} _k  {\mathbf{X}} _2 \\
\mathbf{0} &  {\mathbf{Y}} _2  {\mathbf{T}} _k  {\mathbf{X}} _2
\end{pmatrix} \right)
\]
shows that $p_{\mathcal T} ( {\boldsymbol{\gamma}} )$ factors as in equation \eqref{eq:gmamFactor}. The fact that the polynomial $g( {\boldsymbol{\gamma}} )$ is nonzero is a consequence of the assumption that ${\mathcal T}$ is slice mix invertible, hence $p_{\mathcal T}$ is not identically zero.

To complete the proof, first assume that ${\mathcal T}$ has $ {\mathbb{K}} $-rank $R$ and let ${\mathcal T}={\cp{ {\mathbf{A}} , {\mathbf{B}} , {\mathbf{C}} }}$ be the CPD of ${\mathcal T}$. We have assumed ${\mathcal T}$ is slice mix invertible, so Lemma \ref{lemma:invertslice} shows that $ {\mathbf{A}} $ and $ {\mathbf{B}} $ are invertible. Using $ {\mathbf{T}} _k =  {\mathbf{A}}  D_k ( {\mathbf{C}} )  {\mathbf{B}} ^\mathrm{T} $ shows that
\[
p_{\mathcal T} ( {\boldsymbol{\gamma}} )= \det \left(  {\mathbf{A}}  \sum_{k=1}^K \gamma_k D_k( {\mathbf{C}} )  {\mathbf{B}} ^\mathrm{T} \right) = \det ( {\mathbf{A}} )  \det \left(\sum_{k=1}^K \gamma_k D_k( {\mathbf{C}} ) \right) \det(  {\mathbf{B}} ^\mathrm{T}).
\]
The final term in this equality easily factors as desired.

 Now assume that ${\mathcal T}$ has border $ {\mathbb{K}} $-rank $R$. Then ${\mathcal T}=\lim {\mathcal T}^n$ for some sequence $\{{\mathcal T}^n\}_{n=1}^\infty$ of $ {\mathbb{K}} $-rank $R$ tensors. Since ${\mathcal T}$ is slice mix invertible, we may assume that each ${\mathcal T}^n$ is slice mix invertible. Thus, as we have shown above, for each $n$ we have
\[
p_{{\mathcal T}^n}  ( {\boldsymbol{\gamma}} ) = \alpha_n \Pi_{r=1}^R \langle  {\vec{\lambda}} _{r,n}, \overline{ {\boldsymbol{\gamma}} } \rangle.
\]
Here the constant $\alpha_n \in  {\mathbb{K}} $ has been added so that we may scale to have $|| {\vec{\lambda}} _{r,n}||_2=1$ for each $r=1,\dots, R$ and for all $n \in  {\mathbb N}$.

Since all the $ {\vec{\lambda}} _{r,n}$ have norm $1$, by passing to a common subsequence if necessary, the sequence $\{ {\vec{\lambda}} _{r,n}\}_{n=1}^\infty$ converges to some unit norm $ {\vec{\lambda}} _r \in  {\mathbb{K}} ^K$ for each $r=1,\dots, R$. We next show that the constants $\alpha_n$ are bounded in norm. To this end note that there must exist some unit vector $ {\boldsymbol{\gamma}} _* \in  {\mathbb{K}} ^K$ such that
\[
 \Pi_{r=1}^R \langle   {\vec{\lambda}} _{r}, \overline{ {\boldsymbol{\gamma}} _*} \rangle \neq 0,
\]
as otherwise would imply that the union of the orthogonal complements of the $ {\vec{\lambda}} _r$ is dense in $ {\mathbb{K}} ^K$.

Since the ${\mathcal T}^n$ converge to ${\mathcal T}$ we must have $\lim_{n \to \infty} p_{{\mathcal T}^n} (\gamma) = p_{{\mathcal T}} (\gamma)$ from which it follows that
\[
p_{{\mathcal T}} ( {\boldsymbol{\gamma}} _*) = \lim_{n \to \infty} p_{{\mathcal T}^n} ( {\boldsymbol{\gamma}} _*) = \lim_{n \to \infty} \alpha_n \Pi_{r=1}^R \langle   {\vec{\lambda}} _{r,n}, \overline{ {\boldsymbol{\gamma}} _*} \rangle.
\]
Using this together with
\[
\lim_{n \to \infty} \Pi_{r=1}^R \langle  {\vec{\lambda}} _{r,n}, \overline{ {\boldsymbol{\gamma}} _*} \rangle =  \Pi_{r=1}^R \langle  {\vec{\lambda}} _{r}, \overline{ {\boldsymbol{\gamma}} _*} \rangle \neq 0
\]
shows that the sequence $\{\alpha_n\}_{n=1}^\infty$ converges to some $\alpha \in  {\mathbb{K}} $. It follows that for all $ {\boldsymbol{\gamma}}  \in  {\mathbb{K}} ^K$ we have
\[
\lim_{n\to \infty} p_{{\mathcal T}^n} ( {\boldsymbol{\gamma}} ) = \lim_{n\to \infty} \alpha_n \Pi_{r=1}^R \langle  {\vec{\lambda}} _{n,r}, \overline{ {\boldsymbol{\gamma}} } \rangle =  \alpha \Pi_{r=1}^R \langle  {\vec{\lambda}} _{r}, \overline{ {\boldsymbol{\gamma}} } \rangle = p_{{\mathcal T}} ( {\boldsymbol{\gamma}} ).
\]
Furthermore, since ${\mathcal T}$ is slice mix invertible, we must have $\alpha \neq 0$. We conclude that $p_{{\mathcal T}} ( {\boldsymbol{\gamma}} )$ factors as desired.
\end{proof}

\subsection{Multiplicities for border rank $R$ tensors}

We now examine geometric multiplicities of JGE values for slice mix invertible border $ {\mathbb{K}} $-rank $R$ tensors. Namely we show that JGE values of a slice mix invertible tensor with border $ {\mathbb{K}} $-rank $R$ have nonzero geometric multiplicity. This lemma plays an important role in the proof of Theorem \ref{theorem:rankvsmult} \eqref{it:NonderogatoryImpliesRankR}, as it allows one to conclude that a simple border $ {\mathbb{K}} $-rank $R$ tensor  has a basis of JGE vectors.

\begin{lemma}
\label{lemma:BorderRankGeo}
Let ${\mathcal T} \in  {\mathbb{K}} ^{R \times R \times K}$ be a slice mix invertible tensor and assume ${\mathcal T}$ has multilinear rank $(R,R,K)$ and border $ {\mathbb{K}} $-rank $R$. If ${{\mathscr L}} \subset  {\mathbb{R}} ^K$ is a JGE value of ${\mathcal T}$ with $\mathrm{am}({{\mathscr L}})=m \geq 1$, then $\mathrm{gm}({{\mathscr L}}) \geq 1$.
\end{lemma}

\begin{proof}
Let $\{{\mathcal T}^n\}_{n=1}^\infty$ be a sequence of $ {\mathbb{K}} $-rank $R$ tensors which converges to ${\mathcal T}$ and let $ {\vec{\lambda}}  \in  {\mathbb{K}} ^K$ be a unit vector such that $\mathrm{span}( {\vec{\lambda}} )={{\mathscr L}}$.  Proposition \ref{proposition:alggreatergeo} shows that $p_{\mathcal T} ( {\boldsymbol{\gamma}} )$ and all the $p_{{\mathcal T}^n} ( {\boldsymbol{\gamma}} )$ factor as products of linear terms.  Note that $\lim {\mathcal T}^n = {\mathcal T}$ implies that $\lim p_{{\mathcal T}^n} = p_{\mathcal T}$, so the normalized linear factors of the $p_{{\mathcal T}^n}$ converge to the linear factors of $p_{\mathcal T}$. It follows that there exists a sequence of unit vectors $\{ {\vec{\lambda}} _n\}_{n=1}^\infty \subset  {\mathbb{K}} ^K$ such that $\mathrm{span}( {\vec{\lambda}} _n)$ is a JGE value of ${\mathcal T}^n$ for each $n$ and such that
\[
\lim_{n \to \infty}  {\vec{\lambda}} _n =  {\vec{\lambda}} .
\]

As a consequence of Theorem \ref{theorem:rankvsmult}, each $\mathrm{span}( {\vec{\lambda}} _n)$ has geometric multiplicity at least $1$ as a JGE value of ${\mathcal T}^n$. For each $n$, let $ {\vec{x}} _n \in  {\mathbb{K}} ^R$ be a unit norm JGE vector of ${\mathcal T}^n$ corresponding to $ {\vec{\lambda}} _n$, and let $ {\vec{y}} _n \in  {\mathbb{K}} ^R$ be the nonzero vector such that the tuple $({\mathcal T}^n, {\vec{\lambda}} _n, {\vec{x}} _n, {\vec{y}} _n)$ satisfies equation \eqref{eq:JGEVdef} for each $n$.

By construction the  sequences $\{ {\vec{x}} _n\}$ and $\{ {\vec{\lambda}} _n\}$ lie in a compact sets. Additionally, since the ${\mathcal T}^n$ converge to ${\mathcal T}$, the norms of the frontal slices of ${\mathcal T}^n$ are uniformly bounded, so since the $ {\vec{\lambda}} _n$ all have norm $1$, the sequence $\{ {\vec{y}} _n\}$ must also lie in a compact set. Therefore, by passing to a common subsequence if necessary, the sequences $\{ {\vec{x}} _n\}$ and $\{ {\vec{y}} _n\}$ have limits $ {\vec{x}} , {\vec{y}}  \in  {\mathbb{K}} ^R$. For each $k=1,\dots K$ we then have
\begin{equation}
\label{eq:limJGEV}
 {\mathbf{T}} _k  {\vec{x}}  = \lim_{n \to \infty}  {\mathbf{T}} _k^n  {\vec{x}} _n =  \lim_{n \to \infty}  {\lambda} _{n,k}  {\vec{y}} ^n =  {\lambda} _k  {\vec{y}}
\end{equation}
where $ {\lambda} _{n,k}$ denotes the $k$th entry of $ {\vec{\lambda}} _n$.

By assumption, ${\mathcal T}$ is slice mix invertible, so the frontal slices of ${\mathcal T}$ cannot have a common null space and the vector $ {\vec{y}} $ cannot be equal to zero. We conclude that $({{\mathscr L}}, {\vec{x}} )$ is a joint generalized eigenpair of ${\mathcal T}$ and that $\mathrm{gm}({{\mathscr L}}) \geq 1$.
\end{proof}

\subsection{Proof of Theorem \ref{theorem:rankvsmult}}
\label{sec:rankvsmultproof}

We are now in position to give the proof of our first main result, Theorem \ref{theorem:rankvsmult}.

\begin{proof}

We first prove Item \eqref{it:RankIFFBasis}. Note that an orthogonal compression ${\mathcal T}$ of ${\mathcal T}'$ has $ {\mathbb{K}} $-rank $R$ if and only if ${\mathcal T}'$ has $ {\mathbb{K}} $-rank $R$, so throughout the proof we WLOG assume that ${\mathcal T}'$ has size $R \times R \times K$ and full multilinear rank, hence ${\mathcal T}={\mathcal T}'$.

First assume that ${\mathcal T}$ has $R$ JGE values in $ {\mathbb{K}} ^K$ counting geometric multiplicity. That is, assume that ${\mathcal T}$ has a basis of JGE vectors. Then there exists an invertible matrix $ {\mathbf{X}}  \in  {\mathbb{K}} ^{R \times R}$, a collection of diagonal matrices\footnote{The matrices $\{ {\vec{\Lambda}} _k\}_{k=1}^K$ are closely related to the JGE values of ${\mathcal T}$. In particular the subspace $\mathrm{span}(( {\vec{\Lambda}} _1 (r,r), {\vec{\Lambda}} _2 (r,r),\dots, {\vec{\Lambda}} _K (r,r)) \subset  {\mathbb{K}} ^K$  is a JGE value of ${\mathcal T}$ for each $r=1,\dots,R$. Here $ {\vec{\Lambda}} _k (r,r)$ denotes the $(r,r)$th entry of $ {\vec{\Lambda}} _k$.} $\{ {\vec{\Lambda}} _k\}_{k=1}^K$, and matrix $ {\mathbf{Y}}  \in  {\mathbb{K}} ^{R \times R}$ such that
\[
 {\mathbf{T}} _k  {\mathbf{X}}  =  {\mathbf{Y}}   {\vec{\Lambda}} _k \qquad \mathrm{for \ all \ } k=1, \dots, K.
\]
It follows that
\[
 {\mathbf{T}} _k =  {\mathbf{Y}}   {\vec{\Lambda}} _k  {\mathbf{X}} ^{-1} \qquad \mathrm{for \ all \ } k=1, \dots, K.
\]

Set $ {\mathbf{A}} = {\mathbf{Y}} $ and $ {\mathbf{B}} =  {\mathbf{X}} ^{-\mathrm{T}}$ and let $ {\mathbf{C}}  \in  {\mathbb{K}} ^{R \times K}$ be the matrix whose $k$th row is given by the diagonal entries of $ {\vec{\Lambda}} _k$. Then ${\cp{ {\mathbf{A}} , {\mathbf{B}} , {\mathbf{C}} }}={\mathcal T}$ is an $R$-term CPD for ${\mathcal T}$, hence ${\mathcal T}$ has rank less than or equal to $R$. That ${\mathcal T}$ has rank $R$ then follows from the assumption that ${\mathcal T}$ has multilinear rank $(R,R,K)$ and the fact that the rank of a tensor dominates its multilinear rank.

Now, similar to the proof of Proposition \ref{proposition:alggreatergeo}, assume that ${\mathcal T}$ has rank $R$ with CPD ${\mathcal T}= {\cp{ {\mathbf{A}} , {\mathbf{B}} , {\mathbf{C}} }}$. By assumption ${\mathcal T}$ has full multilinear rank, so Lemma \ref{lemma:invertslice} shows that the factor matrices $ {\mathbf{A}} , {\mathbf{B}}  \in  {\mathbb{K}} ^{R \times R}$ are invertible. Using $ {\mathbf{T}} _k =  {\mathbf{A}}  D_k ( {\mathbf{C}} )  {\mathbf{B}} ^\mathrm{T} $, we then have
\[
 {\mathbf{T}} _k  {\mathbf{B}} ^{-\mathrm{T}} =  {\mathbf{A}}  D_k ( {\mathbf{C}} ) \qquad \mathrm{for \ all \ } k=1, \dots, K.
\]

For $r=1,\dots, R$, let $ {\vec{x}} _r$ be the $r$th column of the matrix $ {\mathbf{B}} ^{-\mathrm{T}}$ and let $ {\vec{y}} _r$ be the $r$th column of $ {\mathbf{A}} $. Then $ {\vec{x}} _r, {\vec{y}} _r \in  {\mathbb{K}} ^K$ are nonzero vectors for all $r$ and for each $r=1,\dots,R$ and $k=1,\dots,K$ we have
\[
 {\mathbf{T}} _k  {\vec{x}} _r = c_{kr}  {\vec{y}} _r
\]
where $c_{kr}$ denotes the $(k,r)$ entry of $ {\mathbf{C}} $. Setting ${{\mathscr L}}_r = \mathrm{span} \big(c_{1r}, \dots, c_{Kr}\big)$, it follows that $({{\mathscr L}}_r,  {\vec{x}} _r)$ is a joint generalized eigenpair of ${\mathcal T}$ for each $r=1, \dots, R$. We conclude that ${\mathcal T}$ has a JGE basis, that the JGE values of ${\mathcal T}$ are given by the spans of the columns of $ {\mathbf{C}} $, and that the JGE vectors of ${\mathcal T}$ are given by the columns of $ {\mathbf{B}} ^{-\mathrm{T}}$, as claimed.

Item \eqref{it:NonderogatoryImpliesRankR} follows immediately from Item \eqref{it:RankIFFBasis} together with Lemma \ref{lemma:BorderRankGeo}. To prove Item \eqref{it:BRankLessRankImpliesDefective}, suppose ${\mathcal T}'$ has border $ {\mathbb{K}} $-rank $R$ and assume towards a contradiction that $\mathrm{am}({{\mathscr L}})=\mathrm{gm}({{\mathscr L}})$ for all JGE values ${{\mathscr L}}$ of ${\mathcal T}$. Proposition \ref{proposition:alggreatergeo} shows that $p_{\mathcal T}$ factors as a product of linear terms, so in this case ${\mathcal T}$ must have a basis of JGE vectors which using Item \eqref{it:RankIFFBasis} shows that ${\mathcal T}'$ has rank $R$, a contradiction. It follows that ${\mathcal T}$ has some JGE value with algebraic multiplicity strictly greater than geometric multiplicity. The proof is completed by Lemma \ref{lemma:BorderRankGeo} which shows that every JGE value of ${\mathcal T}$ has geometric multiplicity at least one.
\end{proof}

\section{Perturbation bounds for joint generalized eigenvalues}
\label{sec:perturbation}

We now examine how joint generalized eigenvalues are affected by a perturbation of the tensor of interest. Our main goal at this point is to provide a way to determine if border $ {\mathbb{K}} $-rank $R$ tensors near a given slice mix invertible $ {\mathbb{K}} $-rank $R$ tensor have distinct JGE values. That is, we want to show that nearby border $ {\mathbb{K}} $-rank $R$ tensors are nondefective hence have rank equal to border rank. To achieve this goal, we provide a Bauer--Fike like bound for JGE eigenvalues which may be converted to a spectral matching distance bound in the case $K=2$.

As in the previous section, it will be helpful to reduce to the case where ${\mathcal T}$ has size $R \times R \times K$ and full multilinear rank.
To make this reduction we consider a ``tensor Procrustes problem" which serves as an analogue of the classical matrix Procrustes problem \cite{S66}.

In particular we will show that given two tensors ${\mathcal W}',\hat{{\mathcal W}'} \in  {\mathbb{K}} ^{I_1 \times I_2 \times I_3}$, both having multilinear rank bounded above by $(R,R,K)$, there exist orthogonal compressions ${\mathcal W},\hat{{\mathcal W}} \in  {\mathbb{K}} ^{R \times R \times K}$ of ${\mathcal W}'$ and $\hat{{\mathcal W}'}$ such that
\begin{equation}
\label{eq:ProcrustesEq}
\|{\mathcal W} - \hat{{\mathcal W}} \|_\mathrm{F} \leq \| {\mathcal W}' - \hat{{\mathcal W}'} \|_\mathrm{F}.
\end{equation}
Though it is straightforward how to compute these orthogonal compressions in the special case that the mode-$j$ fibers of ${\mathcal W}'$ and $\hat{{\mathcal W}'}$ span the same space for each $j=1,2,3$, one should not expect an arbitrary pair of low multilinear rank tensors to have this special property.

In Section \ref{sec:TensorProcrustes} we explain how to compute orthogonal compressions of ${\mathcal W}'$ and $\hat{{\mathcal W}'}$ which satisfy equation \eqref{eq:ProcrustesEq}. Additionally, we explain our perspective on orthogonal compressions which differs slightly from the perspective most commonly taken in applications. While the tensor Procrustes problem plays an important role in our theoretical results, the solution may be unsurprising to readers familiar with the matrix Procrustes problem. Such readers can advance to Section \ref{sec:existenceCorProof}.

\subsection{ Orthogonal compressions of pairs of tensors: The tensor Procrustes problem}
\label{sec:TensorProcrustes}

As previously discussed, the situation most common to application is that one has access to a measurement ${{\mathcal M}}'={\mathcal T}'+ {\mathcal N} '$ of a rank $R$ tensor ${\mathcal T}'$ which is corrupted by measurement error $ {\mathcal N} '$. In order to approximate the factors underlying the tensor ${\mathcal T}'$, the practitioner will compute a rank $R$ approximation of ${{\mathcal M}}'$. Typically the rank of ${\mathcal T}'$ is significantly less than its dimensions, so to reduce the computational complexity of computing a rank $R$ approximation of ${{\mathcal M}}'$ one first computes a low  multilinear rank approximation ${\mathcal W}'$ of ${{\mathcal M}}'$. One then computes a $R$ rank approximation of the low multilinear rank approximation ${\mathcal W}'$.

While this approach is natural to use in applied settings, it can have adverse effects when searching for theoretical guarantees. Letting ${\hat{\cT}}'$ denote a best border rank $R$ approximation of ${{\mathcal M}}'$, to simplify to the $R \times R \times K$ setting, we need orthogonal compressions of both ${\mathcal T}'$ and ${\hat{\cT}}'$. If one independently computes orthogonal compressions ${\mathcal T}$ and ${\hat{\cT}}$ of ${\mathcal T}'$ and ${\hat{\cT}}'$, then it is possible that the distance between ${\mathcal T}$ and ${\hat{\cT}}$ is greater than the distance between ${\mathcal T}'$ and ${\hat{\cT}}'$. This difficulty must be treated with particular care if the span of the mode-$j$ subspaces of ${\mathcal T}'$ is not the same as the span of the mode-$j$ subspace of ${\hat{\cT}}$.

Roughly speaking this difficulty is resolved by first applying an orthogonal transformation to ${\hat{\cT}}$ which both does not increase the distance between ${\mathcal T}$ and ${\hat{\cT}}$ and makes it so that the mode-$j$ subspaces of the orthogonal transformation of ${\hat{\cT}}$ are the same as those of ${\mathcal T}$. From this point it is straightforward how to compute the appropriate orthogonal compressions. As the (border) rank of ${\hat{\cT}}'$ is invariant under orthogonal transformations, for our purposes there is no loss in studying an orthogonal transformation of ${\hat{\cT}}'$ rather than ${\hat{\cT}}'$ itself.\footnote{In practice one rarely has access to both ${\mathcal T}'$ and ${\hat{\cT}}'$, but this is immaterial to our present discussion.}

\begin{theorem}
\label{theorem:TensorProcrustes} Let ${\mathcal W}', \hat{{\mathcal W}'} \in  {\mathbb{K}} ^{I_1 \times \cdots \times I_\ell}$ be tensors having multilinear rank $(R_1, \cdots, R_\ell)$ and $( {\hat{R}}  _1, \cdots,  {\hat{R}}  _\ell)$, respectively. Let $ {\mathbf{V}} :  {\mathbb{K}} ^{R_j} \to  {\mathbb{K}} ^{I_j}$ be column-wise orthonormal with range equal to the subspace spanned by the mode-$j$ fibers of ${\mathcal W}'$. If $ {\hat{R}}  _j \leq R_j$, then there exists a column-wise orthonormal matrix $ {\hat{\mathbf{V}}} :  {\mathbb{R}} ^{R_j} \to  {\mathbb{R}} ^{I_j}$ with range containing the subspace spanned by the mode-$j$ fibers of ${\hat{\cW}}$ such that
\[
\|{\mathcal W}' \cdot_j  {\mathbf{V}} ^\mathrm{H} - \hat{{\mathcal W}'} \cdot_j  {\hat{\mathbf{V}}} ^\mathrm{H} \|_\mathrm{F} \leq \|{\mathcal W}'-\hat{{\mathcal W}'}\|_\mathrm{F}.
\]
As an immediate consequence, letting $R^*_i = \max \{R_i, {\hat{R}}  _i\}$ for each $i=1,\dots,\ell$ there exist orthogonal compressions ${\mathcal W},{\hat{\cW}} \in  {\mathbb{K}} ^{R^*_1 \times \cdots \times R^*_\ell} $ of ${\mathcal W}'$ and $\hat{{\mathcal W}'}$, respectively, such that
\[
\|{\mathcal W}-{\hat{\cW}}\|_\mathrm{F} \leq \|{\mathcal W}' -\hat{{\mathcal W}'}\|_\mathrm{F}.
\]
\end{theorem}

As will become clear, we may  use Theorem \ref{theorem:TensorProcrustes} to without loss of generality restrict to the case of $R \times R \times K$ tensors. Before proving the theorem we give a technical lemma.

\begin{lemma}
\label{lemma:ProcrustesIsometries}
Let $ {\mathbf{M}} , {\mathbf{N}}  \in  {\mathbb{K}} ^{I_1 \times I_2}$ be matrices of rank $R_ {\mathbf{M}} $ and $R_ {\mathbf{N}} $, respectively, where $R_ {\mathbf{M}}  \geq R_ {\mathbf{N}} $. There exists a unitary matrix $  {\mathbf{U}}  \in  {\mathbb{K}} ^{I_1 \times I_1}$ such that $\mathrm{ran \ }  {\mathbf{M}}  \supseteq \mathrm{ran \ }  {\mathbf{U}}   {\mathbf{N}} $ and so that
\[
\| {\mathbf{M}} - {\mathbf{U}}   {\mathbf{N}} \|_\mathrm{F} \leq \| {\mathbf{M}}  - {\mathbf{N}}  \|_\mathrm{F}.
\]
As a consequence, given any column-wise orthonormal matrix $ {\mathbf{U}} _ {\mathbf{M}}  \in  {\mathbb{K}} ^{I_1 \times R_ {\mathbf{M}} }$ mapping onto $\mathrm{ran \ }  {\mathbf{M}} $, there is a column-wise orthonormal matrix $ {\mathbf{U}} _ {\mathbf{N}}   \in  {\mathbb{K}} ^{I_1 \times R_ {\mathbf{N}} }$ such that $\mathrm{ran \ }  {\mathbf{U}} _ {\mathbf{N}}  \supseteq \mathrm{ran \ }  {\mathbf{N}} $ and such that
\[
\|  {\mathbf{U}} _ {\mathbf{M}} ^\mathrm{H}  {\mathbf{M}}  -  {\mathbf{U}} _ {\mathbf{N}} ^\mathrm{H}  {\mathbf{N}}  \|_\mathrm{F} \leq \| {\mathbf{M}} - {\mathbf{N}}  \|_\mathrm{F}.
\]
\end{lemma}
\begin{proof}
The proof of Lemma \ref{lemma:ProcrustesIsometries} is a routine argument using the solution of the orthogonal Procrustes problem \cite{S66}. Details are given in the supplementary materials.
\end{proof}

We now give the proof of Theorem \ref{theorem:TensorProcrustes}.

\begin{proof}
The main strategy of the proof is to apply Lemma \ref{lemma:ProcrustesIsometries} to each mode of the tensors ${\mathcal W}'$ and $\hat{{\mathcal W}'}$. To this end, fix $j$ and let ${\mathcal W}'_{(j)}$ and $\hat{{\mathcal W}'}_{(j)}$ be mode-$j$ unfoldings of ${\mathcal W}'$ and $\hat{{\mathcal W}'}$, respectively, with the same ordering of the fibers.  Recall that ${\mathcal W}_{(j)}$  has rank $R_j$ and $\hat{{\mathcal W}'}_{(j)}$ has rank $ {\hat{R}}  _j$, and WLOG assume $R_j \geq  {\hat{R}}  _j$.

Using Lemma \ref{lemma:ProcrustesIsometries}, given any column-wise orthonormal matrix $ {\mathbf{V}}  \in  {\mathbb{R}} ^{I_j \times R_j}$ with range equal to $\mathrm{ran \ } {\mathcal W}'_{(j)}$, there exists a column-wise orthonormal matrix $ {\hat{\mathbf{V}}}  \in  {\mathbb{K}} ^{I_j \times R_j}$ such that the range of $ {\hat{\mathbf{V}}} $ contains $\mathrm{ran \ } \hat{{\mathcal W}'}_{(j)}$ and such that
\[
\|  {\mathbf{V}} ^\mathrm{H} {\mathcal W}'_{(j)} -  {\hat{\mathbf{V}}} ^\mathrm{H} \hat{{\mathcal W}'}_{(j)} \|_\mathrm{F} \leq \|{\mathcal W}'_{(j)}-\hat{{\mathcal W}'}_{(j)} \|_\mathrm{F}.
\]
Since the ordering of the mode-$j$ fibers is the same in ${\mathcal W}_{(j)}$ and $\hat{{\mathcal W}'}_{(j)}$ we have
\[
\|  {\mathbf{V}} ^\mathrm{H} {\mathcal W}'_{(j)} -  {\hat{\mathbf{V}}} ^\mathrm{H} \hat{{\mathcal W}'}_{(j)}\|_\mathrm{F} = \|{\mathcal W}' \cdot_j  {\mathbf{V}} ^\mathrm{H} - \hat{{\mathcal W}'} \cdot_j  {\hat{\mathbf{V}}} ^\mathrm{H}\|_\mathrm{F}  \quad \quad \mathrm{and} \quad \quad \|{\mathcal W}'_{(j)}-\hat{{\mathcal W}'}_{(j)}\|_\mathrm{F}=\|{\mathcal W}'-\hat{{\mathcal W}'} \|_\mathrm{F}
\]
from which the first conclusion follows.

The second part of the result follows from a standard argument using induction on the mode of the fibers.
\end{proof}

It is worth noting that the Procrustes compressions ${\mathcal W}$ and ${\hat{\cW}}$ one arrives at using the methodology in the proof are not necessarily optimal. That is, there may exist cores ${\mathcal W}^*$ and $\hat{{\mathcal W}^*}$ such that
\[
\|{\mathcal W}^*-\hat{{\mathcal W}^*}\|_\mathrm{F} < \|{\mathcal W} - {\hat{\cW}} \|_\mathrm{F} \leq \|{\mathcal W}'- \hat{{\mathcal W}'}\|_\mathrm{F}.
\]
One may compute locally optimal cores by using an ALS approach for instance where one continues to iterate over the modes of the tensors and solves the matrix Procrustes problem for each corresponding unfolding. Methods for computing globally optimal cores are outside the scope of this article.

A similar issue occurs when computing a best low multilinear rank approximation of a given tensor ${{\mathcal M}}' \in  {\mathbb{K}} ^{I_1 \times I_2 \times I_3}$. Here one wants to minimize $\|{{\mathcal M}}'-{\mathcal W}'\|$ subject to a multilinear rank constraint on ${\mathcal W}'$. A common approach is to compute a MLSVD of ${{\mathcal M}}'$ and truncate at the specified multilinear rank. While this approach is in a sense optimal ``per mode", since the MLSVD treats each mode independently the result is rarely globally optimal when all modes are considered together.

\subsection{Proof of Proposition \ref{proposition:existenceCor}}
\label{sec:existenceCorProof}

We are now in position to prove Proposition \ref{proposition:existenceCor}.

\begin{proof}
The assumption that every border $ {\mathbb{K}} $-rank $R$ tensor in an open $\epsilon$ ball around ${\mathcal T}$ is simple in fact implies that this $\epsilon$ ball cannot contain a tensor with border $ {\mathbb{K}} $-rank strictly less than $R$. We temporarily assume this implication holds and give a detailed proof at the end.

Let ${{\mathcal M}}' \in  {\mathbb{K}} ^{I_1 \times I_2 \times I_3}$ be any tensor such that $\|{{\mathcal M}}'-{\mathcal T}'\| < \epsilon/2$, and let ${\hat{\cT}}'$ be a best border $ {\mathbb{K}} $-rank $R$ approximation of ${{\mathcal M}}'$. Note that ${\mathcal T}'$ has multilinear rank $(R,R,K)$ by assumption, so inequality \eqref{eq:BRankVsMLRank} shows that ${\mathcal T}'$ has border $ {\mathbb{K}} $-rank $R$ in addition to having $ {\mathbb{K}} $-rank $R$. It follows that
\[
\|{\hat{\cT}}' - {{\mathcal M}}' \| \leq \| {\mathcal T}' - {{\mathcal M}}' \|,  \quad \quad \mathrm{hence} \quad \quad \| {\mathcal T}' - {\hat{\cT}}' \| < \epsilon,
\]
as otherwise ${\mathcal T}'$ would be a strictly better border $ {\mathbb{K}} $-rank $R$ approximation of ${{\mathcal M}}'$ than ${\hat{\cT}}'$.  To show that ${{\mathcal M}}'$ has a best $ {\mathbb{K}} $-rank $R$ approximation, it is sufficient to show that ${\hat{\cT}}'$ has $ {\mathbb{K}} $-rank equal to $R$.

Since ${\hat{\cT}}'$ has border $ {\mathbb{K}} $-rank $R$, the rank of each mode-$i$ unfolding of ${\hat{\cT}}'$ is less than or equal to $R$. Therefore,
using Theorem \ref{theorem:TensorProcrustes}, we may simultaneously orthogonally compress ${\mathcal T}'$ and ${\hat{\cT}}'$ to tensors of size $R \times R \times K^*$ for some $K^* \leq R$ without increasing the Frobenius distance between the tensors. Therefore we may WLOG we may assume ${\mathcal T}'={\mathcal T}$ and ${\hat{\cT}}'={\hat{\cT}}$ both have size $R \times R \times K^*$. After this restriction, our assumptions imply that ${\hat{\cT}}$ is simple. An application of Theorem \ref{theorem:rankvsmult} \eqref{it:NonderogatoryImpliesRankR} then shows that ${\hat{\cT}}$ has rank $R$.

We next prove that every rank $R$ tensor in the open ball of Frobenius radius $\epsilon$ centered at ${\mathcal T}$ has a unique CPD. By assumption every border rank $R$ tensor in this ball is simple. As previously mentioned, this assumption implies that the $\epsilon$ ball does not contain a tensor with border rank less than $R$, hence every rank $R$ tensor in the ball is simple. Using Theorem \ref{theorem:rankvsmult} \eqref{it:RankIFFBasis} together with Lemma \ref{lemma:DistinctENsImpliesEigenbasis}, it follows that every tensor inside this ball satisfies Kruskal's condition \cite{K77}, hence every tensor in this ball has a unique CPD.

We now show that the set of tensors contained in this open ball which do not have a unique best $ {\mathbb{K}} $-rank $R$ approximation has Lebesgue measure $0$. For tensors with entries in $ {\mathbb{C}} $ (equipped with Lebesgue measure over the complexes) this is an immediate consequence of \cite{QML19}. For tensors with entries in $ {\mathbb{R}} $ (equipped with Lebesgue measure over the reals) this follows from  \cite[Corollary 18]{FO14} together with the first part of this result.

It remains to show that every tensor in an $\epsilon$ ball centered at ${\mathcal T}'$ has border $ {\mathbb{K}} $-rank at least $R$. To this end, suppose towards a contradiction that there is some tensor ${\mathcal S}'$ such that $\|{\mathcal T}'-{\mathcal S}'\|<\epsilon$ and such that ${\mathcal S}'$ has border $ {\mathbb{K}} $-rank $N < R$. Also fix $\delta_1,\delta_2 > 0$ such that $\|{\mathcal T}'-{\mathcal S}'\|+\delta_1+\delta_2 < \epsilon$. Since ${\mathcal S}'$ has border $ {\mathbb{K}} $-rank $N$, there must exist some tensor ${\mathcal W}'=[\![ {\mathbf{A}} ', {\mathbf{B}} ', {\mathbf{C}} ']\!]$ such that $\|{\mathcal W}'-{\mathcal S}'\|<\delta_1$ and such that $ {\mathbb{K}} $-$\mbox{rank}({\mathcal W}')= {\mathbb{K}} $-$\underline{\mbox{rank}}({\mathcal W}') = N$. Furthermore, a routine argument shows that ${\mathcal W}$ can be chosen so that its factors $ {\mathbf{A}} '$ and $ {\mathbf{B}} '$ have full column rank.

One can then extend the factors $ {\mathbf{A}} ', {\mathbf{B}} ', {\mathbf{C}} '$ to matrices
\[
\tilde{ {\mathbf{A}} }'=[ {\mathbf{A}} ' \ \  {\mathbf{E}} _{ {\mathbf{A}} }'] \qquad \tilde{ {\mathbf{B}} }'=[ {\mathbf{B}} ' \ \  {\mathbf{E}} _{ {\mathbf{B}} }'] \qquad \tilde{ {\mathbf{C}} }'=[ {\mathbf{C}} ' \ \   {\mathbf{E}} _{ {\mathbf{C}} }']
\]
each having $R$ columns where $ {\mathbf{E}} _{ {\mathbf{A}} }', {\mathbf{E}} _{ {\mathbf{B}} }'$ and $ {\mathbf{E}} _{ {\mathbf{C}} }'$ are chosen so that $\tilde{ {\mathbf{A}} }'$ and $ {\tilde{\bB}} '$ have full column rank, $\tilde{ {\mathbf{C}} }'$ has a repeated column, and $\|{\mathcal W}' - [\![ {\tilde{\bA}} ', {\tilde{\bB}} ', {\tilde{\bC}} ']\!]\| < \delta_2$. Using Theorem \ref{theorem:rankvsmult} shows that the tensor $[\![ {\tilde{\bA}} ', {\tilde{\bB}} ', {\tilde{\bC}} ']\!]$ has both $ {\mathbb{K}} $-rank and $ {\mathbb{K}} $-border rank equal to $R$. Furthermore, by construction we have
\[
\|{\mathcal T}'-[\![ {\tilde{\bA}} ', {\tilde{\bB}} ', {\tilde{\bC}} ']\!]\| < \|{\mathcal T}'-{\mathcal S}'\|+\delta_1+\delta_2 < \epsilon.
\]
However, Theorem \ref{theorem:rankvsmult} also shows that this tensor is not simple, since $ {\tilde{\bC}} '$ has repeated columns. This contradicts our assumption that every border rank $R$ tensor in an $\epsilon$ ball around ${\mathcal T}'$ is simple.
\end{proof}

\subsection{Metrics for joint generalized eigenvalues}

We now discuss in more detail the measures of distance we use for joint generalized eigenvalues. Recall the chordal metric between  two one-dimensional subspaces ${{\mathscr L}}_1,{{\mathscr L}}_2 \subseteq  {\mathbb{K}} ^K$ is the sine of the angle between ${{\mathscr L}}_1$ and ${{\mathscr L}}_2$. Letting $ {\vec{\lambda}} _i$ be a unit vector such that $\mathrm{span}( {\vec{\lambda}} _i)={{\mathscr L}}_i$ for $i=1,2$, it is not difficult to show that the chordal metric may be computed as
\begin{equation}
\label{eq:ChordalEqs}
\chi({{\mathscr L}}_1,{{\mathscr L}}_2)=\| {\mathbf{P}} _{({{\mathscr L}}_1)^\perp}  {\vec{\lambda}} _2 \|_2 = \max_{\substack{ {\boldsymbol{\gamma}}  \in ({{\mathscr L}}_1)^\perp \\ \| {\boldsymbol{\gamma}}  \|_2=1}} |\langle  {\vec{\lambda}} _2, {\boldsymbol{\gamma}}  \rangle | { = \sqrt{1- |\langle  {\vec{\lambda}} _1, {\vec{\lambda}} _2 \rangle|^2.}}
\end{equation}
Here $({{\mathscr L}}_1)^\perp \subset  {\mathbb{K}} ^K$ is the orthogonal complement of ${{\mathscr L}}_1 \subset  {\mathbb{K}} ^K$ and $ {\mathbf{P}} _{({{\mathscr L}}_1)^\perp}$ denotes the projection onto $({{\mathscr L}}_1)^\perp$. Note that $ |\langle  {\vec{\lambda}} _1, {\vec{\lambda}} _2 \rangle|$ is the cosine of the angle between the subspaces ${{\mathscr L}}_1$ and ${{\mathscr L}}_2$.

Also recall the matching distance defined in equation \eqref{eq:MatchingDistanceDef} which gives a metric between the spectra of two tensors. In addition to the matching distance we will need a weaker measure of separation between the spectra of two tensors. To this end, for $i=1,2$, let ${\mathcal T}_i \in  {\mathbb{K}} ^{R \times R \times K}$ be a tensor having $R$ JGE values counting algebraic multiplicity, and let $\{{{\mathscr L}}_{i,k}\}_{k=1}^R$ be the spectrum of ${\mathcal T}_i$. The \df{spectral variation} from ${\mathcal T}_1$ to ${\mathcal T}_2$, denoted ${\mathrm{sv}}[{\mathcal T}_1,{\mathcal T}_2]$, is defined by
\[
{\mathrm{sv}}[{\mathcal T}_1,{\mathcal T}_2] = \max_i \min_j \chi({{\mathscr L}}_{1,j}, {{\mathscr L}}_{2,i}).
\]

The key difference between spectral variation and matching distance is that, in matching distance, each JGE value of ${\mathcal T}_2$ must be paired with a unique JGE value of ${\mathcal T}_1$, while in spectral variation the JGE values of ${\mathcal T}_2$ need not be paired with unique JGE values of ${\mathcal T}_1$. Due to the lack of unique pairings in spectral variation, the spectral variation from ${\mathcal T}_1$ to ${\mathcal T}_2$ need not equal the spectral variation from ${\mathcal T}_2$ to ${\mathcal T}_1$. Furthermore, the spectral variation from ${\mathcal T}_1$ to ${\mathcal T}_2$ may be equal to zero while the spectra of ${\mathcal T}_1$ and ${\mathcal T}_2$ are not equal. As a consequence, while matching distance defines a metric on sets of JGE values, spectral variation does not.

\begin{exa}
\label{exa:SVdefective}
Let ${\mathcal T}_1,{\mathcal T}_2 \in  {\mathbb{K}} ^{R \times R \times 2}$. Suppose ${\mathcal T}_1$ has JGE values
\[
{{\mathscr L}}_{1,1}=\mathrm{span} ( {\mathbf{e}} _1) {\quad \quad \mathrm{and} \quad \quad} {{\mathscr L}}_{1,2}= \mathrm{span} ( {\mathbf{e}} _{ 2})
\]
and that ${\mathcal T}_2$ has JGE values
\[
{{\mathscr L}}_{2,1} = {{\mathscr L}}_{2,2} =\mathrm{span} ( {\mathbf{e}} _1).
\]
Here $ {\mathbf{e}} _1$ and $ {\mathbf{e}} _2$ are the standard basis vectors in $ {\mathbb{K}} ^2$. Then one has
\[
0 = {\mathrm{sv}}[{\mathcal T}_1,{\mathcal T}_2] \neq {\mathrm{sv}}[{\mathcal T}_2, {\mathcal T}_1] = 1.
\]
\end{exa}

Of course, matching distance bounds are more desirable to obtain than spectral variation bounds; however, they are also typically much more difficult to obtain. A standard strategy to obtain a matching distance bound is to first find a spectral variation bound which grows linearly with the magnitude of a perturbation.

\subsection{Bauer--Fike Theorem for JGE values}
We now present our main perturbation theoretic bound for the spectral variation between two tensors of $ {\mathbb{K}} $-rank $R$. The bound we present is in essence determined by the spectral norm of (a transformation of) the error tensor.

\begin{theorem}
\label{theorem:TensBauerFike}
Let ${\mathcal T}$ and ${\mathcal W}$ be $R \times R \times K$ tensors with entries in $ {\mathbb{K}} $. Assume that ${\mathcal T}$ has $ {\mathbb{K}} $-rank $R$ with CPD ${\cp{ {\mathbf{A}} , {\mathbf{B}} , {\mathbf{C}} }}$ where the columns of $ {\mathbf{C}} $ are normalized to have unit norm. Let $\{{{\mathscr L}}_k\}_{k=1}^K$ be the spectrum of ${\mathcal T}$ and set $ {\mathcal E}  = {\mathcal T} -{\mathcal W}$. If ${{\mathscr W}}$ is a JGE value of ${\mathcal W}$ with geometric multiplicity at least one, then
\begin{equation}
\label{eq:JGEvaluebound}
\min_k \chi ({{\mathscr L}}_k, {{\mathscr W}}) { \leq \sqrt{R} \Big\| {\mathcal E}  \cdot_1  {\mathbf{A}} ^{-1} \cdot_2  {\mathbf{B}} ^{-1} \cdot_3  {\mathbf{P}} _{{{\mathscr W}}^\perp} \Big\|_{\mathrm{sp}} } \leq   \frac{\sqrt{R} \spnorm{{  {\mathcal E} }}}{\sigma_{\min} ( {\mathbf{A}} ) \sigma_{\min} ( {\mathbf{B}} )}.
\end{equation}
Here $\sigma_{\min} ( {\mathbf{A}} )$ and $\sigma_{\min} ( {\mathbf{B}} )$ denote the smallest\footnote{Lemma \ref{lemma:invertslice} shows that $ {\mathbf{A}} $ and $ {\mathbf{B}} $ are invertible, hence these singular values are nonzero.} singular values of $ {\mathbf{A}} $ and $ {\mathbf{B}} $, respectively and $ {\mathbf{P}} _{{{\mathscr W}}^\perp}$ denotes the projection onto the orthogonal complement of ${{\mathscr W}}$. As a consequence, if ${\mathcal W}$ is a slice mix invertible border $ {\mathbb{K}} $-rank $R$ tensor, then
\begin{equation}
\label{eq:JGEBauerFike}
{\mathrm{sv}} [{\mathcal T}, {\mathcal W}] {  \leq  \sqrt{R} \Big\|  {\mathcal E}  \cdot_1  {\mathbf{A}} ^{-1} \cdot_2  {\mathbf{B}} ^{-1}  \Big\|_{\mathrm{sp}} } \leq \frac{\sqrt{R} \spnorm{{  {\mathcal E} }}}{\sigma_{\min} ( {\mathbf{A}} ) \sigma_{\min} ( {\mathbf{B}} )}.
\end{equation}
In the following special cases the coefficient $\sqrt{R}$ may be dropped from equations \eqref{eq:JGEvaluebound} and \eqref{eq:JGEBauerFike}$\mathrm{:}$
\begin{enumerate}
\item $K=2$.
\item ${\mathcal W}$ is also slice mix invertible and has $ {\mathbb{K}} $-rank $R$  with CPD $[\![ {\mathbf{A}} _{\mathcal W}, {\mathbf{B}} _{\mathcal W}, {\mathbf{C}} _{\mathcal W}]\!]$ where $ {\mathbf{A}} _{{\mathcal W}}= {\mathbf{A}} $ or $ {\mathbf{B}} _{\mathcal W}= {\mathbf{B}} $.
\end{enumerate}
That is, in these special cases one has
\begin{equation}
\label{eq:JGEBauerFikeSpecial}
{\mathrm{sv}} [{\mathcal T}, {\mathcal W}] {  \leq   \Big\|  {\mathcal E}  \cdot_1  {\mathbf{A}} ^{-1} \cdot_2  {\mathbf{B}} ^{-1}  \Big\|_{\mathrm{sp}} } \leq \frac{ \spnorm{{  {\mathcal E} }}}{\sigma_{\min} ( {\mathbf{A}} ) \sigma_{\min} ( {\mathbf{B}} )}.
\end{equation}
\end{theorem}

In the general case, the proof of Theorem \ref{theorem:TensBauerFike} follows a similar strategy to the proof of \cite[Chapter VI, Theorem 2.7]{SS90} with modifications to address issues arising when $K>2$. In particular, \cite[Chapter VI, Theorem 2.7]{SS90} shows
\[
\min_k \chi ({{\mathscr L}}_k, {{\mathscr W}})  \leq   \frac{\tnorm{{  {\mathcal E} _{(1)}}}}{\sigma_{\min} ( {\mathbf{A}} ) \sigma_{\min} ( {\mathbf{B}} )}
\]
in the $K=2$ case. As the proof in the general case is similar to Stewart and Sun's argument, we give it in the supplementary materials. However, if we additionally assume that ${\mathcal W}$ has $ {\mathbb{K}} $-rank $R$, then a new, simplified proof which is more intuitive can be given. This we now describe.

\begin{proof}

Assume that ${\mathcal W} \in  {\mathbb{K}} ^{R \times R \times R}$ has rank $R$ and is slice mix invertible and let $[\![  {\tilde{\bA}} ,  {\tilde{\bB}} ,  {\tilde{\bC}}  ]\!]$ be a CPD of ${\mathcal W} \cdot_1  {\mathbf{A}} ^{-1} \cdot_2  {\mathbf{B}} ^{-1}$ where the columns of $ {\tilde{\bC}} $ have unit norm.\footnote{Up to permutation and scaling one has $ {\tilde{\bA}}  = {\mathbf{A}} ^{-1}  {\mathbf{A}} _{{\mathcal W}}$ and $ {\tilde{\bB}}  =  {\mathbf{B}} ^{-1}  {\mathbf{B}} _{{\mathcal W}}$ and $ {\tilde{\bC}} = {\mathbf{C}} _{\mathcal W}$.} Additionally set $ {\mathcal E}  := {\mathcal W}-{\mathcal T}$.  As a consequence of Theorem \ref{theorem:rankvsmult} \eqref{it:RankIFFBasis}, the spans of the columns of $ {\mathbf{C}} $ and $ {\tilde{\bC}} $ give the JGE values of ${\mathcal T}$ and ${\mathcal W}$, respectively. For each column $ {\tilde{\bc}} _j$ of $ {\tilde{\bC}} $ we will show that the spectral norm of $ {\mathcal E}  \cdot_1  {\mathbf{A}} ^{-1} \cdot_2  {\mathbf{B}} ^{-1} \cdot_3  {\mathbf{P}} _ {\tbc_j^\perp} $ gives an upper bound for
\[
\min_{r} \chi ( {\tilde{\bc}} _j, {\vec{c}} _r) = \min_r \|  {\mathbf{P}} _{ {\tbc_j^\perp} }  {\vec{c}} _r \|_2
\]
from which the result follows. Here $ {\mathbf{P}} _ {\tbc_j^\perp} $ is the orthogonal projection on the orthogonal complement of $\mathrm{span}( {\tilde{\bc}} _j)$.

Expanding out the spectral norm of the transformed error tensor, we have
\begin{equation}
\label{eq:ErrorSpectralProjection}
\spnorm{ {\mathcal E}  \cdot_1  {\mathbf{A}} ^{-1} \cdot_2  {\mathbf{B}} ^{-1} \cdot_3  {\mathbf{P}} _ {\tbc_j^\perp} } = \max_{\| {\vec{u}} \|_2=\| {\vec{v}} \|_2=\| {\vec{z}} \|_2=1} \bigg | \bigg(\sum_{r \neq j}  {\tilde{\ba}} _r \otimes  {\tilde{\bb}} _r \otimes \left(  {\mathbf{P}} _ {\tbc_j^\perp}   {\tilde{\bc}} _r \right) - \sum_{r=1}^R  {\mathbf{e}} _r \otimes  {\mathbf{e}} _r \otimes \left( {\mathbf{P}} _ {\tbc_j^\perp}   {\vec{c}} _r\right) \bigg) \cdot_1  {\vec{u}}  \cdot_2  {\vec{v}}  \cdot_3  {\vec{z}}  \bigg |.
\end{equation}

Note that  since $ {\mathbf{P}} _ {\tbc_j^\perp}   {\tilde{\bc}} _j =  {\vec{0}} $, the tensor ${\tilde{\cT}} \cdot_1  {\mathbf{A}} ^{-1} \cdot_2  {\mathbf{B}} ^{-1} \cdot_3  {\mathbf{P}} _ {\tbc_j^\perp} $ which appears on the left hand side of the difference in equation \eqref{eq:ErrorSpectralProjection} has rank at most $R-1$. It follows that there is some vector $ {\vec{y}}  \in  {\mathbb{K}} ^R$ which is orthogonal to the subspace spanned by the mode$-2$ fibers of this tensor. In particular, let $ {\vec{y}}  \in  {\mathbb{K}} ^R$ be a unit vector such that $\langle  {\vec{y}} , {\tilde{\bb}} _r \rangle =0$ for all $r \neq j$. With this choice of $ {\vec{y}} $ we have
\[
\sum_{r \neq j}\langle  {\vec{y}} , {\tilde{\bb}} _r \rangle  {\tilde{\ba}} _r  \otimes \left(  {\mathbf{P}} _ {\tbc_j^\perp}   {\tilde{\bc}} _r \right) = 0
\]
from which it follows that
\begin{equation}
\begin{array}{rclcl}
\spnorm{ {\mathcal E}  \cdot_1  {\mathbf{A}} ^{-1} \cdot_2  {\mathbf{B}} ^{-1}  {\vec{y}} ^{\mathrm{H}} \cdot_3  {\mathbf{P}} _ {\tbc_j^\perp} } & = & \max_{\| {\vec{u}} \|_2=\| {\vec{z}} \|_2=1} \left |\left(\sum_{r=1}^R \langle  {\vec{y}} , {\mathbf{e}} _r \rangle  {\mathbf{e}} _r  \otimes \left( {\mathbf{P}} _ {\tbc_j^\perp}   {\vec{c}} _r\right) \right) \cdot_1  {\vec{u}}  \cdot_3  {\vec{z}}  \right| \\
 & = & \spnorm{ \sum_{r=1}^R \langle  {\vec{y}} , {\mathbf{e}} _r \rangle  {\mathbf{e}} _r  \otimes \big( {\mathbf{P}} _ {\tbc_j^\perp}   {\vec{c}} _r\big) } \\
 & \geq &  \max_{r} \| \langle  {\vec{y}} , {\mathbf{e}} _r \rangle  {\mathbf{P}} _ {\tbc_j^\perp}   {\vec{c}} _r \|_2.
\end{array}
\end{equation}

Recalling that $ {\vec{y}} $ is a unit vector, there must exist at least one index $r_*$ such that $ | \langle  {\vec{y}} , {\mathbf{e}} _r \rangle | \geq 1/\sqrt{R}$, from which we obtain
\[
\spnorm{ {\mathcal E}  \cdot_1  {\mathbf{A}} ^{-1} \cdot_2  {\mathbf{B}} ^{-1}  {\vec{y}} ^{\mathrm{H}} \cdot_3  {\mathbf{P}} _ {\tbc_j^\perp} } \geq\| \langle  {\vec{y}} , {\mathbf{e}} _{r_*} \rangle  {\mathbf{P}} _ {\tbc_j^\perp}   {\vec{c}} _{r_*} \|_2 \geq  \frac{\chi( {\tilde{\bc}} _j, {\vec{c}} _{r_*})}{\sqrt{R}} \geq \frac{\min_r \chi ( {\tilde{\bc}} _j, {\vec{c}} _{r})}{\sqrt{R}}.
\]
We conclude
\[
\sqrt{R} \tnorm{ {\mathbf{A}} ^{-1}} \tnorm{ {\mathbf{B}} ^{-1}} \spnorm{ {\mathcal E} }\geq \sqrt{R} \spnorm{ {\mathcal E}  \cdot_1  {\mathbf{A}} ^{-1} \cdot_2  {\mathbf{B}} ^{-1}  {\vec{y}} ^{\mathrm{H}} \cdot_3  {\mathbf{P}} _ {\tbc_j^\perp} } \geq \min_r \chi ( {\tilde{\bc}} _j, {\vec{c}} _{r}).
\]
The result in this special case follows from combining the above inequality with Theorem \ref{theorem:rankvsmult} \eqref{it:RankIFFBasis}. If additionally one has $ {\mathbf{B}} _{\mathcal W} =  {\mathbf{B}} $, then one can take $ {\tilde{\bB}}  =  {\mathbf{I}} $. The fact that the coefficient $\sqrt{R}$ may be dropped in this case easily follows. The case $ {\mathbf{A}} _{\mathcal W} = {\mathbf{A}} $ is similar.
\end{proof}

\subsubsection{Scaling of $ {\mathbf{A}} $ and $ {\mathbf{B}} $}

The quality of the bound obtained in Theorem \ref{theorem:TensBauerFike} depends on the scaling of the factor matrices $ {\mathbf{A}} $ and $ {\mathbf{B}} $ in the CPD ${\mathcal T}=[\![ {\mathbf{A}} , {\mathbf{B}} , {\mathbf{C}} ]\!]$. More precisely, for any diagonal matrix $ {\mathbf{D}} $, one has that ${\mathcal T}=[\![ {\mathbf{A}}  {\mathbf{D}} , {\mathbf{B}}  {\mathbf{D}} ^{-1}, {\mathbf{C}} ]\!]$ is a CPD of ${\mathcal T}$, and the quantity of $\sigma_{\min}( {\mathbf{A}}   {\mathbf{D}} ) \sigma_{\min} ( {\mathbf{B}}   {\mathbf{D}} ^{-1})$ is dependent on the choice of $ {\mathbf{D}} $. To obtain an optimal bound, one wants to minimize the right hand side of equation \eqref{eq:JGEBauerFike}, hence one wants to maximize $\sigma_{\min}( {\mathbf{A}}   {\mathbf{D}} ) \sigma_{\min} ( {\mathbf{B}}   {\mathbf{D}} ^{-1})$ over invertible diagonal matrices $ {\mathbf{D}} $.

One strategy to compute the optimal $ {\mathbf{D}} $ is to use an ALS approach in which one sets all but one  of the diagonal entries of $ {\mathbf{D}} $ to be equal to $1$, then maximizes over the remaining diagonal entry of $ {\mathbf{D}} $. By repeating this approach while varying which entry of $ {\mathbf{D}} $ is not fixed, one is able to obtain a (local) optimum for equation \eqref{eq:JGEBauerFike}.

As an alternative strategy, one may simply normalize $ {\mathbf{A}} $ and $ {\mathbf{B}} $ so that the $r$th column of $ {\mathbf{A}} $ and $ {\mathbf{B}} $ have the same norm for each $r=1,\dots,R$. Although this method cannot be expected to be optimal, it appears to perform well for randomly generated $ {\mathbf{A}} $ and $ {\mathbf{B}} $, and the associated computational cost is minimal.

\subsubsection{On the coefficient $\sqrt{R}$ }

Knowing that there are various special cases in which the coefficient $\sqrt{R}$ may be dropped from equation \eqref{eq:JGEBauerFike}, it is natural to wonder if there are examples where this coefficient is necessary, or if it is simply a relic of the proof. The following example shows that some constant is indeed necessary. That is, equation \eqref{eq:JGEBauerFike} does not hold for arbitrary pairs of slice mix invertible $ {\mathbb{K}} $-rank $R$ tensors if the coefficient $\sqrt{R}$ is removed.

\begin{exa}
Let ${\mathcal T}$ have CPD ${\cp{ {\mathbf{A}} , {\mathbf{B}} , {\mathbf{C}} }}$ where
\[
 {\mathbf{A}}  =  {\mathbf{B}}  =  {\mathbf{I}} _3 {\quad \quad \mathrm{and} \quad \quad}  {\mathbf{C}}  = \begin{pmatrix} 1 & 0 & \frac{\sqrt{2}}{2} \\ 0 & 1 & \frac{\sqrt{2}}{2} \\ 0 & 0 & 0 \end{pmatrix}
\]
and let ${\mathcal W}$ have CPD $[\![ {\tilde{\bA}} , {\tilde{\bB}} , {\tilde{\bC}} ]\!]$ where
\[
 {\tilde{\bA}}  = {\tilde{\bB}}   = \begin{pmatrix} \frac{1}{10\sqrt{3}} & \frac{1}{\sqrt{3}} & 0 \\ \frac{1}{10\sqrt{3}} & -\frac{1}{2\sqrt{3}} & -\frac{1}{\sqrt{2}} \\ \frac{1}{10\sqrt{3}} & -\frac{1}{2\sqrt{3}} & \frac{1}{\sqrt{2}}
\end{pmatrix} {\quad \quad \mathrm{and} \quad \quad}  {\tilde{\bC}}  = \begin{pmatrix} 0 & 1 & 0 \\ 0 & 0 & 1 \\ 1 & 0 & 0 \end{pmatrix}
\]
Then ${\mathcal T}$ and ${\mathcal W}$ are both slice mix invertible $ {\mathbb{R}} $-rank $3$ tensors and the spectral variation ${\mathrm{sv}}[{\mathcal T},{\mathcal W}]$ is equal to $1$, as seen from the fact that the first column of $ {\tilde{\bC}} $ is orthogonal to all three columns of $ {\mathbf{C}} $. However, setting $ {\mathcal E} ={\mathcal W}-{\mathcal T}$ one has
\[
\Big\|  {\mathcal E} \cdot_1  {\mathbf{A}} ^{-1} \cdot_2  {\mathbf{B}} ^{-1} \Big\|_{\mathrm{sp}} < 1
\]
illustrating that equation \eqref{eq:JGEBauerFike} does not hold for this pair\footnote{Computation of the spectral norm is NP-hard \cite{HL13}; however, the mode three unfolding $ {\mathcal E} \cdot_1  {\mathbf{A}} ^{-1} \cdot_2  {\mathbf{B}} ^{-1}$ has spectral norm $\approx 0.983$ and this norm is as an upper bound for the tensor spectral norm. Using Tensorlab \cite{VDSBL16} suggests $\Big\|  {\mathcal E} \cdot_1  {\mathbf{A}} ^{-1} \cdot_2  {\mathbf{B}} ^{-1} \Big\|_{\mathrm{sp}} \approx 0.882$.}.
\end{exa}

\subsection{Matching distance bounds and the $K=2$ case}

While the spectral variation bound from Theorem \ref{theorem:TensBauerFike} provides insight as to how a CPD is affected by a perturbation, as illustrated by Example \ref{exa:SVdefective}, a spectral variation bound cannot be used to determine that a perturbation of a given tensor has distinct JGE values. To accomplish this we need to use a matching distance bound such as Theorem \ref{theorem:TensMDBound}.

As briefly mentioned before, a spectral variation bound can often be converted to a matching distance bound using the following idea:  Supposing temporarily that ${\mathcal W}$ is slice mix invertible and has border $ {\mathbb{K}} $-rank $R$ and letting $\{{{\mathscr L}}_r\}_{r=1}^R$ be the  spectrum of ${\mathcal T}$, the bound given in Theorem \ref{theorem:TensBauerFike} gives a collection of neighborhoods $\{\mathfrak{N}_\epsilon ({{\mathscr L}}_r)\}_{r=1}^R$ in which the JGE values of ${\mathcal W}$ must lie. Additionally, it is not difficult to show that JGE values are continuous in a perturbation along the set of slice mix invertible border rank $R$ tensors.\footnote{For example, one may use the fact that the characteristic polynomial of a slice mix invertible border $ {\mathbb{K}} $-rank $R$ tensor factors as a product of $R$ linear terms to show this.}

More precisely, if $\mathfrak{T}:[0,1]\to  {\mathbb{K}} ^{R \times R \times K} $ is a continuous function such that $\mathfrak{T} (0) = {\mathcal T}$ and $\mathfrak{T}(t)$ has border $ {\mathbb{K}} $-rank $R$ and is slice mix invertible for all $t \in [0,1]$, then there exists a continuous function $\mathfrak{L}:[0,1] \to [\text{Gr}(1, {\mathbb{K}} ^K)]^R$ such that $\mathfrak{L}(t)$ is equal to the  spectrum of $\mathfrak{T}(t)$ for all $t \in [0,1]$. Here $\text{Gr}(1, {\mathbb{K}} ^K)$ denotes the set of one-dimensional subspaces of $ {\mathbb{K}} ^K$ and $[\text{Gr}(1, {\mathbb{K}} ^K)]^R$ denotes a collection of $R$ elements of $\text{Gr}(1, {\mathbb{K}} ^K)$ with repetition allowed.\footnote{The collection of one-dimensional subspaces of $ {\mathbb{K}} ^K$ is a metric space under the chordal metric. In general, the collection of $d$-dimensional subpsaces of $ {\mathbb{K}} ^K$ is called a Grassmanian.}

The above discussion implies that if ${\mathcal W}$ is slice mix invertible and has border $ {\mathbb{K}} $-rank $R$ and if there is a path of slice mix invertible border $ {\mathbb{K}} $-rank $R$ tensors from ${\mathcal T}$ to ${\mathcal W}$ lying in the ball of radius $\spnorm{{\mathcal T}-{\mathcal W}}$ (i.e. if $\mathfrak{T}$ can be chosen so $\mathfrak{T}(0)={\mathcal T}$ and $\mathfrak{T}(1) = {\mathcal W}$ and $\spnorm{{\mathcal T}-\mathfrak{T}(t)} \leq  \spnorm{{\mathcal T}-{\mathcal W}}$ for all $t \in [0,1]$) and if the $\mathfrak{N}_\epsilon ({{\mathscr L}}_i) \cap \mathfrak{N}_\epsilon ({{\mathscr L}}_j)=\emptyset$ for each $i \neq j$, then each $\mathfrak{N}_\epsilon ({{\mathscr L}}_r)$ must contain exactly one JGE value of ${\mathcal W}$. As a consequence, if these assumptions hold, then the bound given by \ref{theorem:TensBauerFike} is in fact a matching distance bound and may be used to guarantee the existence of distinct JGE values for ${\mathcal W}$. See Figure \ref{fig:EigPertFig} for an illustration of the functions $\mathfrak{T}$ and $\mathfrak{L}$ and of the neighborhoods $\NepsL{r}$.

\begin{figure}
\label{fig:EigPertFig}
\hspace*{-1cm}
	\scalebox{.8}{
\input{EigPertFigWithNhds.tex}}
\caption{Illustration of continuous functions $\mathfrak{T}$ and $\mathfrak{L}$ defined on $[0,1]$. The function $\mathfrak{T}$ maps to a path of border $ {\mathbb{K}} $-rank $R$ tensors from ${\mathcal T}$ to ${\hat{\cT}}$. The function $\mathfrak{L}$ is defined so that $\mathfrak{L}(t)$ is equal to the joint generalized spectrum of $\mathfrak{T}(t)$ for all $t \in [0,1]$.  All JGE values of $\hat{{\mathcal T}}$ are contained in the union of the regions $\NepsL{1}$ and $\NepsL{2}$. Furthermore the neighborhood $\NepsL{r}$ contains a JGE value of $\hat{{\mathcal T}}$ for each $r=1,2$.}
\end{figure}

The primary issue with using such an argument for general $K$ is that, it is not known if there is a continuous function $\mathfrak{T}$ which satisfies the necessary assumptions. That is, it is not known if there is a path of slice mix invertible border rank $R$ tensors from ${\mathcal T}$ to ${\mathcal W}$ which lies in the ball of radius $\spnorm{{\mathcal T}-{\mathcal W}}$ centered at ${\mathcal T}$. However, in the case that $K=2$, the condition needed to guarantee that the $\mathfrak{N}_\epsilon ({{\mathscr L}}_r)$ do not intersect in fact guarantees the existence of such a path.

As discussed in \cite[Chapter VI]{SS90}, if ${\mathcal T}$ is a slice mix invertible tensor of size $R \times R \times 2$, then the (joint) generalized eigenvalue problem for ${\mathcal T}$ is equivalent to a classical eigenvalue problem. E.g. if $ {\mathbf{T}} _1$ is invertible, then the JGE eigenvectors of ${\mathcal T}$ are exactly equal to the classical eigenvectors of $ {\mathbf{T}} _1^{-1}  {\mathbf{T}} _2$. Using this equivalence, we for free obtain basic facts about multiplicities of JGE values for slice mix invertible $R \times R \times 2$ tensors. Namely, a slice mix invertible $R \times R \times 2$ tensor always has $R$ JGE values in $ {\mathbb{C}} ^2$ counting algebraic multiplicity, and each JGE value of such a tensor must have geometric multiplicity at least one. We will make frequent use of these facts in the upcoming results.

\begin{proposition}
\label{proposition:RegularityBound}
Let ${\mathcal T} \in  {\mathbb{K}} ^{R \times R \times 2}$ be a tensor of $ {\mathbb{K}} $-rank $R$ with CPD ${\mathcal T}={\cp{ {\mathbf{A}} , {\mathbf{B}} , {\mathbf{C}} }}$, where $ {\mathbf{C}} $ has columns of unit norm and let $\{{{\mathscr L}}_r\}_{r=1}^R$ be the set of JGE values of ${\mathcal T}$. Let $ {\mathcal E}  \in  {\mathbb{K}} ^{R \times R \times 2}$ be any tensor which satisfies
\[
\spnorm{ {\mathcal E} } \leq \frac{\sigma_{\min} ( {\mathbf{A}} ) \sigma_{\min} ( {\mathbf{B}} ) \min_{i \neq j} \chi ({{\mathscr L}}_i, {{\mathscr L}}_j)}{2}.
\]
Then ${\mathcal T}+ {\mathcal E} $ is slice mix invertible and has border $ {\mathbb{C}} $-rank equal to $R$.
\end{proposition}
\begin{proof}
To simplify notation, set
\[
\delta = \min_{i \neq j} \frac{\chi ({{\mathscr L}}_i,{{\mathscr L}}_j)}{2}.
\]
By assumption the columns of $ {\mathbf{C}} $ all have unit norm, therefore, by reindexing columns if necessary, Theorem \ref{theorem:rankvsmult} shows that for each $r=1,\dots,R$ the column $ {\vec{c}} _r$ satisfies ${{\mathscr L}}_r=\mathrm{span}( {\vec{c}} _r)$. To show ${\mathcal T}+ {\mathcal E} $ is slice mix invertible, we first show that there is a unit vector $ {{\vec{q}}} \in  {\mathbb{K}} ^2$ such that $|\langle  {\vec{c}} _r,  {{\vec{q}}}  \rangle| \geq \delta$ for all $r$. Using this we will show that the matrix $\bar{q}_1 ( {\mathbf{T}} _1+ {\mathbf{E}} _1)+ \bar{q}_2 ( {\mathbf{T}} _2+ {\mathbf{E}} _2)$ is invertible, hence ${\mathcal T}+ {\mathcal E} $ is slice mix invertible.

To this end let ${{\mathscr Z}} \subset  {\mathbb{K}} ^2$ be a one-dimensional subspace of $ {\mathbb{K}} ^2$ such that $\chi({{\mathscr L}}_1,{{\mathscr Z}})=\delta$. Then using the triangle inequality we must have
\begin{equation}
\label{eq:PsiChordalGap}
\chi({{\mathscr L}}_r,{{\mathscr Z}}) \geq \delta \qquad \mathrm{for \ all \ } r=1,\dots, R,
\end{equation}
as otherwise would imply that there is an index $\ell \neq 1$ such that
\[
\chi({{\mathscr L}}_\ell,{{\mathscr Z}})<\delta \qquad \mathrm{hence} \qquad \chi({{\mathscr L}}_1,{{\mathscr L}}_\ell) \leq \chi({{\mathscr L}}_1,{{\mathscr Z}})+\chi({{\mathscr L}}_\ell,{{\mathscr Z}})< 2\delta= \min_{i \neq j} \chi ({{\mathscr L}}_i,{{\mathscr L}}_j),
\]
which is a contradiction.

Now set ${{\mathscr Q}}={{\mathscr Z}}^\perp \subset  {\mathbb{K}} ^2$ and let $\mathbf{q} \in  {\mathbb{K}} ^2$ be a unit vector such that ${{\mathscr Q}} = \mathrm{span}(\mathbf{q})$. Using equation \eqref{eq:PsiChordalGap}, for each $r$ we have
\[
\delta \leq \chi({{\mathscr L}}_\ell,{{\mathscr Z}})= \sin (\angle ({{\mathscr L}}_r,{{\mathscr Z}}))=\sin(\pi/2 -\angle ({{\mathscr L}}_r,{{\mathscr Q}}))= \cos (\angle ({{\mathscr L}}_r,{{\mathscr Q}}))=|\langle  {\vec{c}} _r, \mathbf{q} \rangle|
\]
where $\angle(\cdot,\cdot)$ denotes the angle between two one-dimensional subspaces of $ {\mathbb{K}} ^2$.

Letting $D_k ( {\mathbf{C}} )  \in  {\mathbb{K}} ^{R \times R}$ be the diagonal matrix whose diagonal entries are given by the $k$th row of $ {\mathbf{C}} $ as usual, we have that $\bar{q}_1 D_1 ( {\mathbf{C}} ) +\bar{q}_2 D_2 ( {\mathbf{C}} )$ is a diagonal matrix with $r$th diagonal entry equal to $\langle  {\vec{c}} _r, \mathbf{q} \rangle$. If follows that
\begin{equation}
\label{eq:OmegaSingVal}
\sigma_{\min} (\bar{q}_1 D_1 ( {\mathbf{C}} ) +\bar{q}_2 D_2 ( {\mathbf{C}} ))= \min (\{|\langle  {\vec{c}} _r, \mathbf{q} \rangle|\}_{r=1}^R) \geq \delta.
\end{equation}
Finally observe that
\[
\begin{array}{rclcl}
\sigma_{\min} \big(\bar{q}_1 ( {\mathbf{T}} _1+ {\mathbf{E}} _1)+ \bar{q}_2 ( {\mathbf{T}} _2+ {\mathbf{E}} _2)\big) &\geq& \sigma_{\min} (\bar{q}_1  {\mathbf{T}} _1 +\bar{q}_2  {\mathbf{T}} _2) - \|\bar{q}_1  {\mathbf{E}} _1 +\bar{q}_2  {\mathbf{E}} _2\|_2 \\
 &\geq & \sigma_{\min} \Big( {\mathbf{A}}  \Big(\bar{q}_1 D_1 \big( {\mathbf{C}} )+\bar{ q}_2 D_2 ( {\mathbf{C}} )\Big)  {\mathbf{B}} ^\mathrm{T}\Big)-\spnorm{ {\mathcal E} } \\
 & \geq & \sigma_{\min} ( {\mathbf{A}} ) \sigma_{\min} ( {\mathbf{B}} ) \sigma_{\min}\big(\bar{q}_1 D_1 \big( {\mathbf{C}} )+\bar{q}_2 D_2 ( {\mathbf{C}} )\big)- \spnorm{ {\mathcal E} }\\
\end{array}
\]

Using inequality \eqref{eq:OmegaSingVal} together with the assumption that $\spnorm{ {\mathcal E} }< \sigma_{\min} ( {\mathbf{A}} ) \sigma_{\min} ( {\mathbf{B}} )  \delta$ shows that ${\mathcal T}+ {\mathcal E} $ is slice mix invertible.

Having shown ${\mathcal T}+ {\mathcal E} $ is slice mix invertible, it follows that ${\mathcal T}+ {\mathcal E} $ has multilinear rank $(R,R,\cdot)$. From equation \eqref{eq:BRankVsMLRank} we conclude that the border rank of ${\mathcal T}+ {\mathcal E} $ is at least $R$. Combining this with the well-known fact that tensors of size $R \times R \times 2$ have border rank less than or equal to $R$, e.g. see \cite[Theorem 5.5.1.1]{Lan12}, shows that ${\mathcal T}+ {\mathcal E} $ has border rank $R$.
\end{proof}

\subsection{Proof of Theorem \ref{theorem:TensMDBound}}

We now give the proof of Theorem \ref{theorem:TensMDBound}.

\begin{proof}

Before proceeding we note that the proof of the $ {\mathbb{K}} = {\mathbb{R}} $ case makes use of the result over $ {\mathbb{C}} $, so we first consider the complex case. Now set
\begin{equation}
\epsilon=\frac{\spnorm{{\mathcal T}-{\mathcal W}}}{\sigma_{\min} ( {\mathbf{A}} ) \sigma_{\min} ( {\mathbf{B}} )} < \frac{ \min_{i \neq j} \chi ({{\mathscr L}}_i, {{\mathscr L}}_j)}{2}
\end{equation}
and for each $r=1,\dots,R$ let $\NepsL{r}$ be the collection of one-dimensional subspaces in $ {\mathbb{C}} ^2$ defined by
\[
\NepsL{r} = \{ {{\mathscr Z}} \subset  {\mathbb{C}} ^2 | \ \dim({{\mathscr Z}})=1 \ \mathrm{and} \ \chi({{\mathscr L}}_r,{{\mathscr Z}}) \leq \epsilon\}.
\]
That is, $\NepsL{r}$ is the ball of chordal radius $\epsilon$ centered at ${{\mathscr L}}_r$. Then our choice of $\epsilon$ together with the triangle inequality gives
\[
\NepsL{i} \cap \NepsL{j} = \emptyset \quad \quad \mathrm{for \ all \ } i \neq j.
\]
To prove inequality \eqref{eq:MatchingDistanceBound}, it is sufficient show that each $\NepsL{r}$ contains exactly one JGE value of ${\mathcal W}$.

To this end let $ {\mathcal E} ={\mathcal W}-{\mathcal T}$ and for $\alpha \in [0,1]$ set ${\mathcal T}^\alpha={\mathcal T}+\alpha  {\mathcal E} $ so that ${\mathcal T}={\mathcal T}^0$ and ${\mathcal W}={\mathcal T}^1$. Note that Proposition \ref{proposition:RegularityBound} shows that ${\mathcal T}^\alpha$ is slice mix invertible and has border $ {\mathbb{C}} $-rank $R$ for each $\alpha \in [0,1]$. It follows from Theorem \ref{theorem:rankvsmult} that each ${\mathcal T}^{\alpha}$ has $R$ JGE values in $ {\mathbb{C}} ^2$ counting algebraic multiplicity.

Now define
\[
\NepsLal{r}= \{ {{\mathscr Z}} \subset  {\mathbb{C}} ^2 | \ \dim({{\mathscr Z}})=1 \ \mathrm{and} \ \chi({{\mathscr L}}_r,{{\mathscr Z}}) \leq \alpha \epsilon\}.
\]
Then $\NepsLal{r} \subset \NepsL{r}$ for each $\alpha \in [0,1]$ and $r$ so $\NepsLal{i} \cap \NepsLal{j} =\emptyset$ whenever $i \neq j$. Furthermore, using Theorem \ref{theorem:TensBauerFike} we have
\[
{\mathrm{sv}}[{\mathcal T},{\mathcal T}^\alpha] \leq \frac{\spnorm{\alpha  {\mathcal E} }}{\sigma_{\min} ( {\mathbf{A}} ) \sigma_{\min} ( {\mathbf{B}} )}  =\alpha \epsilon.
\]
It follows that the spectrum of ${\mathcal T}^\alpha$ is contained in $\cup_{r=1}^R \NepsLal{r}$ for each $\alpha \in [0,1]$.

In particular $\NepsLzero{r}=\{{{\mathscr L}}_r\}$ contains exactly one JGE value of ${\mathcal T}^0={\mathcal T}$ for each $r=1,\dots, R$. Since ${\mathcal T}^\alpha$ is slice mix invertible for all $\alpha \in [0,1]$, \cite[Chapter VI, Theorem 2.1]{SS90} shows that the JGE values of ${\mathcal T}^\alpha$ are continuous in $\alpha$ for $\alpha \in [0,1]$. Therefore, because the $\NepsLal{r}$ remain disjoint as $\alpha$ increases from $0$ to $1$, the continuity of the JGE values of ${\mathcal T}^\alpha$ in $\alpha$ implies that each $\NepsLal{r}$ must contain exactly one JGE value of ${\mathcal T}^\alpha$ for each $\alpha \in [0,1]$ and for each $r=1,\dots,R$. In particular, $\NepsLone{r}=\NepsL{r}$ contains exactly one eigenvalue of ${\mathcal T}^1={\mathcal W}$ for each $r$. By definition, each $\NepsL{r}$ is a ball of chordal radius $\epsilon$ centered ${{\mathscr L}}_r$, so we conclude
\[
{\mathrm{md}}[{\mathcal T},{\mathcal W}] \leq \epsilon = \frac{\spnorm{{\mathcal T}-{\mathcal W}}}{\sigma_{\min} ( {\mathbf{A}} ) \sigma_{\min} ( {\mathbf{B}} )},
\]
and that all the JGE value of ${\mathcal W}$ have algebraic multiplicity one, as claimed.

It remains to show that ${\mathcal W}$ has $ {\mathbb{C}} $-rank $R$.  As previously discussed, all JGE values of a slice mix invertible $R \times R \times 2$ tensor must have geometric multiplicity at least $1$. It follows that ${\mathcal W}$ is simple, so Theorem \ref{theorem:rankvsmult} \eqref{it:NonderogatoryImpliesRankR} shows that ${\mathcal W}$ has $ {\mathbb{C}} $-rank $R$.

Now assume that $ {\mathbb{K}} = {\mathbb{R}} $ and that ${\mathcal T}$ and ${\mathcal W}$ have real entries. Applying the result over the complex field shows that ${\mathcal W}$ must have distinct JGE values over $ {\mathbb{C}} $ so it is sufficient to show that the JGE values of ${\mathcal W}$, hence the JGE vectors of ${\mathcal W}$, must all be real. Here, when we say that ${\mathcal W}$ has real JGE values, we mean that the JGE values of ${\mathcal T}$ can be expressed as the spans of real valued vectors. Equivalently, the classical generalized eigenvalues of the pencil $( {\mathbf{W}} _1, {\mathbf{W}} _2)$ are real valued.

Suppose towards a contradiction that ${\mathcal W}$ has a JGE value which is not real and set
\[
\alpha_*= \inf_{\alpha \in [0,1]} {\mathcal T}+\alpha  {\mathcal E}  {\mathrm{\ has \ a \ JGE \ value \ which \ is \ not \ real }}.
\]
Using the continuity for JGE values in $\alpha$, it is straightforward to show that ${\mathcal T}+\alpha_0  {\mathcal E} $ has a JGE value with algebraic multiplicity at least two. However, ${\mathcal T}+\alpha_0  {\mathcal E} $ satisfies the assumptions of Theorem \ref{theorem:TensMDBound}, so this is a contradiction as we have already shown that all JGE values of such a tensor must have algebraic multiplicity equal to one. We conclude that the JGE values and JGE vectors of ${\mathcal W}$ are {all real}. It follows from Theorem \ref{theorem:rankvsmult} \eqref{it:RankIFFBasis} that ${\mathcal W}$ has $ {\mathbb{R}} $-rank $R$.
\end{proof}

\section{Deterministic bounds for existence of best rank $R$ decompositions}
\label{sec:Bounds}

We now turn to proving the deterministic bounds for existence of best $ {\mathbb{K}} $-rank $R$ approximations presented in the introduction of the article. While we do not have a matching distance bound in the case of general $K$, a key observation is that if the tensor ${\mathcal T}$ is defective in the sense of JGE values, then every subpencil of ${\mathcal T}$ must have a JGE value of algebraic multiplicity greater than one.

\begin{proposition}
\label{proposition:DefectiveSubpencils}
Let ${\mathcal T} \in  {\mathbb{K}} ^{R \times R \times K}$  be a slice mix invertible tensor. Also let $ {\mathbf{U}}  \in  {\mathbb{K}} ^{K \times K}$ be any invertible matrix and set ${\mathcal S} = {\mathcal T} \cdot_3  {\mathbf{U}} $. If ${\mathcal T}$ has a JGE value of algebraic multiplicity $m$, then for any $i \neq j$ the subpencil
\[
\big( {\mathbf{S}} _i, {\mathbf{S}} _{{ j}} \big)
\]
must have a JGE value of algebraic multiplicity at least equal to $m$.
\end{proposition}

Before giving the proof we give a lemma which relates the joint generalized eigenvalues and eigenvectors of ${\mathcal T}$ and ${\mathcal T} \cdot_3  {\mathbf{U}} $.

\begin{lemma}
\label{lemma:EigenvalueModes}
Let ${\mathcal T} \in  {\mathbb{K}} ^{R \times R \times K}$ be a slice mix invertible tensor and let $\mathrm{span}( {\vec{\lambda}} )$ be a JGE value of ${\mathcal T}$. Also let $ {\mathbf{U}}  \in  {\mathbb{K}} ^{K \times K}$ be an invertible matrix and set ${\mathcal S} = {\mathcal T} \cdot_3  {\mathbf{U}} $. Then ${\mathcal S}$ is slice mix invertible  and $\mathrm{span}( {\mathbf{U}}   {\vec{\lambda}} )$ is a JGE value of ${\mathcal S}$ which satisfies $\mathrm{am}( {\vec{\lambda}} )= \mathrm{am}({ {\mathbf{U}} }  {\vec{\lambda}} )$ where $\mathrm{span}( {\vec{\lambda}} )$ and $\mathrm{span}( { {\mathbf{U}} }  {\vec{\lambda}} )$ are treated as JGE values of ${\mathcal T}$ and ${\mathcal S}$, respectively.

Furthermore, if $\mathrm{span}( {\vec{\lambda}} )$ has a corresponding JGE vector $ {\vec{x}} $, then $(\mathrm{span}({ {\mathbf{U}} }  {\vec{\lambda}} ), {\vec{x}} )$ is a JGE pair of ${\mathcal S}$. As a consequence one has $\mathrm{gm}( {\vec{\lambda}} ) = \mathrm{gm}({ {\mathbf{U}} }  {\vec{\lambda}} )$.
\end{lemma}

\begin{proof}
Assume $\mathrm{span}( {\vec{\lambda}} )$ is a JGE value of ${\mathcal T}$ with algebraic multiplicity $m$. View $ {\boldsymbol{\gamma}}  = (\gamma_1, \dots, \gamma_K)$ as a vector in $ {\mathbb{K}} ^K$, and consider the change of variables $ {\boldsymbol{\gamma}}  \to  {\mathbf{U}} ^{{ T}} {\boldsymbol{\gamma}} $. It is a straightforward to verify that
\[
p_{{\mathcal T} \cdot_3  {\mathbf{U}} } ( {\boldsymbol{\gamma}} )= p_{{\mathcal T}} ( {\mathbf{U}} ^{{ T}}  {\boldsymbol{\gamma}} ).
\]
 Since $ {\mathbf{U}} $ is invertible, it follows that $p_{{\mathcal S}} ( {\boldsymbol{\gamma}} ) \not\equiv 0$, hence ${\mathcal S}$ is slice mix invertible. Now set $ {\mathbf{U}} ^{{ T}}  {\boldsymbol{\gamma}}  = \eta$. By assumption $p_{{\mathcal T}} (\eta)=p_{{\mathcal T} \cdot_3  {\mathbf{U}} } ( {\boldsymbol{\gamma}} )$ factors as
\[
p_{{\mathcal T}} (\eta)= \langle  {\vec{\lambda}} , \bar{\eta} \rangle ^m g(\eta) = \langle {  {\mathbf{U}} }  {\vec{\lambda}} , \bar{ {\boldsymbol{\gamma}} } \rangle ^m g( {\mathbf{U}} ^{{ T}}  {\boldsymbol{\gamma}} ),
\]
so $\mathrm{span}({  {\mathbf{U}}   {\vec{\lambda}} })$ is a JGE value of ${\mathcal T} \cdot_3  {\mathbf{U}} $ with algebraic multiplicity $m$.

The proof for geometric multiplicity is straightforward.
\end{proof}

We now give the proof of Proposition \ref{proposition:DefectiveSubpencils}.

\begin{proof}
By assumption ${\mathcal T}$ is slice mix invertible and has a JGE value of algebraic multiplicity $m$, so Lemma \ref{lemma:EigenvalueModes} shows that ${\mathcal S}$ is slice mix invertible and has a JGE value with algebraic multiplicity $m$. Say this JGE value is equal to $\mathrm{span}( {\vec{\lambda}} )$. It follows that the characteristic polynomial of ${\mathcal S}$ factors as
\[
p_{{\mathcal S}}( {\boldsymbol{\gamma}} )=\langle  {\vec{\lambda}} , \bar{ {\boldsymbol{\gamma}} } \rangle^m g( {\boldsymbol{\gamma}} )
\]
where $g( {\boldsymbol{\gamma}} )$ is a nonzero polynomial.

To complete the proof, assume $i < j$ and observe that the characteristic polynomial of the subpencil $( {\mathbf{S}} _i, {\mathbf{S}} _{{ j}})$ is given by
\[
\begin{array}{rclcl}
p_{( {\mathbf{S}} _i, {\mathbf{S}} _{{ j}})}((\gamma_{i},\gamma_{{ j}}))&=&p_{{\mathcal S}} (0,\dots, 0, \gamma_i,0 \dots, 0, \gamma_{{ j}},0 \dots, 0)\\
&=&(\bar{\gamma}_i \lambda_i +\bar{\gamma}_{{ j}} \lambda_{{ j}})^m g(0,\dots, 0, \gamma_i,0 \dots, 0, \gamma_{{ j}},0 \dots, 0)
\end{array}
\]
hence $\mathrm{span}((\lambda_i,\lambda_{{ j}}))$ is a generalized eigenvalue of the subpencil $( {\mathbf{S}} _i, {\mathbf{S}} _{{ j}})$ with multiplicity at least $m$.
\end{proof}

\subsection{Proof of Theorem \ref{theorem:MultiplePencilBound}}
\label{sec:MultiplePencilBoundProof}
We are now in position to prove our multiple pencil based bound for the existence of a best $ {\mathbb{K}} $-rank $R$ approximation.

\begin{proof} Let ${{\mathcal M}}' \in  {\mathbb{K}} ^{I_1 \times I_2 \times I_3}$ be any tensor such that
\[
\|{\mathcal T}'-{{\mathcal M}}'\|_\mathrm{F} {\  < \ } \epsilon/2
\]
and let ${\hat{\cT}}'$ be a best border $ {\mathbb{K}} $-rank $R$ approximation of ${{\mathcal M}}'$ in the Frobenius norm. Since ${\mathcal T}'$ has rank $R$ by assumption, we must have
\[
\|{\hat{\cT}}'-{{\mathcal M}}'\|_\mathrm{F} \leq \|{\mathcal T}'-{{\mathcal M}}'\|_\mathrm{F} \qquad \mathrm{hence} \qquad \|{\hat{\cT}}'-{\mathcal T}'\|_\mathrm{F}  {\  < \ } \epsilon.
\]

Since ${\hat{\cT}}'$ has border $ {\mathbb{K}} $-rank $R$, equation \eqref{eq:BRankVsMLRank} shows that the multilinear rank of ${\hat{\cT}}'$ must be $(R,R,\hat{K})$ for some $\hat{K} \leq R$. It is easy to check that a tensor has $ {\mathbb{K}} $-rank $R$ if and only if an orthogonal compression of the tensor has $ {\mathbb{K}} $-rank $R$, so using Theorem \ref{theorem:TensorProcrustes}, we may WLOG assume that ${\mathcal T}'={\mathcal T}$ and ${\hat{\cT}}'={\hat{\cT}}$ both have size $R \times R \times \max\{K,\hat{K}\}$. There is  little difference in the $K \geq \hat{K}$ and $K< \hat{K}$ cases, so for ease of notation we assume $K \geq \hat{K}$.

Now let $ {\mathbf{U}}  \in  {\mathbb{K}} ^{K \times K}$ be a unitary and set ${\mathcal S}={\mathcal T}\cdot_3  {\mathbf{U}} $ and ${\hat{\cS}}={\hat{\cT}}\cdot_3  {\mathbf{U}} $. Since modal unitary multiplications preserve the Frobenius norm of a tensor we have $\|{\mathcal S} - { {\hat{\cS}}}\|_\mathrm{F}  = \| {\mathcal T} -{\hat{\cT}} \|_\mathrm{F} < \epsilon$. Therefore, for some $i=1,\dots, \lfloor K/2 \rfloor$ we must have that
\[
 \|( {\mathbf{S}} _{2i-1}, {\mathbf{S}} _{2i})-({ \hat{ {\mathbf{S}} }}_{2i-1}, { \hat{ {\mathbf{S}} }}_{2i})\|_\mathrm{F} < \epsilon_i,
\]
as otherwise would imply
\[
\|{\mathcal S} -  {\hat{\cS}}\|_\mathrm{F}^2 \geq  \sum_{k=1}^{\lfloor K/2 \rfloor} \|( {\mathbf{S}} _{2k-1},  {\mathbf{S}} _{2k})-({\hat{ {\mathbf{S}} }}_{2k-1}, {\hat{ {\mathbf{S}} }}_{2k})\|^2_\mathrm{F} \geq   \sum_{k=1}^{\lfloor K/2 \rfloor} \epsilon_k^2 = \epsilon^2.
\]

In particular, from our choice of $\epsilon_i$ the subpencil $({\hat{ {\mathbf{S}} }}_{2i-1},{\hat{ {\mathbf{S}} }}_{2i})$ must be simple. Using Proposition \ref{proposition:DefectiveSubpencils} allows us to conclude that the tensor ${{\hat{\cT}}}$ is simple. The result then follows from Proposition \ref{proposition:existenceCor}.
\end{proof}

One may notice that there is no guarantee that the tensor ${\mathcal T}\cdot_3  {\mathbf{U}} $ has a slice mix invertible subpencil  for a particular choice of $ {\mathbf{U}} $, hence the bound given by Theorem \ref{theorem:MultiplePencilBound} may be equal to $0$ for a particular choice of $ {\mathbf{U}} $. This does not cause serious difficulty: As a consequence of \cite[Section 4.4 (ii)]{DD14}, if ${\mathcal T}$ is slice mix invertible then for generic unitaries every subpencil of ${\mathcal T}$ is slice mix invertible, hence generically the computed bound is nonzero.

This does point toward a greater issue, however. Different choices of unitaries lead to different values for the bound computed by Theorem \ref{theorem:MultiplePencilBound}. Of course, the result holds for each of these values, so one desires to use a unitary which gives the largest possible value of $\epsilon$. Outside of a few special cases, is not known how to optimally pick the unitary $ {\mathbf{U}} $.  Determining how to optimally pick $ {\mathbf{U}} $ is a direction for future research.

A subproblem of picking an optimal unitary is picking the optimal ordering of slices for a fixed unitary $ {\mathbf{U}} $. E.g. if ${\mathcal T} \in  {\mathbb{K}} ^{4 \times 4 \times 4}$ and $ {\mathbf{U}}  \in K^{4 \times 4}$, then for ${\mathcal S}={\mathcal T}\cdot_3  {\mathbf{U}}  \in  {\mathbb{K}} ^{4 \times 4 \times 4}$ one may arrive at different bounds by considering collections of subpencils of the form $\{( {\mathbf{S}} _1,  {\mathbf{S}} _2), ( {\mathbf{S}} _3, {\mathbf{S}} _4)\}$ or $\{( {\mathbf{S}} _1, {\mathbf{S}} _3), ( {\mathbf{S}} _2, {\mathbf{S}} _4)\}$ or $\{( {\mathbf{S}} _1, {\mathbf{S}} _4), ( {\mathbf{S}} _2, {\mathbf{S}} _3)\}$.

Let $ {\mathbf{E}} _ {\mathbf{U}} $ be the matrix with $i,j$ entry equal to the bound given by Theorem \ref{theorem:TensMDBound} for the pencil $( {\mathbf{S}} _i, {\mathbf{S}} _j)$. Then choosing an optimal collection of subpencils (i.e. an optimal ordering for the slices of ${\mathcal S}$) is equivalent to marking the entries of $ {\mathbf{E}} _ {\mathbf{U}} $ so that there is at most one mark in each row and column of $ {\mathbf{E}} _ {\mathbf{U}} $ and so that the norm of the vector of marked entries is maximal.

This problem is  a variant of the longest path problem which is, in a sense, the ``opposite" of the travelling salesman problem (where the total ``distance" is minimized). The longest path problem is known to be NP-hard \cite{UU04}; for our purposes and for small $K$, we may nevertheless consider trying all possible pencil collections/paths.

\section{Numerical experiments}
\label{sec:figures}

In this section, numerical experiments are used to examine the quality of the main perturbation bounds obtained in the paper, i.e. the spectral variation bound obtained in Theorem \ref{theorem:TensBauerFike} and the existence bounds obtained in Theorem \ref{theorem:MultiplePencilBound} and Theorem \ref{theorem:MLSVDExistenceBound}. In addition we illustrate the effect of a tensor Procrustes calculation as in Theorem \ref{theorem:TensorProcrustes} on the error between a rank $R$ tensor and a rank $R$ approximation of a perturbation of the tensor. Throughout the section, we define the signal to noise ratio (SNR) of a signal tensor ${\mathcal T}'$ to a noise tensor $ {\mathcal N} '$ to be $20 \log_{10} \left(\frac{\|{\mathcal T}'\|_\mathrm{F}}{\| {\mathcal N} '\|_\mathrm{F}}\right).$ The experiments discussed in this section all restrict to the real setting. That is, we restrict to considering real valued tensors and real rank.

In all of our experiments, rank $R$ tensors of size $I \times I \times I$ are randomly generated by generating factor matrices $ {\mathbf{A}} ', {\mathbf{B}} ', {\mathbf{C}} '$ of size $I \times R$ with entries independently sampled from a standard normal distribution. We then scale the factor matrices so that each rank $R$ tensor has unit Frobenius norm. The entries of the noise tensors $ {\mathcal N} '$ are all independently sampled from a standard normal distribution and each $ {\mathcal N} '$ is scaled so that the resulting SNR is as specified.

\subsection{Observed spectral variation vs. spectral variation bound}

We first examine the spectral variation bound obtained from the spectral norm\footnote{As previously mentioned, computation of the spectral norm is NP-hard. Using Tensorlab \cite{VDSBL16}, one may compute a (local) optimum for the best rank $1$ approximation of the transformed error. The Frobenius norm of a best rank $1$ approximation to a tensor is equal to the spectral norm of the tensor \cite{FMPS13}.} of (a transformation of) the error tensor in Theorem \ref{theorem:TensBauerFike} and compare it to the observed spectral variation occurring under a perturbation with a specified SNR. Three cases are treated: First, we consider the case of structured perturbations, i.e. perturbations only affecting the factor matrix $ {\mathbf{C}} $. In this special case we may compare to the improved bound presented in equation \eqref{eq:JGEBauerFikeSpecial}. Second, we consider the challenging case of generic perturbations on a rank $R$ tensor of size $R \times R \times R$. Finally we consider the situation most commonly occurring in applications of tensors with rank significantly smaller than dimension. That is, we examine spectral variation bounds for rank $R$ tensors of size $I \times I \times I$ with $I>\!\!>R$.

In each experiment we present  the mean observed spectral variation and either the mean computed spectral variation bound $\sqrt{R} \spnorm{ {\mathcal E}  \cdot_1  {\mathbf{A}} ^{-1} \cdot_2  {\mathbf{B}} ^{-2}}$ or, when appropriate, the mean computed structured spectral variation bound $\spnorm{ {\mathcal E}  \cdot_1  {\mathbf{A}} ^{-1} \cdot_2  {\mathbf{B}} ^{-2}}$.

\subsubsection{Structured Perturbations}

Figure \ref{figure:SVStructPertRcombo} considers the structured perturbation case on tensors of rank $R$ and size $R \times R \times R$ where $R = 4$ or $R = 10$. For various choices of SNR we randomly generate a rank $R$ tensor ${\mathcal T}={\cp{ {\mathbf{A}} , {\mathbf{B}} , {\mathbf{C}} }}$ and an error matrix $ {\mathbf{E}} $ with entries independently sampled from the standard normal distribution. The error matrix $ {\mathbf{E}} $ is then scaled so that the SNR of $ {\mathbf{C}} $ to $ {\mathbf{E}} $ is as specified. A total of 20 trials are performed for each SNR and for each SNR we report the mean observed spectral variation between ${\mathcal T}$ and ${\hat{\cT}}:= [\![  {\mathbf{A}} , {\mathbf{B}} ,  {\mathbf{C}} + {\mathbf{E}} ]\!]$ in addition to the mean computed value of $\| {\mathcal E}  \cdot_1  {\mathbf{A}} ^{-1}  \cdot_2  {\mathbf{B}} ^{-1}  \|_{\mathrm{sp}}$. Here $ {\mathcal E}  := {\hat{\cT}} - {\mathcal T}$. This experiment illustrates the first inequality in equation \eqref{eq:JGEBauerFikeSpecial}.

These figures illustrate that the error bound $\| {\mathcal E}  \cdot_1  {\mathbf{A}} ^{-1} \cdot_2  {\mathbf{B}} ^{-1}  \|_{\mathrm{sp}}$ is highly predictive of the observed spectral variation when considering structured perturbations. We remind the reader that the error bound $\| {\mathcal E}  \cdot_1  {\mathbf{A}} ^{-1}  \cdot_2  {\mathbf{B}} ^{-1}  \|_{\mathrm{sp}}$ is guaranteed to be an upper bound in this special setting.

\begin{figure}[h!]
\hspace*{-0.5cm}
	\input{SVStructPertRcomboSmall.tex}
	\caption{Observed spectral variation and computed bound between tensors of the form ${\cp{ {\mathbf{A}} , {\mathbf{B}} , {\mathbf{C}} }}$ and $[\![ {\mathbf{A}} , {\mathbf{B}} , {\mathbf{C}} + {\mathbf{E}} ]\!]$ with rank $R=4$ or $10$. The structured bound $\| {\mathcal E}  \cdot_1  {\mathbf{A}} ^{-1}  \cdot_2  {\mathbf{B}} ^{-1}  \|_{\mathrm{sp}}$ is observed to be highly predictive of the observed spectral variation.}
	\label{figure:SVStructPertRcombo}
\end{figure}

\subsubsection{Generic Perturbations: Rank equal to dimensions}

Figure \ref{figure:SVGenPertRcombo} illustrates the case of generic perturbations on tensors of rank $R$ and size $R \times R \times R$ where $R = 4$ or $R = 10$. In each trial we randomly generate a rank $R$ tensor ${\mathcal T}$ and noise $ {\mathcal N} $ at the specified SNR, then compute a rank $R$ approximation ${\hat{\cT}}$ of ${\mathcal T}+ {\mathcal N} $ using Tensorlab \cite{VDSBL16}. A total of 50 trials are performed for each SNR, and for each SNR we report the mean observed spectral variation between ${\mathcal T}$ and ${\hat{\cT}}$ as well as the mean computed value for $\sqrt{R} \| {\mathcal E}  \cdot_1  {\mathbf{A}} ^{-1}  \cdot_2  {\mathbf{B}} ^{-1}  \|_{\mathrm{sp}}$ where $ {\mathcal E}  := {\hat{\cT}} - {\mathcal T}$. This experiment illustrates the first inequality in equation \eqref{eq:JGEBauerFike}.

This case is particularly challenging since the perturbation is applied directly to the core of the tensor, which is not typically the case in applications. Furthermore, our bound serves as a worst case scenario bound which is unlikely to be achieved under generic perturbations.

\begin{figure}[h!]
\hspace*{-0.5cm}
	\input{SVGenPertRcomboSmall.tex}
	\caption{Observed spectral variation and computed bound for a rank $R$ tensor ${\mathcal T}$ and a rank $R$ approximation to ${\mathcal T}+ {\mathcal N} $ where $R=4$ or $10$. The bound $\sqrt{R}\| {\mathcal E}  \cdot_1  {\mathbf{A}} ^{-1}  \cdot_2  {\mathbf{B}} ^{-1}  \|_{\mathrm{sp}}$ typically exceeds the observed spectral variation by about one order of magnitude when $R=4$ and by about two orders of magnitude when $R=10$. For perspective we note that $\log_{10} (\sqrt{4}) \approx 0.3$ while $\log_{10} (\sqrt{10}) = 0.5.$}
	\label{figure:SVGenPertRcombo}
\end{figure}

\subsubsection{Generic Perturbations: Rank less than dimensions}

Figures \ref{figure:svscdat4combo} and \ref{figure:svscdat10combo} consider the case which most commonly occurs in applications where the rank of the tensor of interest is significantly less than dimensions of the tensor. Here we work with tensors $ {\mathcal T}'$ of rank $R$ and dimensions $I \times I \times I$ where either $R = 4$ and $I = 10,20,$ or $100$ or $R = 10$ and $I = 20$ or $100$.

Let ${\hat{\cT}}'$ denote a rank $R$ approximation of ${\mathcal T}'+ {\mathcal N} '$. To directly apply the spectral variation bound obtained from the spectral norm in Theorem \ref{theorem:TensBauerFike} to a pair of rank $R$ tensors, the factor matrices of at least one tensor should be invertible. In particular the tensors should have dimensions equal to their rank. Thus it is most appropriate to apply this bound to core tensors ${\mathcal T}$ and ${\hat{\cT}}$ of size $R \times R \times R$ satisfying $\|{\mathcal T} - {\hat{\cT}} \|_\mathrm{F} \leq \|{\mathcal T}' - {\hat{\cT}}'\|_\mathrm{F}$ as described in the tensor Procrustes problem of Theorem \ref{theorem:TensorProcrustes}. Write ${\mathcal T} = [\![ {\mathbf{A}} , {\mathbf{B}} , {\mathbf{C}} ]\!]$. For each SNR we report the mean observed spectral variation between the core tensors ${\mathcal T}$ and ${\hat{\cT}}$ as well as the mean computed value for $\sqrt{R} \| {\mathcal E}  \cdot_1  {\mathbf{A}} ^{-1} \cdot_2  {\mathbf{B}} ^{-1} \|_{\mathrm{sp}}$ where $ {\mathcal E}  := {\hat{\cT}} - {\mathcal T}$. For each choice of $I$ and $R$ the mean is computed over $20$ trials.

In this setting we find that the error bound $\sqrt{R} \| {\mathcal E}  \cdot_1  {\mathbf{A}} ^{-1} \cdot_2  {\mathbf{B}} ^{-1} \|_{\mathrm{sp}}$ from the first inequality in equation \eqref{eq:JGEBauerFike} is highly predictive of the observed spectral variation, with the error bound typically within a half order of magnitude of the observed spectral variation. Furthermore, we find that the observed spectral variation and spectral variation bound become remarkably small as $I$ grows, even for low SNR. It is well known that computation of the core of an observed tensor can significantly increase the effective SNR. Our results further suggest that a significant portion of the error remaining after computing a rank $R$ approximation of an observed tensor ${\mathcal T}'+ {\mathcal N} '$ is in the form of rotational error (i.e. error accounted for by unitary rotations of the mode-$j$ subspaces of the tensor) which disappears after a Tensor Procrustes computation.

While this rotational error is still meaningful from the perspective of quality of the estimated factors, rotational error cannot lead to nonexistence of a best rank $R$ approximation and cannot negatively impact the condition number of computing a CPD (see \cite{V17} for details on this condition number), as these issues are invariant under unitary changes of basis. The observation that, in the setting of independent and identically distributed (i.i.d) noise, a significant portion error observed between the original tensor ${\mathcal T}'$ and a rank $R$ approximation to the observed tensor ${\mathcal T}'+ {\mathcal N} '$ is merely rotational error further explains the success of CPD in many large scale settings. In particular, this observation illustrates that CPD computation in practice is likely to be more stable than one might naively expect. See Section \ref{sec:TensorProcrustesExperiments} for experiments in this direction.

\begin{figure}[h!]
\hspace*{-1cm}
	\input{svscdat4comboSmall.tex}
		\caption{Observed spectral variation and computed bound for cores ${\mathcal T}$ and ${\hat{\cT}}$ where ${\mathcal T}'$ is a rank $4$ tensor of size $I \times I \times I$ with $I = 10,20$ or $100$ and ${\hat{\cT}}'$ is a rank $4$ approximation to ${\mathcal T}'+ {\mathcal N} '$. The bound $\sqrt{R} \| {\mathcal E}  \cdot_1  {\mathbf{A}} ^{-1} \cdot_2  {\mathbf{B}} ^{-1} \|_{\mathrm{sp}}$ provides a strong estimate for the observed spectral variation. The spectral variation bound and observed spectral variation both decrease as $I$ grows.}
	\label{figure:svscdat4combo}
\end{figure}

\begin{figure}[h!]
\hspace*{-0.5cm}
	\input{svscdat10comboSmall.tex}
		\caption{Observed spectral variation and computed bound for cores ${\mathcal T}$ and ${\hat{\cT}}$ where ${\mathcal T}'$ is a rank $10$ tensor of size $I \times I \times I$ with $I = 20$ or $100$ and ${\hat{\cT}}'$ is a rank $10$ approximation to ${\mathcal T}'+ {\mathcal N} '$. The bound $\sqrt{R} \| {\mathcal E}  \cdot_1  {\mathbf{A}} ^{-1} \cdot_2  {\mathbf{B}} ^{-1} \|_{\mathrm{sp}}$ provides a strong estimate for the observed spectral variation. The spectral variation bound and observed spectral variation both decrease as $I$ grows.}
	\label{figure:svscdat10combo}
\end{figure}

\subsection{Bounds for existence of best rank $R$ approximations}

We now present the mean of the bound $\epsilon$ obtained using Theorem \ref{theorem:MultiplePencilBound} for the radius of a neighborhood around a randomly generated rank $R$ tensor of size $R \times R \times R$ of unit Frobenius norm in which best rank $R$ approximations are guaranteed to exist. For each choice of $R$, the bound $\epsilon$ is averaged over $20$ trials to compute the mean bound $ {{\epsilon_{\mathrm{avg}}}}$. The mean bound $ {{\epsilon_{\mathrm{avg}}}}$ is plotted in dB as $-20 \log_{10} ( {{\epsilon_{\mathrm{avg}}}})$. Thus, on average for the tensors in our experiment, a best rank $R$ approximation is guaranteed to exist if one has SNR greater than the reported value. The choices of $R$ considered are $R=2$ through $10$.

As previously discussed, the bound obtained from Theorem \ref{theorem:MultiplePencilBound} depends on the choice of unitary $ {\mathbf{U}} $ chosen when setting ${\mathcal S} = {\mathcal T} \cdot_3  {\mathbf{U}} $. While finding an optimal unitary is computationally challenging, as evidenced by its relationship to the longest path problem, one may simply try many randomly generated unitaries. For each tensor ${\mathcal T}$ we try a total of $1,000$ random unitary matrices  and record the largest bound found after $1,10,100,$ and $1000$ unitaries tried. As a further step, one may find the optimal ordering of tensor slices for each unitary. Bounds computed without reordering slices are presented in the left image of Figure \ref{figure:AvgChordalBoundPlot} while bounds computed with finding the optimal slice ordering are presented in the right image of the same figure. The collection of tensors and unitaries tried in each case are the same.

\begin{figure}[h!]
\hspace*{-0.7cm}
	\input{AvgChordalBoundPlotScalingFixSmall.tex}
	\caption{Mean bound for radius of a ball centered at a slice mix invertible rank $R$ tensor of unit Frobenius norm in which best rank $R$ approximations are guaranteed to exist. The minimum SNR needed for existence of a best low rank approximation to be guaranteed increases as $R$ increases from $2$ to $10$. Increasing the number of unitaries tried can notably improve the bound; however, there are diminishing returns as the number of unitaries tried increases. }
	\label{figure:AvgChordalBoundPlot}
\end{figure}

\subsubsection{Bounds for existence in the case of rank less than dimensions}
\label{sec:existenceImagesBigDim}

We now examine the setting more common in applications where a tensor has rank significantly less than its dimensions. We first note that since the bound obtained from Theorem \ref{theorem:MultiplePencilBound} is a deterministic bound, the bound must take into account a worst case scenario and does not improve as tensor dimensions increase. A worst case scenario achievable for the core of a tensor is always achievable on the full tensor by applying a perturbation directly to the core of the tensor.

Although the deterministic bound given by Theorem \ref{theorem:MultiplePencilBound} does not improve as tensor dimensions grow, {\it the probability of a worst case scenario occurring dramatically decreases.} This is illustrated in Figure \ref{figure:ExistPro} with the following experiment.

For various choices of $R$ and $I$ and SNR, we randomly generate a rank $R$ tensor ${\mathcal T}'$ and a perturbation $ {\mathcal N} '$ of size $I \times I \times I$ such that the SNR of ${\mathcal T}'$ to $ {\mathcal N} '$ is as specified. Then, as is commonly done in practice, we compute an MLSVD of ${\mathcal T}'+ {\mathcal N} '$ and truncate the MLSVD to obtain a multilinear rank $R \times R \times R$ approximation of ${\mathcal T}'+ {\mathcal N} '$. We let ${\mathcal W}$ denote the truncated MLSVD of ${\mathcal T}'+ {\mathcal N} '$.

Theorem \ref{theorem:MLSVDExistenceBound} is used to compute a bound\footnote{For each ${\mathcal W}$, a total of $1000$ randomly generated unitaries are tried when computing this radius. To limit computational complexity, we do not compute the optimal slice orderings. It is expected that either increasing the number of unitaries tried or computing the optimal slice ordering for each unitary would cause the curves in Figure \ref{figure:ExistPro} to shift slightly to the left.} $\epsilon$ such that if there exists a rank $R$ tensor within $\epsilon$ Frobenius distance of ${\mathcal W}$, then ${\mathcal W}$ is guaranteed to have a best rank $R$ approximation. We then compute a rank $R$ approximation ${\hat{\cT}}$ of ${\mathcal W}$ and check if $\|{\mathcal W}-{\hat{\cT}}\|_\mathrm{F} < \epsilon$, in which case we are able to guarantee that ${\mathcal W}$ has a best rank $R$ approximation. For each choice of $I$, $R$, and SNR, a total of $10$ trials are performed and the proportion of truncated MLSVDs guaranteed to have a best rank $R$ approximation is reported.

As expected, we observe that the SNR needed for ${\mathcal W}$ to be guaranteed to have a best rank $R$ approximation greatly decreases as $I$ increases. Additionally, for each choice of $I$ and $R$, we observe a sharp transition from being unable to guarantee the existence of a best rank $R$ approximation of ${\mathcal W}$ to expecting to be able to guarantee the existence of a best rank $R$ approximation of ${\mathcal W}$.

\begin{figure}[h!]
\hspace*{-1cm}
	\input{ExistenceProportionPlotUpdatedSmall.tex}
	\caption{Proportion of $I \times I \times I$ tensors ${\mathcal T}'+ {\mathcal N} '$ for which the truncated MLSVD ${\mathcal W}$ of size $R \times R \times R$ is guaranteed to have a best rank $R$ approximation. The necessary SNR decreases greatly as tensor dimensions grow. For each choice of $I$ and $R$, a sharp transition point is observed. When the SNR exceeds the transition point, the truncated MLSVD is expected to have a best rank $R$ approximation.}
	\label{figure:ExistPro}
\end{figure}

\subsection{0 dB SNR i.i.d noise is not 0 dB SNR noise on the core}
\label{sec:TensorProcrustesExperiments}

In this section we perform an experiment which is intended to help answer the following question: ``Given a tensor ${\mathcal T}'$ and a perturbation $ {\mathcal N} '$, what is the expected distance between a pair of Procrustes cores computed for ${\mathcal T}'$ and a best border rank $R$ approximation to ${\mathcal T}'+ {\mathcal N} '$?" As we now argue, this is a critical measure of distance when considering the existence of a best rank $R$ approximation of ${\mathcal T}'+ {\mathcal N} '$.

Recall the main idea in this article for determining if a tensor ${\mathcal T}'+ {\mathcal N} '$ has a best rank $R$ approximation: The tensor ${\mathcal T}'+ {\mathcal N} '$ is guaranteed to have a best border rank $R$ approximation ${\hat{\cT}}'$, so we need only verify that ${\hat{\cT}}'$ has rank $R$. To determine if ${\hat{\cT}}'$ has rank $R$, we check if ${\hat{\cT}}'$ lies in a neighborhood around ${\mathcal T}'$ whose radius is determined by Theorem \ref{theorem:MultiplePencilBound}.

Moreover, since tensor (border) rank is invariant under modal unitary multiplication, it is sufficient to check that if one computes Procrustes cores ${\mathcal T}$ and ${\hat{\cT}}$ for the pair ${\mathcal T}',{\hat{\cT}}'$, then ${\hat{\cT}}$ lies in the neighborhood around ${\mathcal T}$ given by Theorem \ref{theorem:MultiplePencilBound}. The key point is that \textit{the existence of a best rank $R$ approximation of ${\mathcal T}'+ {\mathcal N} '$ can be determined by checking the distance between ${\mathcal T}$ and ${\hat{\cT}}$.}

The experiment shown in Figure \ref{figure:TPEcombo} is performed as follows.  For various choices of $I$ and $R$ we randomly generate rank $R$ tensors ${\mathcal T}'$ of size $I \times I \times I$ with unit Frobenius norm. For each tensor we generate a noise tensor $ {\mathcal N} '$ and scale $ {\mathcal N} '$ so that the SNR of ${\mathcal T}'$ to $ {\mathcal N} '$ is as specified. After generating ${\mathcal T}'$ and $ {\mathcal N} '$, we compute a rank $R$ approximation\footnote{For the purpose of this experiment, it is not necessary to compute a best rank $R$ approximation.} ${\hat{\cT}}'$ of ${\mathcal T}'+ {\mathcal N} '$ and  the distance $\|{\mathcal T}'-{\hat{\cT}}'\|_\mathrm{F}$ between the tensors $ {\hat{\cT}}'$ and ${\mathcal T}'$. We then compute Procrustes cores ${\mathcal T}$ and ${\hat{\cT}}$ of ${\mathcal T}'$ and ${\hat{\cT}}'$, respectively, for which it is guaranteed that $\|{\mathcal T} -{\hat{\cT}}\|_\mathrm{F} \leq \|{\mathcal T}' - {\hat{\cT}}' \|_\mathrm{F}$ and record the distance $\|{\mathcal T} -{\hat{\cT}}\|_\mathrm{F}$ between the cores. For each choice of $I$ and $R$, a total of $10$ trials are performed and the mean distance between the relevant pairs at each step is reported.

Figure \ref{figure:TPEcombo} illustrates that the distance between the core tensors ${\mathcal T}$ and ${\hat{\cT}}$ is significantly smaller than the distance between ${\mathcal T}'$ and the perturbed tensor ${\mathcal T}'+ {\mathcal N} '$. We observe that the distance between ${\mathcal T}'$ and ${\hat{\cT}}'$ is significantly less than that between ${\mathcal T}'$ and ${\mathcal T}'+ {\mathcal N} '$, and that the distance between ${\mathcal T}$ and ${\hat{\cT}}$ is even smaller still. Furthermore, consistent with Figure \ref{figure:ExistPro}, this experiment shows that when $R <\!\!< I$, issues of nonexistence of a best rank $R$ approximation are unlikely to occur (given our experimental setup of i.i.d factors and i.i.d noise) unless the SNR between ${\mathcal T}'$ and $ {\mathcal N} '$ is significantly lower than suggested by Figure \ref{figure:AvgChordalBoundPlot}.

\begin{figure}[h!]
	\scalebox{.8}{\input{TPEcomboLogSmall.tex}}
		\label{figure:TPEcombo}
	\caption{Distance between ${\mathcal T}'$ and ${\mathcal T}'+ {\mathcal N} '$, between ${\mathcal T}'$ and ${\hat{\cT}}'$, and between ${\mathcal T}$ and ${\hat{\cT}}$. The mean distance between ${\mathcal T}$ and ${\hat{\cT}}$ is observed to be significantly smaller than that between ${\mathcal T}'$ and ${\mathcal T}'+ {\mathcal N} '$, and the difference between these distances increases as the ratio of $I$ to $R$ increases. The mean distance between ${\mathcal T}$ and ${\hat{\cT}}$ is smaller than that between ${\mathcal T}'$ and ${\hat{\cT}}'$; however, the order of magnitude of the gap depends on the initial SNR of ${\mathcal T}'$ to $ {\mathcal N} '$ with smaller gaps when the initial SNR is lower.}
\end{figure}

\section{Conclusion}
We have given a deterministic bound on the radius of an open Frobenius norm ball centered at a simple rank $R$ tensor in which best rank $R$ approximations are guaranteed to exist and in which every rank $R$ tensor has a unique CPD. This illustrates the existence of a neighborhood of ``mathematical truth" around the signal portion of tensors occurring in application and gives a quantitative explanation for the numerical success of the CPD in practical settings. In addition we solved a tensor Procrustes problem which provides a method for computing orthogonal compressions for a pair of tensors.

The bound for existence of best rank $R$ approximations was computed by connecting the CPD to a JGE value decomposition for a tuple of matrices. Several basic results were established for JGE values, including an examination of algebraic and geometric multiplicities. We showed that a slice mix invertible tensor of size $R \times R \times K$ which has  border $ {\mathbb{K}} $-rank $R$ and $ {\mathbb{K}} $-rank strictly greater than $R$ must be defective in the sense of algebraic and geometric multiplicities for JGE values. Furthermore, we developed perturbation theoretic bounds which may be used to guarantee that tensors in a neighborhood of a given simple rank $R$ tensor are simple, thus giving a neighborhood in which every tensor has a best rank $R$ approximation.

\section{Acknowledgements} The authors thank the anonymous referees for their comments which helped improve clarity and exposition in the article.

\renewcommand\thesection{\ifcase\value{section}S\ or S.1\else S\fi}

\section{Supplementary Proofs}

The supplementary materials give the proofs of Theorem \ref{theorem:MLSVDExistenceBound}, Lemma \ref{lemma:ProcrustesIsometries}, and Theorem \ref{theorem:TensBauerFike} in the general case.

\subsection{Proof of Theorem \ref{theorem:MLSVDExistenceBound}}

As previously mentioned, the proof of Theorem \ref{theorem:MLSVDExistenceBound} is nearly identical to that of Theorem \ref{theorem:MultiplePencilBound}, with all notable differences occurring at the start of the proof.

\begin{proof} Let ${{\mathcal M}}',{\mathcal W}',{\mathcal W}$ and $R$ and $K$ be as in the statement of the theorem. Assume there exists some $ {\mathbb{K}} $-rank $R$ tensor ${\tilde{\cT}}'$ such that
\[
\|{{\mathcal M}}-{\tilde{\cT}}'\|_\mathrm{F} < \epsilon - \|{{\mathcal M}}' - {\mathcal W}'\|_\mathrm{F}
\]
and let ${\hat{\cT}}'$ be a best border $ {\mathbb{K}} $-rank $R$ approximation of ${{\mathcal M}}'$ in the Frobenius norm. Since ${\tilde{\cT}}'$ has rank $R$ by assumption, we must have
\[
\|{{\mathcal M}}'-{\hat{\cT}}' \|_\mathrm{F} \leq \|{{\mathcal M}}'-{\tilde{\cT}}'\|_\mathrm{F}  < \epsilon - \|{{\mathcal M}}' - {\mathcal W}'\|_\mathrm{F},
\]
hence
\[
\|{\mathcal W}'-{\hat{\cT}}'\|_\mathrm{F} < \epsilon.
\]

 Note that projecting the mode three fibers of ${\hat{\cT}}'$ to the subspace spanned by the mode three fibers of ${{\mathcal M}}'$ neither increases the Frobenius distance between ${\hat{\cT}}'$ and ${{\mathcal M}}'$ nor the rank of ${\hat{\cT}}'$, thus we can assume $R_3({\hat{\cT}}') \leq R_3 ({{\mathcal M}}')$. Recalling that $K=\min\{({{\mathcal M}}'),R\}$, since ${\hat{\cT}}'$ has border $ {\mathbb{K}} $-rank $R$, the multilinear rank of ${\hat{\cT}}'$ is at most $(R,R,K)$. As consequence of Theorem \ref{theorem:TensorProcrustes}, there exists an orthogonal compression ${\hat{\cT}} \in  {\mathbb{K}} ^{R \times R \times K}$ of ${\hat{\cT}}'$ such that
\[
\|{\mathcal W}-{\hat{\cT}}\|_\mathrm{F} \leq \|{\mathcal W}'-{\hat{\cT}}'\|_\mathrm{F}  < \epsilon
\]
To show ${\hat{\cT}}'$ has $ {\mathbb{K}} $-rank $R$ is sufficient to show that ${\hat{\cT}}$ has $ {\mathbb{K}} $-rank $R$.

Now let $ {\mathbf{U}}  \in  {\mathbb{K}} ^{K \times K}$ be a unitary and set ${\mathcal S}={\mathcal W}\cdot_3  {\mathbf{U}} $ and ${\hat{\cS}}={\hat{\cT}}\cdot_3  {\mathbf{U}} $. Since modal unitary multiplications preserve the Frobenius norm of a tensor we have $\|{\mathcal S}- {\hat{\cS}}\|_\mathrm{F}  = \| {\mathcal W} -{\hat{\cT}} \|_\mathrm{F} < \epsilon $. Therefore, for some $i=1,\dots, \lfloor K/2 \rfloor$ we must have that
\[
\|( {\mathbf{S}} _{2i-1},  {\mathbf{S}} _{2i})-({ \hat{ {\mathbf{S}} }}_{2i-1}, { \hat{ {\mathbf{S}} }}_{2i})\|_\mathrm{F} < \epsilon_i,
\]
as otherwise would imply
\[
\| {\mathcal S} -  {\hat{\cS}}\|_\mathrm{F}^2 \geq  \sum_{k=1}^{\lfloor K/2 \rfloor} \|( {\mathbf{S}} _{2k-1}, {\mathbf{S}} _{2k})-({\hat{ {\mathbf{S}} }}_{2k-1},{\hat{ {\mathbf{S}} }}_{2k})\|^2_\mathrm{F} \geq   \sum_{k=1}^{\lfloor K/2 \rfloor} \epsilon_k ^2=  \epsilon^2 .
\]

In particular, from our choice of $\epsilon_i$ the subpencil $({\hat{ {\mathbf{S}} }}_{2i-1},{\hat{ {\mathbf{S}} }}_{2i})$ must be  simple. Using Proposition \ref{proposition:DefectiveSubpencils} allows us to conclude that the tensor ${{\hat{\cT}}}$ is simple. It follows from Theorem \ref{theorem:rankvsmult} \eqref{it:NonderogatoryImpliesRankR} that ${\hat{\cT}}$ has rank $R$.

The fact that ${\hat{\cT}}$ has a unique CPD is then a straightforward consequence of Kruskal's condition \cite{K77}. In addition, the fact that the set of tensors ${{\mathcal M}}'$ which satisfy the assumptions of the theorem but do not have a unique best rank $R$ approximation has measure zero is an immediate consequence of \cite{QML19} when working over $ {\mathbb{C}} $ and of \cite[Corollary 18]{FO14} when working over $ {\mathbb{R}} $.
\end{proof}

\subsection{Proof of Lemma \ref{lemma:ProcrustesIsometries}}

\begin{proof}
Consider the orthogonal Procrustes problem:
\begin{equation}
\label{eq:Procrustes}
\min_{\Psi \in O(I_1)} \| {\mathbf{M}} - \Psi  {\mathbf{N}} \|_\mathrm{F}.
\end{equation}
where $O(I_1)$ denotes the group of $I_1 \times I_1$ unitary matrices. \cite{S66} shows that $ {\mathbf{U}}   \in O(I_1)$ is a minimizer of \eqref{eq:Procrustes} if and only if $ {\mathbf{U}} $ can be written $ {\mathbf{U}} = {\mathbf{V}}   {\mathbf{Z}} ^\mathrm{H}$ where $ {\mathbf{V}}  , {\mathbf{Z}}  \in O( I_1)$ are unitary matrices such that
\begin{equation}
\label{eq:ProcrustesSVD}
 {\mathbf{M}}   {\mathbf{N}} ^\mathrm{H}= {\mathbf{V}}  \mathbf{\Sigma}  {\mathbf{Z}} ^\mathrm{H}
\end{equation}
gives a singular value decomposition for $ {\mathbf{M}}   {\mathbf{N}} ^\mathrm{H}$ for an appropriate choice of $\mathbf{\Sigma}$.

It is straightforward to show that the unitary matrices $ {\mathbf{V}} $ and $ {\mathbf{Z}} $ in equation \eqref{eq:ProcrustesSVD} may be chosen so that
\[
 {\mathbf{V}}  = ( {\mathbf{V}} _1 \  {\mathbf{V}} _2)  {\quad \quad \mathrm{and} \quad \quad}  {\mathbf{Z}}  = ( {\mathbf{Z}} _1 \  {\mathbf{Z}} _2)
\]
where $ {\mathbf{V}} _1,{  {\mathbf{Z}} _1} \in  {\mathbb{K}} ^{I_1 \times R_ {\mathbf{M}} }$ and $ {\mathbf{V}} _2, {\mathbf{Z}} _2 \in  {\mathbb{K}} ^{I_1 \times (I_1-R_ {\mathbf{M}} )}$ are column-wise orthonormal matrices satisfying
\[
\mathrm{ran \ }  {\mathbf{V}} _1 = \mathrm{ran \ }  {\mathbf{M}}  {\quad \quad \mathrm{and} \quad \quad} \mathrm{ran \ }  {\mathbf{Z}} _1 \supseteq  \mathrm{ran \ } {\mathbf{N}}
\]
\[
\mathrm{ran \ }  {\mathbf{V}} _2  = (\mathrm{ran \ }  {\mathbf{M}} )^\perp {\quad \quad \mathrm{and} \quad \quad} \mathrm{ran \ }  {\mathbf{Z}} _2 \subseteq (\mathrm{ran \ }  {\mathbf{N}} )^\perp.
\]
Using this choice of $ {\mathbf{V}} $ and $ {\mathbf{Z}} $ we have $ {\mathbf{Z}} _2^\mathrm{H}  {\mathbf{N}} = {\vec{0}} $, hence
\[
 {\mathbf{V}}   {\mathbf{Z}} ^\mathrm{H}  {\mathbf{N}} = ( {\mathbf{V}} _1 \  {\mathbf{V}} _2  ) ( {\mathbf{Z}} _1 \  {\mathbf{Z}} _2)^\mathrm{H}  {\mathbf{N}} =  {\mathbf{V}} _1  {\mathbf{Z}} _1^\mathrm{H}  {\mathbf{N}} .
\]
It follows that $\mathrm{ran \ }  {\mathbf{M}}  \supseteq \mathrm{ran \ }  {\mathbf{V}}   {\mathbf{Z}} ^\mathrm{H}  {\mathbf{N}} $ which completes the first part of the proof.

To prove the second claim let $ {\mathbf{U}} _ {\mathbf{M}}  \in  {\mathbb{R}} ^{I_1 \times R_ {\mathbf{M}} }$ be any column-wise orthonormal matrix mapping onto $\mathrm{ran \ }  {\mathbf{M}} $, and set $ {\mathbf{U}} _ {\mathbf{N}}  =  {\mathbf{Z}}   {\mathbf{V}} ^\mathrm{H}  {\mathbf{U}} _ {\mathbf{M}} $ with $ {\mathbf{V}} $ and $ {\mathbf{Z}} $ as above. Then we have
\[
\| {\mathbf{U}} _ {\mathbf{M}} ^\mathrm{H} {  {\mathbf{M}} }-  {\mathbf{U}} _ {\mathbf{N}} ^\mathrm{H}  {\mathbf{N}}  \|_\mathrm{F} \leq \|  {\mathbf{U}} _ {\mathbf{M}} ^\mathrm{H} \|_2 \|  {\mathbf{M}} -   {\mathbf{V}}  {\mathbf{Z}} ^\mathrm{H}  {\mathbf{N}} \|_\mathrm{F}  \leq \| {\mathbf{M}} - {\mathbf{N}} \|_\mathrm{F},
\]
To see that $\mathrm{ran \ }  {\mathbf{U}} _ {\mathbf{N}}  \supseteq \mathrm{ran \ }  {\mathbf{N}} $, we use $ {\mathbf{V}} _2^\mathrm{H}  {\mathbf{M}} = {\vec{0}} $ to find
\[
 {\mathbf{Z}}   {\mathbf{V}} ^\mathrm{H}  {\mathbf{U}} _ {\mathbf{M}}  = ( {\mathbf{Z}} _1 \  {\mathbf{Z}} _2 ) ( {\mathbf{V}} _1 \  {\mathbf{V}} _2)^\mathrm{H}  {\mathbf{U}} _ {\mathbf{M}} = {\mathbf{Z}} _1  {\mathbf{V}} _1^\mathrm{H}  {\mathbf{U}} _ {\mathbf{M}} .
\]
Since $ {\mathbf{U}} _ {\mathbf{M}}  \in  {\mathbb{K}} ^{I_1 \times R_ {\mathbf{M}} }$ and  $ {\mathbf{V}} _1 \in  {\mathbb{R}} ^{I_1 \times R_ {\mathbf{M}} }$ are column-wise orthonormal matrices with the same range, it follows that $ {\mathbf{V}} _1^{\mathrm{H}}  {\mathbf{U}} _ {\mathbf{M}} $ is an invertible map on $ {\mathbb{K}} ^{R_ {\mathbf{M}} }$. Therefore
\[
\mathrm{ran \ }  {\mathbf{U}} _ {\mathbf{N}} = \mathrm{ran \ }  {\mathbf{Z}} _1  {\mathbf{V}} _{1}^\mathrm{H}  {\mathbf{U}} _ {\mathbf{M}}  = \mathrm{ran \ }   {\mathbf{Z}} _1 \supseteq \mathrm{ran \ }  {\mathbf{N}}
\]
as claimed. \ \end{proof}

\subsection{Proof of Theorem \ref{theorem:TensBauerFike} in the general case}

We now give the proof of Theorem \ref{theorem:TensBauerFike} in the general setting. For clarity, the difference between the proof we now give and the proof in the main article is that we now only assume that ${\mathcal W}$ has some JGE value with nonzero geometric multiplicity rather than assuming that ${\mathcal W}$ has $ {\mathbb{K}} $-rank $R$ and is slice mix invertible. This more general setting for example  be applies to slice mix invertible $R \times R \times K$ tensors which have border $ {\mathbb{K}} $-rank $R$ but rank strictly greater than $R$.

\begin{proof}

We begin by proving the first inequality in equation \eqref{eq:JGEvaluebound}. To this end, let ${{\mathscr W}} \subset  {\mathbb{K}} ^K$ be a JGE value of ${\mathcal W}$ and let $ { \vec{\omega}}  \in {{\mathscr W}}$ and $ {\boldsymbol{\gamma}}  \in {{\mathscr W}}^\perp$ be vectors with unit two norm. Let $ {\vec{x}}  \in  {\mathbb{K}} ^R$ be a JGE vector of ${\mathcal W}$ corresponding to ${{\mathscr W}}$ and let $ {\vec{y}}  \in  {\mathbb{K}} ^R$ be a vector such that the tuple $({\mathcal W}, { \vec{\omega}} , {\vec{x}} , {\vec{y}} )$ satisfies equation (2.8). Set $ {\vec{v}}  = \frac{ {\mathbf{B}} ^{\mathrm{T}} {\vec{x}} }{\| {\mathbf{B}} ^{\mathrm{T}} {\vec{x}} \|_2}$ so that $\| {\vec{v}} \|_2=1$. Additionally, for $k=1,\dots,K$ let $ {\mathbf{E}} _k$ denote the $k$th frontal slice of $ {\mathcal E} $.

With this setup we have
\[
 {\vec{0}} = \langle  { \vec{\omega}} , {\boldsymbol{\gamma}}  \rangle  {\mathbf{A}} ^{-1}  {\vec{y}}  = \sum_{k=1}^K \bar{\gamma}_k  {\mathbf{A}} ^{-1}  {\mathbf{W}} _k  {\mathbf{B}} ^{-\mathrm{T}}  {\vec{v}}  = \sum_{k=1}^K \bar{\gamma}_k D_k ( {\mathbf{C}} )  {\vec{v}}  + \sum_{k=1}^K \bar{\gamma}_k  {\mathbf{A}} ^{-1}  {\mathbf{E}} _k  {\mathbf{B}} ^{-\mathrm{T}}  {\vec{v}} .
\]
It follows that
\[
\sum_{k=1}^K \bar{\gamma}_k D_k ( {\mathbf{C}} )  {\vec{v}}  = - \sum_{k=1}^K \bar{\gamma}_k  {\mathbf{A}} ^{-1}  {\mathbf{E}} _k  {\mathbf{B}} ^{-\mathrm{T}}  {\vec{v}} .
\]
Furthermore we have
\begin{align}
\label{eq:JGEspNormBound} \max_{\substack{ {\boldsymbol{\gamma}}  \in {{\mathscr W}}^\perp \\ \| {\boldsymbol{\gamma}}  \|_2=1}} \Big\|\sum_{k=1}^K \bar{\gamma}_k D_k ( {\mathbf{C}} )  {\vec{v}} \Big\|_2
&\leq
\max_{\substack{\| {\vec{u}} \|_2=1 \\ \| {\vec{z}} \|_2=1 }} \max_{\substack{ {\boldsymbol{\gamma}}  \in {{\mathscr W}}^\perp \\ \| {\boldsymbol{\gamma}}  \|_2=1 }} \Big| \sum_{k=1}^K \bar{\gamma}_k   {\vec{u}} ^\mathrm{T} {\mathbf{A}} ^{-1}  {\mathbf{E}} _k  {\mathbf{B}} ^{-\mathrm{T}}  {\vec{z}}  \Big|
\\ \nonumber
& =  \max_{\substack{\| {\vec{w}} \|_2=1 \\ \| {\vec{z}} \|_2=1  \\ \| {\boldsymbol{\gamma}} \|_2=1 }}  \Big| \Big(  {\mathcal E}  \cdot_1  {\mathbf{A}} ^{-1} \cdot_2  {\mathbf{B}} ^{-1} \cdot_3  {\mathbf{P}} _{{{\mathscr W}}^\perp} \Big) \cdot_1  {\vec{u}} ^\mathrm{T} \cdot_2  {\vec{z}} ^\mathrm{T} \cdot_3  {\boldsymbol{\gamma}} ^\mathrm{T} \Big|
\\ \nonumber
& = \Big\|  {\mathcal E}  \cdot_1  {\mathbf{A}} ^{-1} \cdot_2  {\mathbf{B}} ^{-1} \cdot_3  {\mathbf{P}} _{{{\mathscr W}}^\perp} \Big\|_{\mathrm{sp}}.
\end{align}
In the $K=2$ case, one may then follow the proof of \cite[Chapter VI, Theorem 2.7]{SS90} to obtain
\[
\min_k \chi ({{\mathscr L}}_k, {{\mathscr W}}) \leq \Big\|  {\mathcal E}  \cdot_1  {\mathbf{A}} ^{-1} \cdot_2  {\mathbf{B}} ^{-1} \cdot_3  {\mathbf{P}} _{{{\mathscr W}}^\perp} \Big\|_{\mathrm{sp}}.
\]

We next show that there exists some JGE value ${{\mathscr L}}$ of ${\mathcal T}$ such that
\begin{equation}
\label{eq:JGEvalueIntermediateBound}
\frac{\chi({{\mathscr L}},{{\mathscr W}})}{\sqrt{R}} \leq \max_{\substack{ {\boldsymbol{\gamma}}  \in {{\mathscr W}}^\perp \\ \| {\boldsymbol{\gamma}}  \|_2=1}}  \bigg\|\sum_{k=1}^K  {\bar{\boldsymbol{\gamma}}} _k D_k ( {\mathbf{C}} )  {\vec{v}}  \bigg\|_2,
\end{equation}
from which the first inequality in equation \eqref{eq:JGEvaluebound} will follow. For $r=1, \dots, R$, let $ {\vec{c}} _r= (c_{1 r},\dots, c_{K r})$ be the $r$th column of $ {\mathbf{C}} $. Then $c_{k r}$ is the $r$th diagonal entry of $D_k ( {\mathbf{C}} )$ for each $k$, and, as shown by Theorem \ref{theorem:rankvsmult}, the set $\{\mathrm{span} ( {\vec{c}} _r)\}_{r=1}^R$ is the spectrum of ${\mathcal T}$. By way of equation \eqref{eq:ChordalEqs}, to  prove that equation \eqref{eq:JGEvalueIntermediateBound} holds it is sufficient to show that there is an index $r_*$ such that
\begin{equation}
\label{eq:svBoundGoal}
 \max_{\substack{ {\boldsymbol{\gamma}}  \in {{\mathscr W}}^\perp \\ \| {\boldsymbol{\gamma}}  \|_2=1}}\frac{ | \langle  {\vec{c}} _{r_*},  {\boldsymbol{\gamma}}  \rangle | }{\sqrt{R}} \leq \max_{\substack{ {\boldsymbol{\gamma}}  \in {{\mathscr W}}^\perp \\ \| {\boldsymbol{\gamma}}  \|_2=1}} \bigg\|\sum_{k=1}^K  {\bar{\boldsymbol{\gamma}}} _k D_k ( {\mathbf{C}} )  {\vec{v}}  \bigg\|_2 .
\end{equation}

To this end let $ {\mathbf{e}} _r$ be the $r$th standard basis vector in $ {\mathbb{K}} ^R$ and observe
\begin{equation}
\label{eq:BoundOnBasisVec}
\sum_{k=1}^K \bar{\gamma}_k D_k ( {\mathbf{C}} )  {\mathbf{e}} _r   =\sum_{k=1}^K \bar{\gamma}_k c_{kr}  {\mathbf{e}} _r =   \left(\sum_{k=1}^K \bar{\gamma}_k c_{kr}\right)  {\mathbf{e}} _r  = \langle   {\vec{c}} _r,  {\boldsymbol{\gamma}}  \rangle  {\mathbf{e}} _r.
\end{equation}
Decompose $ {\vec{v}} =\sum_{r=1}^R v_r  {\mathbf{e}} _r$ and note that $\| {\vec{v}} \|_2=1$ implies $\sum_{r=1}^R |v_r|^2=1$. Then for each fixed $ {\boldsymbol{\gamma}}  \in {{\mathscr W}}^\perp$ with $\| {\boldsymbol{\gamma}} \|_2=1$ we have
\[
\begin{array}{rcccc}

\|\sum_{k=1}^K \bar{\gamma}_k D_k ( {\mathbf{C}} )  {\vec{v}} \|_2 &=&  \left\| \sum_{k=1}^K \bar{\gamma}_k D_k ( {\mathbf{C}} ) \left(\sum_{r=1}^R v_r  {\mathbf{e}} _r \right) \right\|_2 & = & \left\| \sum_{r=1}^R v_r \left( \sum_{k=1}^K \bar{\gamma}_k D_k ( {\mathbf{C}} )  {\mathbf{e}} _r\right) \right\|_2 \\
&=& \|\sum_{r=1}^R v_r \langle  {\vec{c}} _r,  {\boldsymbol{\gamma}}  \rangle  {\mathbf{e}} _r \|_2 &=& \sqrt{\sum_{r=1}^R |v_r|^2 | \langle   {\vec{c}} _r,  {\boldsymbol{\gamma}}  \rangle|^2}.
\end{array}
\]

Now let $r_*$ be an index such that $|v_{r_*}|\geq |v_r|$ for all $r=1,\dots, R$. Since $\|v\|_2=1$, we must have $|v_{r_*}|\geq \frac{1}{\sqrt{R}}$. Thus we have
\[
\bigg\|\sum_{k=1}^K \bar{\gamma}_k D_k ( {\mathbf{C}} )  {\vec{v}} \bigg\|_2=\sqrt{\sum_{r=1}^R |v_r|^2 | \langle   {\vec{c}} _r,  {\boldsymbol{\gamma}}  \rangle|^2}\geq |v_{r_*} | |\langle  {\vec{c}} _{r_*}, {\boldsymbol{\gamma}}  \rangle | \geq \frac{|\langle  {\vec{c}} _{r_*},  {\boldsymbol{\gamma}}  \rangle |}{\sqrt{R}}.
\]
Noting that the above inequality holds for all $ {\boldsymbol{\gamma}}  \in {{\mathscr W}}^\perp$ with $\| {\boldsymbol{\gamma}} \|_2=1$ gives equation \eqref{eq:svBoundGoal} from which equation \eqref{eq:JGEvalueIntermediateBound} follows. The desired inequality then follows by combining inequalities \eqref{eq:JGEspNormBound} and \eqref{eq:JGEvalueIntermediateBound}.

To prove the second inequality in equation \eqref{eq:JGEvaluebound} observe that
\[
\spnorm{ {\mathcal E}  \cdot_1  {\mathbf{A}} ^{-1} \cdot_2  {\mathbf{B}} ^{-1} \cdot_3  {\mathbf{P}} _{{{\mathscr W}}^\perp}}
 \leq   \| {\mathbf{A}} ^{-1}\|_2 \| {\mathbf{B}} ^{-1}\|_2 \max_{\substack{ {\boldsymbol{\gamma}}  \in {{\mathscr W}}^\perp \\ \| {\boldsymbol{\gamma}}  \|_2=1 }} \left\|\sum _{k=1}^K \bar{\gamma}_k  {\mathbf{E}} _k \right\|_2    \leq  \frac{ \spnorm{{  {\mathcal E} }}}{\sigma_{\min} ( {\mathbf{A}} ) \sigma_{\min} ( {\mathbf{B}} )}
\]
from which the result follows.

It remains to show that the coefficient $\sqrt{R}$ is not necessary in the special case where ${\mathcal W}$ is slice mix invertible with $ {\mathbb{K}} $-rank $R$ and CPD $[\![ {\mathbf{A}} _{\mathcal W}, {\mathbf{B}} _{\mathcal W}, {\mathbf{C}} _{\mathcal W}]\!]$ where $ {\mathbf{A}} _{\mathcal W}= {\mathbf{A}} $ or $ {\mathbf{B}} _{\mathcal W}= {\mathbf{B}} $. It is straightforward to check that the spectral variation ${\mathrm{sv}} [{\mathcal T},{\mathcal W}]$ is invariant under a simultaneous permutation of the first and second mode of ${\mathcal T}$ and ${\mathcal W}$, so it is sufficient to consider the case $ {\mathbf{B}} _{\mathcal W}= {\mathbf{B}} $.

In this case, Theorem \ref{theorem:rankvsmult} shows that the JGE vectors of both ${\mathcal T}$ and ${\mathcal W}$ are given by the columns of $ {\mathbf{B}} ^{-\mathrm{T}}$, hence the vector $ {\vec{v}} $ appearing the proof above is simply equal to $ {\mathbf{e}} _r$ for some index $r \in 1,\dots,R$.  Using equation \eqref{eq:JGEspNormBound} together with equation \eqref{eq:BoundOnBasisVec} then shows that we have
\[
\chi({{\mathscr L}}_r,{{\mathscr W}})=\max_{\substack{ {\boldsymbol{\gamma}}  \in {{\mathscr W}}^\perp \\ \| {\boldsymbol{\gamma}}  \|_2=1}} |\langle  {\vec{c}} _r,  {\boldsymbol{\gamma}} \rangle| \leq  \Big\|  {\mathcal E}  \cdot_1  {\mathbf{A}} ^{-1} \cdot_2  {\mathbf{B}} ^{-1} \cdot_3  {\mathbf{P}} _{{{\mathscr W}}^\perp} \Big\|_{\mathrm{sp}}  \leq   \frac{ \spnorm{{  {\mathcal E} }}}{\sigma_{\min} ( {\mathbf{A}} ) \sigma_{\min} ( {\mathbf{B}} )},
\]
as claimed.
\end{proof}

\bibliographystyle{siamplain}

Numerical Analysis (math.NA); Algebraic Geometry (math.AG)