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