[yt-users] length of octree and smoothed arrays

[Desika Narayanan]

Hi Matt,

I just did yt update very recently (like an hour ago).  an hg id gives:

yt-hg desika$ hg id

ea182fc68fae (yt-3.0) tip





On Thu, Jun 12, 2014 at 2:28 PM, Matthew Turk <matthewturk at gmail.com> wrote:

> Hi Desika,
>
> I actually can't get your script to run.  Which changeset are you on?
>
> -Matt
>
> On Thu, Jun 12, 2014 at 12:13 PM, Desika Narayanan
> <dnarayan at haverford.edu> wrote:
> > Hi Matt,
> >
> > Thanks - the numbers are the same as the length of the numbers of ones in
> > saved['octree'], total sizes of octree and size of density array:
> >
> > saved["octree"].sum()  = 4681
> > saved["octree"].size  = 37499
> > ad["deposit", "all_density"].size = 262144
> >
> >
> >
> >
> > On Thu, Jun 12, 2014 at 2:06 PM, Matthew Turk <matthewturk at gmail.com>
> wrote:
> >>
> >> Hi Desika,
> >>
> >> On Thu, Jun 12, 2014 at 12:04 PM, Desika Narayanan
> >> <dnarayan at haverford.edu> wrote:
> >> > Hi All,
> >> >
> >> > First, I'm on an updated yt 3.0 updated about 5 minutes ago. (hg id:
> >> > ea182fc68fae (yt-3.0) tip)
> >> >
> >> > This might be related to Matt's update on octree on June 8th (in
> >> > response to
> >> > the CIC stuff) that I don't quite understand, but I've run into the
> >> > following issue highlighted in the example code:
> >> >
> >> > http://paste.yt-project.org/show/4781/
> >> >
> >> > Here, I create a fake data set (motivated by Nathan's shared notebook
> >> > earlier this week), construct the octree, and then smooth the density
> >> > onto a
> >> > grid with:
> >> >
> >> >>density = ad["deposit","all_density"]
> >> >
> >> > What I find is that the number of elements in the density array is not
> >> > the
> >> > same as the number of 0's (i.e Falses) in the generated octree
> (revealed
> >> > by
> >> >
> >> >>saved["octree"]
> >> >
> >> > in the above script).   Nominally, I thought the smoothing happened
> onto
> >> > the
> >> > same octree mesh that I generate with
> >> > ds.index.oct_handler.save_octree().
> >> > Is this not the case?  Or, is there a way to make it the case in the
> >> > order
> >> > that I call things?
> >>
> >> No, it should all be the same.  What're the values of:
> >>
> >> saved["octree"].sum()
> >> saved["octree"].size
> >> ad["deposit", "all_density"].size
> >>
> >> >
> >> > Thanks,
> >> > desika
> >> >
> >> >
> >> >
> >> > _______________________________________________
> >> > yt-users mailing list
> >> > yt-users at lists.spacepope.org
> >> > http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
> >> >
> >> _______________________________________________
> >> yt-users mailing list
> >> yt-users at lists.spacepope.org
> >> http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
> >
> >
> >
> > _______________________________________________
> > yt-users mailing list
> > yt-users at lists.spacepope.org
> > http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
> >
> _______________________________________________
> 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/20140612/26f711f2/attachment.html>