def METHOD_NAME(self):
    m = pyo.ConcreteModel()
    m.x = pyo.Var()
    m.y = pyo.Var()
    m.z = pyo.Var()
    m.c1 = pyo.Constraint(expr=2 * m.x >= m.y**2)
    m.c2 = pyo.Constraint(expr=m.x**2 >= m.y**2 + m.z**2)
    m.c3 = pyo.Constraint(expr=m.z >= 0)
    m.c4 = pyo.Constraint(expr=m.x + m.y >= 0)
    opt = pyo.SolverFactory('mosek_persistent')
    opt.set_instance(m)
    self.assertEqual(opt._solver_model.getnumcon(), 4)
    opt.remove_constraint(m.c1)
    self.assertEqual(opt._solver_model.getnumcon(), 3)
    opt.remove_constraints(m.c2, m.c3)
    self.assertEqual(opt._solver_model.getnumcon(), 1)
    opt.add_constraint(m.c1)
    self.assertEqual(opt._solver_model.getnumcon(), 2)
    self.assertRaises(ValueError, opt.remove_constraint, m.c2)