However, when the condition number \( \kappa \)  is large, the convergence factor \( \rho \)  could be very close to 1.
Consequently, the convergence rate could be painfully slow.
To accelerate it, [1]} proposed the heavy-ball method, which is also called the momentum method [2]}.
Specifically, in this work, we investigate the gradient descent with the classical momentum, which updates the estimates as
\(\widehat{\theta }^{(t)} = \widehat{\theta }^{(t-1)} - v^{(t)},\\v^{(t)} = \gamma v^{(t-1)} + \alpha \dot{\mathcal {L}}(\widehat{\theta }^{(t-1)}),\) 
