diff options
Diffstat (limited to 'modules/umfpack/examples/bench.txt')
-rwxr-xr-x | modules/umfpack/examples/bench.txt | 101 |
1 files changed, 101 insertions, 0 deletions
diff --git a/modules/umfpack/examples/bench.txt b/modules/umfpack/examples/bench.txt new file mode 100755 index 000000000..e5ba9f7d4 --- /dev/null +++ b/modules/umfpack/examples/bench.txt @@ -0,0 +1,101 @@ + + A test on few matrices between umfpack V3 and Sparse 1.3 + via their Scilab interface. The test consists in loading + some Harwell-Boeing sparse matrices and solve a linear + system with a random rhs via umfpack and sparse 1.3 : + + t1 = 0 + for i = 1:nb_rep + hrtimer(); x = umfpack(A,"\",b); t1 = t1 + hrtimer(); + end + t1 = t1 / nb_rep + + t2 = 0 + for i = 1:nb_rep + hrtimer(); x = A\b; t2 = t2 + hrtimer(); + end + t2 = t2 / nb_rep + + then display t1 and t2 and also the norm 2 residual ; i have + used nb_rep = 10. + + The computer is a 933 M Hz Pentium III / 1 Go Memory / 256 Mo cache + running Gnu-Linux and the ATLAS BLAS is used. + + n : order of the matrix + nnz : number of non zeros elements + t1 : mean time for solving Ax=b with umfpack + t2 : mean time for solving Ax=b with Sparse 1.3 (current + package to solve sparse systems in scilab) + + + The matrices come from the Florida University Sparse Collection + (www.cise.ufl.edu/research/sparse/matrices) maintained by Tim + Davis the author of UMFPACK + + The hrtimer is a scilab extension package from Christoph L. Spiel + & Lydia E. van Dijk (www.hammersmith-consulting.com). + + ----------------------------------------------------------------- + + arc130.rua : n = 130 , nnz = 1037 + + mean time umfpack V3 : t1 = 0.002783 , ||Ax-b|| = 6.538563e-11 + mean time sparse 1.3 : t2 = 0.015559 , ||Ax-b|| = 1.026040e-10 + speed up : t2/t1 = 5.591677 + ----------------------------------------------------------------- + + bcsstk24.rsa : n = 3562 , nnz = 159910 + + mean time umfpack V3 : t1 = 1.540691 , ||Ax-b|| = 1.298751e-07 + mean time sparse 1.3 : t2 = 52.546995 , ||Ax-b|| = 1.299314e-07 + speed up : t2/t1 = 34.106113 + ----------------------------------------------------------------- + + bcsstk34.rsa : n = 588 , nnz = 21418 + + mean time umfpack V3 : t1 = 0.068137 , ||Ax-b|| = 2.921253e-14 + mean time sparse 1.3 : t2 = 4.162029 , ||Ax-b|| = 6.665581e-14 + speed up : t2/t1 = 61.083241 + ----------------------------------------------------------------- + + bcsstk38.rsa : n = 8032 , nnz = 355460 + + mean time umfpack V3 : t1 = 4.671143 , ||Ax-b|| = 1.465322e-09 + mean time sparse 1.3 : t2 = 295.272744 , ||Ax-b|| = 5.456355e-09 + speed up : t2/t1 = 63.212092 + ----------------------------------------------------------------- + + cavity16.rua : n = 4562 , nnz = 137887 + + mean time umfpack V3 : t1 = 0.915590 , ||Ax-b|| = 3.816630e-11 + mean time sparse 1.3 : t2 = 21.034642 , ||Ax-b|| = 1.278231e-10 + speed up : t2/t1 = 22.973871 + ----------------------------------------------------------------- + + epb2.rua : n = 25228 , nnz = 175027 + + mean time umfpack V3 : t1 = 5.606093 , ||Ax-b|| = 1.141705e-12 + mean time sparse 1.3 : t2 = 951.516539 , ||Ax-b|| = 5.832767e-12 + speed up : t2/t1 = 169.729009 + ----------------------------------------------------------------- + + ex14.rua : n = 3251 , nnz = 65875 + + mean time umfpack V3 : t1 = 1.014328 , ||Ax-b|| = 7.475711e-07 + mean time sparse 1.3 : t2 = 159.702604 , ||Ax-b|| = 2.502100e-06 + speed up : t2/t1 = 157.446646 + ----------------------------------------------------------------- + + utm5940.rua : n = 5940 , nnz = 83842 + + mean time umfpack V3 : t1 = 1.037315 , ||Ax-b|| = 7.392123e-08 + mean time sparse 1.3 : t2 = 288.563648 , ||Ax-b|| = 1.583351e-07 + speed up : t2/t1 = 278.183294 + ----------------------------------------------------------------- + + wang1.rua : n = 2903 , nnz = 19093 + + mean time umfpack V3 : t1 = 0.654956 , ||Ax-b|| = 1.286372e-12 + mean time sparse 1.3 : t2 = 29.773115 , ||Ax-b|| = 6.243728e-12 + speed up : t2/t1 = 45.458204 |