Method Model

IO/IOSS/vtkIOSSModel.cxx:1423–1509

  }

  void Model(Ioss::Region& region) const override
  {
    for (const auto& element : this->ElementCounts)
    {
      const int64_t elementCount = element.second;
      const unsigned char vtk_cell_type = element.first;

      const Ioss::ElementTopology* elementTopology = nullptr;
      try
      {
        elementTopology = vtkIOSSUtilities::GetElementTopology(vtk_cell_type);
      }
      catch (std::runtime_error&)
      {
        continue;
      }
      const auto& elementType = elementTopology->name();
      const int nodeCount = elementTopology->number_nodes();
      const auto blockName = this->GetSubElementBlockInfo(vtk_cell_type, elementType).second;

      auto* entityBlock = this->GetEntity(region, blockName);

      std::vector<int> orderingTransformation;
      const bool needsIdsTransformation =
        NeedsIdsTransformation(vtk_cell_type, orderingTransformation);
      // populate ids.
      std::vector<int32_t> elementIds; // these are global IDs.
      elementIds.reserve(elementCount);

      std::vector<int32_t> connectivity;
      connectivity.reserve(elementCount * nodeCount);

      const int32_t gidOffset = this->Writer->GetOffsetGlobalIds() ? 1 : 0;
      const bool removeGhosts = this->Writer->GetRemoveGhosts();
      for (auto& ds : this->DataSets)
      {
        vtkUnsignedCharArray* ghost = ds->GetCellGhostArray();
        auto* gids = vtkIdTypeArray::SafeDownCast(ds->GetCellData()->GetGlobalIds());
        auto* pointGIDs = vtkIdTypeArray::SafeDownCast(ds->GetPointData()->GetGlobalIds());

        vtkNew<vtkIdList> tempCellPointIds;
        for (vtkIdType cc = 0, max = ds->GetNumberOfCells(); cc < max; ++cc)
        {
          const bool process = !removeGhosts || !ghost || ghost->GetValue(cc) == 0;
          if (process && ds->GetCellType(cc) == vtk_cell_type)
          {
            elementIds.push_back(gidOffset + gids->GetValue(cc));

            vtkIdType numPts;
            vtkIdType const* cellPoints;
            ds->GetCellPoints(cc, numPts, cellPoints, tempCellPointIds);
            assert(numPts == nodeCount);

            if (!needsIdsTransformation)
            {
              // map cell's point to global IDs for those points.
              std::transform(cellPoints, cellPoints + numPts, std::back_inserter(connectivity),
                [&](vtkIdType ptid) { return gidOffset + pointGIDs->GetValue(ptid); });