All matrices are $n \times n$. $C$ is real symmetric positive definite.
How to solve $C = X^\top C X$ for $X$?
I am interested in characterizing both the set of real matrices satisfying the equation and the (possibly larger) set of complex solutions.