Method ComputeAMRBlocksToLoad

Filters/AMR/vtkAMRResampleFilter.cxx:856–890

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

//------------------------------------------------------------------------------
void vtkAMRResampleFilter::ComputeAMRBlocksToLoad(vtkOverlappingAMR* metadata)
{
  assert("pre: metadata is nullptr" && (metadata != nullptr));

  this->BlocksToLoad.clear();

  unsigned int maxLevelToLoad = 0;
  if (this->LevelOfResolution < static_cast<int>(metadata->GetNumberOfLevels()))
  {
    maxLevelToLoad = this->LevelOfResolution + 1;
  }
  else
  {
    maxLevelToLoad = metadata->GetNumberOfLevels();
  }

  unsigned int level = 0;
  for (; level < maxLevelToLoad; ++level)
  {
    unsigned int dataIdx = 0;
    for (; dataIdx < metadata->GetNumberOfBlocks(level); ++dataIdx)
    {
      double grd[6];
      metadata->GetBounds(level, dataIdx, grd);
      if (this->IsBlockWithinBounds(grd))
      {
        this->BlocksToLoad.push_back(metadata->GetAbsoluteBlockIndex(level, dataIdx));
      } // END check if the block is within the bounds of the ROI
    }   // END for all data
  }     // END for all levels

  std::sort(this->BlocksToLoad.begin(), this->BlocksToLoad.end());
  std::cerr << "Number Levels Loaded = " << maxLevelToLoad
            << " Number of Blocks = " << BlocksToLoad.size() << "\n";
}

//------------------------------------------------------------------------------
void vtkAMRResampleFilter::GetDomainParameters(vtkOverlappingAMR* amr, double domainMin[3],