def test_zgbtrf(self):
        M,N = shape(self.comp_mat)
        lu_symm_band, ipiv, info = zgbtrf(self.bandmat_comp, self.KL, self.KU)

        # extract matrix u from lu_symm_band
        u = diag(lu_symm_band[2*self.KL,:])
        for i in xrange(self.KL + self.KU):
            u += diag(lu_symm_band[2*self.KL-1-i,i+1:N], i+1)

        p_lin, l_lin, u_lin = lu(self.comp_mat, permute_l=0)
        assert_array_almost_equal(u, u_lin)
Compare zgbtrf  LU factorisation with the LU factorisation result
dataset/ETHPy150Open scipy/scipy/scipy/linalg/tests/test_decomp.py/TestEigBanded.test_zgbtrf