V>    while (i <= j)
V>    {
V>        while (*matrix(i) < p && i <= last) ++i;
V>        while (*matrix(j) > p && j >= first) --j;
V>        if (i <= j)
V>        {
V>            temp = *matrix(i);
V>            idxTemp = idx[i - 1];
V>            *matrix(i) = *matrix(j);
V>            *matrix(j) = temp;
V>            idx[i - 1] = idx[j - 1];
V>            idx[j - 1] = idxTemp;
V>            ++i; --j;
V>        }
V>    }
V>