def METHOD_NAME(self, codim, subentity_codim):
    assert 0 <= codim <= 1, 'Invalid codimension'
    assert codim <= subentity_codim <= self.dim, 'Invalid subentity codimension'
    if codim == 0:
        if subentity_codim == 0:
            return np.arange(self.size(0), dtype='int32')[:, np.newaxis]
        else:
            return self.__subentities.T
    else:
        return super().METHOD_NAME(codim, subentity_codim)