[yt-users] length of octree and smoothed arrays

Desika Narayanan dnarayan at haverford.edu
Thu Jun 12 11:13:52 PDT 2014


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
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.spacepope.org/pipermail/yt-users-spacepope.org/attachments/20140612/9bd3f9df/attachment.html>


More information about the yt-users mailing list