[yt-dev] Hop in yt-3

Nathan Goldbaum nathan12343 at gmail.com
Mon May 26 23:55:54 PDT 2014


It ended up being due to the halo finding code using a feature of ndarray I
didn't know existed and therefore wasn't reflected in the design of
YTArray.  This defect should hopefully be fixed (and tested!) here:

https://bitbucket.org/yt_analysis/yt/pull-request/926/improved-compatibility-for-ytarray/diff

If you could try that script after applying the PR and then leave a comment
saying it it worked or not, I'd appreciate it.

Thanks!

-Nathan


On Mon, May 26, 2014 at 7:31 PM, Nathan Goldbaum <nathan12343 at gmail.com>wrote:

> Hey Brian,
>
> Sorry for the trouble - this is due to a recent change to make YTArray
> subclasses behave more naturally.
>
> I've managed to reproduce this locally using the following script:
>
> http://paste.yt-project.org/show/4684/
>
> I will hopefully have a PR for you to test soon.
>
> -Nathan
>
>
> On Mon, May 26, 2014 at 6:22 PM, Brian O'Shea <bwoshea at gmail.com> wrote:
>
>> Hi folks,
>>
>> I just tried to do some halo-finding with yt-3 using Hop (the current
>> tip; rev 3cfb38afe34d), and it died somewhere after calling 'regroup' with
>> this error:
>>
>> Traceback (most recent call last):
>>   File "hop.py", line 11, in <module>
>>     halo_list = HaloFinder(pf, subvolume=my_subregion)
>>   File "/mnt/a/u/sciteam/bwoshea/yt-3.0/src/yt-hg/yt/analysis_modules/halo_finding/halo_objects.py", line 2440, in __init__
>>     self._join_halolists()
>>   File "/mnt/a/u/sciteam/bwoshea/yt-3.0/src/yt-hg/yt/analysis_modules/halo_finding/halo_objects.py", line 1889, in _join_halolists
>>     self._groups.sort(haloCmp)
>>   File "/mnt/a/u/sciteam/bwoshea/yt-3.0/src/yt-hg/yt/analysis_modules/halo_finding/halo_objects.py", line 1888, in haloCmp
>>     return cmp(h1.center_of_mass()[0], h2.center_of_mass()[0])
>>   File "/mnt/a/u/sciteam/bwoshea/yt-3.0/src/yt-hg/yt/utilities/parallel_tools/parallel_analysis_interface.py", line 191, in single_proc_results
>>     retval = func(self, *args, **kwargs)
>>   File "/mnt/a/u/sciteam/bwoshea/yt-3.0/src/yt-hg/yt/analysis_modules/halo_finding/halo_objects.py", line 158, in center_of_mass
>>     c = (com * pm).sum(axis=1) / pm.sum()
>>   File "/mnt/a/u/sciteam/bwoshea/yt-3.0/src/yt-hg/yt/units/yt_array.py", line 648, in __rmul__
>>     return YTArray(super(YTArray, self).__rmul__(lo))
>>   File "/mnt/a/u/sciteam/bwoshea/yt-3.0/src/yt-hg/yt/units/yt_array.py", line 910, in __array_wrap__
>>     ret_class = get_binary_op_return_class(cls1, cls2)
>>   File "/mnt/a/u/sciteam/bwoshea/yt-3.0/src/yt-hg/yt/units/yt_array.py", line 1045, in get_binary_op_return_class
>>     raise RuntimeError("Operations are only defined on pairs of objects"
>> RuntimeError: Operations are only defined on pairs of objectsin which one is a subclass of the other
>>
>>
>> It was functional roughly two weeks ago (on May 10), when I successfully
>> did halo-finding on the same machine with yt-3.  Looking at the errors, it
>> seems to be something more general than Hop functionality, but I'm not
>> entirely sure.  Does anybody have any thoughts on this?
>>
>> Also, for the sake of completeness, you can see my batch and python
>> scripts:
>>
>>   http://paste.yt-project.org/show/4681/
>>    http://paste.yt-project.org/show/4682/
>>
>> And the full output log from the run:
>>
>>   http://paste.yt-project.org/show/4683/
>>
>> And, again, this is using the current tip of yt-3.
>>
>> Thanks!
>>
>> --Brian
>>
>> _______________________________________________
>> yt-dev mailing list
>> yt-dev at lists.spacepope.org
>> http://lists.spacepope.org/listinfo.cgi/yt-dev-spacepope.org
>>
>>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.spacepope.org/pipermail/yt-dev-spacepope.org/attachments/20140526/247fb878/attachment.html>


More information about the yt-dev mailing list