Fix of coarse operator. Additive block now also works.

kyrillh · Oct 14, 2024 · c689de7 · c689de7
1 parent 933eee1
commit c689de7
Show file tree

Hide file tree

Showing 2 changed files with 2 additions and 2 deletions.
diff --git a/feddlib/problems/Solver/NonLinearSchwarzSolver/CoarseNonLinearSchwarzOperator_def.hpp b/feddlib/problems/Solver/NonLinearSchwarzSolver/CoarseNonLinearSchwarzOperator_def.hpp
@@ -473,7 +473,7 @@ void CoarseNonLinearSchwarzOperator<SC, LO, GO, NO>::apply(const XMultiVector &x
     feddY->split();
     apply(feddX, feddY, alpha, beta);
     feddY->merge();
-    y.update(alpha, *feddY->getMergedVector()->getXpetraMultiVector(), beta);
+    y.update(ST::one(), *feddY->getMergedVector()->getXpetraMultiVector(), ST::zero());
 }
 
 template <class SC, class LO, class GO, class NO>

diff --git a/feddlib/problems/Solver/NonLinearSchwarzSolver/NonLinearSchwarzOperator_def.hpp b/feddlib/problems/Solver/NonLinearSchwarzSolver/NonLinearSchwarzOperator_def.hpp
@@ -476,7 +476,7 @@ void NonLinearSchwarzOperator<SC, LO, GO, NO>::apply(const XMultiVector &x, XMul
     feddY->split();
     apply(feddX, feddY, alpha, beta);
     feddY->merge();
-    y.update(alpha, *feddY->getMergedVector()->getXpetraMultiVector(), beta);
+    y.update(ST::one(), *feddY->getMergedVector()->getXpetraMultiVector(), ST::zero());
 }
 
 template <class SC, class LO, class GO, class NO>