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