def METHOD_NAME(self):
    loc = pipek.PipekMezey(h2o, mf_h2o.mo_coeff[:,3:8])
    loc.pop_method = 'becke'
    mo = loc.kernel()
    pop = pipek.atomic_pops(h2o, mo)
    z = numpy.einsum('xii,xii->', pop, pop)
    self.assertAlmostEqual(z, 3.548139685217463, 4)