<div dir="ltr">I'll take a look at analysis_modules.  I didn't want to bring in another dependency either, but looking at pyFFTW, it seems that its interface is identical to the NumPy one, so it could be an easy option for someone in the future.</div>
<div class="gmail_extra"><br><br><div class="gmail_quote">On Thu, Apr 17, 2014 at 7:15 AM, Matthew Turk <span dir="ltr"><<a href="mailto:matthewturk@gmail.com" target="_blank">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 Mike,<br>
<br>
On Wed, Apr 16, 2014 at 9:43 PM, Michael Zingale<br>
<<a href="mailto:michael.zingale@stonybrook.edu">michael.zingale@stonybrook.edu</a>> wrote:<br>
> Dear yters,<br>
<br>
Oh, don't say that like you're not one of us.  :-)<br>
<div class=""><br>
><br>
> I am interested in using yt to do power spectra for analyzing turbulence in<br>
> our simulations.  I quickly coded up a script to do this (warning: I haven't<br>
> verified this, so some of the normalization may be wrong, but it doesn't<br>
> assume a cubic domain, and it looked ok by eye):<br>
><br>
> <a href="http://paste.yt-project.org/show/4528/" target="_blank">http://paste.yt-project.org/show/4528/</a><br>
><br>
> The caveat is that we need to put this onto a uniform grid to do the FFT,<br>
> but then we can use the numpy FFT routines.<br>
<br>
</div>This is quite nice.<br>
<div class=""><br>
><br>
> I think that it would be nice to have some sort of power spectrum<br>
> functionality built into yt.  One thing that I don't know is what sort of<br>
> things do people in other fields look at?  For our work, we look at the<br>
> power spectrum of a density-weighted velocity field (rho**1/3 U for<br>
> stratified flows gives the expected Kolmogorov scaling).<br>
<br>
</div>I would strongly support this.  I know Jeff Oishi has looked a little<br>
bit at this, and I think Sam Skillman has lately been doing things<br>
with developing 2PCFs.<br>
<div class=""><br>
><br>
> The other thing is what sort of interface and where should this live in yt?<br>
> Ideally we'd compute some other properties (integral scales, etc.), but it<br>
> would be nice to have some discussion about what people use and then we can<br>
> start woth some basic functionality.<br>
<br>
</div>I think analysis_modules for now would be perfect.  And I tend to<br>
agree that we could do the optional dependency on FFTW, like Britton<br>
said, but I also seem to recall Jeff has a new set of FFTW bindings<br>
that work with MPI.<br>
<br>
-Matt<br>
<div class="im HOEnZb"><br>
><br>
> Mike<br>
><br>
> --<br>
> Michael Zingale<br>
> Associate Professor<br>
><br>
> Dept. of Physics & Astronomy • Stony Brook University • Stony Brook, NY<br>
> 11794-3800<br>
> phone:  <a href="tel:631-632-8225" value="+16316328225">631-632-8225</a><br>
> e-mail: <a href="mailto:Michael.Zingale@stonybrook.edu">Michael.Zingale@stonybrook.edu</a><br>
> web: <a href="http://www.astro.sunysb.edu/mzingale" target="_blank">http://www.astro.sunysb.edu/mzingale</a><br>
><br>
</div><div class="HOEnZb"><div class="h5">> _______________________________________________<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>
><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>
</div></div></blockquote></div><br><br clear="all"><div><br></div>-- <br><div>Michael Zingale</div><div>Associate Professor</div><div><br></div><div>Dept. of Physics & Astronomy • Stony Brook University • Stony Brook, NY 11794-3800</div>
<div><i>phone</i>:  631-632-8225</div><div><i>e-mail</i>: <a href="mailto:Michael.Zingale@stonybrook.edu" target="_blank">Michael.Zingale@stonybrook.edu</a></div><div><i>web</i>: <a href="http://www.astro.sunysb.edu/mzingale" target="_blank">http://www.astro.sunysb.edu/mzingale</a></div>

</div>