[yt-users] Variance in phase/profiles

[Britton Smith]

Hi Dave,

The best I have been able to do for this is to manually bin up a 3d data
object with the cut_region function and then do the calculation myself bin
by bin.
For example, if you have a sphere, you can do
new_sphere = sphere.cut_region(['grid["Temperature"] > 1e5',
'grid["Temperature"] < 1e6'])
That'll give you access to all the cells within that bin for a variance
calculation.  Note the quotes around those expressions.  The cuts are
applied in yt using eval functions.

This will work, but it's slow and requires you to do a lot by hand.

Britton

On Wed, Oct 5, 2011 at 11:51 AM, david collins <antpuncher at gmail.com> wrote:

> Hi, Everybody!
>
> Does anyone out there have a technique for getting the variance out of
> a profile object?  A profile object is good at getting <X> vs. B, I'd
> then like to get < (X - <X>)^2 > vs B.  Matt and I had spittballed the
> possibility some time ago, but I was wondering if anyone out there had
> successfully done it.
>
> Thanks,
> d.
>
> --
> Sent from my computer.
> _______________________________________________
> yt-users mailing list
> yt-users at lists.spacepope.org
> http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.spacepope.org/pipermail/yt-users-spacepope.org/attachments/20111005/810d9d96/attachment.html>