Mid function

Geodesic midpoint between HPD matrices