Function GetShape

IO/IOSS/vtkIOSSCellGridUtilities.cxx:369–416

}

bool GetShape(Ioss::Region* region, const Ioss::GroupingEntity* group_entity,
  vtkCellAttribute::CellTypeInfo& cellShapeInfo, int timestep, vtkDGCell* meta, vtkCellGrid* grid,
  vtkIOSSUtilities::Cache* cache)
{
  (void)timestep;
  (void)group_entity; // TODO: If we ever squeeze points on a per-block basis, we must cache
                      // the nodal coords on group_entity (not nodeblock_entity) and
                      // use the vtkDGCell's cellSpec.Connectivity to subset points when
                      // generating the cache entry.

  auto nodeblock_entity = region->get_entity("nodeblock_1", Ioss::EntityType::NODEBLOCK);
  if (!nodeblock_entity)
  {
    return false;
  }
  vtkSmartPointer<vtkDataArray> cached = (cache
      ? vtkDataArray::SafeDownCast(cache->Find(nodeblock_entity, "__vtk_mesh_model_coordinates__"))
      : nullptr);
  if (cached)
  {
    vtkLogF(TRACE, "using cached mesh_model_coordinates");
  }
  else
  {
    cached = vtkIOSSUtilities::GetData(nodeblock_entity, "mesh_model_coordinates");
    cached = ChangeComponents(cached, 3);

    if (cache)
    {
      cache->Insert(group_entity, "__vtk_mesh_model_coordinates__", cached.GetPointer());
    }
  }

  grid->GetAttributes("coordinates")->AddArray(cached);

  vtkNew<vtkCellAttribute> attribute;
  attribute->Initialize("shape", "ℝ³", 3);
  cellShapeInfo.DOFSharing = "coordinates"_token; // Required for the shape attribute.
  cellShapeInfo.FunctionSpace = meta->GetCellSpec().SourceShape == vtkDGCell::Shape::Vertex
    ? "constant"_token
    : "HGRAD"_token; // Required for the shape attribute.
  cellShapeInfo.ArraysByRole["connectivity"] = meta->GetCellSpec().Connectivity;
  cellShapeInfo.ArraysByRole["values"] = cached;
  attribute->SetCellTypeInfo(meta->GetClassName(), cellShapeInfo);
  grid->AddCellAttribute(attribute);
  grid->SetShapeAttribute(attribute);
  return true;
}

VTK_ABI_NAMESPACE_END
} // namespace vtkIOSSCellGridUtilities