[yt-users] add_ray problem

[Elizabeth Tasker]

Hi Matt,

I did think about that, but it looks from the graph 
DD0010_LineQueryPlot_0_t_Density.png that the peak is a occuring at 
about 0.5, when it should be at ~0.9 in a box that runs from 0-->1.

I've attached the script.

Elizabeth



Matthew Turk wrote:
> 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
>>
>>
>>     
> _______________________________________________
> yt-users mailing list
> yt-users at lists.spacepope.org
> http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
>   

-------------- next part --------------
A non-text attachment was scrubbed...
Name: sedov_yt.py
Type: text/x-python
Size: 1225 bytes
Desc: not available
URL: <http://lists.spacepope.org/pipermail/yt-users-spacepope.org/attachments/20101028/40f6ff95/attachment.py>