Function GetAMRDataSet

Examples/AMR/Cxx/Generate3DAMRDataSetWithPulse.cxx:115–189

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

}
//------------------------------------------------------------------------------
vtkOverlappingAMR* GetAMRDataSet()
{
  vtkNew<vtkOverlappingAMR> data;
  int blocksPerLevel[2] = { 1, 3 };
  double globalOrigin[3] = { -2.0, -2.0, -2.0 };
  data->Initialize(2, blocksPerLevel);
  data->SetOrigin(globalOrigin);
  data->SetGridDescription(vtkStructuredData::VTK_STRUCTURED_XYZ_GRID);

  double origin[3];
  double h[3];
  int ndim[3];

  // Root Block -- Block 0
  ndim[0] = 6;
  ndim[1] = ndim[2] = 5;
  h[0] = h[1] = h[2] = 1.0;
  origin[0] = origin[1] = origin[2] = -2.0;
  int blockId = 0;
  int level = 0;
  vtkSmartPointer<vtkUniformGrid> root;
  root.TakeReference(AMRCommon::GetGrid(origin, h, ndim));
  vtkAMRBox box(origin, ndim, h, data->GetOrigin(), data->GetGridDescription());
  AttachPulseToGrid(root);
  data->SetAMRBox(level, blockId, box);
  data->SetDataSet(level, blockId, root);

  // Block 1
  ndim[0] = 3;
  ndim[1] = ndim[2] = 5;
  h[0] = h[1] = h[2] = 0.5;
  origin[0] = origin[1] = origin[2] = -2.0;
  blockId = 0;
  level = 1;
  vtkSmartPointer<vtkUniformGrid> grid1;
  root.TakeReference(AMRCommon::GetGrid(origin, h, ndim));
  AttachPulseToGrid(grid1);
  vtkAMRBox box1(origin, ndim, h, data->GetOrigin(), data->GetGridDescription());
  data->SetAMRBox(level, blockId, box1);
  data->SetDataSet(level, blockId, grid1);

  // Block 2
  ndim[0] = 3;
  ndim[1] = ndim[2] = 5;
  h[0] = h[1] = h[2] = 0.5;
  origin[0] = 0.0;
  origin[1] = origin[2] = -1.0;
  blockId = 1;
  level = 1;
  vtkSmartPointer<vtkUniformGrid> grid2;
  grid2.TakeReference(AMRCommon::GetGrid(origin, h, ndim));
  AttachPulseToGrid(grid2);
  vtkAMRBox box2(origin, ndim, h, data->GetOrigin(), data->GetGridDescription());
  data->SetAMRBox(level, blockId, box2);
  data->SetDataSet(level, blockId, grid2);

  // Block 3
  ndim[0] = 3;