def METHOD_NAME():
    global cell, kmdf, auxcell, dfobj, kpts
    cell = pgto.M(
        a = np.eye(3) * 3.5,
        mesh = [11] * 3,
        atom = '''He    3.    2.       3.
                  He    1.    1.       1.''',
        basis = 'ccpvdz',
        verbose = 0,
        max_memory = 1000,
    )
    kpts = cell.make_kpts([3,5,6])[[0, 2, 3, 4, 6, 12, 20]]
    kpts[3] = kpts[0]-kpts[1]+kpts[2]
    kmdf = rsdf.RSDF(cell)
    kmdf.linear_dep_threshold = 1e-7
    kmdf.auxbasis = 'weigend'
    kmdf.kpts = kpts
    dfobj = rsdf.RSGDF(cell, kpts)
    dfobj.auxbasis = 'weigend'
    dfobj._rs_build()
    auxcell = dfobj.auxcell