Changeset dfed1c for test/interfaces/interface_sinus.cpp

Timestamp:

Nov 22, 2011, 9:22:10 PM (14 years ago)

Author:

Julian Iseringhausen <isering@…>

Children:

facba0

Parents:

66f24d

Message:

Major vmg update.

git-svn-id: ​https://svn.version.fz-juelich.de/scafacos/trunk@1136 5161e1c8-67bf-11de-9fd5-51895aff932f

File:

• test/interfaces/interface_sinus.cpp (modified) (2 diffs)

test/interfaces/interface_sinus.cpp

@@ -15 +15 @@
 
 #include <cmath>
+#include <iostream>
 
 #include "grid/grid.hpp"
@@ -34 +35 @@
   const double factor = 2.0 * M_PI * grid_l.Extent().Size().Max();
 
+  Index begin, end, begin_local;
+
+  for (int j=0; j<3; ++j) {
+
+    begin[j] = (grid_l.Local().HaloEnd1()[j] - grid_l.Local().HaloBegin1()[j] > 0 ? grid_l.Local().Begin()[j] : 0);
+    end[j] = (grid_l.Local().HaloEnd2()[j] - grid_l.Local().HaloBegin2()[j] > 0 ? grid_l.Local().End()[j] : grid_l.Local().SizeTotal()[j]);
+    begin_local[j] = grid_l.Global().BeginLocal()[j] - (grid_l.Local().HaloEnd1()[j] - grid_l.Local().HaloBegin1()[j] > 0 ? 1 : 0);
+  }
+
   grid_l.ClearBoundary();
 
-  for (i.X()=0; i.X()<grid_l.Local().SizeTotal().X(); ++i.X())
-    for (i.Y()=0; i.Y()<grid_l.Local().SizeTotal().Y(); ++i.Y())
-      for (i.Z()=0; i.Z()<grid_l.Local().SizeTotal().Z(); ++i.Z()) {
+  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()) {
 
-        pos = grid_l.Extent().Begin() + sp * static_cast<Vector>(grid_l.Global().Begin() + i);
+        pos = grid_l.Extent().Begin() + sp * static_cast<Vector>(begin_local + i);
 
         grid_l(i) = 3.0 * factor * factor * sin(factor * pos.X()) * sin(factor * pos.Y()) * sin(factor * pos.Z());