[yt-users] yt volume render/AMRKDtree help

[Renyue Cen]

Hi John,

I am testing Sam's transpose fix now to see if that changes things.
I will let you know very shortly when I am done that and you may then
upgrade to the latest version.

Thanks,
Renyue

On Jan 8, 2013, at 11:56 AM, John Wise wrote:

> Hi Renyue,
> 
> Yes, I reverted it back to v2.4 a few weeks ago.  I can update yt to the latest version to see if it helps with these problems.
> 
> John
> 
> On 01/08/2013 11:38 AM, Renyue Cen wrote:
>> Hi Sam,
>> 
>> I am using the version that John installed on Pleiades.
>> I got some other unrelated problems with profiling before and John was
>> kind enough to
>> revert back to an older version (I think). Anyway, this is what I got:
>> 
>> /u/jhwise/local/lib/python2.7/site-packages/yt-2.4-py2.7-linux-x86_64.egg
>> 
>> Separately, I will try your transpose fix and will let you know.
>> 
>> But the KDTree was really smoking fast (50-100 times compared to
>> without) that
>> makes things very efficient. I would like to thank you for that.
>> 
>> Thanks very much,
>> Renyue
>> 
>> 
>> On Jan 8, 2013, at 11:25 AM, Sam Skillman wrote:
>> 
>>> Hi Renyue,
>>> 
>>> I have actually seen this on one other system but I could never
>>> reproduce it. Could you tell me the version of yt you are using? If
>>> you run "yt instinfo" it should tell you. I haven't seen this bug for
>>> the last few months so I thought it was gone.
>>> 
>>> Additionally, can you try the following at the end of your script:
>>> 
>>> im=cam.snapshot('test1.png')
>>> 
>>> write_bitmap(im, 'test2.png', transpose=False)
>>> 
>>> I think there may have been something wrong with the transposing of
>>> the image for a short time period a few months ago.  I have very
>>> sporadic and unreliable internet at AAS but I'll try to respond as
>>> quickly as possible.
>>> 
>>> Sorry for the trouble,
>>> 
>>> Sam
>>> 
>>> 
>>> On Jan 7, 2013 10:11 AM, "Renyue Cen" <cen at astro.princeton.edu
>>> <mailto:cen at astro.princeton.edu>> wrote:
>>> 
>>>    Hi,
>>> 
>>>    I did some simple volume rendering with the following script:
>>> 
>>>        volume2 = AMRKDTree(pf, fields=["Dark_Matter_Density"],
>>>                    no_ghost=False, tree_type="domain",
>>>                    le=c-0.5*WW, re=c+0.5*WW)
>>>        cam = pf.h.camera(c, L, W, N, tf, volume=volume2, no_ghost=False,
>>>                    north_vector=L, steady_north=True)
>>>        cam.snapshot(fn="%s_iso-DMdensity-%3.3d.png" % (filenameTHIS, j))
>>> 
>>> 
>>>    I got rather strange results in that the pictures look symmetric,
>>>    which I am pretty
>>>    sure can not be true.
>>>    I attach the obtained plot.
>>>    Note that I am using KD tree and using 32 cores.
>>> 
>>>    Your help at your earliest convenience is appreciated.
>>>    Best,
>>>    Renyue
>>> 
>>> 
>>>    <C15z1600.0043_iso-DMdensity-000.png>
>>> 
>> 
> 
> -- 
> John Wise
> Assistant Professor of Physics
> Center for Relativistic Astrophysics, Georgia Tech