Method SplitCell

Filters/Geometry/vtkProjectSphereFilter.cxx:447–481

------------------------------------------------------------------------------

//------------------------------------------------------------------------------
void vtkProjectSphereFilter::SplitCell(vtkPointSet* input, vtkPointSet* output,
  vtkIdType inputCellId, vtkIncrementalPointLocator* locator, vtkCellArray* connectivity,
  int splitSide)
{
  // i screw up the canonical ordering of the cell but apparently this
  // gets fixed by vtkCell::Clip().
  vtkCell* cell = input->GetCell(inputCellId);
  vtkNew<vtkDoubleArray> cellScalars;
  cellScalars->SetNumberOfTuples(cell->GetNumberOfPoints());
  double coord[3];
  for (vtkIdType pt = 0; pt < cell->GetNumberOfPoints(); pt++)
  {
    output->GetPoint(cell->GetPointId(pt), coord);
    if (splitSide == 0 && coord[0] > this->SplitLongitude + 180.)
    {
      coord[0] -= 360.;
    }
    else if (splitSide == 1 && coord[0] < this->SplitLongitude + 180.)
    {
      coord[0] += 360.;
    }
    cellScalars->SetValue(pt, coord[0]);
    cell->GetPoints()->SetPoint(pt, coord);
  }
  vtkIdType numberOfCells = output->GetNumberOfCells();
  double splitLocation = (splitSide == 0 ? -180. : 180.);
  cell->Clip(splitLocation, cellScalars.GetPointer(), locator, connectivity, output->GetPointData(),
    output->GetPointData(), input->GetCellData(), inputCellId, output->GetCellData(), splitSide);
  // if the grid was an unstructured grid we have to update the cell
  // types and locations for the created cells.
  if (vtkUnstructuredGrid* ugrid = vtkUnstructuredGrid::SafeDownCast(output))
  {
    this->SetCellInformation(ugrid, cell, output->GetNumberOfCells() - numberOfCells);
  }
}

//------------------------------------------------------------------------------
void vtkProjectSphereFilter::SetCellInformation(