We present and analyze an abstract step size selection algorithm which ensures asymptotic stability of numerical approximations to asymptotically stable ODEs. A particular implementation of this algorithm is proposed and tested with two numerical examples. The application to ODEs solving nonlinear optimization problems on manifolds is explained and illustrated by means of the Rayleigh flow for computing eigenvalues of symmetric matrices.