def METHOD_NAME(self):
    pmol = mol.copy()
    pmol.charge = 1
    pmol.spin = 1
    pmol.build(False, False)
    mf = scf.density_fit(scf.ROHF(pmol), auxbasis='weigend')
    self.assertAlmostEqual(mf.scf(), -75.626515724371814, 8)