Matt,<div>  OK, thanks for clearing that up.  The deposition of that quantity (or others like it) could be useful in a number of ways.  It is a very specific type of need, but I think there could be others.  For instance, one could use this to post-process simulations with rad transport. One could have a collection of points and construct a radiation field quantity on the grid for instance in yt based on distance from the points and assumed luminosity.  That is, if one did not do active rad transport in a run.  </div>
<div><br></div><div>Anyway, interesting.</div><div><br></div><div>Thanks.</div><div><br></div><div>Eric<br><br><div class="gmail_quote">On Tue, Oct 25, 2011 at 11:30 AM, Matthew Turk <span dir="ltr"><<a href="mailto:matthewturk@gmail.com">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 Eric,<br>
<div class="im"><br>
On Tue, Oct 25, 2011 at 11:26 AM, Eric Hallman <<a href="mailto:hallman13@gmail.com">hallman13@gmail.com</a>> wrote:<br>
> Matt and everyone,<br>
>   I'm also interested in a tool of this type.  Stephen  has a particular use<br>
> in mind, but also, one could imagine calculating an "all-sky" light-cone<br>
> type of thing using the similar code tool.  Does the HEALPIX stuff do that?<br>
<br>
</div>Yes, the existing HEALpix will do this; what I believe Stephen is<br>
addressing is the deposition of column density back in the grids.  The<br>
existing HEALpix will calculate column density (and can also calculate<br>
star particle contributions) between two radii and return that to the<br>
user.<br>
<div class="im"><br>
>  We had discussed in the past doing a projection outward from a point, then<br>
> stacking the higher z outputs in shells around the z=0 set, to make a<br>
> spherical lightcone projection.<br>
<br>
</div>I believe this is possible with existing infrastructure, but keep in<br>
mind that it will only currently work with a static-resolution HEALpix<br>
pixelization.  Stephen's will deposit column density back on the grid<br>
so that it can then be treated as any other field in yt.<br>
<br>
-Matt<br>
<div><div></div><div class="h5"><br>
> Not sure if the current tools are capable of the individual pieces of such<br>
> an operation, but what Stephen's describing would certainly do part of it.<br>
> Eric<br>
><br>
> On Tue, Oct 25, 2011 at 11:22 AM, Matthew Turk <<a href="mailto:matthewturk@gmail.com">matthewturk@gmail.com</a>><br>
> wrote:<br>
>><br>
>> Hi Stephen,<br>
>><br>
>> For what it's worth, I also agree with Cameron that calling it<br>
>> "ColumnDensity" is a bit too broad, and maybe it should be called<br>
>> something like "RadialColumnDensity" or something similarly indicative<br>
>> of its nature to indicate it's not the same as a projection.<br>
>><br>
>> Can you also maybe take a minute to outline a couple use cases?<br>
>><br>
>> -Matt<br>
>><br>
>> On Tue, Oct 25, 2011 at 9:23 AM, Stephen Skory <<a href="mailto:s@skory.us">s@skory.us</a>> wrote:<br>
>> > Cameron & Matt,<br>
>> ><br>
>> > What it's called isn't a big concern to me, but I see what you're<br>
>> > saying, Cameron.<br>
>> ><br>
>> >> The issue about field detector is an interesting one.  I guess I don't<br>
>> >> understand why asking for 'x' is a problem.  Your solution sounds good<br>
>> >> to me.<br>
>> ><br>
>> > When I wasn't looking for the field detector as I am now, what was<br>
>> > happening is data['x' or 'y' or 'z'] would return some values that<br>
>> > weren't cell centers, which when passed to the interpolation stuff<br>
>> > would not work. So asking for 'x' wasn't a problem, but the values it<br>
>> > returned were not what I wanted.<br>
>> ><br>
>> >> How fast is the code?  It looks to me like it probably does quite a<br>
>> >> few expensive operations<br>
>> ><br>
>> > Running on a 40^3 dataset on my 2.0Ghz i7 lappy on battery power gives<br>
>> > about 300,000 cells/second for the whole process (HEALPix with 5<br>
>> > surfaces + interpolation). I think I'm close to making it parallel,<br>
>> > but some weird stuff is popping up that I don't quite understand just<br>
>> > yet.<br>
>> ><br>
>> >>... would you be willing at some point in the<br>
>> >> future to explore replacing it with an actual adaptive healpix<br>
>> >> operation?<br>
>> ><br>
>> > Perhaps. It seems to me before you said that that would be quite a bit<br>
>> > more difficult, which looking at the source is true. Everything in<br>
>> > this current attempt is using numpy vectorization, so I don't know how<br>
>> > much more speed can be squeezed out of this method.<br>
>> ><br>
>> > --<br>
>> > Stephen Skory<br>
>> > <a href="mailto:s@skory.us">s@skory.us</a><br>
>> > <a href="http://stephenskory.com/" target="_blank">http://stephenskory.com/</a><br>
>> > <a href="tel:510.621.3687" value="+15106213687">510.621.3687</a> (google voice)<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>
><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>
_______________________________________________<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>