  In this work, we are presenting an efficient way to compute the geometric
mean of two positive definite matrices times a vector. For this purpose, we are
inspecting the application of methods based on Krylov spaces to compute the
square root of a matrix. These methods, using only matrix-vector products, are
capable of producing a good approximation of the result with a small
computational cost.
