./cs_demo1 < ../Matrix/t1 T: CSparse Version 2.0.7, Dec 12, 2006. Copyright (c) Timothy A. Davis, 2006 triplet: 4-by-4, nzmax: 16 nnz: 10 2 2 : 3 1 0 : 3.1 3 3 : 1 0 2 : 3.2 1 1 : 2.9 3 0 : 3.5 3 1 : 0.4 1 3 : 0.9 0 0 : 4.5 2 1 : 1.7 A: CSparse Version 2.0.7, Dec 12, 2006. Copyright (c) Timothy A. Davis, 2006 4-by-4, nzmax: 10 nnz: 10, 1-norm: 11.1 col 0 : locations 0 to 2 1 : 3.1 3 : 3.5 0 : 4.5 col 1 : locations 3 to 5 1 : 2.9 3 : 0.4 2 : 1.7 col 2 : locations 6 to 7 2 : 3 0 : 3.2 col 3 : locations 8 to 9 3 : 1 1 : 0.9 AT: CSparse Version 2.0.7, Dec 12, 2006. Copyright (c) Timothy A. Davis, 2006 4-by-4, nzmax: 10 nnz: 10, 1-norm: 7.7 col 0 : locations 0 to 1 0 : 4.5 2 : 3.2 col 1 : locations 2 to 4 0 : 3.1 1 : 2.9 3 : 0.9 col 2 : locations 5 to 6 1 : 1.7 2 : 3 col 3 : locations 7 to 9 0 : 3.5 1 : 0.4 3 : 1 D: CSparse Version 2.0.7, Dec 12, 2006. Copyright (c) Timothy A. Davis, 2006 4-by-4, nzmax: 16 nnz: 16, 1-norm: 139.58 col 0 : locations 0 to 3 1 : 13.95 3 : 15.75 0 : 100.28 2 : 9.6 col 1 : locations 4 to 7 1 : 88.62 3 : 12.91 0 : 13.95 2 : 4.93 col 2 : locations 8 to 11 1 : 4.93 3 : 0.68 2 : 81.68 0 : 9.6 col 3 : locations 12 to 15 1 : 12.91 3 : 83.2 0 : 15.75 2 : 0.68 ./cs_demo2 < ../Matrix/t1 --- Matrix: 4-by-4, nnz: 10 (sym: 0: nnz 0), norm: 1.11e+01 blocks: 1 singletons: 0 structural rank: 4 QR natural time: 0.00 resid: 1.53e-17 QR amd(A'*A) time: 0.00 resid: 1.53e-17 LU natural time: 0.00 resid: 1.04e-17 LU amd(A+A') time: 0.00 resid: 4.94e-18 LU amd(S'*S) time: 0.00 resid: 4.94e-18 LU amd(A'*A) time: 0.00 resid: 4.94e-18 ./cs_demo2 < ../Matrix/ash219 --- Matrix: 219-by-85, nnz: 438 (sym: 0: nnz 0), norm: 9.00e+00 blocks: 1 singletons: 0 structural rank: 85 QR natural time: 0.00 resid: 1.61e-02 QR amd(A'*A) time: 0.00 resid: 1.61e-02 ./cs_demo2 < ../Matrix/bcsstk01 --- Matrix: 48-by-48, nnz: 224 (sym: -1: nnz 400), norm: 3.57e+09 blocks: 1 singletons: 0 structural rank: 48 QR natural time: 0.00 resid: 2.83e-19 QR amd(A'*A) time: 0.00 resid: 5.19e-19 LU natural time: 0.00 resid: 2.63e-19 LU amd(A+A') time: 0.00 resid: 8.63e-20 LU amd(S'*S) time: 0.00 resid: 2.04e-19 LU amd(A'*A) time: 0.00 resid: 2.04e-19 Chol natural time: 0.00 resid: 1.90e-19 Chol amd(A+A') time: 0.00 resid: 2.01e-19 ./cs_demo2 < ../Matrix/fs_183_1 --- Matrix: 183-by-183, nnz: 988 (sym: 0: nnz 0), norm: 1.70e+09 zero entries dropped: 71 tiny entries dropped: 10 blocks: 38 singletons: 37 structural rank: 183 QR natural time: 0.01 resid: 1.09e-27 QR amd(A'*A) time: 0.00 resid: 5.34e-28 LU natural time: 0.01 resid: 3.08e-28 LU amd(A+A') time: 0.00 resid: 1.42e-27 LU amd(S'*S) time: 0.00 resid: 7.11e-28 LU amd(A'*A) time: 0.00 resid: 7.11e-28 ./cs_demo2 < ../Matrix/mbeacxc --- Matrix: 492-by-490, nnz: 49920 (sym: 0: nnz 0), norm: 9.29e-01 blocks: 10 singletons: 8 structural rank: 448 QR natural time: 0.26 resid: nan QR amd(A'*A) time: 0.30 resid: nan ./cs_demo2 < ../Matrix/west0067 --- Matrix: 67-by-67, nnz: 294 (sym: 0: nnz 0), norm: 6.14e+00 blocks: 2 singletons: 1 structural rank: 67 QR natural time: 0.00 resid: 3.42e-17 QR amd(A'*A) time: 0.00 resid: 1.95e-17 LU natural time: 0.00 resid: 3.85e-17 LU amd(A+A') time: 0.00 resid: 1.95e-17 LU amd(S'*S) time: 0.00 resid: 2.60e-17 LU amd(A'*A) time: 0.00 resid: 2.60e-17 ./cs_demo2 < ../Matrix/lp_afiro --- Matrix: 27-by-51, nnz: 102 (sym: 0: nnz 0), norm: 3.43e+00 blocks: 1 singletons: 0 structural rank: 27 QR natural time: 0.00 resid: 9.54e-17 QR amd(A'*A) time: 0.00 resid: 1.89e-16 ./cs_demo2 < ../Matrix/bcsstk16 --- Matrix: 4884-by-4884, nnz: 147631 (sym: -1: nnz 290378), norm: 7.01e+09 blocks: 75 singletons: 74 structural rank: 4884 QR amd(A'*A) time: 2.43 resid: 2.02e-22 LU amd(A+A') time: 1.44 resid: 1.06e-22 LU amd(S'*S) time: 1.41 resid: 1.21e-22 LU amd(A'*A) time: 1.46 resid: 1.70e-22 Chol amd(A+A') time: 0.67 resid: 9.28e-23 ./cs_demo3 < ../Matrix/bcsstk01 --- Matrix: 48-by-48, nnz: 224 (sym: -1: nnz 400), norm: 3.57e+09 chol then update/downdate amd(A+A') symbolic chol time 0.00 numeric chol time 0.00 solve chol time 0.00 original: resid: 2.01e-19 update: time: 0.00 update: time: 0.00 (incl solve) resid: 8.39e-19 rechol: time: 0.00 (incl solve) resid: 1.20e-18 downdate: time: 0.00 downdate: time: 0.00 (incl solve) resid: 3.30e-17 ./cs_demo3 < ../Matrix/bcsstk16 --- Matrix: 4884-by-4884, nnz: 147631 (sym: -1: nnz 290378), norm: 7.01e+09 chol then update/downdate amd(A+A') symbolic chol time 0.03 numeric chol time 0.60 solve chol time 0.02 original: resid: 9.28e-23 update: time: 0.00 update: time: 0.01 (incl solve) resid: 8.62e-24 rechol: time: 0.62 (incl solve) resid: 8.72e-24 downdate: time: 0.00 downdate: time: 0.01 (incl solve) resid: 3.60e-22