Method ClipBoxInOut1D

Filters/General/vtkBoxClipDataSet.cxx:5753–5970

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

void vtkBoxClipDataSet::ClipBoxInOut1D(vtkPoints* newPoints, vtkGenericCell* cell,
  vtkIncrementalPointLocator* locator, vtkCellArray** lines, vtkPointData* inPD,
  vtkPointData** outPD, vtkCellData* inCD, vtkIdType cellId, vtkCellData** outCD)
{
  vtkIdType cellType = cell->GetCellType();
  vtkIdList* cellIds = cell->GetPointIds();
  vtkCellArray* arrayline = vtkCellArray::New();
  vtkPoints* cellPts = cell->GetPoints();
  vtkIdType npts = cellPts->GetNumberOfPoints();
  std::vector<vtkIdType> cellptId(npts);
  vtkIdType iid[2];
  const vtkIdType* v_id = nullptr;
  vtkIdType ptId;
  vtkIdType tab_id[2];
  vtkIdType ptsline = 2;

  int i, j;
  unsigned int allInside;
  unsigned int planes;
  unsigned int cutInd;

  double value;
  double t;
  double v[3], x[3];
  double v_line[2][3];

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

  // Convert all 1d cells to single line.
  this->CellGrid(cellType, npts, cellptId.data(), arrayline);

  unsigned int totalnewline = arrayline->GetNumberOfCells();
  for (unsigned int idlinenew = 0; idlinenew < totalnewline; idlinenew++)
  {
    arrayline->GetNextCell(ptsline, v_id);

    for (allInside = 1, i = 0; i < 2; i++)
    {
      ptId = cellIds->GetId(v_id[i]);
      cellPts->GetPoint(v_id[i], v);

      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])))
      {
        // outside
        allInside = 0;
      }

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