[Yt-svn] commit/yt: MatthewTurk: Major fix for RAMSES. It took a sadly very long time to identify this typo,

[Bitbucket]

1 new changeset in yt:

http://bitbucket.org/yt_analysis/yt/changeset/afaef03af0ba/
changeset:   afaef03af0ba
branch:      yt
user:        MatthewTurk
date:        2011-08-23 13:14:10
summary:     Major fix for RAMSES.  It took a sadly very long time to identify this typo,
but this prevents duplicate creation of ride-side-split grids and reduces
substantially the overlap between grids; I can no longer identify any overlap
at all.  However, projecting Ones still does not give uniform 1.0 everywhere,
which may be a different problem.
affected #:  2 files (33 bytes)

--- a/yt/frontends/ramses/_ramses_reader.pyx	Sun Aug 21 10:22:51 2011 -0400
+++ b/yt/frontends/ramses/_ramses_reader.pyx	Tue Aug 23 07:14:10 2011 -0400
@@ -797,9 +797,9 @@
         self.grid_file_locations = gfl
         for gi in range(ng):
             if mask[gi] == 1:
-                grid_file_locations[gi,3] = left_edges[gi, 0]
-                grid_file_locations[gi,4] = left_edges[gi, 1]
-                grid_file_locations[gi,5] = left_edges[gi, 2]
+                grid_file_locations[gi,3] = left_edges[gi,0]
+                grid_file_locations[gi,4] = left_edges[gi,1]
+                grid_file_locations[gi,5] = left_edges[gi,2]
                 for i in range(6):
                     gfl[used, i] = grid_file_locations[gi,i]
                 used += 1
@@ -1047,7 +1047,7 @@
     if L.efficiency > 1.0: raise RuntimeError
     if L.efficiency <= 0.0: rv_l = []
     elif L.efficiency < min_eff:
-        rv_l = recursive_patch_splitting(L, dims_r, li_r,
+        rv_l = recursive_patch_splitting(L, dims_l, li_l,
                 left_index, fl, num_deep + 1)
     else:
         rv_l = [L]


--- a/yt/frontends/ramses/data_structures.py	Sun Aug 21 10:22:51 2011 -0400
+++ b/yt/frontends/ramses/data_structures.py	Tue Aug 23 07:14:10 2011 -0400
@@ -227,6 +227,7 @@
                 self.grid_levels[gi,:] = level
                 grids.append(self.grid(gi, self, level, fl, props[0,:]))
                 gi += 1
+        self.proto_grids = []
         self.grids = na.array(grids, dtype='object')
 
     def _populate_grid_objects(self):

Repository URL: https://bitbucket.org/yt_analysis/yt/

--

This is a commit notification from bitbucket.org. You are receiving
this because you have the service enabled, addressing the recipient of
this email.