[yt-users] Streamlines error
Mark Krumholz
mkrumhol at ucsc.edu
Wed Jul 1 06:01:33 PDT 2015
Hi Caroline,
I can’t help with the crash, but the huge slowdown in streamlining in 3.x versus 2.x is a known issue <https://bitbucket.org/yt_analysis/yt/pull-request/1514/speed-up-amrkdtree-if-grids-exist/diff>, which is not yet been fixed unfortunately.
> On Jul 1, 2015, at 8:35 AM, Caroline Van Borm <borm at astro.rug.nl> wrote:
>
> Hi yt users,
>
> I'm having a problem generating streamlines. I'm following the cookbook script and it works fine for the IsolatedGalaxy dataset, but not for my own data (also Enzo). The error I'm getting is the following:
>
> ---------------------------------------------------------------------------
> AttributeError Traceback (most recent call last)
> <ipython-input-20-d3a0fa7307d4> in <module>()
> 5 pos = c+pos_dx
> 6
> ----> 7 streamlines = Streamlines(pf,pos,'x-velocity', 'y-velocity', 'z-velocity', length=1.0
> )
> 8 #streamlines = Streamlines(pf,pos,'x-velocity', 'y-velocity', 'z-velocity', length=pf.domain_width[0])
> 9 streamlines.integrate_through_volume()
>
> /data/users/borm/YT/yt-x86_64/src/yt-hg/yt/visualization/streamlines.pyc in __init__(self, pf, positions, xfield, yfield, zfield, volume, dx, length, direction, get_magnitude)
> 102 volume = AMRKDTree(self.pf, fields=[self.xfield,self.yfield,self.zfield],
> 103 log_fields=[False,False,False])
> --> 104 volume.join_parallel_trees()
> 105 self.volume = volume
> 106 if dx is None:
>
> /data/users/borm/YT/yt-x86_64/src/yt-hg/yt/utilities/amr_kdtree/amr_kdtree.pyc in join_parallel_trees(self)
> 456 nid = np.array(nid)
> 457 new_tree = self.rebuild_tree_from_array(nid, pid, lid,
> --> 458 rid, les, res, gid, splitdims, splitposs)
> 459
> 460 def get_node_arrays(self):
>
> /data/users/borm/YT/yt-x86_64/src/yt-hg/yt/utilities/amr_kdtree/amr_kdtree.pyc in rebuild_tree_from_array(self, nids, pids, lids, rids, les, res, gids, splitdims, splitposs)
> 505 N = nids.shape[0]
> 506 for i in xrange(N):
> --> 507 n = self.get_node(nids[i])
> 508 n.set_left_edge(les[i])
> 509 n.set_right_edge(res[i])
>
> /data/users/borm/YT/yt-x86_64/src/yt-hg/yt/utilities/amr_kdtree/amr_kdtree.pyc in get_node(self, nodeid)
> 212 for depth in range(1,len(path)):
> 213 if path[depth] == '0':
> --> 214 temp = temp.left
> 215 else:
> 216 temp = temp.right
>
> AttributeError: 'NoneType' object has no attribute 'left'
> Though sometimes with a different dataset I also get
> AttributeError: 'NoneType' object has no attribute 'right'
> I'm using yt 2.6, but I've also tried with yt 3.1 and I get the same error there (only it takes much, much longer before it crashes?) My datasets do not have particles, not sure if that may be relevant.
>
> Any help is much appreciated!
>
> Cheers,
> Caroline
>
>
>
> _______________________________________________
> yt-users mailing list
> yt-users at lists.spacepope.org
> http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
--
Mark Krumholz
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.spacepope.org/pipermail/yt-users-spacepope.org/attachments/20150701/9c3a21f3/attachment.html>
More information about the yt-users
mailing list