def METHOD_NAME(self):
    beta, eps, w_max = 1.337, 1e-12, 10.
    m = MeshDLR(beta, 'Fermion', w_max, eps)
    rf = np.array([ p.value for p in m ])
    g = Gf(mesh=m, target_shape=[])
    g.data[:] = np.random.randn(len(m))
    # block gf test
    Bgdlr = BlockGf(block_list=[g, g], make_copies=True)
    for tau in np.linspace(0, beta, num=10):
        ref = np.sum(g.data * -np.exp(-rf * tau/beta) / ( 1 + np.exp(-rf) ))
        np.testing.assert_almost_equal(g(tau), ref)
        for block, gdlr in Bgdlr:
            np.testing.assert_almost_equal(gdlr(tau), ref)