Method GetNodeSet

IO/IOSS/vtkIOSSCellGridReaderInternal.cxx:325–407

}

std::vector<vtkSmartPointer<vtkCellGrid>> vtkIOSSCellGridReaderInternal::GetNodeSet(
  const std::string& blockName, vtkIOSSReader::EntityType vtk_entity_type,
  const DatabaseHandle& handle, int timestep, vtkIOSSCellGridReader* self)
{
  (void)self;
  (void)vtk_entity_type;
  auto region = this->GetRegion(handle);
  auto group_entity = region->get_entity(blockName, Ioss::EntityType::NODESET);
  if (!group_entity)
  {
    vtkErrorWithObjectMacro(self, "No group entity for node-set \"" << blockName << "\".");
    return {};
  }
  auto grid = vtkSmartPointer<vtkCellGrid>::New();
  auto meta = vtkCellMetadata::NewInstance("vtkDGVert"_token, grid);
  auto* dg = vtkDGCell::SafeDownCast(meta);
  if (!meta || !dg)
  {
    vtkErrorWithObjectMacro(
      self, "Could not create metadata for node-set \"" << blockName << "\".");
    return {};
  }
  if (!grid->AddCellMetadata(meta))
  {
    vtkErrorWithObjectMacro(
      self, "Could not add metadata for node-set \"" << blockName << "\" to grid.");
    return {};
  }
  // Fetch the IDs of the file-global points included in the node-set,
  // offsetting by -1 so they are 0-indexed:
  auto transform = Ioss::TransformFactory::create("offset");
  transform->set_property("offset", -1);
  auto ids_raw = vtkIOSSUtilities::GetData(group_entity, "ids_raw", transform);
  ids_raw->SetNumberOfComponents(1);

  // Add the ID array to a vtkDataSetAttributes instance corresponding to
  // the number of cells of the node-set. Since a separate cell-grid holds
  // each node-set, we use the name of the cell type ("vtkDGVert") for the
  // array group:
  auto* cellGroup = grid->GetAttributes(dg->GetClassName());
  cellGroup->AddArray(ids_raw);
  dg->GetCellSpec().Connectivity = ids_raw;
  dg->GetCellSpec().SourceShape = vtkDGCell::Shape::Vertex;
  dg->GetCellSpec().Blanked = false;

  // From the shape of cells in the block, the connectivity size, and the order,
  // we need to infer vtkDGCell::CellTypeInfo data (FunctionSpace, Basis, Order).
  vtkCellAttribute::CellTypeInfo cellShapeInfo;
  cellShapeInfo.FunctionSpace = "constant"_token;
  cellShapeInfo.Basis = "C"_token;
  cellShapeInfo.Order = 0;

  // Read node coordinates as the shape attribute.
  // This must always be a "CG" (continuous) attribute.
  vtkIOSSCellGridUtilities::GetShape(
    region, group_entity, cellShapeInfo, timestep, dg, grid, &this->Cache);
  // Apply displacements before reading other cell-attributes as
  // computing the range of HDIV/HCURL attributes **must** use