2
matrices = {'Asmall', 'Amiddle', 'Alarge'} ;
5
for i = 1:length(matrices)
7
fprintf ('\n\nLoading matrix: %s\n', matrices {i}) ;
12
% solve A*x = y with KLU + BTF + AMD
14
% symtol: partial pivoting tolerance. Use 1e-3 (for example) to prefer diagonal
17
fprintf ('calling klu\n') ;
19
[x, Info] = klus (A, y, symtol, 'print info') ;
24
% Info (2): nz in off diagonal part
25
% Info (3): # of blocks
26
% Info (4): max nz in diagonal blocks of A
27
% Info (5): dimension of largest block
28
% Info (6): estimated nz in L, incl. diagonal, excludes off-diagonal entries
29
% Info (7): estimated nz in U, incl. diagonal, excludes off-diagonal entries
31
% Info (8): nz in L, including diagonal, excludes off-diagonal entries
32
% Info (9): nz in U, including diagonal, excludes off-diagonal entries
33
% Info (10): analyze cputime
34
% Info (11): factor cputime
35
% Info (12): solve cputime
36
% Info (13): refactorize cputime (if computed)
37
% Info (14): # off-diagonal pivots chosen
39
% b may be n-by-m with m > 1. It must be dense.
42
fprintf ('total wallclock time %g\n', t) ;
43
fprintf ('Info (2): %g nz in offdiag part\n', Info (2)) ;
44
fprintf ('Info (3): %g # of blocks\n', Info (3)) ;
45
fprintf ('Info (4): %g max nz in diagonal blocks of A\n', Info (4)) ;
46
fprintf ('Info (5): %g dimension of largest block\n', Info (5)) ;
47
fprintf ('Info (6): %g estimated nz in L, incl. diagonal, excludes off-diagonal entries\n', Info (6)) ;
48
fprintf ('Info (7): %g estimated nz in U, incl. diagonal, excludes off-diagonal entries\n', Info (7)) ;
49
fprintf ('Info (8): %g nz in L, including diagonal, excludes off-diagonal entries\n', Info (8)) ;
50
fprintf ('Info (9): %g nz in U, including diagonal, excludes off-diagonal entries\n', Info (9)) ;
51
fprintf ('Info (10): %g analyze cputime\n', Info (10)) ;
52
fprintf ('Info (11): %g factor cputime\n', Info (11)) ;
53
fprintf ('Info (12): %g solve cputime\n', Info (12)) ;
54
fprintf ('Info (13): %g refactorize cputime (if computed)\n', Info (13)) ;
55
fprintf ('Info (14): %g # off-diagonal pivots chosen\n', Info (14)) ;