Logm function

Riemannian HPD logarithmic map