Method UnstructuredExecute

Filters/General/vtkDataSetTriangleFilter.cxx:191–380

3D cells use the ordered triangulator. The ordered triangulator is used to create templates on the fly. Once the templates are created then they are used to produce the final triangulation.

// are used to produce the final triangulation.
//
void vtkDataSetTriangleFilter::UnstructuredExecute(
  vtkDataSet* dataSetInput, vtkUnstructuredGrid* output)
{
  vtkPointSet* input = static_cast<vtkPointSet*>(dataSetInput); // has to be
  vtkIdType numCells = input->GetNumberOfCells();
  vtkGenericCell* cell;
  vtkIdType newCellId, j;
  int k;
  vtkCellData* inCD = input->GetCellData();
  vtkCellData* outCD = output->GetCellData();
  vtkIdList* cellPtIds;
  vtkIdType ptId, numTets, ncells;
  int numPts, type;
  int numSimplices, dim;
  vtkIdType pts[4];
  double x[3];

  if (numCells == 0)
  {
    outCD->CopyAllocate(inCD, 0);
    output->GetPointData()->CopyAllocate(input->GetPointData(), 0);
    return;
  }

  vtkUnstructuredGrid* inUgrid = vtkUnstructuredGrid::SafeDownCast(dataSetInput);
  if (inUgrid)
  {
    int allSimplices = 1;
    for (vtkIdType cellId = 0; cellId < inUgrid->GetNumberOfCells() && allSimplices; cellId++)
    {
      switch (inUgrid->GetCellType(cellId))
      {
        case VTK_TETRA:
          break;
        case VTK_VERTEX:
        case VTK_LINE:
        case VTK_TRIANGLE:
          if (this->TetrahedraOnly)
          {
            allSimplices = 0; // don't shallowcopy need to strip non tets
          }
          break;
        default:
          allSimplices = 0;
          break;
      }
    }
    if (allSimplices)
    {
      output->ShallowCopy(input);
      return;
    }
  }

  cell = vtkGenericCell::New();
  cellPtIds = vtkIdList::New();

  // Create an array of points