additive_schwarz = richardson + additive_schwarz_correction

Author: fverdugo

extensions/PartitionedSolvers/src/PartitionedSolvers.jl

@@ -22,6 +22,7 @@ export jacobi_correction
 export richardson
 export jacobi
 export gauss_seidel
+export additive_schwarz_correction
 export additive_schwarz
 include("smoothers.jl")
 

extensions/PartitionedSolvers/src/amg.jl

@@ -312,7 +312,7 @@ function smoothed_aggregation(;
 end
 
 function amg_level_params(;
-    pre_smoother = richardson(additive_schwarz(gauss_seidel(;iters=1));iters=1),
+    pre_smoother = additive_schwarz(gauss_seidel(;iters=1);iters=1),
     coarsening = smoothed_aggregation(;),
     cycle = v_cycle,
     pos_smoother = pre_smoother,

extensions/PartitionedSolvers/src/smoothers.jl

@@ -102,8 +102,11 @@ function gauss_seidel(;iters=1,sweep=:symmetric)
     linear_solver(;setup,setup!,solve!)
 end
 
+function additive_schwarz(local_solver;iters=1)
+    richardson(additive_schwarz_correction(local_solver);iters)
+end
 
-function additive_schwarz(local_solver)
+function additive_schwarz_correction(local_solver)
     function build_local_operators(O::MatrixWithNullspace)
         A = matrix(O)
         ns = map(i->own_values(i),nullspace(O))