Dada una matriz diagonal $D$ con elementos diagonales dados por el vector $\mathbf{d}$ . Representando esto en notación de Einstein se obtiene
$$ D_{ij} = \delta_{ijk} d_k $$
donde
$$ \delta_{ijk} = \begin{cases} 1 & \text{if } i = j = k \\ 0 & \text{ otherwise} \end{cases} $$
Si ahora aplico esto en una multiplicación de matrices, por ejemplo
$$ (AD)_{lj} = A_{li} D_{ij} = A_{li} \delta_{ijk} d_k = A_{li} d_i $$ o $$ (ADA)_{ml} = A_{mi} D_{ij} A_{jl} = A_{mi} \delta_{ijk} d_k A_{jl} = A_{mi} d_i A_{il} $$
El primer ejemplo tiene sentido desde el punto de vista de la entrada, pero sólo si no hay suma sobre $i$ . Además, los índices del LHS y del RHS ya no coinciden. Lo mismo ocurre con el segundo ejemplo, pero sólo si se suman los tres $i$ 's.
Esto obviamente viola la notación de Einstein, pero no veo en qué paso se hace una suposición falsa. Por lo tanto, mis preguntas son:
- ¿Dónde me equivoco en mi razonamiento?
- ¿Hay alguna otra manera de explotar el hecho de que $D$ es diagonal (en notación de índices)?