[yt-users] add_ray problem

[Matthew Turk]

Hi Elizabeth,

One possible issue is that the non-ortho ray is normalized to having a
length in the "t" direction of 1.0; it's solved as a simple line
integrated from point A to point B, and the way that's done is just
with normalization..  Can you show us the entirety of the plotting
script?

-Matt

On Thu, Oct 28, 2010 at 11:47 AM, Elizabeth Tasker <taskere at mcmaster.ca> wrote:
> Hi,
>
> I've got some (different) problems with add_ray. I have output from the
> Sedov blast test (so unitless sytems) which is run with a boxsize of 10. I
> can analyse the results using:
>
> ray = pf.h.ray([5., 5., 5.], [10.,10.,10.], "Density")
> p = pc.add_ray([5.,5.,5.], [10.,10.,10.], "Density")
> oray = pf.h.ortho_ray(0, [5., 5.], "Density")
>
> but only ortho_ray gives the correct result, and it's limited to going along
> the axes.
>
> I've attached the results from the 3 techniques. fig0010.png shows the
> analytical solution (red line), the ortho_ray solution (green x) and the
> pf.h.ray solution (blue circles). DD0010_LineQueryPlot shows the results
> from pc.add_ray.
>
> Does anyone know why this isn't behaving? A slice of the shock confirms it
> looks nice and spherical.
>
> Elizabeth
>
> _______________________________________________
> yt-users mailing list
> yt-users at lists.spacepope.org
> http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
>
>