Changeset 9e3e0a

Timestamp:

Apr 19, 2013, 8:39:37 PM (13 years ago)

Author:

Julian Iseringhausen <isering@…>

Children:

6f7f90

Parents:

540c9f8

Message:

Fixed offset for GSRB smoother.

File:

• src/smoother/gsrb.cpp (modified) (1 diff)

src/smoother/gsrb.cpp

@@ -82 +82 @@
   const vmg_float prefactor_inv = 1.0 / MG::GetDiscretization()->OperatorPrefactor(sol);
   const vmg_float diag_inv = 1.0 / mat.GetDiag();
-  const int off = rhs.Global().LocalBegin().Sum() - rhs.Global().GlobalBegin().Sum() - rhs.Local().HaloSize1().Sum();
+  const int off = rhs.Global().LocalBegin().Sum() - rhs.Global().GlobalBegin().Sum() + rhs.Local().HaloSize1().Sum();
   const LocalIndices& local = rhs.Local();
   Comm& comm = *MG::GetComm();