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)
Compare zgbtrs  solutions for linear equation system  A*x = b
dataset/ETHPy150Open scipy/scipy/scipy/linalg/tests/test_decomp.py/TestEigBanded.test_zgbtrs