Changeset 14d38c for src/grid/grid_index_translations.cpp

Timestamp:

Apr 13, 2013, 4:28:08 PM (13 years ago)

Author:

Julian Iseringhausen <isering@…>

Children:

13b272

Parents:

1e45786

Message:

Renamed folder samples to discretization and worked on inner boundaries for open boundary conditions.

File:

• src/grid/grid_index_translations.cpp (modified) (1 diff)

src/grid/grid_index_translations.cpp

@@ -156 +156 @@
 
 void GridIndexTranslations::GetInnerBoundary(const Grid& grid_c, Index& b1_c, Index& b2_c,
-                                             const Grid& grid_f,   Index& b1_f, Index& b2_f)
+                                             const Grid& grid_f, Index& b1_f, Index& b2_f)
 {
   if (grid_f.Global().BoundaryType() == LocallyRefined) {
 
+    const Boundary& bc = MG::GetComm()->BoundaryConditions();
+
     const int global_mult = Helper::intpow(2, MG::GetSol()->MaxLevel() - grid_c.Level());
+
+    for (int i=0; i<3; ++i) {
+      if (bc[i] == Open) {
+        int begin_f_global = grid_f.Global().LocalBegin()[i] + (grid_f.Global().LocalBegin()[i] % 2 == 0 ? 0 : 1);
+        int end_f_global = grid_f.Global().LocalEnd()[i] + (grid_f.Global().LocalEnd()[i] % 2 == 0 ? 1 : 0);
+      } else {
+
+      }
+    }
 
     Index begin_finest = GlobalToGlobalFinest(grid_c, grid_c.Global().LocalBegin())