2
function growth = pivgrowth(A)
4
[p,q,r,lnz,Info1] = klua (A) ;
5
[l,u,off,pnum,rs,Info2] = kluf (A, p,q,r,lnz,Info1) ;
7
A = rs(pnum, pnum) \ A(pnum,q) ;
8
growth = ones(Info1(2), 1);
17
growth(col) = max(abs(A(start:en,col))) / max(abs(u(:,col))) ;
20
growth = min(growth) ;