[yt-dev] power spectra in yt

Michael Zingale michael.zingale at stonybrook.edu
Thu Apr 17 05:11:02 PDT 2014


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.


On Thu, Apr 17, 2014 at 7:15 AM, Matthew Turk <matthewturk at gmail.com> wrote:

> Hi Mike,
>
> On Wed, Apr 16, 2014 at 9:43 PM, Michael Zingale
> <michael.zingale at stonybrook.edu> wrote:
> > Dear yters,
>
> Oh, don't say that like you're not one of us.  :-)
>
> >
> > I am interested in using yt to do power spectra for analyzing turbulence
> in
> > our simulations.  I quickly coded up a script to do this (warning: I
> haven't
> > verified this, so some of the normalization may be wrong, but it doesn't
> > assume a cubic domain, and it looked ok by eye):
> >
> > http://paste.yt-project.org/show/4528/
> >
> > The caveat is that we need to put this onto a uniform grid to do the FFT,
> > but then we can use the numpy FFT routines.
>
> This is quite nice.
>
> >
> > I think that it would be nice to have some sort of power spectrum
> > functionality built into yt.  One thing that I don't know is what sort of
> > things do people in other fields look at?  For our work, we look at the
> > power spectrum of a density-weighted velocity field (rho**1/3 U for
> > stratified flows gives the expected Kolmogorov scaling).
>
> I would strongly support this.  I know Jeff Oishi has looked a little
> bit at this, and I think Sam Skillman has lately been doing things
> with developing 2PCFs.
>
> >
> > The other thing is what sort of interface and where should this live in
> yt?
> > Ideally we'd compute some other properties (integral scales, etc.), but
> it
> > would be nice to have some discussion about what people use and then we
> can
> > start woth some basic functionality.
>
> I think analysis_modules for now would be perfect.  And I tend to
> agree that we could do the optional dependency on FFTW, like Britton
> said, but I also seem to recall Jeff has a new set of FFTW bindings
> that work with MPI.
>
> -Matt
>
> >
> > Mike
> >
> > --
> > Michael Zingale
> > Associate Professor
> >
> > Dept. of Physics & Astronomy * Stony Brook University * Stony Brook, NY
> > 11794-3800
> > phone:  631-632-8225
> > e-mail: Michael.Zingale at stonybrook.edu
> > web: http://www.astro.sunysb.edu/mzingale
> >
> > _______________________________________________
> > yt-dev mailing list
> > yt-dev at lists.spacepope.org
> > http://lists.spacepope.org/listinfo.cgi/yt-dev-spacepope.org
> >
> _______________________________________________
> yt-dev mailing list
> yt-dev at lists.spacepope.org
> http://lists.spacepope.org/listinfo.cgi/yt-dev-spacepope.org
>



-- 
Michael Zingale
Associate Professor

Dept. of Physics & Astronomy * Stony Brook University * Stony Brook, NY
11794-3800
*phone*:  631-632-8225
*e-mail*: Michael.Zingale at stonybrook.edu
*web*: http://www.astro.sunysb.edu/mzingale
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.spacepope.org/pipermail/yt-dev-spacepope.org/attachments/20140417/ec52dded/attachment.html>


More information about the yt-dev mailing list