'Compare zhbevx eigenvalues and eigenvectors
with the result of linalg.eig.'

def test_zhbevx(self):

    (N, N) = shape(self.herm_mat)
    (w, evec, num, ifail, info) = zhbevx(self.bandmat_herm, 0.0, 0.0, 1, N, compute_v=1, range=2)
    evec_ = evec[:, argsort(w)]
    assert_array_almost_equal(sort(w), self.w_herm_lin)
    assert_array_almost_equal(abs(evec_), abs(self.evec_herm_lin))
