  We study the convergence of the Riemannian steepest descent algorithm on the
Grassmann manifold for minimizing the block version of the Rayleigh quotient of
a symmetric matrix. Even though this problem is non-convex in the Euclidean
sense and only very locally convex in the Riemannian sense, we discover a
structure for this problem that is similar to geodesic strong convexity,
namely, weak-strong convexity. This allows us to apply similar arguments from
convex optimization when studying the convergence of the steepest descent
algorithm but with initialization conditions that do not depend on the eigengap
$\delta$. When $\delta>0$, we prove exponential convergence rates, while
otherwise the convergence is algebraic. Additionally, we prove that this
problem is geodesically convex in a neighbourhood of the global minimizer of
radius $\mathcal{O}(\sqrt{\delta})$.
