Hi Stephen,<br><br>I'm having a problem with parallelHOP where one (and seemingly only one) of the processors is having an array out of bounds error with the densest_in_chain_real_index array. I printed out the size of this array and the densest_in_chain array and in the one processor where it fails the densest_in_chain_real_index array has one less element. Below is the traceback:<br>
<br>The three numbers printed out at the top are the sizes of densest_in_chain, densest_in_chain_real_index, and sort, in that order.<br><br>P005 yt INFO 2009-12-01 10:12:28,438 Locally sorting chains...<br>
3794<br>3793<br>3794<br>Traceback (most recent call last):<br> File "do_hop.py", line 9, in <module><br> fancy_padding=True, rearrange=True)<br> File "/duanestorage/home/student/brittons/local/x86_64/lib/python2.6/site-packages/yt-1.6dev-py2.6-linux-x86_64.egg/yt/lagos/HaloFinding.py", line 958, in __init__<br>
threshold=threshold, dm_only=dm_only, rearrange=rearrange)<br> File "/duanestorage/home/student/brittons/local/x86_64/lib/python2.6/site-packages/yt-1.6dev-py2.6-linux-x86_64.egg/yt/lagos/HaloFinding.py", line 598, in __init__<br>
HaloList.__init__(self, data_source, dm_only)<br> File "/duanestorage/home/student/brittons/local/x86_64/lib/python2.6/site-packages/yt-1.6dev-py2.6-linux-x86_64.egg/yt/lagos/HaloFinding.py", line 339, in __init__<br>
self._run_finder()<br> File "/duanestorage/home/student/brittons/local/x86_64/lib/python2.6/site-packages/yt-1.6dev-py2.6-linux-x86_64.egg/yt/lagos/HaloFinding.py", line 609, in _run_finder<br> self.threshold, rearrange=self.rearrange)<br>
File "/duanestorage/home/student/brittons/local/x86_64/lib/python2.6/site-packages/yt-1.6dev-py2.6-linux-x86_64.egg/yt/lagos/parallelHOP/parallelHOP.py",line 56, in __init__<br> self._chain_hop()<br> File "/duanestorage/home/student/brittons/local/x86_64/lib/python2.6/site-packages/yt-1.6dev-py2.6-linux-x86_64.egg/yt/lagos/parallelHOP/parallelHOP.py",line 1423, in _chain_hop<br>
chain_count = self._preconnect_chains(chain_count)<br> File "/duanestorage/home/student/brittons/local/x86_64/lib/python2.6/site-packages/yt-1.6dev-py2.6-linux-x86_64.egg/yt/lagos/parallelHOP/parallelHOP.py",line 538, in _preconnect_chains<br>
self.densest_in_chain_real_index = self.densest_in_chain_real_index[sort]<br>IndexError: index 3793 out of bounds 0<=index<3793<br><br>Do you have any idea what's going on here?<br><br>Thanks,<br>Britton<br>
<br><div class="gmail_quote">On Wed, Nov 25, 2009 at 4:45 PM, Stephen Skory <span dir="ltr"><<a href="mailto:stephenskory@yahoo.com">stephenskory@yahoo.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
J. S.,<br>
<div class="im"><br>
<br>
> > 2. Install Forthon in the standard python way.<br>
> > <a href="http://hifweb.lbl.gov/Forthon/" target="_blank">http://hifweb.lbl.gov/Forthon/</a><br>
> Can I ask what made you use Forthon over f2py? I have been using f2py,<br>
> and have found it quite useful for wrapping fortran code. Given that<br>
> f2py comes with numpy, I'm wondering what the specific advantage for<br>
> Forthon is, if any.<br>
<br>
</div>I was unable to convince f2py to build static objects, which is required on Kraken. That is it's big advantage.<br>
<div class="im"><br>
_______________________________________________________<br>
<a href="mailto:sskory@physics.ucsd.edu">sskory@physics.ucsd.edu</a> o__ Stephen Skory<br>
<a href="http://physics.ucsd.edu/%7Esskory/" target="_blank">http://physics.ucsd.edu/~sskory/</a> _.>/ _Graduate Student<br>
________________________________(_)_\(_)_______________<br>
_______________________________________________<br>
</div><div><div></div><div class="h5">Yt-dev mailing list<br>
<a href="mailto:Yt-dev@lists.spacepope.org">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>
</div></div></blockquote></div><br>