human
I wrote the below code

    (mu, nu, rho, sigma) = tensor_indices('mu,    nu,    rho,    sigma', LorentzIndex)
    (a1, a2, a3, a4, a5, a6) = tensor_indices('a1:7', LorentzIndex)
    (mu11, mu12, mu21, mu31, mu32, mu41, mu51, mu52) = tensor_indices('mu11,    mu12,    mu21,    mu31,    mu32,    mu41,    mu51,    mu52', LorentzIndex)
    (mu61, mu71, mu72) = tensor_indices('mu61,    mu71,    mu72', LorentzIndex)
    (m0, m1, m2, m3, m4, m5, m6) = tensor_indices('m0:7', LorentzIndex)
    def g(xx, yy):
        return (((G(xx) * G(yy)) + (G(yy) * G(xx))) / 2)
    if (D == 4):
        t = ((((((G(a1) * G(mu11)) * G(a2)) * G(mu21)) * G((- a1))) * G(mu31)) * G((- a2)))
        assert _is_tensor_eq(tfunc(t), (((((-4) * G(mu11)) * G(mu31)) * G(mu21)) - (((4 * G(mu31)) * G(mu11)) * G(mu21))))
        t = ((((((((((((((((((G(a1) * G(mu11)) * G(mu12)) * G(a2)) * G(mu21)) * G(a3)) * G(mu31)) * G(mu32)) * G(a4)) * G(mu41)) * G((- a2))) * G(mu51)) * G(mu52)) * G((- a1))) * G(mu61)) * G((- a3))) * G(mu71)) * G(mu72)) * G((- a4)))
        assert _is_tensor_eq(tfunc(t), ((((((((((((((16 * G(mu31)) * G(mu32)) * G(mu72)) * G(mu71)) * G(mu11)) * G(mu52)) * G(mu51)) * G(mu12)) * G(mu61)) * G(mu21)) * G(mu41)) + (((((((((((16 * G(mu31)) * G(mu32)) * G(mu72)) * G(mu71)) * G(mu12)) * G(mu51)) * G(mu52)) * G(mu11)) * G(mu61)) * G(mu21)) * G(mu41))) + (((((((((((16 * G(mu71)) * G(mu72)) * G(mu32)) * G(mu31)) * G(mu11)) * G(mu52)) * G(mu51)) * G(mu12)) * G(mu61)) * G(mu21)) * G(mu41))) + (((((((((((16 * G(mu71)) * G(mu72)) * G(mu32)) * G(mu31)) * G(mu12)) * G(mu51)) * G(mu52)) * G(mu11)) * G(mu61)) * G(mu21)) * G(mu41))))
    def add_delta(ne):
        return (ne * eye(4))
    t = (G(mu) * G((- mu)))
    ts = add_delta(D)
    assert _is_tensor_eq(tfunc(t), ts)
    t = (((G(mu) * G(nu)) * G((- mu))) * G((- nu)))
    ts = add_delta(((2 * D) - (D ** 2)))
    assert _is_tensor_eq(tfunc(t), ts)
    t = (((G(mu) * G(nu)) * G((- nu))) * G((- mu)))
    ts = add_delta((D ** 2))
    assert _is_tensor_eq(tfunc(t), ts)
    t = (((((G(mu) * G(nu)) * G((- rho))) * G((-...

What's a good function header?

Here's a good function header?
assistant
#Here's a good function header.

def execute_gamma_simplify_tests_for_function(tfunc, D):