This is regarding the LU decomposition at page 52 of NR3.
It seems that only one index is tracked after exchange rows. Should we track both to have
indx[k]=imax;
index[imax]=k;
instead of only indx[k]=imax? Thanks!