if pos < nr_sites - var_sites:
                    # We always do this, because we don't do the last site again.
                    self.canonicalize(right=pos_end)
                    lvecs[pos + 1] = None
                    rvecs[pos] = _adapt_to_add_r(rvecs[pos + 1], self._lt[pos_end],