Md. Asif Uddin
    Problem I.7.B03

    A problem no learning rate solves

    counterexample▲▲▲

    Six decimal places on multipliers; step counts exact.

    STATEMENT

    Exhibit a quadratic on which every learning rate fails — not one badly chosen rate, but the whole stable range — and identify what is responsible. Then show that remembering the previous step fixes it, and by how much.

    GIVEN

    H=diag(100, 0.01)\mat{H} = \mathrm{diag}(100,\ 0.01), so κ=104\kappa = 10^4. Start at θ0=(1,1)\vec{\theta}_0 = (1, 1); the target is an error norm below 10−310^{-3}.

    FIND

    The stability limit; the step count for every η\eta across the stable range; the count at the optimal η∗\eta^*; and the same for heavy-ball momentum at its own optimum, with the speed-up.

    STRATEGY

    The two directions decouple, as in I.7.B02, so the difficulty can be stated before anything is run: the steep direction demands η<0.02\eta < 0.02, and at any such η\eta the flat direction’s multiplier is 1−0.01η≥0.99981 - 0.01\eta \ge 0.9998. One rate cannot serve both. Sweeping confirms there is no exception rather than merely no obvious one.

    SOLUTION

    Step 1 — the bind, stated in advance. Stability requires η<2/λmax⁡=0.02\eta < 2/\lambda_{\max} = 0.02. At the largest permissible rate the flat direction is multiplied by

    1−ηλmin⁡>1−0.02(0.01)=0.99981 - \eta\lambda_{\min} > 1 - 0.02(0.01) = 0.9998

    so it loses at most two parts in ten thousand per step, however η\eta is chosen. That is the counterexample in one line; the sweep is the demonstration.

    Step 2 — the sweep.

    η\etasteep multiplierflat multipliersteps to 10−310^{-3}
    0.0200000.0200001.0000001.0000000.9998000.999800diverges
    0.0100000.0100000.0000000.0000000.9999000.99990069 07569\,075
    0.0198000.0198000.9800000.9800000.9998020.99980234 88534\,885
    0.0198010.0198010.9801000.9801000.9998020.99980234 88334\,883
    0.0199000.0199000.9900000.9900000.9998010.99980134 70934\,709
    0.0200100.0200101.0010001.0010000.9998000.999800diverges
    0.0500000.0500004.0000004.0000000.9995000.999500diverges

    Nothing in the stable range does better than about 34 70034\,700 steps for a thousandfold reduction on a two-dimensional problem with an exactly known Hessian.

    Step 3 — the most instructive row. At η=0.01\eta = 0.01 the steep multiplier is exactly zero: 1−0.01(100)=01 - 0.01(100) = 0. That direction is solved perfectly, in one step, and the run still takes 69 07569\,075 steps — twice as many as the rates that solve it imperfectly. Perfect handling of the direction that constrains stability buys nothing, because the cost is entirely in the other direction. This is why “tune the learning rate” is not a complete answer to slow training.

    Step 4 — the optimum, and its transient. η∗=0.0199980\eta^* = 0.0199980 gives ρ=0.99980002\rho = 0.99980002 and 36 27236\,272 steps — slightly worse than η=0.0199\eta = 0.0199 at 34 70934\,709. As in I.7.B01, η∗\eta^* minimises the asymptotic rate while the run is measured over a finite budget. The gap is small here because κ\kappa is large enough that the asymptote arrives early.

    Step 5 — momentum. Heavy ball keeps a running average of gradients, so the recurrence is second-order and its rate depends on κ\sqrt{\kappa}. At the classical optimum

    β=(κ−1κ+1)2=0.96078816,η=4(λmax⁡+λmin⁡)2=0.03921184\beta = \left(\frac{\sqrt{\kappa} - 1}{\sqrt{\kappa} + 1}\right)^2 = 0.96078816, \qquad \eta = \frac{4}{(\sqrt{\lambda_{\max}} + \sqrt{\lambda_{\min}})^2} = 0.03921184

    with ρ=(κ−1)/(κ+1)=0.98019802\rho = (\sqrt{\kappa} - 1)/(\sqrt{\kappa} + 1) = 0.98019802 and 708 steps, against 36 27236\,272: a speed-up of 51.2×51.2\times.

    Step 6 — why it works. Momentum does not make the steep direction stable at a larger rate; note that η=0.0392\eta = 0.0392 is above the plain-descent threshold of 0.020.02. It works because the steep direction’s gradient alternates in sign when overshooting, so consecutive contributions to the moving average partly cancel, while the flat direction’s gradient keeps the same sign and accumulates. The average suppresses oscillation and amplifies drift, which is exactly the asymmetry the problem needs.

    Answer

    Stability requires η<0.02\eta < 0.02, and every rate in that range needs at least about 34 70034\,700 steps; η∗=0.0199980\eta^* = 0.0199980 takes 36 27236\,272. The obstruction is κ=104\kappa = 10^4 itself, not the choice of rate.

    Heavy-ball momentum at β=0.96078816\beta = 0.96078816, η=0.03921184\eta = 0.03921184 reaches the same target in 708 steps, a 51.2×51.2\times speed-up.

    Check — sanity

    The measured momentum rate matches the theory. ρ=0.98019802=99/101=(κ−1)/(κ+1)\rho = 0.98019802 = 99/101 = (\sqrt{\kappa} - 1)/(\sqrt{\kappa} + 1) for κ=100\sqrt{\kappa} = 100, exactly. And it equals the plain-descent rate for a problem of condition number 100 — momentum turns κ=104\kappa = 10^4 into an effective 10210^2, which is what “the square root” means concretely.

    The speed-up is 51×51\times, not the 100×100\times the rates suggest. Both methods are measured to a fixed target from a fixed start, so both pay a transient the asymptotic rates ignore, and momentum’s is proportionally larger because it needs several steps before its average carries information. Quoting κ\sqrt{\kappa} as the speed-up would be quoting a limit as a measurement.

    Divergence appears immediately above 0.020.02. η=0.02001\eta = 0.02001 gives a steep multiplier of 1.0011.001 and diverges; η=0.0199\eta = 0.0199 gives 0.990.99 and converges. Five parts in ten thousand separate them, which is the sharp edge of I.7.B01 in a stiffer problem.

    Where this breaks

    The momentum parameters here were computed from λmax⁡\lambda_{\max} and λmin⁡\lambda_{\min}, which nobody has for a real network. In practice β\beta is set to 0.90.9 and left there, which is optimal for κ≈361\kappa \approx 361 and merely serviceable elsewhere — the tuned value above is 0.960.96, and using 0.90.9 on this problem gives up much of the advantage. The theory says what is available; it does not hand over the setting.

    Momentum also does not help with the part of the difficulty that is noise. Averaging gradients reduces variance, but the accelerated rate of I.7.T2 assumes exact gradients, and with minibatch noise the improvement degrades towards the noise floor. On a stochastic problem momentum is worth having and the factor is smaller than κ\sqrt{\kappa}.

    Finally, this is a fixed quadratic and a neural network is not. Its curvature changes as training proceeds, so the tuned (η,β)(\eta, \beta) for step 1 is wrong by step 10 000 — which is one honest reading of why schedules exist, and why adaptive methods that estimate curvature per coordinate as they go, rather than being told it once, are the ones that ended up in use.

    Variation

    Set β=0.9\beta = 0.9 — the default nobody tunes — and keep η\eta at its value above. Predict whether the run gets faster or slower before measuring, then find the κ\kappa for which β=0.9\beta = 0.9 is optimal. That number explains where the default came from.

    Draws on