Closed-form geodesics and trust-region method to calculate Riemannian logarithms on Stiefel and its quotient manifolds