“…First, it exhibits the important role of the reduction to the limiting case ( μ , μ c )=(1,0)—it is more than a routine standardisation; rather, it allows us, by removing the skew‐symmetric part (since μ c =0), to use the much simpler, negative‐curvature geometry and rich structure of the space of positive definite matrices with the Riemannian metric d . In connection with geodesics and exponentials, we should recall that the ordinary matrix exponential from GL( n ) to GL + ( n ) is not surjective and in fact not identical to the Riemannian exponential, , p1535 while for the positive definite matrices, surjectivity holds, as we already noted, and the two exponentials are the same. Second, the preceding computation, which is based on the fact that α I commutes with every other matrix for any α >0, suggests that the “degree of commutativity” of the factor U from the paragraph before , which is related to the sub‐matrices (or subspaces) on which commutativity does hold, should play a major role in the coordinate system that we wish to set up for problem .…”