<p dir="ltr"><br>
On Oct 7, 2014 10:41 PM, "Nathan Goldbaum" <<a href="mailto:nathan12343@gmail.com">nathan12343@gmail.com</a>> wrote:<br>
><br>
><br>
><br>
> On Tue, Oct 7, 2014 at 7:45 PM, Matthew Turk <<a href="mailto:matthewturk@gmail.com">matthewturk@gmail.com</a>> wrote:<br>
>><br>
>> On Tue, Oct 7, 2014 at 7:55 PM, Andrew Myers <<a href="mailto:atmyers2@gmail.com">atmyers2@gmail.com</a>> wrote:<br>
>> > Hi Folks,<br>
>> ><br>
>> > I'm trying to understand how the particle deposition fields work for<br>
>> > grid-based codes. Specifically, how does yt handle particles that are close<br>
>> > to the grid boundaries, so that the CIC interpolation stencil involves cells<br>
>> > on grid(s) other than the one that contains the particle position?. Like,<br>
>> > say I did:<br>
>> ><br>
>> > import yt<br>
>> > ds = yt.load(fn)<br>
>> > g = ds.index.grids[0]<br>
>> > print g['deposit', 'all_cic']<br>
>> ><br>
>> > How does yt handle the fact that particles other than those living on "g"<br>
>> > can contribute to the total CIC density? Does it grab a layer of "ghost<br>
>> > particles" somehow? I can't seem to find where that happens.<br>
>> ><br>
>> > Sorry if I'm just missing something obvious here...<br>
><br>
><br>
> Andrew's original e-mail prompted me to make the following notebook, which might prove useful for constructing a unit test:<br>
><br>
> <a href="http://nbviewer.ipython.org/gist/ngoldbaum/eb40dfe74cd58a4e16ba">http://nbviewer.ipython.org/gist/ngoldbaum/eb40dfe74cd58a4e16ba</a><br>
>  <br>
>><br>
>><br>
>> It doesn't, and you're not.  This is something it should do, but does<br>
>> not currently.  If you're using a "smoothing" operation, it will, as<br>
>> then it handles non-local particles, but deposit operations don't do<br>
>> ghost zones.<br>
><br>
><br>
> Smoothing is also only implemented for octrees.<br>
>  <br>
>><br>
>><br>
>> This is something we should fix, and that we probably could by simply<br>
>> making all the CIC fields ValidateSpatial(1).  Can you give that a<br>
>> shot?<br>
><br>
><br>
> Will this have interactions with <a href="https://trello.com/c/SFuSSoIo">https://trello.com/c/SFuSSoIo</a> ?<br>
><br>
> I think the answer is no, since we're querying particle fields and selection doesn't require cascading interpolation, but I just want to double-check.</p>
<p dir="ltr">It should be pretty fast.</p>
<p dir="ltr">><br>
>><br>
>> ><br>
>> > Best,<br>
>> > Andrew<br>
>> ><br>
>> > _______________________________________________<br>
>> > yt-dev mailing list<br>
>> > <a href="mailto:yt-dev@lists.spacepope.org">yt-dev@lists.spacepope.org</a><br>
>> > <a href="http://lists.spacepope.org/listinfo.cgi/yt-dev-spacepope.org">http://lists.spacepope.org/listinfo.cgi/yt-dev-spacepope.org</a><br>
>> ><br>
>> _______________________________________________<br>
>> yt-dev mailing list<br>
>> <a href="mailto:yt-dev@lists.spacepope.org">yt-dev@lists.spacepope.org</a><br>
>> <a href="http://lists.spacepope.org/listinfo.cgi/yt-dev-spacepope.org">http://lists.spacepope.org/listinfo.cgi/yt-dev-spacepope.org</a><br>
><br>
><br>
><br>
> _______________________________________________<br>
> yt-dev mailing list<br>
> <a href="mailto:yt-dev@lists.spacepope.org">yt-dev@lists.spacepope.org</a><br>
> <a href="http://lists.spacepope.org/listinfo.cgi/yt-dev-spacepope.org">http://lists.spacepope.org/listinfo.cgi/yt-dev-spacepope.org</a><br>
><br>
</p>