Index: vtkPointSet.cxx
===================================================================
RCS file: /cvsroot/ParaView3/ParaView3/VTK/Filtering/vtkPointSet.cxx,v
retrieving revision 1.8
diff -u -r1.8 vtkPointSet.cxx
--- vtkPointSet.cxx	27 Oct 2008 18:25:46 -0000	1.8
+++ vtkPointSet.cxx	11 Mar 2009 23:00:22 -0000
@@ -195,7 +195,32 @@
       return (-1); //if point completely outside of data
       }
 
-    this->GetPointCells(ptId, cellIds);
+    // It is possible that the toplogy is not fully connected.  That is, two
+    // points in the point list could be coincident.  Handle that by looking
+    // at every point within the tolerance and consider all cells connected.
+    double ptCoord[3];
+    this->GetPoint(ptId, ptCoord);
+    vtkIdList *coincidentPtIds = vtkIdList::New();
+    this->Locator->FindPointsWithinRadius(tol2, ptCoord, coincidentPtIds);
+    if (coincidentPtIds->GetNumberOfIds() > 0)
+      {
+      this->GetPointCells(coincidentPtIds->GetId(0), cellIds);
+      if (coincidentPtIds->GetNumberOfIds() > 1)
+        {
+        vtkIdList *tempCellIds = vtkIdList::New();
+        for (vtkIdType i = 1; i < coincidentPtIds->GetNumberOfIds(); i++)
+          {
+          this->GetPointCells(coincidentPtIds->GetId(1), tempCellIds);
+          for (vtkIdType j = 0; j < tempCellIds->GetNumberOfIds(); j++)
+            {
+            cellIds->InsertUniqueId(tempCellIds->GetId(j));
+            }
+          }
+        tempCellIds->Delete();
+        }
+      }
+    coincidentPtIds->Delete();
+
     curIdListPosn = 0;
     if ( cellIds->GetNumberOfIds() > 0 )
       {