def METHOD_NAME(self):
    sampler = HigherOrderComposite(self.tracker)
    J = {'abc': -1, 'ab': 1}
    sampleset = sampler.sample_hising({}, J,
                                      initial_state={'a': 1, 'b': 1, 'c': -1})
    bqm = self.tracker.input['bqm']
    initial_state = self.tracker.input['initial_state']
    samples = dimod.ExactSolver().sample(bqm).samples()
    # make sure that the initial-state is minimzed over the product/aux
    mask = (samples[:, ['a', 'b', 'c']] == [1, 1, -1]).all(axis=1)
    for v, val in initial_state.items():
        if v in ['a', 'b', 'c']:
            continue
        self.assertTrue(samples[mask, [v]][0, 0], val)