10
$\begingroup$

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?

$\endgroup$
7
  • $\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$ Aug 1, 2011 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, 2011 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, 2011 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$ Aug 1, 2011 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, 2011 at 12:42

1 Answer 1

13
$\begingroup$

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

$\endgroup$
10
  • 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$ Aug 1, 2011 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, 2011 at 17:22
  • 1
    $\begingroup$ Sorry for the sloppiness, by "invariant" I mean that $S(A\sharp B)S^*=(SAS^*) \sharp (SBS^*)$. $\endgroup$ Aug 1, 2011 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$ Aug 2, 2011 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, 2011 at 18:19

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge you have read our privacy policy.

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