Let $A,B$ be positive definite (Hermitian) matrices. Define the Arithmetic-geometric means of positive matrices by $A_0=A, G_0=B$, $A_{n+1}=\frac{A_n+G_n}{2}, G_{n+1}=A_n\natural G_n$, where $A_n\natural G_n$ means the geometric mean defined in http://www.isid.ac.in/~statmath/eprints/2011/isid201102.pdf

Will $\{A_n\}$ and $\{G_n\}$ converge to the same matrix?

  • $\begingroup$ I computed the geometric mean of $A=[1,-1/2;-1/2,1]$ and $B=[1,1/2;1/2,1]$ (both matrices are symmetric positive definite) using Eq.(1) from the link. The result is a matrix whose diagonal entries are not real (in particular, it isn't Hermitian). You have to find a proper set of matrices (or a better geometric mean) to have your question well defined. $\endgroup$ – Wadim Zudilin Aug 1 '11 at 11:56
  • $\begingroup$ @Wadim: Are you sure with your calculation? $A\natural B$ must be Hermitian provided it is defined (i.e., $A,B$ are psd). $\endgroup$ – Sunni Aug 1 '11 at 13:25
  • $\begingroup$ @Wadim: Did you use the definition? $a\sharp b := a^{1/2}(a^{-1/2}ba^{-1/2})^{1/2}a^{1/2}$ $\endgroup$ – Suvrit Aug 1 '11 at 17:09
  • $\begingroup$ Sounds fishy, since that matrix is clearly Hermitian from its definition. Is the loss of Hermitianity small? It could be due to numerical errors and/or a poor sqrtm() implementation. (Incidentally, you're not using Matlab's sqrt() by chance, are you?) $\endgroup$ – Federico Poloni Aug 1 '11 at 19:29
  • 1
    $\begingroup$ Hi Wadim:Inspired from this paper lab.rockefeller.edu/cohenje/PDFs/… Instead of entrywise definition, I would like to see a natural definition. $\endgroup$ – Russel Sep 12 '11 at 12:42

Recall that since A and B are (hermitian) positive definite, we can without loss of generality (see below for proof) assume that $A=I$ and $B=D$, where $D$ is some positive diagonal matrix. With this observation, merely recall the convergence theory for the scalar case to conclude that the sequences $\{A_n\}$ and $\{B_n\}$ converge to the arithmetic-geometric mean of $I$ and $D$.

Note: (Added to improve clarity)

For positive definite $A$ and $B$, let $A=Q\Lambda Q^T$, $S=\Lambda^{-1/2}Q$, and let $U$ diagonalize $S^TQ^TBQS$ to $D$. Then, with $P=QSU$, we have $$P^TAP = U^TS^TQ^TQ\Lambda Q^TQSU = U^TQ^T\Lambda^{-1/2}\Lambda \Lambda^{-1/2}QU = I,$$ and by construction, $$P^TBP = U^TS^TQ^TBQSU=D.$$

| cite | improve this answer | |
  • 1
    $\begingroup$ I don't understand the "no loss of generality" here. Of course if $A,B$ commute, then we can reduce to the diagonal case. $\endgroup$ – Gerald Edgar Aug 1 '11 at 1:13
  • $\begingroup$ @Gerald: No commutability was assumed; I had first put a comment, but now I have edited it into the answer. $\endgroup$ – Suvrit Aug 1 '11 at 17:22
  • 1
    $\begingroup$ Sorry for the sloppiness, by "invariant" I mean that $S(A\sharp B)S^*=(SAS^*) \sharp (SBS^*)$. $\endgroup$ – Federico Poloni Aug 1 '11 at 19:32
  • 2
    $\begingroup$ @Gerald: this result is quite well-known; it is easier to prove using other characterizations of the geometric mean. For instance, using the fact that $Af(BA)=f(AB)A$ for all holomorphic $f$ (and in particular for the square root function), you can rewrite the definition as $A \sharp B = A(A^{-1}B)^{1/2}$, and from this form congruence invariance follows easily. $\endgroup$ – Federico Poloni Aug 2 '11 at 7:34
  • 2
    $\begingroup$ Maybe the easiest way to prove this is to note that the geometric mean is the unique positive definite solution to the Riccati equation: $XA^{-1}X=B$. $\endgroup$ – Suvrit Aug 2 '11 at 18:19

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Not the answer you're looking for? Browse other questions tagged or ask your own question.