NEM::MSH::geoMeshBase::GeoMesh Struct Reference

Detailed Description

Definition at line 418 of file geoMeshBase.H.

Public Member Functions
void	findSide2OrigCell ()
	For each edge/polygon in sideSet, find the corresponding mesh element in mesh. More...
int	getDimension () const
	Get dimension of mesh. More...

Public Attributes
vtkSmartPointer< vtkUnstructuredGrid >	mesh
std::string	geo
std::string	link
SideSet	sideSet

Member Function Documentation

findSide2OrigCell()

void NEM::MSH::geoMeshBase::GeoMesh::findSide2OrigCell

(

)

Does nothing if "OrigCellIds"/"CellFaceIds" arrays already present. Sets the sideSet's "OrigCellIds"/"CellFaceIds" arrays. Orientation of side only checked if multiple candidate original cells. Both "OrigCellIds" and "CellFaceIds" are set to -1 for any cells in the sideSet that can't be matched.

Definition at line 314 of file geoMeshBase.C.

References mesh.

                                           {
  if (sideSet.sides && !sideSet.getOrigCellArr() && !sideSet.getCellFaceArr()) {
    vtkIdType numSides = sideSet.sides->GetNumberOfCells();
    auto entArr = sideSet.getGeoEntArr();
    vtkNew<vtkIdTypeArray> origCellArr;
    origCellArr->SetNumberOfComponents(2);
    origCellArr->SetNumberOfTuples(numSides);
    origCellArr->FillTypedComponent(0, -1);
    origCellArr->FillTypedComponent(1, -1);
    vtkNew<vtkIntArray> cellFaceArr;
    cellFaceArr->SetNumberOfComponents(2);
    cellFaceArr->SetNumberOfTuples(numSides);
    cellFaceArr->FillTypedComponent(0, -1);
    cellFaceArr->FillTypedComponent(1, -1);
    // Using points->cell, intersect the set of cells that have the same points
    // as the side cell
    mesh->BuildLinks();
    auto pt2Cell = mesh->GetCellLinks();
    for (vtkIdType i = 0; i < sideSet.sides->GetNumberOfCells(); ++i) {
      auto sidePoints = sideSet.sides->GetCell(i)->GetPointIds();
      auto cellList = pt2Cell->GetLink(sidePoints->GetId(0));
      std::vector<vtkIdType> candidates{cellList.cells,
                                        cellList.cells + cellList.ncells};
      std::sort(candidates.begin(), candidates.end());
      for (auto j = 1; j < sidePoints->GetNumberOfIds(); ++j) {
        std::vector<vtkIdType> newCandidates{};
        cellList = pt2Cell->GetLink(sidePoints->GetId(j));
        for (auto iterCell = cellList.cells;
             iterCell < cellList.cells + cellList.ncells; ++iterCell) {
          auto findIter =
              std::lower_bound(candidates.begin(), candidates.end(), *iterCell);
          if (findIter != candidates.end() && *iterCell == *findIter) {
            newCandidates.emplace_back(*iterCell);
          }
        }
        candidates = std::move(newCandidates);
        std::sort(candidates.begin(), candidates.end());
      }
      // For each candidate cell, check each edge/face for a match
      bool foundMatch = false;
      if (candidates.size() == 1) {
        // If there is only one candidate (side is on boundary, not material
        // interface), check by sorting points - note no check for orientation
        std::vector<vtkIdType> sortedSidePoints(sidePoints->GetNumberOfIds());
        auto pointsBegin = sidePoints->GetPointer(0);
        std::partial_sort_copy(
            pointsBegin, pointsBegin + sidePoints->GetNumberOfIds(),
            sortedSidePoints.begin(), sortedSidePoints.end());
        auto candOrigCell = mesh->GetCell(candidates.at(0));
        auto dim = candOrigCell->GetCellDimension();
        for (int j = 0; j < (dim == 3 ? candOrigCell->GetNumberOfFaces()
                                      : candOrigCell->GetNumberOfEdges());
             ++j) {
          auto candSide =
              dim == 3 ? candOrigCell->GetFace(j) : candOrigCell->GetEdge(j);
          std::vector<vtkIdType> sortedCandPoints(sortedSidePoints.size());
          auto candPointsBegin = candSide->GetPointIds()->GetPointer(0);
          // Note use of sidePoints->GetNumberOfIds() in case, eg, candSide is
          // a vtkQuadraticTriangle, but side is a vtkTriangle
          auto candPointsEnd = candPointsBegin + sidePoints->GetNumberOfIds();
          std::partial_sort_copy(candPointsBegin, candPointsEnd,
                                 sortedCandPoints.begin(),
                                 sortedCandPoints.end());
          if (sortedSidePoints == sortedCandPoints) {
            foundMatch = true;
            origCellArr->SetTypedComponent(i, 0, candidates.at(0));
            cellFaceArr->SetTypedComponent(i, 0, j);
            break;
          }
        }
      } else {
        // If there are two candidates, match the orientation of the side
        assert(candidates.size() == 2);
        for (const auto &candIdx : candidates) {
          auto candOrigCell = mesh->GetCell(candIdx);
          auto dim = candOrigCell->GetCellDimension();
          bool foundSide = false;
          for (int j = 0; j < (dim == 3 ? candOrigCell->GetNumberOfFaces()
                                        : candOrigCell->GetNumberOfEdges()) &&
                          !foundSide;
               ++j) {
            auto candSide =
                dim == 3 ? candOrigCell->GetFace(j) : candOrigCell->GetEdge(j);
            auto candPointsBegin = candSide->GetPointIds()->GetPointer(0);
            // Note use of sidePoints->GetNumberOfIds() in case, eg, candSide is
            // a vtkQuadraticTriangle, but side is a vtkTriangle
            auto candPointsEnd = candPointsBegin + sidePoints->GetNumberOfIds();
            auto findFirstPoint =
                std::find(candPointsBegin, candPointsEnd, sidePoints->GetId(0));
            if (findFirstPoint != candPointsEnd) {
              std::vector<vtkIdType> rotatedCandPoints(
                  sidePoints->GetNumberOfIds());
              std::rotate_copy(candPointsBegin, findFirstPoint, candPointsEnd,
                               rotatedCandPoints.begin());
              if (std::equal(rotatedCandPoints.begin(), rotatedCandPoints.end(),
                             sidePoints->GetPointer(0))) {
                foundMatch = true;
                foundSide = true;
                origCellArr->SetTypedComponent(i, 0, candIdx);
                cellFaceArr->SetTypedComponent(i, 0, j);
              } else if (std::equal(rotatedCandPoints.rbegin(),
                                    rotatedCandPoints.rend(),
                                    sidePoints->GetPointer(0))) {
                foundSide = true;
                origCellArr->SetTypedComponent(i, 1, candIdx);
                cellFaceArr->SetTypedComponent(i, 1, j);
              }
            }
          }
          if (foundMatch) { break; }
        }
      }
      if (!foundMatch) {
        std::cerr << "No cell in mesh found to correspond to side set cell "
                  << i << ".\n";
      }
    }
    sideSet.setOrigCellArr(origCellArr);
    sideSet.setCellFaceArr(cellFaceArr);
  }
}

