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