<div dir="ltr"><div>Hi Cameron,</div><div><br></div>No, and I just checked, it's not in my draft of the 3.0 paper yet.  It might be in a YTEP somewhere, but here's the outline.<div><br></div><div>Here's what it does:</div><div><br></div><div> * Build "mesh" octree from all the particles in a simulation</div><div> * When requested, a "particle" octree is created for particle searching</div><div> * Each cell in the mesh octree is iterated over; neighbor leaves in the "particle" octree are identified, and the N nearest particles are found.  (This is a "gather" operation, but because the particle octree is not balanced, it is not guaranteed to find enough neighbors.)</div><div> * Using the chosen SPH kernel (thanks to Bili Dong, we now have cubic, quartic, quintic, wendland2, wendland4, and wendland6), the smoothing is conducted on each cell by looking at all the nearest particles.  This is implemented in the process function on the class VolumeWeightedSmooth in yt/geometry/particle_smooth.pyx</div><div><br></div><div>I should note that Meagan Lang has implemented a particle kD-tree that she's integrating to replace the "particle" octree, which will then ensure a balanced search.  She's verified that when using this, there are no longer any zero values, and it's also about the same speed, which is pretty rad.</div><div><br></div><div>-Matt</div></div><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Sep 8, 2016 at 12:42 PM, Cameron Hummels <span dir="ltr"><<a href="mailto:chummels@gmail.com" target="_blank">chummels@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 dir="ltr">Hey everyone,<div><br></div><div>Britton, Devin, and I are writing up the Trident method paper, and we wanted to be explicit about how exactly yt deposits SPH field quantities on to an octree for particle-based fields.  Are there any formal equations for this written down anywhere?  I realize SPH support was added since the first yt method paper, so it isn't in there, and I cannot seem to find anything in the docs.  I wonder if anyone has anything written down for this or if I'm missing it somewhere?</div><div><br></div><div>Thanks!</div><span class="HOEnZb"><font color="#888888"><div><br></div><div>Cameron</div><div><br></div><div><br clear="all"><div><br></div>-- <br><div data-smartmail="gmail_signature"><div dir="ltr"><div>Cameron Hummels<div><span style="font-size:12.8000001907349px">NSF Postdoctoral Fellow</span></div><div><span style="font-size:12.8000001907349px">Department of Astronomy</span></div><div><span style="font-size:12.8000001907349px">California Institute of Technology</span><br></div><div><a href="http://chummels.org" style="font-size:12.8000001907349px" target="_blank">http://chummels.org</a><br></div></div></div></div>
</div></font></span></div>
<br>______________________________<wbr>_________________<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" rel="noreferrer" target="_blank">http://lists.spacepope.org/<wbr>listinfo.cgi/yt-dev-spacepope.<wbr>org</a><br>
<br></blockquote></div><br></div>