<div dir="ltr">Hi Matt,<div><br></div><div>Just making all the CIC fields ValidateSpatial(1) doesn't quite work - it looks like it kicks off an endless recursion, where in order to generate ('deposit', 'io_cic'), you have to generate ('deposit', 'io_cic'), etc, etc. I'll try to look into it some more.</div><div><br></div><div>Nathan, thanks for that notebook - that will be very useful. I can't quite get it to work, <a href="http://nbviewer.ipython.org/gist/atmyers/f1e31ccbe18268bfb466">though</a>. Ever see that before?</div><div><br></div><div>-Andrew</div></div><div class="gmail_extra"><br><div class="gmail_quote">On Tue, Oct 7, 2014 at 7:45 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"><div class="HOEnZb"><div class="h5">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>
</div></div>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>
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>
> 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" target="_blank">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" target="_blank">http://lists.spacepope.org/listinfo.cgi/yt-dev-spacepope.org</a><br>
</blockquote></div><br></div>