Mid: Geodesic midpoint between HPD matrices

In JorisChau/pdSpecEst: An Analysis Toolbox for Hermitian Positive Definite Matrices

Description

Mid calculates the geodesic midpoint between two HPD matrices under the affine-invariant Riemannian metric as in \insertCiteB09pdSpecEst[Chapter 6].

Usage

Mid(A, B)

Arguments

A, B	Hermitian positive definite matrices (of equal dimension).

References

\insertAllCited

See Also

pdMean

Examples

 ## Generate two random HPD matrices
 a <- matrix(complex(real = rnorm(9), imaginary = rnorm(9)), nrow = 3)
 A <- t(Conj(a)) %*% a
 b <- matrix(complex(real = rnorm(9), imaginary = rnorm(9)), nrow = 3)
 B <- t(Conj(b)) %*% b
 ## Compute midpoint
 Mid(A, B)
 ## Midpoint coincides with two-point intrinsic Karcher mean
 all.equal(pdMean(array(c(A, B), dim = c(3, 3, 2))), Mid(A, B))

JorisChau/pdSpecEst documentation built on Jan. 13, 2020, 6:08 p.m.