def METHOD_NAME(self):
    bqm = oj.BinaryQuadraticModel(self.h, self.J, 'SPIN')
    ising_matrix = np.array([
        [0, -1,  0,  0, 1],
        [0, 0, -3, 0, -2],
        [0, 0, 0, 0.5, 0],
        [0, 0, 0, 0, 0],
        [0, 0, 0, 0, 1]
    ])
    np.testing.assert_array_equal(
        bqm.interaction_matrix(), ising_matrix
    )