I have n-by-n correlation matrix $C_{ij}$, it's symmetric positive definite, and $C_{ii}=1$ and $C_{ij}\le1$ etc.

I want to generate a bunch of matrices $M_{ij}$ which are on the same distance from this initial matrix. The distance measure is Riemannian, like in [this answer][1]:
$$d(C,M)=\sqrt{\sum_{k=1}^n \log^2 \lambda_k}$$
where $\lambda_k$ are eigenvalues of $CM^{-1}$

I'm looking for an algorithm that can find me the matrices $M$.


  [1]: https://mathoverflow.net/questions/108280/left-invariant-metric-on-sln/109951#109951