[Paraview] Issue with structured grids and extraction of boundaries with cell data

[Mickael Philit]

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,

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://public.kitware.com/pipermail/paraview/attachments/20150325/82b46514/attachment.html>