def METHOD_NAME(self):
    numpy.random.seed(1)
    nao = mol.nao_nr(cart=True)
    dm = numpy.random.random((nao,nao))
    vhfopt = _vhf.VHFOpt(mol, 'int2e_cart', 'CVHFnrs8_prescreen',
                         'CVHFsetnr_direct_scf',
                         'CVHFsetnr_direct_scf_dm')
    vj0, vk0 = _vhf.direct(dm, mol._atm, mol._bas, mol._env,
                           vhfopt=vhfopt, hermi=0, cart=True)
    vj = _vhf.direct_mapdm('int2e_cart', 's1', 'kl->s1ij', dm, 1,
                           mol._atm, mol._bas, mol._env, vhfopt)
    self.assertTrue(numpy.allclose(vj0, vj))
    vk = _vhf.direct_mapdm('int2e_cart', 's1', 'jk->s1il', [dm]*2, 1,
                           mol._atm, mol._bas, mol._env, vhfopt)
    self.assertTrue(numpy.allclose(vk0, vk[0]))
    self.assertTrue(numpy.allclose(vk0, vk[1]))