Issue with contours and Exodus output
I have a simplified example of a 2x2x2 cuboid mesh with constant values on each element, stored in the variable group_0 in Exodus format (this is all done with MOOSE). My environment is Yt 3.4.1 and Python 3.6.5. If I do a SlicePlot and try to use annotate_contour, I end up with an indexing error. I've tried to follow the traceback but frankly I'm not well versed enough with Yt to make any sense of it. Interestingly, if I instead attempt to plot a 2x2x1 mesh cut in the z-plane, the contours will plot without issue. I've uploaded the full output, the simple (4 line) example, and my Exodus file: http://use.yt/upload/078c018a Thank you in advance for the assistance.
Hi Logan,
I'm able to reproduce this issue, thank you for the detailed report.
When yt explodes like this we generally consider it a bug. When you see it
with your own scripts please feel free to open an issue on github.
I think there's an issue with how yt is masking the data here. Internally
when we ask for the widths of the elements in the slice, we get back data
for 4 elements as expected. However when we ask for the field data for the
group_0 variable, we get back an array with 8 elements, with four of the
array entries set to zero. I think what we should be getting back is an
array of shape (4,). Unfortunately I don't know as much about the internals
of how yt supports unstructured meshes so I don't have a satisfactory fix
yet, I'd like to hear back from one or two other yt developers.
I've gone ahead and filed an issue which you can subscribe to for
notifications:
https://github.com/yt-project/yt/issues/1916
-Nathan
On Tue, Jul 24, 2018 at 9:51 AM, Logan Harbour
I have a simplified example of a 2x2x2 cuboid mesh with constant values on each element, stored in the variable group_0 in Exodus format (this is all done with MOOSE). My environment is Yt 3.4.1 and Python 3.6.5.
If I do a SlicePlot and try to use annotate_contour, I end up with an indexing error. I've tried to follow the traceback but frankly I'm not well versed enough with Yt to make any sense of it. Interestingly, if I instead attempt to plot a 2x2x1 mesh cut in the z-plane, the contours will plot without issue.
I've uploaded the full output, the simple (4 line) example, and my Exodus file: http://use.yt/upload/078c018a
Thank you in advance for the assistance. _______________________________________________ yt-users mailing list -- yt-users@python.org To unsubscribe send an email to yt-users-leave@python.org
participants (2)
-
Logan Harbour
-
Nathan Goldbaum