Hello,
my previous mail was wrong. The present problem is not related to those
other bugs.
Here, when using ExtractSubset the extent for cell data in miscalculated
if i==imax or j==jmax or k == kmax.
in vtkExtractGrid.cxx: we have the call
"this->Internal->CopyCellData(inExt,outExt,cd,outCD);" that use point extent
and then in vtkExtractStructuredGridHelper.cxx the extent is computed
for cells
"vtkStructuredData::GetCellExtentFromPointExtent(outExt,outCellExt);"
but this later function does not know about the whole extent of cell
data. Thus you need to pass down a valid extent for cells before calling it.
If you have inExt = {0, 30, 0, 20, 0, 10} and you pass outExt =
{0,30,0,20, 10,10} the result outCellExt will be {0,29,0,19,*10,10*} and
then the copy of cell data will be done outside of the array leading to
segfault.
The outExt extent should have been corrected to be {0,30,0,20,9,10} or
{0,30,0,20,9,9}...
I hope it is clearer now.
Best Regards
Mickael
On 20/03/2015 16:26, Etienne Tang wrote:
Hello,
I'm working with a 3D structured dataset containing cell data. I tried
to extract a 2D boundary of this dataset corresponding to the maximum
index of a dimension with the Extract Subset filter and ParaView 4.3.1,
and I found an issue: the resulting dataset is incorrect. It holds cell
data from the opposite boundary (eg. if I try to extract the i = imax
boundary, I get values from i = 0 boundary), and these values are
shifted of one index in the other dimensions.
This case worked with ParaView 4.1.0, maybe it is related to this other
bug about structured data :
http://markmail.org/message/fvatgidqbvnufgi6#query:+page:1+mid:fvatgidqbvnufgi6+state:results
?
Here are some pictures showing the problem on a simple dataset which can
be used to reproduce it:
- The simple .vts file : http://z1.gb2n.org/pv/test.vts
- The original dataset : http://z1.gb2n.org/pv/original_dataset.png
- The result of the ExtractSubset filter to extract the i = imax
boundary : http://z1.gb2n.org/pv/extract_subset.png
Best regards,
_______________________________________________
Powered by www.kitware.com
Visit other Kitware open-source projects at
http://www.kitware.com/opensource/opensource.html
Please keep messages on-topic and check the ParaView Wiki at:
http://paraview.org/Wiki/ParaView
Search the list archives at: http://markmail.org/search/?q=ParaView
Follow this link to subscribe/unsubscribe:
http://public.kitware.com/mailman/listinfo/paraview
