[yt-dev] Issue #1017: Streamlines crashes with python 3 (yt_analysis/yt)

[Jill Naiman]

New issue 1017: Streamlines crashes with python 3
https://bitbucket.org/yt_analysis/yt/issue/1017/streamlines-crashes-with-python-3

Jill Naiman:

I'm on version:

Version = 3.2-dev
Changeset = 39f9c8e0e97c+ yt

And I built thinks with python 3.  

Code is here: http://paste.yt-project.org/show/5529/

Error is:

```
RuntimeError                              Traceback (most recent call last)
/Users/jillnaiman1/astroblend-dev/exampleScripts/streamlines6.py in <module>()
     23 pos = c+pos_dx
     24 
---> 25 streamlines = Streamlines(ds,pos,'velocity_x', 'velocity_y', 'velocity_z', length=1.0)
     26 
     27 

/Users/jillnaiman1/yt-x86_64/lib/python3.4/site-packages/yt-3.2dev-py3.4-macosx-10.10-x86_64.egg/yt/visualization/streamlines.py in __init__(self, ds, positions, xfield, yfield, zfield, volume, dx, length, direction, get_magnitude)
    118         if volume is None:
    119             volume = AMRKDTree(self.ds)
--> 120             volume.set_fields([self.xfield,self.yfield,self.zfield],[False,False,False], False)
    121             volume.join_parallel_trees()
    122         self.volume = volume

/Users/jillnaiman1/yt-x86_64/lib/python3.4/site-packages/yt-3.2dev-py3.4-macosx-10.10-x86_64.egg/yt/utilities/amr_kdtree/amr_kdtree.py in set_fields(self, fields, log_fields, no_ghost)
    168         self.brick_dimensions = []
    169         bricks = []
--> 170         for b in self.traverse():
    171             bricks.append(b)
    172         self.bricks = np.array(bricks)

/Users/jillnaiman1/yt-x86_64/lib/python3.4/site-packages/yt-3.2dev-py3.4-macosx-10.10-x86_64.egg/yt/utilities/amr_kdtree/amr_kdtree.py in traverse(self, viewpoint)
    182     def traverse(self, viewpoint=None):
    183         for node in kd_traverse(self.tree.trunk, viewpoint=viewpoint):
--> 184             yield self.get_brick_data(node)
    185 
    186     def slice_traverse(self, viewpoint = None):

/Users/jillnaiman1/yt-x86_64/lib/python3.4/site-packages/yt-3.2dev-py3.4-macosx-10.10-x86_64.egg/yt/utilities/amr_kdtree/amr_kdtree.py in get_brick_data(self, node)
    269             dds = []
    270             for i, field in enumerate(self.fields):
--> 271                 vcd = grid.get_vertex_centered_data(field, smoothed=True, no_ghost=self.no_ghost).astype('float64')
    272                 if self.log_fields[i]: vcd = np.log10(vcd)
    273                 dds.append(vcd)

/Users/jillnaiman1/yt-x86_64/lib/python3.4/site-packages/yt-3.2dev-py3.4-macosx-10.10-x86_64.egg/yt/data_objects/grid_patch.py in get_vertex_centered_data(self, field, smoothed, no_ghost)
    277                                    new_field, output_left)
    278         else:
--> 279             cg = self.retrieve_ghost_zones(1, field, smoothed=smoothed)
    280             np.add(new_field, cg[field][1: ,1: ,1: ], new_field)
    281             np.add(new_field, cg[field][:-1,1: ,1: ], new_field)

/Users/jillnaiman1/yt-x86_64/lib/python3.4/site-packages/yt-3.2dev-py3.4-macosx-10.10-x86_64.egg/yt/data_objects/grid_patch.py in retrieve_ghost_zones(self, n_zones, fields, all_levels, smoothed)
    254                 level, new_left_edge,
    255                 field_parameters = field_parameters,
--> 256                 **kwargs)
    257         else:
    258             cube = self.ds.covering_grid(level, new_left_edge,

/Users/jillnaiman1/yt-x86_64/lib/python3.4/site-packages/yt-3.2dev-py3.4-macosx-10.10-x86_64.egg/yt/data_objects/construction_data_containers.py in __init__(self, *args, **kwargs)
    811                self.ds.domain_dimensions.astype("float64"))
    812         self.global_endindex = None
--> 813         YTCoveringGridBase.__init__(self, *args, **kwargs)
    814         self._final_start_index = self.global_startindex
    815 

/Users/jillnaiman1/yt-x86_64/lib/python3.4/site-packages/yt-3.2dev-py3.4-macosx-10.10-x86_64.egg/yt/data_objects/construction_data_containers.py in __init__(self, level, left_edge, dims, fields, ds, num_ghost_zones, use_pbar, field_parameters)
    487                     self.ds.domain_left_edge)/self.dds).astype('int64')
    488         self._setup_data_source()
--> 489         self.get_data(fields)
    490 
    491     @property

/Users/jillnaiman1/yt-x86_64/lib/python3.4/site-packages/yt-3.2dev-py3.4-macosx-10.10-x86_64.egg/yt/data_objects/construction_data_containers.py in get_data(self, fields)
    570         fill, gen, part, alias = self._split_fields(fields_to_get)
    571         if len(part) > 0: self._fill_particles(part)
--> 572         if len(fill) > 0: self._fill_fields(fill)
    573         for a, f in sorted(alias.items()):
    574             self[a] = f(self)

/Users/jillnaiman1/yt-x86_64/lib/python3.4/site-packages/yt-3.2dev-py3.4-macosx-10.10-x86_64.egg/yt/data_objects/construction_data_containers.py in _fill_fields(self, fields)
    849                             chunk.ires, domain_dims, self.ds.refine_by)
    850             if level == 0 and tot != 0:
--> 851                 raise RuntimeError
    852             self._update_level_state(ls)
    853         for name, v in zip(fields, ls.fields):

RuntimeError: 
```


_______________________________________________
yt-dev mailing list
yt-dev at lists.spacepope.org
http://lists.spacepope.org/listinfo.cgi/yt-dev-spacepope.org