Method StripGhostLayersFromGrid

Common/DataModel/vtkAMRUtilities.cxx:215–260

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

//------------------------------------------------------------------------------
vtkUniformGrid* vtkAMRUtilities::StripGhostLayersFromGrid(vtkUniformGrid* grid, int ghost[6])
{
  assert("pre: input grid is nullptr" && (grid != nullptr));

  // STEP 0: Get the grid properties, i.e., origin, dims, extent, etc.
  double origin[3];
  double spacing[3];
  int dims[3];
  int ghostedGridDims[3];

  grid->GetOrigin(origin);
  grid->GetSpacing(spacing);
  grid->GetDimensions(ghostedGridDims);
  grid->GetDimensions(dims);

  int copyExtent[6];
  grid->GetExtent(copyExtent);

  // STEP 1: Adjust origin, copyExtent, dims according to the supplied ghost
  // vector.
  for (int i = 0; i < 3; ++i)
  {
    if (ghost[i * 2] > 0)
    {
      copyExtent[i * 2] += ghost[i * 2];
      dims[i] -= ghost[i * 2];
      origin[i] += ghost[i * 2] * spacing[i];
    }
    if (ghost[i * 2 + 1] > 0)
    {
      dims[i] -= ghost[i * 2 + 1];
      copyExtent[i * 2 + 1] -= ghost[i * 2 + 1];
    }
  } // END for all dimensions

  // STEP 2: Initialize the unghosted grid
  vtkUniformGrid* myGrid = vtkUniformGrid::New();
  myGrid->Initialize();
  myGrid->SetOrigin(origin);
  myGrid->SetSpacing(spacing);
  myGrid->SetDimensions(dims);

  // STEP 3: Copy the field data within the real extent
  vtkAMRUtilities::CopyFieldsWithinRealExtent(copyExtent, grid, myGrid);
  return (myGrid);
}

//------------------------------------------------------------------------------
void vtkAMRUtilities::StripGhostLayers(