If you check out the IsBound function in derived quantities, the work is done by a c function in PointCombine.c that accepts a list of x, y, z, and mass. It should be simple to write something to throw particle lists instead of cell lists. It's not the fastest thing in the world, as there are (N^2 - N)/2 operations. Matt has also implemented a version that uses the CUDA libraries and works on the gpu and is significantly faster.<br>
<br>I haven't seen any objections to removing the extra factor of 2 in front of the potential energy. I'm going to do that now.<br><br>Britton<br><br><div class="gmail_quote">On Sat, Jan 17, 2009 at 9:59 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;"><div class="Ih2E3d">> I (reluctantly?) agree with Britton that there's an error. He's going<br>
to fix it, unless somebody speaks now. (Else, holding peace forever is<br>
expected.)<br>
<br>
</div>I give up none of my rights!<br>
<br>
IsBound interests me. What would be the best way to make this work with HOP haloes? Put a kind of wrapper in the HopGroup class that cuts out a sphere centered on .center_of_mass() with radius .maximum_radius()? I'd prefer to only operate on the DM particles, not the gas mass. Any suggestions?<br>
<br>
Thanks!<br>
<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>
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>
</blockquote></div><br>