Function testCellIterator

Common/DataModel/Testing/Cxx/TestCellIterators.cxx:24–80

------------------------------------------------------------------------------ Compare the cell type, point ids, and points in 'grid' with those returned in 'iter'.

// Compare the cell type, point ids, and points in 'grid' with those returned
// in 'iter'.
bool testCellIterator(vtkCellIterator* iter, vtkUnstructuredGrid* grid)
{
  vtkIdType cellId = 0;
  vtkNew<vtkGenericCell> cell;
  iter->InitTraversal();
  while (!iter->IsDoneWithTraversal())
  {
    grid->GetCell(cellId, cell);

    if (iter->GetCellType() != cell->GetCellType())
    {
      std::cerr << "Type mismatch for cell " << cellId << std::endl;
      return false;
    }

    vtkIdType numPoints = iter->GetNumberOfPoints();
    if (numPoints != cell->GetNumberOfPoints())
    {
      std::cerr << "Number of points mismatch for cell " << cellId << std::endl;
      return false;
    }

    for (vtkIdType pointInd = 0; pointInd < numPoints; ++pointInd)
    {
      if (iter->GetPointIds()->GetId(pointInd) != cell->PointIds->GetId(pointInd))
      {
        std::cerr << "Point id mismatch in cell " << cellId << std::endl;
        return false;
      }

      double iterPoint[3];
      double cellPoint[3];
      iter->GetPoints()->GetPoint(pointInd, iterPoint);
      cell->Points->GetPoint(pointInd, cellPoint);
      if (iterPoint[0] != cellPoint[0] || iterPoint[1] != cellPoint[1] ||
        iterPoint[2] != cellPoint[2])
      {
        std::cerr << "Point mismatch in cell " << cellId << std::endl;
        return false;
      }
    }

    iter->GoToNextCell();
    ++cellId;
  }

  // ensure that we checked all of the cells
  if (cellId != grid->GetNumberOfCells())
  {
    std::cerr << "Iterator did not cover all cells in the dataset!" << std::endl;
    return false;
  }

  //  std::cout << "Verified " << cellId << " cells with a " << iter->GetClassName()
  //       << "." << std::endl;
  return true;
}