Solution (source code)

= Solution

Put $h_{ab}=\delta g_{ab}$. Varying $g^{ac}g_{cb}=\delta^a_b$ gives
$$
\boxed{\delta g^{ab}=-g^{ac}g^{bd}h_{cd}=-h^{ab}}.
$$
The determinant identity $\delta\log|\det g|=g^{ab}h_{ab}$ gives
$$
\boxed{\delta(d\operatorname{vol}_g)
=\frac12g^{ab}h_{ab}\,d\operatorname{vol}_g}.
$$
Finally, varying the formula for the <Christoffel symbols> and rewriting partial derivatives covariantly gives the tensor
$$
\boxed{\delta\Gamma^a{}_{bc}
=\frac12g^{ad}(\nabla_bh_{cd}+\nabla_ch_{bd}-\nabla_dh_{bc})}.
$$
These are the basic <metric variation> identities.