Optimizing ideas

+ Currently we scan ALL elements for the one that contains one or more particles in an lapply loop.  Can we reduce the search effort by row-binding particles as we go? 

+ Can we reduce the search area by cropping the region around the point?  Like cropping the mesh?

+ if user selects a fixed z then skip retrieval of `ww`?
