[yt-users] cutting plane not working or changed

Nathan Goldbaum nathan12343 at gmail.com
Fri Mar 10 14:30:16 PST 2017


On Fri, Mar 10, 2017 at 1:50 PM, Slavin, Jonathan <jslavin at cfa.harvard.edu>
wrote:

> ​Hi Nathan,
>
> It turns out to be a bit different from what I said.  I tried with
> somewhat different parameters so as to make the frb square and it worked.
> So the issue seems to be with making non-square cuts, or maybe in making
> the frb from the​ cut.  I wonder if this is a recurrence of issue #1023,
> but it doesn't look like it.  I don't have time at the moment to post my
> dataset, but I will soon.  It is fairly large so I hope that's not a
> problem..
>

No worries, I think the limit on the curldrop is 10 GB.

In my test script I was making a non-square cutting plane and wasn't seeing
NaNs, but it's also possible that I wasn't triggering the failure mode
you're seeing. Unfortunately due to the nature of yt working with opaque
data on-disk usually it's sometimes very difficult to figure out how to
reproduce an issue without a test dataset to work with (even if the issue
isn't specific to that dataset).


>
> Jon
>
> On Fri, Mar 10, 2017 at 4:05 PM, <yt-users-request at lists.spacepope.org>
> wrote:
>
>> Date: Fri, 10 Mar 2017 12:14:55 -0800
>> From: Nathan Goldbaum <nathan12343 at gmail.com>
>> To: Discussion of the yt analysis package
>>         <yt-users at lists.spacepope.org>
>> Subject: Re: [yt-users] cutting plane not working or changed behavior
>> Message-ID:
>>         <CAJXewOmWwY1CU4nrafG5_=GrAScGBnq=VMNxxh3s0OUzDzaGVw at mail.gm
>> ail.com>
>> Content-Type: text/plain; charset="utf-8"
>>
>> Hi Jonathan,
>>
>> I'm not able to reproduce the behavior you're seeing with a test dataset
>> (e.g. see http://paste.yt-project.org/show/7080/, this uses a dataset
>> from
>> yt-project.org/data).
>>
>> Is there any chance you can share a dataset that triggers this?
>> Alternatively, can you reproduce this using one of the test datasets on
>> yt-project.org/data?
>>
>> The yt curldrop might be useful to share a dataset that triggers the issue
>> you're seeing, see docs.hub.yt. If you're not comfortable sharing the
>> dataset publicly, feel free to share it with me off-list.
>>
>> -Nathan
>>
>> On Fri, Mar 10, 2017 at 12:02 PM, Slavin, Jonathan <
>> jslavin at cfa.harvard.edu>
>> wrote:
>>
>> > Hi all,
>> >
>> > I have some code that uses a cutting plane to create a non-axis aligned
>> > slice plot:
>> > ds = yt.load(file)
>> > ctr = YTArray([15.,15.,15.],'pc')
>> > cut = ds.cutting(normal=(1.,-1.,0.), center=ctr,
>> north_vector=(0.,0.,1.))
>> > frb = cut.to_frb(width=(42.427,'pc'),height=(30.,'pc'),resolution=
>> > (512,724))
>> > image = frb['density']/2.126E-24
>> >
>> > I then view the image using the matplotlib imshow command.
>> > The data is a 30 pc**3 cartesian FLASH AMR dataset with maximum
>> resolution
>> > of
>> > ?512**3.  I used this code to produce an image, which I have, so I know
>> it
>> > used to work (last April).  Now however it does not seem to project in
>> the
>> > same direction and seems to be shaped incorrectly with a top portion of
>> the
>> > image filled with NaNs.  Has anyone else experienced this behavior?  I
>> > tried various things, for example reducing the resolution or switching
>> the
>> > order of those values (that leads to a core dump!).  The goal here is to
>> > create a slice along the plane aligned with the z axis and 45 degrees
>> away
>> > from the x and y axes (the data goes from 0 to 30 pc in each direction).
>> > Any help would be appreciated.?
>> >
>> > Jon
>> >
>> > --
>> > ________________________________________________________
>> > Jonathan D. Slavin                 Harvard-Smithsonian CfA
>> > jslavin at cfa.harvard.edu       60 Garden Street, MS 83
>> > phone: (617) 496-7981       Cambridge, MA 02138-1516
>> > cell: (781) 363-0035             USA
>> > ________________________________________________________
>>
>
>
>
>
> --
> ________________________________________________________
> Jonathan D. Slavin                 Harvard-Smithsonian CfA
> jslavin at cfa.harvard.edu       60 Garden Street, MS 83
> phone: (617) 496-7981       Cambridge, MA 02138-1516
> cell: (781) 363-0035             USA
> ________________________________________________________
>
>
> _______________________________________________
> yt-users mailing list
> yt-users at lists.spacepope.org
> http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.spacepope.org/pipermail/yt-users-spacepope.org/attachments/20170310/a293a361/attachment.htm>


More information about the yt-users mailing list