<div dir="ltr">Hi Matt,<div><br></div><div>Thanks - it's the factor of 8 that was throwing me off. Once I account for that, the numbers make sense.</div><div><br></div><div>thanks</div><div>-desika</div></div><div class="gmail_extra">
<br><br><div class="gmail_quote">On Thu, Jun 12, 2014 at 2:38 PM, 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 Desika,<br>
<br>
Your particles are pretty randomly distributed, so it's refining<br>
everywhere up to a given level. Try changing to np.random.normal, but<br>
you might need to clip. There's also a factor of eight for leaf<br>
nodes, because a leaf node by default becomes eight cells.<br>
<br>
<a href="http://paste.yt-project.org/show/4782/" target="_blank">http://paste.yt-project.org/show/4782/</a><br>
<br>
This gives the results I was expecting. Do they look okay to you?<br>
<br>
-Matt<br>
<br>
On Thu, Jun 12, 2014 at 12:30 PM, Desika Narayanan<br>
<div class="HOEnZb"><div class="h5"><<a href="mailto:dnarayan@haverford.edu">dnarayan@haverford.edu</a>> wrote:<br>
> Hi Matt,<br>
><br>
> I just did yt update very recently (like an hour ago). an hg id gives:<br>
><br>
> yt-hg desika$ hg id<br>
><br>
> ea182fc68fae (yt-3.0) tip<br>
><br>
><br>
><br>
><br>
><br>
> On Thu, Jun 12, 2014 at 2:28 PM, Matthew Turk <<a href="mailto:matthewturk@gmail.com">matthewturk@gmail.com</a>> wrote:<br>
>><br>
>> Hi Desika,<br>
>><br>
>> I actually can't get your script to run. Which changeset are you on?<br>
>><br>
>> -Matt<br>
>><br>
>> On Thu, Jun 12, 2014 at 12:13 PM, Desika Narayanan<br>
>> <<a href="mailto:dnarayan@haverford.edu">dnarayan@haverford.edu</a>> wrote:<br>
>> > Hi Matt,<br>
>> ><br>
>> > Thanks - the numbers are the same as the length of the numbers of ones<br>
>> > in<br>
>> > saved['octree'], total sizes of octree and size of density array:<br>
>> ><br>
>> > saved["octree"].sum() = 4681<br>
>> > saved["octree"].size = 37499<br>
>> > ad["deposit", "all_density"].size = 262144<br>
>> ><br>
>> ><br>
>> ><br>
>> ><br>
>> > On Thu, Jun 12, 2014 at 2:06 PM, Matthew Turk <<a href="mailto:matthewturk@gmail.com">matthewturk@gmail.com</a>><br>
>> > wrote:<br>
>> >><br>
>> >> Hi Desika,<br>
>> >><br>
>> >> On Thu, Jun 12, 2014 at 12:04 PM, Desika Narayanan<br>
>> >> <<a href="mailto:dnarayan@haverford.edu">dnarayan@haverford.edu</a>> wrote:<br>
>> >> > Hi All,<br>
>> >> ><br>
>> >> > First, I'm on an updated yt 3.0 updated about 5 minutes ago. (hg id:<br>
>> >> > ea182fc68fae (yt-3.0) tip)<br>
>> >> ><br>
>> >> > This might be related to Matt's update on octree on June 8th (in<br>
>> >> > response to<br>
>> >> > the CIC stuff) that I don't quite understand, but I've run into the<br>
>> >> > following issue highlighted in the example code:<br>
>> >> ><br>
>> >> > <a href="http://paste.yt-project.org/show/4781/" target="_blank">http://paste.yt-project.org/show/4781/</a><br>
>> >> ><br>
>> >> > Here, I create a fake data set (motivated by Nathan's shared notebook<br>
>> >> > earlier this week), construct the octree, and then smooth the density<br>
>> >> > onto a<br>
>> >> > grid with:<br>
>> >> ><br>
>> >> >>density = ad["deposit","all_density"]<br>
>> >> ><br>
>> >> > What I find is that the number of elements in the density array is<br>
>> >> > not<br>
>> >> > the<br>
>> >> > same as the number of 0's (i.e Falses) in the generated octree<br>
>> >> > (revealed<br>
>> >> > by<br>
>> >> ><br>
>> >> >>saved["octree"]<br>
>> >> ><br>
>> >> > in the above script). Nominally, I thought the smoothing happened<br>
>> >> > onto<br>
>> >> > the<br>
>> >> > same octree mesh that I generate with<br>
>> >> > ds.index.oct_handler.save_octree().<br>
>> >> > Is this not the case? Or, is there a way to make it the case in the<br>
>> >> > order<br>
>> >> > that I call things?<br>
>> >><br>
>> >> No, it should all be the same. What're the values of:<br>
>> >><br>
>> >> saved["octree"].sum()<br>
>> >> saved["octree"].size<br>
>> >> ad["deposit", "all_density"].size<br>
>> >><br>
>> >> ><br>
>> >> > Thanks,<br>
>> >> > desika<br>
>> >> ><br>
>> >> ><br>
>> >> ><br>
>> >> > _______________________________________________<br>
>> >> > yt-users mailing list<br>
>> >> > <a href="mailto:yt-users@lists.spacepope.org">yt-users@lists.spacepope.org</a><br>
>> >> > <a href="http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org" target="_blank">http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org</a><br>
>> >> ><br>
>> >> _______________________________________________<br>
>> >> yt-users mailing list<br>
>> >> <a href="mailto:yt-users@lists.spacepope.org">yt-users@lists.spacepope.org</a><br>
>> >> <a href="http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org" target="_blank">http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org</a><br>
>> ><br>
>> ><br>
>> ><br>
>> > _______________________________________________<br>
>> > yt-users mailing list<br>
>> > <a href="mailto:yt-users@lists.spacepope.org">yt-users@lists.spacepope.org</a><br>
>> > <a href="http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org" target="_blank">http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org</a><br>
>> ><br>
>> _______________________________________________<br>
>> yt-users mailing list<br>
>> <a href="mailto:yt-users@lists.spacepope.org">yt-users@lists.spacepope.org</a><br>
>> <a href="http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org" target="_blank">http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org</a><br>
><br>
><br>
><br>
> _______________________________________________<br>
> yt-users mailing list<br>
> <a href="mailto:yt-users@lists.spacepope.org">yt-users@lists.spacepope.org</a><br>
> <a href="http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org" target="_blank">http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org</a><br>
><br>
_______________________________________________<br>
yt-users mailing list<br>
<a href="mailto:yt-users@lists.spacepope.org">yt-users@lists.spacepope.org</a><br>
<a href="http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org" target="_blank">http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org</a><br>
</div></div></blockquote></div><br></div>