Let $p<n$ and
-$H\in\mathbb{R}^{n\times n}$ be symmetric with eigendecompoistion being $H=U\Lambda U^{\text{T}}$,
-$A\in\mathbb{R}^{n\times p}$,
-$D\in\mathbb{R}^{p\times p}$ be a diagonal matrix.
We have $H=U\Lambda U^{\text{T}}=A D A^{\text{T}}$. I want to find $A$ and $D$ as functions of $U$ and $\Lambda$.
For the above equation to be true only $p$ eigenvalues of $H$ in $\Lambda$ are nonzero (RHS is of rank $p$). Forming $A$ by the corresponding $p$ eigenvectors (columns of $U$) and $D$ with those nonezero eigenvalues is a solution.
My guess is that is the only solution but don't know how to prove it.