Method StripGhostLayers

Common/DataModel/vtkAMRUtilities.cxx:263–342

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

//------------------------------------------------------------------------------
void vtkAMRUtilities::StripGhostLayers(
  vtkOverlappingAMR* ghostedAMRData, vtkOverlappingAMR* strippedAMRData)
{
  assert("pre: input AMR data is nullptr" && (ghostedAMRData != nullptr));
  assert("pre: outputAMR data is nullptr" && (strippedAMRData != nullptr));
  double spacing[3];

  if (!vtkAMRUtilities::HasPartiallyOverlappingGhostCells(ghostedAMRData))
  {
    strippedAMRData->CompositeShallowCopy(ghostedAMRData);
    return;
  }

  // TODO: At some point we should check for overlapping cells within the
  // same level, e.g., consider a level 0 with 2 abutting blocks that is
  // ghosted by N !!!!
  std::vector<unsigned int> blocksPerLevel(ghostedAMRData->GetNumberOfLevels());
  for (unsigned int i = 0; i < blocksPerLevel.size(); i++)
  {
    blocksPerLevel[i] = ghostedAMRData->GetNumberOfBlocks(i);
  }
  strippedAMRData->Initialize(blocksPerLevel);
  strippedAMRData->SetOrigin(ghostedAMRData->GetOrigin());
  strippedAMRData->SetGridDescription(ghostedAMRData->GetGridDescription());

  ghostedAMRData->GetSpacing(0, spacing);
  strippedAMRData->SetSpacing(0, spacing);
  unsigned int dataIdx = 0;
  for (; dataIdx < ghostedAMRData->GetNumberOfBlocks(0); ++dataIdx)
  {
    vtkCartesianGrid* cg = ghostedAMRData->GetDataSetAsCartesianGrid(0, dataIdx);
    // XXX: The vtkCartesianGrid could be used instead of using spacing, which would provide
    // an alternative implementation when spacing is not available, which can happen when not using
    // vtkUniformGrid.
    vtkUniformGrid* grid = vtkUniformGrid::SafeDownCast(cg);
    if (cg && !grid)
    {
      // this method supports nullptr data but not cartesian grid.
      vtkErrorWithObjectMacro(nullptr, "Cannot StripGhostLayers for an AMR of non vtkUniformGrid");
      return;
    }
    const vtkAMRBox& box = ghostedAMRData->GetAMRBox(0, dataIdx);
    strippedAMRData->SetAMRBox(0, dataIdx, box);
    strippedAMRData->SetDataSet(0, dataIdx, grid);
  } // END for all data at level 0

  int ghost[6];
  unsigned int levelIdx = 1;
  for (; levelIdx < ghostedAMRData->GetNumberOfLevels(); ++levelIdx)
  {
    dataIdx = 0;
    ghostedAMRData->GetSpacing(levelIdx, spacing);
    strippedAMRData->SetSpacing(levelIdx, spacing);
    for (; dataIdx < ghostedAMRData->GetNumberOfBlocks(levelIdx); ++dataIdx)
    {
      vtkUniformGrid* grid =
        vtkUniformGrid::SafeDownCast(ghostedAMRData->GetDataSetAsCartesianGrid(levelIdx, dataIdx));
      int r = ghostedAMRData->GetRefinementRatio(levelIdx);