Given $(M,g)$ Riemannian manifold, $f \in \mathcal{C}^{\infty}(M)$, I want to prove that the following holds (denoting by $D\,df $ the Hessian of the function $f$ and $\text{grad} f$ its gradient):
$$2D\,df (\text{grad} f,\text{grad} f) = \langle \text{grad} \,|\text{grad} f|^2, \text{grad} f\rangle$$
I tried doing this:
Start from $$D\,df (\text{grad} f,\text{grad} f) = D_{\text{grad} f}D_{\text{grad} f}(f)$$ Then since we can write in coordinates $\text{grad} f = g^{ij}\frac{\partial f}{\partial x^i} \, \frac{\partial f}{\partial x^j}$ then we have:
$$D_{\text{grad} f} \, (g^{ij}\frac{\partial f}{\partial x^i} D_{\frac{\partial }{\partial x^j}}f)$$
Now, since $f$ is a function, the covaiant derivative of $f$ along $\frac{\partial }{\partial x^j} $ is $ \frac{\partial f}{\partial x^j}$ so we get:
$$D_{g^{ab}\frac{\partial f}{\partial x^a}\frac{\partial f}{\partial x^b}}g^{ij}\frac{\partial f}{\partial x^i}\frac{\partial f}{\partial x^j}$$ getting $$g^{ab}\frac{\partial f}{\partial x^a}\left[\frac{\partial}{\partial x^b}\left(g^{ij}\frac{\partial f}{\partial x^i}\frac{\partial f}{\partial x^j}\right)\right]$$
In an analogous fashion, if I develop in coordinates this term $\langle \text{grad} f \,|\text{grad} f|^2, \text{grad} f\rangle$ I get: $$\langle \text{grad} \left( g^{i \alpha}\frac{\partial f}{\partial x^i}\frac{\partial f}{\partial x^\alpha} \right), \text{grad} f\rangle$$ and tracing with the metric we end up with:
$$g^{aA}\frac{\partial f}{\partial x^A}\left[\frac{\partial}{\partial x^a}\left(g^{i\alpha}\frac{\partial f}{\partial x^i}\frac{\partial f}{\partial x^\alpha}\right)\right]$$
First I ask if I am doing the computations correctly, and then why I do not get a $2$, I think I misssed something. Thanks!