<div dir="ltr">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:<div>
<br></div><div><a href="https://bitbucket.org/yt_analysis/yt/pull-request/926/improved-compatibility-for-ytarray/diff">https://bitbucket.org/yt_analysis/yt/pull-request/926/improved-compatibility-for-ytarray/diff</a><br>
</div>
<div><br></div><div>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.</div><div><br></div><div>Thanks!</div><div><br></div><div>-Nathan</div></div>
<div class="gmail_extra"><br><br><div class="gmail_quote">On Mon, May 26, 2014 at 7:31 PM, Nathan Goldbaum <span dir="ltr"><<a href="mailto:nathan12343@gmail.com" target="_blank">nathan12343@gmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hey Brian,<div><br></div><div>Sorry for the trouble - this is due to a recent change to make YTArray subclasses behave more naturally.</div>
<div><br></div><div>I've managed to reproduce this locally using the following script:</div>
<div><br></div><div><a href="http://paste.yt-project.org/show/4684/" target="_blank">http://paste.yt-project.org/show/4684/</a><br></div><div><br></div><div>I will hopefully have a PR for you to test soon.</div><div><br>
</div><div>-Nathan</div>
</div><div class="gmail_extra"><br><br><div class="gmail_quote"><div><div class="h5">On Mon, May 26, 2014 at 6:22 PM, Brian O'Shea <span dir="ltr"><<a href="mailto:bwoshea@gmail.com" target="_blank">bwoshea@gmail.com</a>></span> wrote:<br>
</div></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div class="h5"><div dir="ltr">Hi folks,<div><br></div><div>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:</div>
<div><br>
</div><div><pre style="margin-top:0px;margin-bottom:0px;padding:5px 0px;font-family:'Bitstream Vera Sans Mono',monospace;font-size:13px;color:rgb(0,0,0)">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</pre><pre style="margin-top:0px;margin-bottom:0px;padding:5px 0px;font-family:'Bitstream Vera Sans Mono',monospace;font-size:13px;color:rgb(0,0,0)">
<br></pre></div><div>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?</div>
<div><br></div><div>Also, for the sake of completeness, you can see my batch and python scripts:</div><div><br></div><div> <a href="http://paste.yt-project.org/show/4681/" target="_blank">http://paste.yt-project.org/show/4681/</a><br>
</div>
<div> <a href="http://paste.yt-project.org/show/4682/" target="_blank">http://paste.yt-project.org/show/4682/</a><br></div><div><br></div><div>And the full output log from the run:</div><div><br></div><div><div> <a href="http://paste.yt-project.org/show/4683/" target="_blank">http://paste.yt-project.org/show/4683/</a><br>
</div></div><div><br></div><div>And, again, this is using the current tip of yt-3.</div><div><br></div><div>Thanks!</div><span><font color="#888888"><div><br></div><div>--Brian</div></font></span></div>
<br></div></div>_______________________________________________<br>
yt-dev mailing list<br>
<a href="mailto:yt-dev@lists.spacepope.org" target="_blank">yt-dev@lists.spacepope.org</a><br>
<a href="http://lists.spacepope.org/listinfo.cgi/yt-dev-spacepope.org" target="_blank">http://lists.spacepope.org/listinfo.cgi/yt-dev-spacepope.org</a><br>
<br></blockquote></div><br></div>
</blockquote></div><br></div>