i_max  := argmax (i = 1 ... k , abs(A[i, k]))
if A[i_max, k] = 0
     k = k - 1
     continue
swap rows(k, i_max)