function c = cond1est (A) %COND1EST 1-norm condition estimate. % Example: % c = cond1est(A) % See also: cs_demo % Copyright 2006-2007, Timothy A. Davis. % http://www.cise.ufl.edu/research/sparse [m n] = size (A) ; if (m ~= n | ~isreal (A)) %#ok error ('A must be square and real') ; end if isempty(A) c = 0 ; return ; end [L,U,P,Q] = lu (A) ; if (~isempty (find (abs (diag (U)) == 0))) %#ok c = Inf ; else c = norm (A,1) * norm1est (L,U,P,Q) ; end