<div>Hi Stephen,</div><div><br></div>I've been comparing results of the halo mass function with the old ones on the small dataset, it seems the new fit function is systematically higher than the old ones.  I've also ran it on a 20 Mpc 800 cube unigrid, and the fit is also higher than the simulated data by almost 4x.  I'll send the pictures in another private email so I don't spam the list.  <div>
<br></div><div>Just want to let people know it no longer hangs due to not enough memory during calculation of the HMF.  Does anyone have a halo profile data at hand they can just check the halo mass function against the Warren fit in particular?  Calculating the halo mass function is quick if the haloes are already profiled and saved, took less than a minute on a single processor for the 800 cube.<br>
<div><br><div>From</div><div>G.S.</div><div><br><div class="gmail_quote">On Wed, Sep 21, 2011 at 8:00 PM, Stephen Skory <span dir="ltr"><<a href="mailto:s@skory.us">s@skory.us</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
Hi Geoffrey,<br>
<div class="im"><br>
> Thinking memory is the issue (mine or<br>
> my laptop's), I then tried to run the script on triton with the TCC node, it<br>
> still hangs and runs out of memory, I then tried the PDAF node and got the<br>
> same results even after updating to the latest dev (22ac529a1dab (yt) tip).<br>
<br>
</div>I see the same issue as you. I turned down the precision of the<br>
integrator in the halo mass function which should make these kind of<br>
memory overruns less likely. Please update your install and let me<br>
know if A) you still see problems or B) if the fit looks fishy (due to<br>
the lowering of the integrator's precision).<br>
<br>
Good luck!<br>
<br>
--<br>
<font color="#888888">Stephen Skory<br>
<a href="mailto:s@skory.us">s@skory.us</a><br>
<a href="http://stephenskory.com/" target="_blank">http://stephenskory.com/</a><br>
<a href="tel:510.621.3687" value="+15106213687">510.621.3687</a> (google voice)<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>
</font></blockquote></div><br></div></div></div>