piv pivsel mat;nrow;ncol;irow;jcol;v;w;i;m
LK'';''
nrowK(+/1^R(LKmKmat)),0RncolK(+/1URmat)
mat[v;w]K-m[v;wK((jcol=Incol)/Incol)]Pmat[irow;jcol]K`m[irow;jcol]
mat[v;]Kmat[v;]-mat[vK((irowBInrow)/Inrow);jcol]J.Pmat[irow;]Kmat[irow;]`mat[irowKpiv[1];jcolKpiv[2]]
mat
