Sea $\mathbf A$ ser un $n\times n$ matriz, y $\mathbf B$ ser un $nr\times nm$ matriz diagonal de bloques, siendo cada bloque una matriz de tamaño $r\times m$ . Además $\otimes$ denotan el producto de Kronecker.
¿Cómo puedo realizar la siguiente operación?
$$ \mathbf B (\mathbf A \otimes \mathbf I_m ) \mathbf B ^\top?$$
Me gustaría reescribir este término en el siguiente formato
$$ \sum_{i=1}^n\sum_{j=1}^n F(\mathbf A)G(\mathbf B_i \mathbf B_j^\top),$$
donde $\mathbf B_i$ corresponde al $i$ -en la diagonal de $\mathbf B$ y $F$ y $G$ son operadores lineales.