Function CloneUnstructuredGrid

Parallel/DIY/vtkDIYGhostUtilities.cxx:3810–3880

---------------------------------------------------------------------------- * We're doing a homebrewed shallow copy because we do not want to share any pointer with the input, * which is the case for unstructured grid cell connectivity information. */

 * which is the case for unstructured grid cell connectivity information.
 */
void CloneUnstructuredGrid(
  vtkUnstructuredGrid* ug, vtkUnstructuredGrid* clone, ::UnstructuredGridInformation& info)
{
  ::CloneDataObject(ug, clone);
  ::ClonePointData(ug, clone, info);
  ::ClonePoints(ug, clone, info);
  ::CloneCellData(ug, clone, info);

  if (!ug->GetPolyhedronFaces() && clone->GetPolyhedronFaces() &&
    clone->GetPolyhedronFaceLocations())
  {
    // initialize empty polyhedron when none are present but they are required for output.

    // use generic cell array API since speed is not that important for this rare case.
    auto polyhedronFaceLocations = clone->GetPolyhedronFaceLocations();
    for (vtkIdType pos = 0; pos < polyhedronFaceLocations->GetNumberOfCells() + 1; ++pos)
    {
      polyhedronFaceLocations->SetOffset(pos, 0);
    }
    // use generic cell array API since speed is not that important for this rare case.
    auto polyhedronFaces = clone->GetPolyhedronFaces();
    for (vtkIdType pos = 0; pos < polyhedronFaces->GetNumberOfCells() + 1; ++pos)
    {
      polyhedronFaces->SetOffset(pos, 0);
    }
  }

  if (vtkIdList* redirectionMap = info.OutputToInputCellIdRedirectionMap)
  {
    ::DeepCopyCells(
      ug->GetCells(), clone->GetCells(), redirectionMap, info.InputToOutputPointIdRedirectionMap);
    ug->GetCellTypes()->GetTuples(redirectionMap, clone->GetCellTypes());

    vtkCellArray* ugFaceLocations = ug->GetPolyhedronFaceLocations();
    if (clone->GetPolyhedronFaceLocations() && ugFaceLocations &&
      ugFaceLocations->GetNumberOfCells())
    {
      ::DeepCopyPolyhedrons(ug, clone, info);
    }
  }
  else
  {
    vtkCellArray* ugCellArray = ug->GetCells();
    vtkCellArray* cloneCellArray = clone->GetCells();
    vtkDataArray* ugConnectivity = ugCellArray->GetConnectivityArray();
    vtkDataArray* ugOffsets = ugCellArray->GetOffsetsArray();

    ugConnectivity->GetTuples(
      0, ugConnectivity->GetNumberOfTuples() - 1, cloneCellArray->GetConnectivityArray());
    ugOffsets->GetTuples(0, ugOffsets->GetNumberOfTuples() - 1, cloneCellArray->GetOffsetsArray());
    ug->GetCellTypes()->GetTuples(0, ug->GetNumberOfCells() - 1, clone->GetCellTypes());

    vtkCellArray* ugFaces = ug->GetPolyhedronFaces();

    if (clone->GetPolyhedronFaces() && ugFaces && ugFaces->GetNumberOfCells())
    {
      vtkDataArray* ugFaceConnectivity = ug->GetPolyhedronFaces()->GetConnectivityArray();
      vtkDataArray* ugFaceOffsets = ug->GetPolyhedronFaces()->GetOffsetsArray();