Consider massless free scalar theory $$S = \int d^4x \sqrt{-g}L = \int d^4x \sqrt{-g} \;g^{cd}\nabla_c\phi \nabla_d \phi.\tag{1}$$ To compute the Hilbert stress-energy tensor we require $\sim\delta S/\delta g_{ab}$. I want to focus on $\delta L/ \delta g_{ab}$. $$\frac{\delta L}{\delta g_{ab}} = \frac{\delta g^{cd}}{\delta g_{ab}} \nabla_c\phi \nabla_d \phi = - g^{ac} g^{db}\nabla_c\phi \nabla_d \phi= -\nabla^a\phi \nabla^b \phi,\tag{2}$$ where I have used $\delta (g_{ab}g^{bc})=0$ to derive the implied expression for $\delta g^{ab}$ in terms of $\delta g_{ab}$. This is correct according to H. Reall's notes: (eq. 10.31) http://www.damtp.cam.ac.uk/user/hsr1000/part3_gr_lectures_2017.pdf .
However $$L=g^{cd}\nabla_c\phi \nabla_d \phi=g_{cd}\nabla^c\phi \nabla^d \phi\tag{3}$$ so I could also have $$\frac{\delta L}{\delta g_{ab}} = \frac{\delta g_{cd}}{\delta g_{ab}} \nabla^c\phi \nabla^d \phi.\tag{4}$$ For the two results to match, it would seem that we need $$\frac{\delta g_{cd}}{\delta g_{ab}} \sim - \frac{1}{2}(\delta^a_c \delta^b_d +\delta^a_d \delta^b_c )\;?\tag{5}$$ but I'm pretty sure there should be no minus sign!