  Quantum Optimal Control (QOC) is the field devoted to the production of
external control protocols that actively guide quantum dynamics. Solutions to
QOC problems were shown to constitute continuous submanifolds of control space.
A solution navigation method exploiting this property to achieve secondary
features in the control protocols was proposed [Larocca \textit{et al},
arXiv:1911.07105]. Originally, the technique involved the computation of the
exact Hessian matrix. In this paper, we show that the navigation can be
alternatively performed with a finite-difference approximation scheme, thus
enabling the application of this procedure in systems where computing the exact
Hessian is out of reach.
