Geometric means of quasi-Toeplitz matrices

We study means of geometric type of quasi-Toeplitz matrices, that are semi-infinite matrices A = (a(i,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 sigma(infinity)(k=-infinity) a(k)(eikt,) in the sense that (T(a))(i, j) = a(j-i) .If a is real valued and essentially bounded, then these matrices represent bounded self-adjoint operators on L-2. We prove that if a(1), ... , a (p) 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 a(1), ... , a(p), are quasi-Toeplitz matrices associated with the function (a(1) middot middot middot a(p)) (1/p) , 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.