MathOverflow is a question and answer site for professional mathematicians. Join them; it only takes a minute:

Sign up
Here's how it works:
  1. Anybody can ask a question
  2. Anybody can answer
  3. The best answers are voted up and rise to the top

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

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

share|cite|improve this question
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. – Wadim Zudilin Aug 1 '11 at 11:56
@Wadim: Are you sure with your calculation? $A\natural B$ must be Hermitian provided it is defined (i.e., $A,B$ are psd). – Sunni Aug 1 '11 at 13:25
@Wadim: Did you use the definition? $a\sharp b := a^{1/2}(a^{-1/2}ba^{-1/2})^{1/2}a^{1/2}$ – Suvrit Aug 1 '11 at 17:09
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?) – Federico Poloni Aug 1 '11 at 19:29
Hi Wadim:Inspired from this paper… Instead of entrywise definition, I would like to see a natural definition. – Russel Sep 12 '11 at 12:42
up vote 12 down vote accepted

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.$$

share|cite|improve this answer
I don't understand the "no loss of generality" here. Of course if $A,B$ commute, then we can reduce to the diagonal case. – Gerald Edgar Aug 1 '11 at 1:13
@Gerald: No commutability was assumed; I had first put a comment, but now I have edited it into the answer. – Suvrit Aug 1 '11 at 17:22
Sorry for the sloppiness, by "invariant" I mean that $S(A\sharp B)S^*=(SAS^*) \sharp (SBS^*)$. – Federico Poloni Aug 1 '11 at 19:32
@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. – Federico Poloni Aug 2 '11 at 7:34
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$. – Suvrit Aug 2 '11 at 18:19

Your Answer


By posting your answer, you agree to the privacy policy and terms of service.

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