1function err = test_all_cod 2%TEST_ALL_COD test the COD factorization 3% 4% Example 5% test_all_cod 6% 7% See also test_cod 8 9% Copyright 2011-2012, Timothy A. Davis, http://www.suitesparse.com 10 11err = 0 ; 12for m = 0:10 13 for n = 0:10 14 for sp = 0:1 15 for im = 0:1 16 A = sprand (m, n, 0.2) ; 17 if (im) 18 A = A + 1i*sprand (m, n, 0.1) ; 19 end 20 if (~sp) 21 A = full (A) ; 22 end 23 err = max (err, test_cod (A)) ; 24 err = max (err, test_cod (A, 1e-12)) ; 25 end 26 end 27 end 28end 29 30fprintf ('test COD, COD_SPARSE, and RQ: error %g\n', err) ; 31