<div dir="ltr">Hi Matt,<div style>+1 for subclassing a DomainSubset common to ARTIO, ART and RAMSES.</div><div style><br></div><div style>chris</div></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Wed, Apr 17, 2013 at 1:35 PM, Douglas Harvey Rudd <span dir="ltr"><<a href="mailto:drudd@uchicago.edu" target="_blank">drudd@uchicago.edu</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi Matt,<br>
<br>
I'm +1 on specifying the DomainSubset interface though subclasses, and it makes<br>
sense to integrate it with selection.<br>
<br>
One question I have relating to deposition in particular is whether the spatial fields<br>
are expected to impose any sort of order on the octs they return.  If we select a region<br>
including both particles and octs, but the two are completely unordered, matching<br>
one set to the other may be very expensive.<br>
<br>
Do the deposition routines already use some organization scheme (oct tree, say) to<br>
match particles to grids?<br>
<br>
Doug<br>
<br>
<br>
On Apr 17, 2013, at 3:24 PM, Matthew Turk <<a href="mailto:matthewturk@gmail.com">matthewturk@gmail.com</a>><br>
<div class="HOEnZb"><div class="h5"> wrote:<br>
<br>
> Hi all,<br>
><br>
> Today in IRC, Chris and I went back and forth on the particle<br>
> deposition stuff in advance of the meeting next week.  My task right<br>
> now is to come up with spatial field chunking for Octs in a consistent<br>
> way; of course I'll write this up into the extant YTEP, which will<br>
> allow for feedback particularly from ARTIO, ART and RAMSES folks, but<br>
> as a quick first pass I was wondering what everyone thought of:<br>
><br>
> 1) Unifying DomainSubsets into a single object with subclasses<br>
> 2) Making that new base class a subclass, itself, of YTSelectionContainer<br>
><br>
> This would mean that this new base class would be very similar to the<br>
> GridPatch object -- which it already is -- and would implement, most<br>
> crucially, the __getitem__ accessor pattern.  This would return a<br>
> 2x2x2xNoct field.<br>
><br>
> [+-][01]?  Specifically, Doug, Chris, and SamL?<br>
><br>
> Sidenote: I discovered that we will not have to rewrite spatial fields<br>
> that currently work for grids; my understanding of the numpy<br>
> broadcasting system was not complete, and I didn't realize that this<br>
> script does *exactly what I want it to.<br>
><br>
> <a href="http://paste.yt-project.org/show/3391/" target="_blank">http://paste.yt-project.org/show/3391/</a><br>
><br>
> This means that our Oct spatial fields will be (2+NGZ*2, 2+NGZ*2,<br>
> 2+NGZ*2, NOCT) shaped.<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>
</div></div></blockquote></div><br></div>