MantisBT - ParaView
View Issue Details
0008316ParaView(No Category)public2008-12-18 09:352011-01-13 17:00
Jean M. Favre 
Berk Geveci 
normalcrashalways
closedfixed 
 
3.8 
0008316: paraview crashes during initialization of the slicer with a dataset with very large bounds
The dataset below crashes paraview immediately after creating a Slice filter

Slicing itself works very well, since the following python code (look at plane.Origin) is able to cut it. The problem is not the slicer, but the slice widget initialization:
=======================================================================
Program received signal SIGSEGV, Segmentation fault.
[Switching to Thread 0x7fa971e20740 (LWP 12072)]
0x00007fa960295856 in vtkMergePoints::InsertUniquePoint (this=0x2682920, x=0x7fff79f89e50, id=@0x7fff79f89e20)
    at /local/apps/ParaView-3.4.0/VTK/Filtering/vtkMergePoints.cxx:136
136 bucket = this->HashTable[idx];
(gdb) where
#0 0x00007fa960295856 in vtkMergePoints::InsertUniquePoint (this=0x2682920, x=0x7fff79f89e50, id=@0x7fff79f89e20)
    at /local/apps/ParaView-3.4.0/VTK/Filtering/vtkMergePoints.cxx:136
#1 0x00007fa9620dbb4f in vtkFeatureEdges::RequestData (this=0x277e690, inputVector=0x277e1d0, outputVector=0x277d630)
    at /local/apps/ParaView-3.4.0/VTK/Graphics/vtkFeatureEdges.cxx:316
=======================================================================

# this python script slices the data without crashing, since it does not
# create a slice widget

cutter = servermanager.filters.Cut()
cutter.Input = servermanager.OutputPort(reader, 0)

plane = servermanager.implicit_functions.Plane()
plane.Normal = [0, 0, 1]
plane.Origin = [25e18, 25e18, 25e18]
cutter.CutFunction = plane

cutterRep = servermanager.CreateRepresentation(cutter, view)
pm.RegisterProxy('sources', 'cutter', cutter)
pm.RegisterProxy('representations', 'cutterRep', cutterRep)
<?xml version="1.0"?>
<VTKFile type="ImageData" version="0.1" byte_order="LittleEndian" compressor="vtkZLibDataCompressor">
  <ImageData WholeExtent="0 50 0 50 0 50" Origin="0 0 0" Spacing="1e+18 1e+18 1e+18">
  <Piece Extent="0 50 0 50 0 50">
    <PointData>
    </PointData>
    <CellData>
    </CellData>
  </Piece>
  </ImageData>
</VTKFile>
No tags attached.
txt patch.txt (1,950) 2009-05-15 09:16
https://www.vtk.org/Bug/file/7251/patch.txt
Issue History
2008-12-18 09:35Jean M. FavreNew Issue
2008-12-29 13:36Berk GeveciStatusbacklog => tabled
2008-12-29 13:36Berk GeveciAssigned To => Berk Geveci
2008-12-29 13:36Berk GeveciCategory3.4 => 3.6
2009-02-17 15:59Utkarsh AyachitCategory3.6 => 3.8
2009-05-13 13:40Utkarsh AyachitTarget Version => 3.8
2009-05-15 09:15Jean M. FavreNote Added: 0016471
2009-05-15 09:16Jean M. FavreFile Added: patch.txt
2009-08-26 13:47Berk GeveciNote Added: 0017219
2009-08-26 13:47Berk GeveciStatustabled => closed
2009-08-26 13:47Berk GeveciResolutionopen => fixed
2011-01-13 17:00Source_changeset_attached => VTK master a2bd8391
2011-01-13 17:00Source_changeset_attached => VTK master 020ef709
2011-06-16 13:10Zack GalbreathCategory => (No Category)

Notes
(0016471)
Jean M. Favre   
2009-05-15 09:15   
I have found a simple fix after looking again at the source code of ParaView3/VTK/Filtering/vtkMergePoints.cxx

the problem is an overflow because of a mix of int and vtkIdType.
I replaced the ints by vtkIdTypes and the segmentation fault goes away

here is a patch file. Apply with patch -p0
(0017219)
Berk Geveci   
2009-08-26 13:47   
Committed Jean's patch.