Overview
{% \vec{r} = \textbf{D} \textbf{H} + \vec{\epsilon} %}
where {% \textbf{Var}(\textbf{H}) = \mathbb{I} %}
Diagonalization
Diagonalization