search

JuliaLinearAlgebra / AlgebraicMultigrid.jl

Algebraic Multigrid in Julia
Other
116 stars 23 forks source link

SOR Smoother #112

Open Fastaxx opened 3 months ago

Fastaxx commented 3 months ago

Add a simple SOR Smoother.

A = poisson(1000)
sor_smoother = SOR(1.2; iter=2)
ml = smoothed_aggregation(A; presmoother=sor_smoother, postsmoother=sor_smoother)
x = AlgebraicMultigrid._solve(ml, A * ones(1000))