# square-n
# cg-iter sans precond: asymp en O(n^3)=O(N^1.5)
# mumps en O(n^2.25)=O(N^1.125)
# cg-iter avec precond mic0: asymp en O(n^2.5)=O(N^1.25) et plus petite cte ?
# => meth iter devient interessante en 2d ?
# n cg-iter mumps_fac mumps_sol umfpack_tot cholmod_tot
40 0.00282979 0.0146101 0.000816822 0.00716496 0.0154388
80 0.0149441 0.073535 0.00297904 0.0468049 0.104896
160 0.108872 0.417163 0.013027 0.353471 1.43207
320 0.860483 2.48603 0.0615819 2.87632 8.60481
640 8.78795 15.3644 0.274312 25.7161 76.3135
1280 75.0865 110.212 1.22439 226.163 612.257
|