[yt-dev] Halo Catalog unit error

[Rasmi Elasmar]

Hi all,

I'm running into more issues with Halo Catalogs and units. I'm not sure if
the last PR <https://bitbucket.org/yt_analysis/yt/pull-requests/2208/> is
causing this or not. I'm running the latest code from the repo as of right
now (fd8796c8e06d). Here's the error, generated on this data
<https://drive.google.com/open?id=0BwK-7Z3S5X_yMDNyZ1RTcG56SFU> with this
script <https://gist.github.com/rasmi/6cf6017b50277390735dd484c120855f>.

/work/03330/tg826294/applications/pythonenv/lib/python2.7/site-packages/yt-3.3.dev0-py2.7-linux-x86_64.egg/yt/units/yt_array.py:976:
> RuntimeWarning: divide by zero encountered in divide
>   return super(YTArray, self).__div__(ro)
> Traceback (most recent call last):
>   File "/work/03330/tg826294/applications/scripts/findhalos.py", line 65,
> in <module>
>     hc.create()
>   File
> "/work/03330/tg826294/applications/pythonenv/lib/python2.7/site-packages/yt-3.3.dev0-py2.7-linux-x86_64.egg/yt/analysis_modules/halo_analysis/halo_catalog.py",
> line 335, in create
>     self._run(save_halos, save_catalog, njobs=njobs, dynamic=dynamic)
>   File
> "/work/03330/tg826294/applications/pythonenv/lib/python2.7/site-packages/yt-3.3.dev0-py2.7-linux-x86_64.egg/yt/utilities/parallel_tools/parallel_analysis_interface.py",
> line 302, in barrierize
>     return func(*args, **kwargs)
>   File
> "/work/03330/tg826294/applications/pythonenv/lib/python2.7/site-packages/yt-3.3.dev0-py2.7-linux-x86_64.egg/yt/analysis_modules/halo_analysis/halo_catalog.py",
> line 427, in _run
>     action(new_halo)
>   File
> "/work/03330/tg826294/applications/pythonenv/lib/python2.7/site-packages/yt-3.3.dev0-py2.7-linux-x86_64.egg/yt/analysis_modules/halo_analysis/halo_callbacks.py",
> line 60, in __call__
>     self.function(halo, *self.args, **self.kwargs)
>   File
> "/work/03330/tg826294/applications/pythonenv/lib/python2.7/site-packages/yt-3.3.dev0-py2.7-linux-x86_64.egg/yt/analysis_modules/halo_analysis/halo_callbacks.py",
> line 571, in iterative_center_of_mass
>     sphere = halo.halo_catalog.data_ds.sphere(center_orig,
> halo.quantities[radius_field])
>   File
> "/work/03330/tg826294/applications/pythonenv/lib/python2.7/site-packages/yt-3.3.dev0-py2.7-linux-x86_64.egg/yt/data_objects/selection_data_containers.py",
> line 649, in __init__
>     if radius < self.index.get_smallest_dx():
>   File
> "/work/03330/tg826294/applications/pythonenv/lib/python2.7/site-packages/yt-3.3.dev0-py2.7-linux-x86_64.egg/yt/units/yt_array.py",
> line 1095, in __lt__
>     return super(YTArray, self).__lt__(oth)
>   File
> "/work/03330/tg826294/applications/pythonenv/lib/python2.7/site-packages/yt-3.3.dev0-py2.7-linux-x86_64.egg/yt/units/yt_array.py",
> line 1225, in __array_wrap__
>     raise YTUfuncUnitError(context[0], unit1, unit2)
> yt.utilities.exceptions.YTUfuncUnitError: The NumPy <ufunc 'less'>
> operation is only allowed on objects with identical units. Convert one of
> the arrays to the other's units first. Received units (code_length) and
> (code_length).


Has anyone run into this before? yt seems to think these two units aren't
the same -- is it possible the HaloCatalog unit import is being done
incorrectly? At this point, I haven't written anything to the disk, so I'm
not sure what the issue might be.


Thanks,

Rasmi
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.spacepope.org/pipermail/yt-dev-spacepope.org/attachments/20160612/0ab4f3cd/attachment.htm>