1function test13
2%TEST13 test cholmod2 and MATLAB on large tridiagonal matrices
3% Example:
4%   test13
5% See also cholmod_test
6
7% Copyright 2007, Timothy A. Davis, http://www.suitesparse.com
8
9fprintf ('=================================================================\n');
10fprintf ('test13: test cholmod2 and MATLAB on large tridiagonal matrices\n') ;
11
12for n = [10000 1e4 1e5 1e6]
13
14    e = ones (n,1) ;
15    A = spdiags ([e 4*e e], -1:1, n, n) ;
16    clear e
17    b = rand (n,1) ;
18
19    tic ;
20    x = cholmod2 (A,b) ;
21    t2 = toc ;
22    e = norm (A*x-b,1) ;
23    fprintf ('n %9d   cholmod2 %8.2f  err %6.1e\n', n, t2, e) ;
24
25    tic ;
26    x = A\b ;
27    t1 = toc ;
28    e = norm (A*x-b,1) ;
29    fprintf ('n %9d   matlab  %8.2f  err %6.1e\n', n, t1, e) ;
30
31    clear A b
32
33end
34
35
36