I think you should have $\frac{1}{2}$ in front of $t^THt$. But the reason is that the gradient of the function $f(t) = t^THt - t^Tv$
$$\nabla_t(t^THt - t^Tv) = (H+H^T)t-v = 2Ht-v$$
This is by linearity and for the bi-linear part see for example here.
For a critical point we get $\frac{1}{2}H^{-1}v$. This is a local minimum because the Hessian is actually $H$ and it's positive definite. To see that it's a global minimum notice that for any direction $t$ (a unit vector) the function $f(st), s\in{[0, \infty)}$ tends to infinity. This is because
$$f(st) = s^2 t^THt - s t^Tv = s(s t^THt - t^Tv) > s,$$
when $s>\frac{ t^Tv}{t^THt}$. By C-S we have $t^Tv \leq |t||v| = |v|$ and $|t^THt|$ is bounded by the largest eigenvalue of $H$, so outside a big enough ball $f$ will be larger than in the critical point.
I don't have the access to the original paper so I won't find out who made the initial mistake :P
– Jack Yoon Jul 12 '22 at 11:00