Metric on an open subset of $\mathbb{R}^d$ and Christoffel symbol of the second kind I need to prove the following identity. $\Omega \subset \mathbb{R}^d$ is open and $g$ is a metric field on $\Omega$. Further $\Gamma_{\,kl}^j$ denotes the Christoffel symbol of second kind. $g^{ij}$ is the inverse matrix to $g_{ij}$ and $g$ denotes the determinant.
$$
 \Gamma_{\,jl}^j = \frac{\partial}{\partial x^l} \log \sqrt{\vert g \vert}
$$
Here is what I have so far:
$$
\Gamma^j_{\;jl} = \frac{1}{2} g^{jk} \left( \frac{\partial }{ \partial x^l} g_{kj} + \frac{\partial}{ \partial x^j} g_{kl} - \frac{\partial}{ \partial x^k} g_{jl} \right) = \frac{1}{2} g^{jk}  \frac{\partial }{ \partial x^l} g_{kj}
$$
The second and the third term cancel each other out. Now as stated here (it's about the sixth equation, after the line starting with "The contracting relations..."
$$ \frac{1}{2} g^{jk}  \frac{\partial }{ \partial x^l} g_{kj} =
\frac{1}{2g} \frac{\partial}{ \partial x^l} g = \frac{\partial}{ \partial x^l} \log \sqrt{\vert g \vert}
$$
I don't get the second last step. I suspect it has got something to do with the Laplace formula and the Cramer Rule. However i seem not to be able to figure it out, the derivation is bugging me.
Any hints are appreciated!
 A: First the slick method: write $\bar g_{ik}$ for the value of $g_{ik}$ at the specific $x^l = \bar x^l$ (so $\bar g_{ik}$ is constant).  Let $x^l$ be a function of $\epsilon$ that is a perturbation of $\bar x^l$, so that $x^l = \bar x^l$ at $\epsilon=0$.
Let $$h_{ij} = \bar g^{ik} g_{kj} . \tag1$$  Then $h_{ij}$ is a perturbation of the identity: $$h_{ij} = \delta_{ij} + \epsilon \frac{\partial h_{ij}}{\partial \epsilon} + o(\epsilon).$$  Then it is easily computed that
$$ \frac{\partial h}{\partial \epsilon} = \frac{\partial h_{ii}}{\partial \epsilon} + o(\epsilon)$$
Differentiate $(1)$ with respect to $\epsilon$
$$ \frac{\partial h_{ii}}{\partial\epsilon} = \bar g^{ik} \frac{\partial g_{ki}}{\partial \epsilon} $$
Also $h = \bar g^{-1} g$.  Hence
$$ \bar g^{-1} \frac{\partial g}{\partial \epsilon} = \bar g^{ik} \frac{\partial g_{ki}}{\partial \epsilon} + o(\epsilon). $$
In particular, at $\epsilon = 0$ we can remove the overlines, and we get:
$$ g^{-1} \frac{\partial g}{\partial \epsilon} = g^{ik} \frac{\partial g_{ki}}{\partial \epsilon}. $$
Next the direct way: so $$g = \sum_\pi\sigma(\pi) \prod_i g_{i\pi(i)}.$$  Here the sum is over permutations $\pi$ of $\{1,\dots,n\}$.  Then
$$\frac{\partial g}{\partial x^l} = \sum_\pi \sigma(\pi) \sum_i \frac{\partial g_{i\pi(i)}}{\partial x^l} \prod_{j \ne i} g_{j\pi(j)} \\
=\sum_{ik} \frac{\partial g_{ik}}{\partial x^l}\sum_{\pi:\pi(i)=k} \prod_{j\ne i}\sigma(\pi) g_{j\pi(j)} \\
=\sum_{ik} \frac{\partial g_{ik}}{\partial x^l} g \ g^{ki}. $$