getDimension()

int NEM::MSH::geoMeshBase::GeoMesh::getDimension

(

)

const

All cells in mesh assumed to have same dimension.

Returns

dimension of mesh

Definition at line 436 of file geoMeshBase.C.

References mesh.

Referenced by NEM::MSH::exoGeoMesh::stitch().

                                           {
  if (this->mesh->GetNumberOfCells() > 0) {
    return this->mesh->GetCell(0)->GetCellDimension();
  } else {
    return -1;
  }
}

Member Data Documentation

geo

std::string NEM::MSH::geoMeshBase::GeoMesh::geo

Definition at line 421 of file geoMeshBase.H.

Referenced by NEM::MSH::exoGeoMesh::addCellsToBlock(), NEM::MSH::exoGeoMesh::addSideSet(), NEM::MSH::gmshGeoMesh::GM2gmsh(), NEM::MSH::exoGeoMesh::reconstructGeo(), NEM::MSH::geoMeshBase::reconstructGeo(), NEM::MSH::exoGeoMesh::scaleNodes(), and NEM::MSH::exoGeoMesh::stitch().

link

std::string NEM::MSH::geoMeshBase::GeoMesh::link

Definition at line 424 of file geoMeshBase.H.

Referenced by NEM::MSH::exoGeoMesh::addCellsToBlock(), NEM::MSH::exoGeoMesh::addSideSet(), NEM::MSH::gmshGeoMesh::GM2gmsh(), NEM::MSH::oshGeoMesh::GM2osh(), NEM::MSH::smeshGeoMesh::GMToSMESH(), NEM::MSH::oshGeoMesh::reconstructGeo(), NEM::MSH::exoGeoMesh::reconstructGeo(), NEM::MSH::geoMeshBase::reconstructGeo(), NEM::MSH::exoGeoMesh::resetNative(), NEM::MSH::exoGeoMesh::scaleNodes(), NEM::MSH::exoGeoMesh::setPhysGrpPropertyName(), and NEM::MSH::exoGeoMesh::stitch().

