Method ExtractRegion

Filters/AMR/vtkAMRResampleFilter.cxx:817–853

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

//------------------------------------------------------------------------------
void vtkAMRResampleFilter::ExtractRegion(
  vtkOverlappingAMR* amrds, vtkMultiBlockDataSet* mbds, vtkOverlappingAMR* vtkNotUsed(metadata))
{

  assert("pre: input AMR data-structure is nullptr" && (amrds != nullptr));
  assert("pre: resampled grid should not be nullptr" && (mbds != nullptr));

  //  std::cout << "NumBlocks: " << this->ROI->GetNumberOfBlocks() << std::endl;
  //  std::cout << "NumProcs: "  << this->Controller->GetNumberOfProcesses() << std::endl;
  //  std::cout.flush();

  assert("pre: NumProcs must be less than or equal to NumBlocks" &&
    (!this->Controller ||
      (static_cast<int>(this->ROI->GetNumberOfBlocks()) <=
        this->Controller->GetNumberOfProcesses())));

  mbds->SetNumberOfBlocks(this->ROI->GetNumberOfBlocks());
  for (unsigned int block = 0; block < this->ROI->GetNumberOfBlocks(); ++block)
  {
    if (this->CheckAbort())
    {
      break;
    }
    if (this->IsRegionMine(block))
    {
      vtkUniformGrid* grid = vtkUniformGrid::New();
      grid->ShallowCopy(this->ROI->GetBlock(block));
      this->TransferSolution(grid, amrds);
      mbds->SetBlock(block, grid);
      grid->Delete();
    }
    else
    {
      mbds->SetBlock(block, nullptr);
    }
  } // END for all blocks
}

//------------------------------------------------------------------------------
void vtkAMRResampleFilter::ComputeAMRBlocksToLoad(vtkOverlappingAMR* metadata)