JuliaLinearAlgebra / AlgebraicMultigrid.jl

Algebraic Multigrid in Julia
Other
117 stars 22 forks source link

Not cutting off at the right number of levels #18

Closed ranjanan closed 6 years ago

ranjanan commented 6 years ago

The hierarchy adds one level extra.

A = float.(poisson(10^6))
smoothed_aggregation(A)

Multilevel Solver
-----------------
Operator Complexity: 1.5
Grid Complexity: 1.5
No. of Levels: 11
Coarse Solver: AMG.Pinv()
Level     Unknowns     NonZeros
-----     --------     --------
    1      1000000      2999998 [66.67%]
    2       333334      1000000 [22.22%]
    3       111112       333334 [ 7.41%]
    4        37038       111112 [ 2.47%]
    5        12346        37036 [ 0.82%]
    6         4116        12346 [ 0.27%]
    7         1372         4114 [ 0.09%]
    8          458         1372 [ 0.03%]
    9          153          457 [ 0.01%]
   10           51          151 [ 0.00%]
   11           17           49 [ 0.00%]