mesh

vtkSmartPointer<vtkUnstructuredGrid> NEM::MSH::geoMeshBase::GeoMesh::mesh

Definition at line 419 of file geoMeshBase.H.

Referenced by NEM::MSH::exoGeoMesh::addCellsToBlock(), NEM::MSH::exoGeoMesh::addElemBlock(), NEM::MSH::exoGeoMesh::addNodeSet(), NEM::MSH::exoGeoMesh::addSideSet(), NEM::MSH::exoGeoMesh::exoGeoMesh(), NEM::MSH::geoMeshBase::getCellDataArrayCopy(), NEM::MSH::exoGeoMesh::getElemBlock(), NEM::MSH::exoGeoMesh::getElemBlockId(), NEM::MSH::geoMeshBase::getFieldDataArrayCopy(), NEM::MSH::geoMeshBase::getPointDataArrayCopy(), NEM::MSH::foamGeoMesh::GM2foam(), NEM::MSH::gmshGeoMesh::GM2gmsh(), NEM::MSH::oshGeoMesh::GM2osh(), NEM::MSH::smeshGeoMesh::GMToSMESH(), NEM::MSH::geoMeshBase::mergeGeoMesh(), NEM::MSH::exoGeoMesh::reassignCells(), NEM::MSH::oshGeoMesh::reconstructGeo(), NEM::MSH::exoGeoMesh::reconstructGeo(), NEM::MSH::geoMeshBase::reconstructGeo(), NEM::MSH::geoMeshBase::report(), NEM::MSH::exoGeoMesh::resetElemBlocks(), NEM::MSH::gmshGeoMesh::resetNative(), NEM::MSH::inpGeoMesh::resetNative(), NEM::MSH::exoGeoMesh::resetNative(), NEM::MSH::geoMeshBase::setGeoMesh(), NEM::MSH::exoGeoMesh::setPhysGrpPropertyName(), NEM::MSH::exoGeoMesh::stitch(), and NEM::MSH::vtkGeoMesh::write().

sideSet

SideSet NEM::MSH::geoMeshBase::GeoMesh::sideSet

Definition at line 426 of file geoMeshBase.H.

Referenced by NEM::MSH::exoGeoMesh::addCellsToBlock(), NEM::MSH::exoGeoMesh::addSideSet(), NEM::MSH::foamGeoMesh::GM2foam(), NEM::MSH::gmshGeoMesh::GM2gmsh(), NEM::MSH::oshGeoMesh::GM2osh(), NEM::MSH::smeshGeoMesh::GMToSMESH(), NEM::MSH::exoGeoMesh::reassignCells(), NEM::MSH::oshGeoMesh::reconstructGeo(), NEM::MSH::exoGeoMesh::reconstructGeo(), NEM::MSH::geoMeshBase::reconstructGeo(), NEM::MSH::inpGeoMesh::resetNative(), NEM::MSH::exoGeoMesh::scaleNodes(), NEM::MSH::geoMeshBase::setGeoMesh(), and NEM::MSH::exoGeoMesh::stitch().

---

The documentation for this struct was generated from the following files:

• geoMeshBase.H
• geoMeshBase.C