[Yt-svn] yt-commit r592 - trunk/yt/lagos
mturk at wrangler.dreamhost.com
mturk at wrangler.dreamhost.com
Sat Jun 21 10:04:00 PDT 2008
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
More information about the yt-svn
mailing list