Matt and everyone,<div>  I'm also interested in a tool of this type.  Stephen  has a particular use in mind, but also, one could imagine calculating an "all-sky" light-cone type of thing using the similar code tool.  Does the HEALPIX stuff do that?  We had discussed in the past doing a projection outward from a point, then stacking the higher z outputs in shells around the z=0 set, to make a spherical lightcone projection.  </div>
<div><br></div><div>Not sure if the current tools are capable of the individual pieces of such an operation, but what Stephen's describing would certainly do part of it.</div><div><br></div><div>Eric<br><br><div class="gmail_quote">
On Tue, Oct 25, 2011 at 11:22 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 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>
<div class="im"><br>
On Tue, Oct 25, 2011 at 9:23 AM, Stephen Skory <<a href="mailto:s@skory.us">s@skory.us</a>> wrote:<br>
</div><div><div></div><div class="h5">> 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>
</div></div></blockquote></div><br></div>