def METHOD_NAME(self):
    mol = gto.Mole()
    mol.atom = ''' O                  0.00000000    0.00000000   -0.11081188
                   H                 -0.00000000   -0.84695236    0.59109389
                   H                 -0.00000000    0.89830571    0.52404783 '''
    mol.verbose = 0
    mol.basis = '6-31g'
    mol.build()
    coords = [(0.5,0.6,0.1)]
    charges = [-0.1]
    mf = itrf.add_mm_charges(scf.RHF(mol), coords, charges).run(conv_tol=1e-10)
    mc = mcscf.CASCI(mf, 4, 4).run()
    self.assertAlmostEqual(mc.e_tot, -75.98156095286714, 8)
    mf = scf.RHF(mol).run(conv_tol=1e-10)
    mc = itrf.add_mm_charges(mcscf.CASCI(mf, 4, 4), coords, charges).run()
    self.assertAlmostEqual(mc.e_tot, -75.98156095286714, 8)