Flat FLRW Einstein-tensor divergence (source code)

= Flat FLRW Einstein-tensor divergence

For $ds^2=-dt^2+a(t)^2d\mathbf x^2$, the nonzero connection coefficients are
$$
\Gamma^t_{ij}=a\dot a\,\delta_{ij},
\qquad
\Gamma^i_{tj}=\Gamma^i_{jt}=H\delta^i_j.
$$
With $G_t^t=-3H^2$ and
$G_i^j=-(2\ddot a/a+H^2)\delta_i^j$, direct substitution into the mixed-tensor covariant derivative gives zero.