[yt-dev] problem with disk data container and related fields in yt-3.0

Matthew Turk matthewturk at gmail.com
Wed Aug 13 08:37:45 PDT 2014


On Wed, Aug 13, 2014 at 10:32 AM, Douglas Harvey Rudd
<drudd at uchicago.edu> wrote:
>>
>>>
>>> There is definitely something strange going on that eats away that portion of the lower part of the cylinder, which must be something to do with the way the code computes the periodic radial distance, then subtracts the portion projected onto the cylinder normal to obtain the cylindrical radius.
>>
>> My understanding from looking into it is that the only problem with
>> the disk bbox selection is in selecting bboxes that have no corners
>> within the cylinder.  All cell/point selection is correct, I believe.
>
> I'm not sure how that's possible, given the shape of the carved out region.  That appears to have a refined region, so most of those grids should have all 8 corners inside the cylinder.  I've reproduced the problem and will probably throw an hour or two at it today to see if I can find a simple fix.  Either way, I agree with Britton that the no-op bbox selector should be swapped in immediately.

Thank you, Doug -- I really appreciate it.

>
> Doug
>
>
>
>>>
>>>
>>> On Aug 12, 2014, at 7:32 AM, Matthew Turk <matthewturk at gmail.com> wrote:
>>>
>>>> Hi Britton,
>>>>
>>>> On Tue, Aug 12, 2014 at 7:30 AM, Britton Smith <brittonsmith at gmail.com> wrote:
>>>>> Hi all,
>>>>>
>>>>> Some issues seem to have cropped up in the disk data container in yt-3.0.  I
>>>>> made a disk object in a periodic simulation in the following way:
>>>>> http://paste.yt-project.org/show/5021/
>>>>>
>>>>> I then made a project of the "cylindrical_z" field, which yielded the
>>>>> following:
>>>>> http://i.imgur.com/0XSv5Zs.png
>>>>>
>>>>> For reference, the domain is ~45 Mpc in length on a side, so the sphere
>>>>> should be just the object at the bottom, extending a little less than
>>>>> halfway up.  I don't know what that top part is about.  My guess is that
>>>>> it's an issue with periodicity.
>>>>
>>>> Yes, looks like it to me.
>>>>
>>>>>
>>>>> Secondly, I made a similar projection of the "cylindrical_r" field, which
>>>>> gave this:
>>>>> http://i.imgur.com/h8SEJp8.png
>>>>>
>>>>> As you can see, it is not symmetric as one would expect.  For reference,
>>>>> this is what I get from yt-2.x:
>>>>> http://i.imgur.com/O1xT9MI.png
>>>>> This is what I would expect.
>>>>
>>>> Are you doing a MIP projection here?  Could you do a MIP projection so
>>>> we can compare?
>>>>
>>>> It almost looks to me like the center is set incorrectly.
>>>>
>>>>>
>>>>> Does anyone have any insight on this?
>>>>>
>>>>> Britton
>>>>>
>>>>> _______________________________________________
>>>>> yt-dev mailing list
>>>>> yt-dev at lists.spacepope.org
>>>>> http://lists.spacepope.org/listinfo.cgi/yt-dev-spacepope.org
>>>>>
>>>> _______________________________________________
>>>> yt-dev mailing list
>>>> yt-dev at lists.spacepope.org
>>>> http://lists.spacepope.org/listinfo.cgi/yt-dev-spacepope.org
>>>
>>> _______________________________________________
>>> yt-dev mailing list
>>> yt-dev at lists.spacepope.org
>>> http://lists.spacepope.org/listinfo.cgi/yt-dev-spacepope.org
>> _______________________________________________
>> yt-dev mailing list
>> yt-dev at lists.spacepope.org
>> http://lists.spacepope.org/listinfo.cgi/yt-dev-spacepope.org
>
> _______________________________________________
> yt-dev mailing list
> yt-dev at lists.spacepope.org
> http://lists.spacepope.org/listinfo.cgi/yt-dev-spacepope.org



More information about the yt-dev mailing list