Changeset bd61b5c

Timestamp:

Apr 9, 2013, 1:39:16 PM (13 years ago)

Author:

Julian Iseringhausen <isering@…>

Children:

61edd6

Parents:

0bd47e

Message:

Fixed offset calculation in GSRB4.

File:

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

src/smoother/gsrb_poisson_4.cpp

@@ -77 +77 @@
 {
   const vmg_float prefactor_inv = Helper::pow_2(sol.Extent().MeshWidth().Max()) / 4.0;
-  const int off = rhs.Global().LocalBegin().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();