1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
|
function c = cond1est (A)
%COND1EST 1-norm condition estimate.
% Example:
% c = cond1est(A)
% See also: cs_demo
% CXSparse, Copyright (c) 2006-2022, Timothy A. Davis. All Rights Reserved.
% SPDX-License-Identifier: LGPL-2.1+
[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
|
