The Riemannian Barzilai-Borwein method with nonmonotone line search and the matrix geometric mean computation