def METHOD_NAME(self):
    linear = {0: -0.5, 1: -0.3, 2: -0.8}
    quadratic = {(0, 1, 2): -1.7}
    sampler = HigherOrderComposite(ExactSolver())
    response = sampler.sample_hising(linear, quadratic, penalty_strength=10,
                                     keep_penalty_variables=True,
                                     discard_unsatisfied=True)
    self.assertEqual(len(response.first.sample), 5)
    self.assertAlmostEqual(response.first.energy, -3.3)
    self.assertTrue(response.first.penalty_satisfaction)