The cone ${\bf SPD}\_n({\mathbb R})$ of symmetric positive definite matrices is endowed with a nice geometrical structure. The midpoint of the (unique) geodesic between $A$ and $B$ is the so-called *geometric mean*
$$A\sharp B=A^{1/2}(A^{-1/2}BA^{-1/2})^{1/2}A^{1/2}.$$



> 
> What is the associated distance ?
> 
> 
> 



metric on ${\bf SPD}_n({\mathbb R})$
If I'm not mistaken, this midpoint comes from the structure of Riemannian symmetric space on the space of symmetric positive definite matrices. The Riemannian metric is defined by $g(S)(h\_1,h\_2)=trace(S^{-1}h\_2S^{-1}h\_2)$, where $S$ is in $SDP\_n$ and $h\_1,h\_2$ are symmetric matrices (tangent vectors at $S$).


This metric has a lot of nice properties: 


* $Gl(n)$ acts transitively and isometrically on $SDP\_n$ by $g\cdot S=gSg^{T}$, making it a homogeneous Riemannian manifold. (In particular it is complete.)
* $S\mapsto S^{-1}$ is an isometric involution whose differential at $Id\in SDP\_n$ is minus the identity. This makes $SDP\_n$ a Riemannian symmetric space of non compact type (in particular it has nonpositive curvature). In the classification of Riemannian symmetric spaces it is often described as the quotient $GL(n)/O(n)$.
* This metric also arises as the Fisher information metric on the space of centered Gaussian random vectors on $\mathbb{R}^n$ (which are characterized by their covariance matrix). This makes it useful in some applications, as a mean of treating covariance matrices in a basis independent manner.


All these properties makes it a very nice object, and it turns out that most of the objects of Riemannian geometry (connection, geodesics, curvature) can be described explicitly in terms of matrix computations. The distance is given by: $d(A,B)=\|\log(AB^{-1})\|\_F$ , where $\|\ \|\_F$ is the Frobenius norm. I found the formula in this article: <http://arxiv.org/pdf/0807.4462v2.pdf>.


