def METHOD_NAME(self):
    dcm = DipoleChargeModifier(
        os.path.join(modifier_datapath, "dipole.pb"),
        [-1, -3],
        [1, 1, 1, 1, 1],
        1,
        0.25,
    )
    ve0, vf0, vv0 = dcm.eval(self.coords0, self.box0, self.atom_types0)
    ve1, vf1, vv1 = dcm.eval(self.coords1, self.box1, self.atom_types1)
    vf01 = vf0[:, self.idx_map, :]
    np.testing.assert_almost_equal(ve0, ve1, err_msg="energy should match")
    np.testing.assert_almost_equal(vv0, vv1, err_msg="virial should match")
    np.testing.assert_almost_equal(vf01, vf1, err_msg="force dose not match")