def METHOD_NAME(self):
    nao = mol.nao_nr()
    naoaux = auxmol.nao_nr()
    eri0 = numpy.empty((nao,nao,naoaux))
    pi = 0
    for i in range(mol.nbas):
        pj = 0
        for j in range(mol.nbas):
            pk = 0
            for k in range(mol.nbas, mol.nbas+auxmol.nbas):
                shls = (i, j, k)
                buf = gto.moleintor.getints_by_shell('int3c2e_sph',
                                                     shls, atm, bas, env)
                di, dj, dk = buf.shape
                eri0[pi:pi+di,pj:pj+dj,pk:pk+dk] = buf
                pk += dk
            pj += dj
        pi += di
    j3c = df.incore.aux_e2(mol, auxmol, intor='int3c2e_sph', aosym='s1')
    self.assertTrue(numpy.allclose(eri0, j3c.reshape(nao,nao,naoaux)))
    self.assertAlmostEqual(lib.finger(j3c), 45.27912877994409, 9)
    idx = numpy.tril_indices(nao)
    j3c = df.incore.aux_e2(mol, auxmol, intor='int3c2e_sph', aosym='s2ij')
    self.assertTrue(numpy.allclose(eri0[idx], j3c))
    self.assertAlmostEqual(lib.finger(j3c), 12.407403711205063, 9)