Square-root barycenter of covariance operators (source code)

= Square-root barycenter of covariance operators

For covariance operators $C_1,\ldots,C_n$, their Fréchet mean under the square-root distance is
$$
\overline C_R=\left(\frac1n\sum_{i=1}^nC_i^{1/2}\right)^2.
$$
This follows because taking positive square roots embeds the operators into the Hilbert space of self-adjoint Hilbert-Schmidt operators.