Hi Matt,<div><br></div><div>Looks like some good speedup:</div><div><div>Old: 5.53041e+01</div><div>New: 2.57021e+01</div><div><br></div><div>This was on a 6 level AMR everywhere.  I will note that as you mentioned, the disk cache does do a lot, as here was the results from the first run:</div>

<div><div>Old: 5.51695e+01</div><div>New: 2.23177e+02</div></div><div><br></div><div>Cheers,</div><div>Sam</div><div><br></div><br><div class="gmail_quote">On Fri, Aug 27, 2010 at 1:06 PM, Matthew Turk <span dir="ltr"><<a href="mailto:matthewturk@gmail.com">matthewturk@gmail.com</a>></span> wrote:<br>

<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">Hi everyone,<br>
<br>
I added the quad tree projection to the tip of the hg repository.<br>
<br>
If you could please put it through a few paces, that'd be great.  You<br>
can access it as the attribute "quad_proj" on a hierarchy.  I've done<br>
some speed testing and the generality necessary to get it to work with<br>
sources and whatnot makes it a bit slower than if it was just fed the<br>
data, but it's still about twice as fast on my datasets as the old<br>
projection method.<br>
<br>
It doesn't work in parallel yet, but give it a shot on some big<br>
datasets anyway -- I tested it on a 1024^3 L7 and it worked fine in<br>
serial.  But, I know how to parallelize it, and the parallelization<br>
process should also speed up the serial projections.  (It will allow<br>
us to better choose the batching of IO.)  I think this projection<br>
method should also scale *much* better for parallelism, because we<br>
shouldn't have to do domain decomposition to parallelize, which was<br>
always a killer on nested grid simulations.<br>
<br>
I'd like to push forward with this, so if you get a chance to try it,<br>
report back with any issues.  I'm sure there will be several...<br>
<br>
I've placed a comparison test here:<br>
<br>
<a href="http://paste.enzotools.org/show/1134/" target="_blank">http://paste.enzotools.org/show/1134/</a><br>
<br>
(But note that disk cache will favor the old projection method the<br>
first time you run, so run at least twice.  :)<br>
<br>
-Matt<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><br clear="all"><br>-- <br>Samuel W. Skillman<br>DOE Computational Science Graduate Fellow<br>Center for Astrophysics and Space Astronomy<br>University of Colorado at Boulder<br>samuel.skillman[at]<a href="http://colorado.edu">colorado.edu</a><br>


</div>