Geometric means of quasi-Toeplitz matrices

We study means of geometric type of quasi-Toeplitz matrices, that are semi-infinite matrices A=(ai,j)i,j=1,2,… of the form A= T(a) + E, where E represents a compact operator, and T(a) is a semi-infinite Toeplitz matrix associated with the function a, with Fourier series ∑k=-∞∞akeikt, in the sense that (T(a)) i,j= aj-i. If a is real valued and essentially bounded, then these matrices represent bounded self-adjoint operators on ℓ2. We prove that if a1, … , ap are continuous and positive functions, or are in the Wiener algebra with some further conditions, then matrix geometric means, such as the ALM, the NBMP and the weighted mean of quasi-Toeplitz positive definite matrices associated with a1, … , ap, are quasi-Toeplitz matrices associated with the function (a1⋯ap)1p, which differ only by the compact correction. We introduce numerical algorithms for their computation and show by numerical tests that these operator means can be effectively approximated numerically.