'Compare zgbtrs  solutions for linear equation system  A*x = b
with solutions of linalg.solve.'

def test_zgbtrs(self):

    (lu_symm_band, ipiv, info) = zgbtrf(self.bandmat_comp, self.KL, self.KU)
    (y, info) = zgbtrs(lu_symm_band, self.KL, self.KU, self.bc, ipiv)
    y_lin = linalg.solve(self.comp_mat, self.bc)
    assert_array_almost_equal(y, y_lin)
