Method ProcessTrees

Filters/HyperTree/vtkHyperTreeGridExtractGhostCells.cxx:31–111

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

//------------------------------------------------------------------------------
int vtkHyperTreeGridExtractGhostCells::ProcessTrees(
  vtkHyperTreeGrid* input, vtkDataObject* outputDO)
{
  // Downcast output data object to hyper tree grid
  vtkHyperTreeGrid* output = vtkHyperTreeGrid::SafeDownCast(outputDO);
  if (!output)
  {
    vtkErrorMacro("Incorrect type of output: " << outputDO->GetClassName());
    return 0;
  }

  output->ShallowCopy(input);

  // Retrieve and copy input mask if it exists
  this->InMask = input->HasMask() ? input->GetMask() : nullptr;
  if (this->InMask)
  {
    this->OutMask->DeepCopy(this->InMask);
  }
  else
  {
    this->OutMask->SetNumberOfTuples(output->GetNumberOfCells());
  }
  output->SetMask(this->OutMask);

  // Retrieve ghost array
  if (input->HasAnyGhostCells())
  {
    this->InGhost = input->GetGhostCells();
  }
  else
  {
    vtkWarningMacro(<< "Input does not have a ghost array. Output HTG will be empty.");
  }

  // Iterate over
  vtkIdType inIndex = 0;
  vtkHyperTreeGrid::vtkHyperTreeGridIterator it;
  output->InitializeTreeIterator(it);
  vtkNew<vtkHyperTreeGridNonOrientedCursor> outCursor;
  while (it.GetNextTree(inIndex))
  {
    if (this->CheckAbort())
    {
      break;
    }

    output->InitializeNonOrientedCursor(outCursor, inIndex, true);

    if (!this->InGhost)
    {
      // Input has no ghost cell, mask the whole tree
      this->OutMask->InsertValue(outCursor->GetGlobalNodeIndex(), 1);
    }
    else
    {
      this->RecursivelyMaskNonGhost(outCursor);
    }