View Issue Details [ Jump to Notes ]

[ Print ]

Project

Category

View Status

Date Submitted

Last Update

0008316

ParaView

(No Category)

public

2008-12-18 09:35

2011-01-13 17:00

Reporter

Jean M. Favre

Assigned To

Berk Geveci

Priority

normal

Severity

crash

Reproducibility

always

Status

closed

Resolution

fixed

Platform

OS Version

Product Version

Target Version

3.8

Fixed in Version

Summary

0008316: paraview crashes during initialization of the slicer with a dataset with very large bounds

Description

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)

Additional Information

<?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>

Tags

No tags attached.

Project

Topic Name

Type

Attached Files

patch.txt [^] (1,950 bytes) 2009-05-15 09:16 [Show Content] [Hide Content]

--- ParaView-3.4.0/VTK/Filtering/vtkMergePoints.cxx	2007-11-06 18:13:33.000000000 +0100
+++ ParaView3/VTK/Filtering/vtkMergePoints.cxx	2009-05-15 14:13:40.000000000 +0200
@@ -28,21 +28,21 @@
 // -1.
 vtkIdType vtkMergePoints::IsInsertedPoint(const double x[3])
 {
-  int i, ijk0, ijk1, ijk2;
+  vtkIdType i, ijk0, ijk1, ijk2;
   vtkIdType idx;
   vtkIdList *bucket;
   //
   //  Locate bucket that point is in.
   //
-  ijk0 = static_cast<int>(
+  ijk0 = static_cast<vtkIdType>(
     static_cast<double> ((x[0] - this->Bounds[0]) / 
                          (this->Bounds[1] - this->Bounds[0]))
     * (this->Divisions[0] - 1));
-  ijk1 = static_cast<int>(
+  ijk1 = static_cast<vtkIdType>(
     static_cast<double> ((x[1] - this->Bounds[2]) / 
                          (this->Bounds[3] - this->Bounds[2]))
     * (this->Divisions[1] - 1));
-  ijk2 = static_cast<int>(
+  ijk2 = static_cast<vtkIdType>(
     static_cast<double> ((x[2] - this->Bounds[4]) / 
                          (this->Bounds[5] - this->Bounds[4]))
     * (this->Divisions[2] - 1));
@@ -110,22 +110,22 @@
 
 int vtkMergePoints::InsertUniquePoint(const double x[3], vtkIdType &id)
 {
-  int i, ijk0, ijk1, ijk2;
+  vtkIdType i, ijk0, ijk1, ijk2;
   vtkIdType idx;
   vtkIdList *bucket;
 
   //
   //  Locate bucket that point is in.
   //
-  ijk0 = static_cast<int>(
+  ijk0 = static_cast<vtkIdType>(
     static_cast<double> ((x[0] - this->Bounds[0]) / 
                          (this->Bounds[1] - this->Bounds[0]))
     * (this->Divisions[0] - 1));
-  ijk1 = static_cast<int>(
+  ijk1 = static_cast<vtkIdType>(
     static_cast<double> ((x[1] - this->Bounds[2]) / 
                          (this->Bounds[3] - this->Bounds[2]))
     * (this->Divisions[1] - 1));
-  ijk2 = static_cast<int>(
+  ijk2 = static_cast<vtkIdType>(
     static_cast<double> ((x[2] - this->Bounds[4]) / 
                          (this->Bounds[5] - this->Bounds[4]))
     * (this->Divisions[2] - 1));

Relationships

Notes
(0016471) Jean M. Favre (reporter) 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 (administrator) 2009-08-26 13:47	Committed Jean's patch.

Issue History
Date Modified	Username	Field	Change
2008-12-18 09:35	Jean M. Favre	New Issue
2008-12-29 13:36	Berk Geveci	Status	backlog => tabled
2008-12-29 13:36	Berk Geveci	Assigned To	=> Berk Geveci
2008-12-29 13:36	Berk Geveci	Category	3.4 => 3.6
2009-02-17 15:59	Utkarsh Ayachit	Category	3.6 => 3.8
2009-05-13 13:40	Utkarsh Ayachit	Target Version	=> 3.8
2009-05-15 09:15	Jean M. Favre	Note Added: 0016471
2009-05-15 09:16	Jean M. Favre	File Added: patch.txt
2009-08-26 13:47	Berk Geveci	Note Added: 0017219
2009-08-26 13:47	Berk Geveci	Status	tabled => closed
2009-08-26 13:47	Berk Geveci	Resolution	open => fixed
2011-01-13 17:00		Source_changeset_attached	=> VTK master a2bd8391
2011-01-13 17:00		Source_changeset_attached	=> VTK master 020ef709
2011-06-16 13:10	Zack Galbreath	Category	=> (No Category)