[Yt-svn] yt-commit r592 - trunk/yt/lagos

[mturk at wrangler.dreamhost.com]

Author: mturk
Date: Sat Jun 21 10:03:58 2008
New Revision: 592
URL: http://yt.spacepope.org/changeset/592

Log:
Fixed for periodic boundary conditions.  Now it works on everything I have tested.  (Thanks for the pointer to the problem, DC!)



Modified:
   trunk/yt/lagos/BaseDataTypes.py
   trunk/yt/lagos/HierarchyType.py

Modified: trunk/yt/lagos/BaseDataTypes.py
==============================================================================
--- trunk/yt/lagos/BaseDataTypes.py	(original)
+++ trunk/yt/lagos/BaseDataTypes.py	Sat Jun 21 10:03:58 2008
@@ -1460,7 +1460,7 @@
         if na.any(self.left_edge < self.pf["DomainLeftEdge"]) or \
            na.any(self.right_edge > self.pf["DomainRightEdge"]):
             grids,ind = self.pf.hierarchy.get_periodic_box_grids(
-                            self.left_edge.copy(), self.right_edge.copy())
+                            self.left_edge, self.right_edge)
         else:
             grids,ind = self.pf.hierarchy.get_box_grids(
                             self.left_edge, self.right_edge)

Modified: trunk/yt/lagos/HierarchyType.py
==============================================================================
--- trunk/yt/lagos/HierarchyType.py	(original)
+++ trunk/yt/lagos/HierarchyType.py	Sat Jun 21 10:03:58 2008
@@ -629,13 +629,13 @@
         dr = self.parameters["DomainRightEdge"]
         db = right_edge - left_edge
         for off_x in [-1, 0, 1]:
-            nle = left_edge[:]
-            nre = left_edge[:]
-            nle[0] = dl[0] + (dr[0]-dl[0])*off_x
+            nle = left_edge.copy()
+            nre = left_edge.copy()
+            nle[0] = dl[0] + (dr[0]-dl[0])*off_x + left_edge[0]
             for off_y in [-1, 0, 1]:
-                nle[1] = dl[1] + (dr[1]-dl[1])*off_y
+                nle[1] = dl[1] + (dr[1]-dl[1])*off_y + left_edge[1]
                 for off_z in [-1, 0, 1]:
-                    nle[2] = dl[2] + (dr[2]-dl[2])*off_z
+                    nle[2] = dl[2] + (dr[2]-dl[2])*off_z + left_edge[2]
                     nre = nle + db
                     g, gi = self.get_box_grids(nle, nre)
                     mask[gi] = True