1function test16 2%TEST16 test cholmod2 on a large matrix 3% Example: 4% test16 5% See also cholmod_test 6 7% Copyright 2007, Timothy A. Davis, http://www.suitesparse.com 8 9fprintf ('=================================================================\n'); 10fprintf ('test16: test cholmod2 on a large matrix\n') ; 11 12rand ('state',1) ; 13randn ('state',1) ; 14 15Prob = ssget (936) %#ok 16A = Prob.A ; 17% tic 18% [L,s,p] = lchol (A) ; 19% toc 20% norm (L,1) 21 22n = size (A,1) ; 23b = rand (n,1) ; 24tic 25x = cholmod2(A,b) ; 26t = toc ; 27fprintf ('time %g\n', t) ; 28err = norm (A*x-b) ; 29 30if (err > 1e-5) 31 error ('!') ; 32end 33 34fprintf ('test16 passed\n') ; 35