Changeset 43f24e5

Timestamp:

Apr 15, 2013, 1:43:31 PM (13 years ago)

Author:

Julian Iseringhausen <isering@…>

Children:

0174593

Parents:

d23483

Message:

Fixed CommMPI::PrintGrid in case of HAVE_VTK.

File:

• src/comm/comm_mpi.cpp (modified) (3 diffs)

src/comm/comm_mpi.cpp

@@ -344 +344 @@
 }
 
+//TODO: Variable length buffer
 void CommMPI::Print(const OutputLevel level, const char* format, ...)
 {
@@ -605 +606 @@
 
     Index i;
-    const Index offset = grid.Global().LocalBegin() - grid.Local().Begin();
-
-    for (i.X()=grid.Local().Begin().X(); i.X()<grid.Local().End().X(); ++i.X())
-      for (i.Y()=grid.Local().Begin().Y(); i.Y()<grid.Local().End().Y(); ++i.Y())
-        for (i.Z()=grid.Local().Begin().Z(); i.Z()<grid.Local().End().Z(); ++i.Z())
+    Index begin, end, offset;
+
+    for (int j=0; j<3; ++j) {
+      if (grid.Local().BoundarySize1()[j] > 0) {
+        begin[j] = grid.Local().BoundaryBegin1()[j];
+        offset[j] = grid.Global().LocalBegin()[j];
+      } else {
+        begin[j] = grid.Local().Begin()[j];
+        offset[j] =  grid.Global().LocalBegin()[j] - grid.Local().Begin()[j];
+      }
+      end[j] = (grid.Local().BoundarySize2()[j] > 0 ? grid.Local().BoundaryEnd2()[j] : grid.Local().End()[j]);
+    }
+
+    for (i.X()=begin.X(); i.X()<end.X(); ++i.X())
+      for (i.Y()=begin.Y(); i.Y()<end.Y(); ++i.Y())
+        for (i.Z()=begin.Z(); i.Z()<end.Z(); ++i.Z())
           image->SetScalarComponentFromDouble(i.X() + offset.X(),
                                               i.Y() + offset.Y(),
@@ -634 +646 @@
   }
 
+//FIXME
 #else /* HAVE_VTK */
   Index i;