Method ExchangeHTGMetadata

Filters/ParallelMPI/vtkHyperTreeGridRedistribute.cxx:462–605

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

//------------------------------------------------------------------------------
bool vtkHyperTreeGridRedistribute::ExchangeHTGMetadata()
{
  // Make sure all ranks share the same HTG metadata.
  // Metadata mismatch can happen, for instance, when we read a HTG from a .htg file in parallel;
  // all ranks have unconfigured HTGs except one rank.

  // Get minimum rank id with an initialized input HTG (correct bounds).
  // This rank will broadcast its metadata.
  int metadataSourceProcess = 0;
  double* bounds = this->InputHTG->GetBounds();
  int processInit = bounds[0] <= bounds[1] ? this->Controller->GetLocalProcessId()
                                           : std::numeric_limits<int>::max();
  this->Controller->AllReduce(&processInit, &metadataSourceProcess, 1, vtkCommunicator::MIN_OP);

  if (metadataSourceProcess == std::numeric_limits<int>::max())
  {
    vtkDebugMacro("No valid HTG found in any process");
    return false;
  }

  vtkDebugMacro("Metadata source process is " << metadataSourceProcess);

  this->OutputHTG->Initialize();

  // Exchange BranchFactor
  int branchFactor = this->InputHTG->GetBranchFactor();
  this->Controller->Broadcast(&branchFactor, 1, metadataSourceProcess);
  this->OutputHTG->SetBranchFactor(branchFactor);

  // Exchange DepthLimiter
  int depth = this->InputHTG->GetDepthLimiter();
  this->Controller->Broadcast(&depth, 1, metadataSourceProcess);
  this->OutputHTG->SetDepthLimiter(depth);

  // Exchange mask info
  int hasMask = this->InputHTG->HasMask();
  this->Controller->Broadcast(&hasMask, 1, metadataSourceProcess);
  if (hasMask)
  {
    this->OutMask = vtkSmartPointer<vtkBitArray>::New();
    this->HasMask = true;
  }

  // Exchange TransposedRootIndexing
  int transposedRoot = this->InputHTG->GetTransposedRootIndexing();
  this->Controller->Broadcast(&transposedRoot, 1, metadataSourceProcess);
  this->OutputHTG->SetTransposedRootIndexing(transposedRoot);

  // Exchange Dimensions
  int dims[3];
  this->InputHTG->GetDimensions(dims);
  this->Controller->Broadcast(dims, 3, metadataSourceProcess);
  this->OutputHTG->SetDimensions(dims);

  // Exchange Interface
  int hasInterface = this->InputHTG->GetHasInterface();
  this->Controller->Broadcast(&hasInterface, 1, metadataSourceProcess);
  this->OutputHTG->SetHasInterface(hasInterface);