Un corto de prueba puede ser determinado utilizando la norma de Frobenius: Como ya se señaló en la pregunta que está conectado, tenemos
$$
\mathrm{TD}(a,b,c) = a^T \cdot \mathrm{diag}(b_1, \dotsc, b_n) \cdot c
$$
donde $\mathrm{diag}(b_1, \dotsc, b_n)$ denota la diagonal de la matriz diagonal con entradas de $b_1, \dotsc, b_n$. El uso de la submultiplicativity de la norma de Frobenius nos encontramos con que
\begin{align*}
\mathrm{TD}(a,b,c)
&\leq |\mathrm{TD}(a,b,c)|
= \|\mathrm{TD}(a,b,c)\|_F
= \|a^T \cdot \mathrm{diag}(b_1, \dotsc, b_n) \cdot c\|_F \\
&\leq \|a^T\|_F \cdot \|\mathrm{diag}(b_1, \dotsc, b_n)\|_F \cdot \|c\|_F
= \|a\| \cdot \|b\| \cdot \|c\|,
\end{align*}
donde $\|\cdot\|_F$ denota la norma de Frobenius.
Otra prueba puede ser dada por el uso de la de Cauchy-Schwarz desigualdad: podemos suponer w.l.o.g. que $a_i, b_i, c_i \geq 0$ por cada $1 \leq i \leq n$. Debido a $b_i \geq 0$ todos los $1 \leq i \leq n$ nos encontramos con que la forma bilineal $\langle \cdot, \cdot \rangle_b$ definidos a través de
$$
\langle x,y \rangle_b
= x^T \cdot \mathrm{diag}(b_1, \dotsc, b_n) \cdot y
= \mathrm{TD}(x,b,y)
$$
es simétrica y positiva semidefinite con
$$
\|x\|_b
= \sqrt{\langle x,x \rangle_b}
= \sqrt{ \sum_{i=1}^n b_i x_i^2 }.
$$
Observe que para todos los $1 \leq i,k \leq n$ tenemos $b_i b_k \leq \sum_{j=1}^n b_j^2$: Si $i = k$ esto es claro y si $i \neq k$
$$
b_i b_k \leq 2 b_i b_k \leq b_i^2 + b_k^2 \leq \sum_{j=1}^n b_j^2.
$$
Así que a partir de Cauchy-Schwarz desigualdad se sigue que
\begin{align*}
\mathrm{TD}(a,b,c)
&= \langle a, c \rangle_b
\leq \|a\|_b \cdot \|c\|_b
= \sqrt{ \sum_{i=1}^n b_i a_i^2} \cdot \sqrt{\sum_{k=1}^n b_k c_k^2 } \\
&= \sqrt{ \sum_{i,k=1}^n a_i^2 b_i b_k c_k^2}
\leq \sqrt{ \sum_{i,j,k=1}^n a_i^2 b_j^2 c_k^2 } \\
&= \sqrt{\sum_{i=1}^n a_i^2} \cdot \sqrt{\sum_{j=1}^n b_j^2}
\cdot \sqrt{\sum_{k=1}^n c_k^2}
= \|a\| \cdot \|b\| \cdot \|c\|.
\end{align*}
Es importante hacer notar que el uso de la aproximación a través de la norma de Frobenius también podemos directamente generalizar nuestros resultados arbitrarios $x^1, \dotsc, x^m \in \mathbb{C}^n$, en el sentido de que
\begin{align*}
\left|\sum_{i=1}^n x^1_i \dotsm x^m_i\right|
&= \left\|
(x^1)^T
\cdot \mathrm{diag}(x^2_1, \dotsc, x^2_n)
\dotsm \mathrm{diag}(x^{m-1}_1, \dotsc, x^{m-1}_n)
\cdot x^m
\right\|_F \\
&\leq
\|(x^1)^T\|_F
\cdot \|\mathrm{diag}(x^2_1, \dotsc, x^2_n)\|_F
\dotsm \|\mathrm{diag}(x^{m-1}_1, \dotsc, x^{m-1}_n)\|_F
\cdot \|x^m\|_F \\
&= \|x^1\| \dotsm \|x^m\|.
\end{align*}