Commits

Matthew Turk committed 707a83d

Removing sort-by-oct in NMSU-ART. This fixes most issues.

Level 0 seems confused still.

Comments (0)

Files changed (2)

yt/frontends/art/data_structures.py

                           LL[1]:RL[1]:NX[1]*1j,
                           LL[2]:RL[2]:NX[2]*1j ]
         root_fc= np.vstack([p.ravel() for p in root_fc]).T
-        nocts_check = oct_handler.add(1, 0, root_octs_side**3, root_fc, \
-                                      self.domain_id)
+        nocts_check = oct_handler.add(1, 0, root_octs_side**3,
+                                      root_fc, self.domain_id)
         assert(oct_handler.nocts == root_fc.shape[0])
         nocts_added = root_fc.shape[0]
         mylog.debug("Added %07i octs on level %02i, cumulative is %07i",

yt/frontends/art/io.py

 
     assert np.unique(iocts).shape[0] == nLevel
     
-    #ioct tries to access arrays much larger than le & fl
-    #just make sure they appear in the right order, skipping
-    #the empty space in between
-    idx = np.argsort(iocts)
-    
-    #now rearrange le & fl in order of the ioct
-    le = le[idx]
-    fl = fl[idx]
-
     #left edges are expressed as if they were on 
     #level 15, so no matter what level max(le)=2**15 
     #correct to the yt convention