Changeset da7c98

Timestamp:

Apr 20, 2013, 4:15:44 PM (13 years ago)

Author:

Julian Iseringhausen <isering@…>

Children:

252f5c

Parents:

881e4b

Message:

Experimental fix.

File:

• src/discretization/discretization_poisson_fv.cpp (modified) (2 diffs)

src/discretization/discretization_poisson_fv.cpp

@@ -106 +106 @@
   const Boundary& bc = comm.BoundaryConditions();
   const Index off((GridIndexTranslations::LocalToGlobal(sol_f, sol_f.Local().Begin())[0] % 2 == 0 ? 0 : 1),
-      (GridIndexTranslations::LocalToGlobal(sol_f, sol_f.Local().Begin())[1] % 2 == 0 ? 0 : 1),
-      (GridIndexTranslations::LocalToGlobal(sol_f, sol_f.Local().Begin())[2] % 2 == 0 ? 0 : 1));
+                  (GridIndexTranslations::LocalToGlobal(sol_f, sol_f.Local().Begin())[1] % 2 == 0 ? 0 : 1),
+                  (GridIndexTranslations::LocalToGlobal(sol_f, sol_f.Local().Begin())[2] % 2 == 0 ? 0 : 1));
 
   const Index begin_f = sol_f.Local().Begin() - off;
@@ -114 +114 @@
   const Index begin_c = GridIndexTranslations::GlobalFinestToLocal(sol_c, GridIndexTranslations::LocalToGlobalFinest(sol_f, begin_f));
 
-  const Index b1_finest = GridIndexTranslations::LocalToGlobalFinest(sol_f, sol_f.Local().Begin());
-  const Index b2_finest = GridIndexTranslations::LocalToGlobalFinest(sol_f, sol_f.Local().End()-1);
-
   const Index b1_f = sol_f.Local().BoundaryBegin1();
   const Index b2_f = sol_f.Local().BoundaryBegin2();
+
+  const Index b1_finest = GridIndexTranslations::LocalToGlobalFinest(sol_f, sol_f.Local().Begin() + off);
+  const Index b2_finest = GridIndexTranslations::LocalToGlobalFinest(sol_f, sol_f.Local().End()-1);
 
   const Index b1_c = GridIndexTranslations::GlobalFinestToLocal(sol_c, b1_finest) - 1;