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

def test_dgbtrs(self):

    (lu_symm_band, ipiv, info) = dgbtrf(self.bandmat_real, self.KL, self.KU)
    (y, info) = dgbtrs(lu_symm_band, self.KL, self.KU, self.b, ipiv)
    y_lin = linalg.solve(self.real_mat, self.b)
    assert_array_almost_equal(y, y_lin)
