Method ClipBox2D

Filters/General/vtkBoxClipDataSet.cxx:4015–4361

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

//-------------------------------------------------------
void vtkBoxClipDataSet::ClipBox2D(vtkPoints* newPoints, vtkGenericCell* cell,
  vtkIncrementalPointLocator* locator, vtkCellArray* tets, vtkPointData* inPD, vtkPointData* outPD,
  vtkCellData* inCD, vtkIdType cellId, vtkCellData* outCD)
{
  vtkIdType cellType = cell->GetCellType();
  vtkIdList* cellIds = cell->GetPointIds();
  vtkCellArray* arraytriangle = vtkCellArray::New();
  vtkPoints* cellPts = cell->GetPoints();
  vtkIdType npts = cellPts->GetNumberOfPoints();
  std::vector<vtkIdType> cellptId(npts);
  vtkIdType iid[3];
  const vtkIdType* v_id = nullptr;
  vtkIdType *verts, v1, v2;
  vtkIdType ptId;
  vtkIdType tab_id[3];
  vtkIdType ptstriangle = 3;

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

  /* Edges Triangle*/
  vtkIdType edges[3][2] = {
    { 0, 1 },
    { 1, 2 },
    { 2, 0 },
  };
  double value, deltaScalar;
  double t, *p1, *p2;
  double v[3], x[3];
  double v_triangle[3][3];

  // To avoid false positive warning.
  tab_id[0] = 0;
  tab_id[1] = 0;
  tab_id[2] = 0;

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

  // Convert all 2d cells to triangle
  this->CellGrid(cellType, npts, cellptId.data(), arraytriangle);

  unsigned int totalnewtriangle = arraytriangle->GetNumberOfCells();
  unsigned int idtrianglenew;

  for (idtrianglenew = 0; idtrianglenew < totalnewtriangle; idtrianglenew++)
  {
    arraytriangle->GetNextCell(ptstriangle, v_id);

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

      if (!(((v[0] >= this->BoundBoxClip[0][0]) && (v[0] <= this->BoundBoxClip[0][1]) &&