Method ClipBoxInOut0D

Filters/General/vtkBoxClipDataSet.cxx:6431–6490

//------------------------------------------------------------------------------

void vtkBoxClipDataSet::ClipBoxInOut0D(vtkGenericCell* cell, vtkIncrementalPointLocator* locator,
  vtkCellArray** verts, vtkPointData* inPD, vtkPointData** outPD, vtkCellData* inCD,
  vtkIdType cellId, vtkCellData** outCD)
{
  vtkIdType cellType = cell->GetCellType();
  vtkIdList* cellIds = cell->GetPointIds();
  vtkCellArray* arrayvert = vtkCellArray::New();
  vtkPoints* cellPts = cell->GetPoints();
  vtkIdType npts = cellPts->GetNumberOfPoints();
  std::vector<vtkIdType> cellptId(npts);
  vtkIdType iid;
  const vtkIdType* v_id = nullptr;
  vtkIdType ptId;
  vtkIdType ptsvert = 1;

  int i;

  double v[3];

  for (i = 0; i < npts; i++)
  {
    cellptId[i] = cellIds->GetId(i);
  }

  // Convert all 0d cells to single vert.
  this->CellGrid(cellType, npts, cellptId.data(), arrayvert);

  unsigned int totalnewvert = arrayvert->GetNumberOfCells();
  for (unsigned int idlinenew = 0; idlinenew < totalnewvert; idlinenew++)
  {
    arrayvert->GetNextCell(ptsvert, v_id);

    // One way or another, we are adding the point.
    ptId = cellIds->GetId(v_id[0]);
    cellPts->GetPoint(v_id[0], v);

    if (locator->InsertUniquePoint(v, iid))
    {
      outPD[0]->CopyData(inPD, ptId, iid);
      outPD[1]->CopyData(inPD, ptId, iid);
    }

    // Clipping verts is easy.  Either it is inside the box or it isn't.
    if ((v[0] >= this->BoundBoxClip[0][0]) && (v[0] <= this->BoundBoxClip[0][1]) &&
      (v[1] >= this->BoundBoxClip[1][0]) && (v[1] <= this->BoundBoxClip[1][1]) &&
      (v[2] >= this->BoundBoxClip[2][0]) && (v[2] <= this->BoundBoxClip[2][1]))
    {
      // Vert is inside.
      int newCellId = verts[0]->InsertNextCell(1, &iid);
      outCD[0]->CopyData(inCD, cellId, newCellId);
    }
    else
    {
      // Vert is outside.
      int newCellId = verts[1]->InsertNextCell(1, &iid);
      outCD[1]->CopyData(inCD, cellId, newCellId);
    }
  }