[yt-dev] Gridding Particles

[Christopher Moody]

Hi everyone,
I just asked the IRC list this, but I guess no one was there. I'm trying to
grid a few million particles onto the yt grids, but checking if the left
and right edges of my particles are inside every cell is slow:

#find particles inside cell
idx = na.logical_and(na.all(le < pos,axis=1), na.all(re > pos,axis=1))
#assign them to a grid
...
#remove particles from future searches
pos = pos[~idx]

I start with the finest grids at the highest levels, and once the particles
are gridded they're removed -- so the list of particles is always
decreasing, and the coarsest mesh tends to have only a few particles. I'm
guessing that somewhere this has been done before - I'd rather not reinvent
the wheel. Has anyone used any relevant kdtree (or cython) code in the yt
utilities that can select a subvolume of particles quickly?

Thanks ahead of time!
chris
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.spacepope.org/pipermail/yt-dev-spacepope.org/attachments/20120415/43efebf1/attachment.htm>