Put on $\mathbb{R}^2$ the standard riemannian metric $g=dx^2+dy^2$ and consider a differentiable function $f:\mathbb{R}^2\rightarrow \mathbb{R}$. I'll write $\nabla f_p$ for the gradient of $f$ computed in $p\in \mathbb{R}^2$. I want to understand if $g(\nabla f_p,\nabla f_p)$ is independent of the coordinates (the writing $g(\nabla f_p,\nabla f_p)$ seems independent of the coordinates. I'll explain myself better:
$g(\nabla f_p,\nabla f_p)$ is clearly equal to $(\partial_xf_p)^2+(\partial_y f_p)^2$. Now, consider $V$, $W$ vector fields which are orthonormal in every point of a neightborhood $U$ of $0\in \mathbb{R}^2$ with respect to $g$. Is it true $(\partial_xf_p)^2+(\partial_y f_p)^2=(Vf_p)^2+(Wf_p)^2$?
Intuitively this seems true, but I don't know how to do the right calculations:
suppose $V=v_x\partial_x+v_y\partial_y$ and $W=w_x\partial_x+w_y\partial_y$, with $g(V,V)=v_x^2+v_y^2=1$,$g(W,W)=w_x^2+w_y^2=1$ and $g(V,W)=v_xw_x+v_yw_y=0$.
Then $(Vf_p)^2+(Wf_p)^2=(\partial_xf_p)^2(v_x^2+w_x^2)+(\partial_yf_p)^2(v_y^2+w_y^2)+2\partial_xf\partial_yf(v_xv_y+w_xw_y)$ which doesn't seem equal to $(\partial_xf)^2+(\partial_yf)^2$... What am I doing wrong?
What you're forgetting is that if the columns of a matrix are orthonormal, then so are the rows! So there's nothing wrong with your calculation. :)