Skip to content
Snippets Groups Projects

ENH: Adds basic information for faceZone and cellZones during checkMesh (fixes #560)

Merged Prashant Sonakar requested to merge feature-zoneInformation into develop
3 unresolved threads

zone information

  • name
  • points
  • boundingBox

@andy : Please confirm if parallel OK

Edited by Prashant Sonakar

Merge request reports

Loading
Loading

Activity

Filter activity
  • Approvals
  • Assignees & reviewers
  • Comments (from bots)
  • Comments (from users)
  • Commits & branches
  • Edits
  • Labels
  • Lock status
  • Mentions
  • Merge request status
  • Tracking
589 Pout.setf(ios_base::left);
590
591 const faceZoneMesh& faceZones = mesh.faceZones();
592
593 if (faceZones.size())
594 {
595 Info<< " "
596 << setw(20) << "FaceZone"
597 << setw(9) << "Faces"
598 << setw(9) << "Points"
599 << setw(13) << "BoundingBox" <<endl;
600
601 const pointField& points = mesh.points();
602 const faceList& faces = mesh.faces();
603
604 DynamicList<point> localPoints;
  • 614 const pointField facePoints(f.points(points));
    615
    616 forAll(facePoints, pointI)
    617 {
    618 const point& pt = facePoints[pointI];
    619 localPoints.append(pt);
    620 }
    621 }
    622 boundBox bb(localPoints, true);
    623
    624 Info<< " "
    625 << setw(20) << fZone.name()
    626 << setw(9) << returnReduce(fZone.size(), sumOp<label>())
    627 << setw(9)
    628 << returnReduce(localPoints.size(), sumOp<label>())
    629 << setw(3) << bb << endl;
  • 662 localPoints.setCapacity(10*faces.size());
    663
    664 forAll(cZone, i)
    665 {
    666 const label cellI = cZone[i];
    667 const cell& c = cells[cellI];
    668 const pointField cellPoints(c.points(faces, points));
    669
    670 forAll(cellPoints, pointI)
    671 {
    672 const point& pt = cellPoints[pointI];
    673 localPoints.append(pt);
    674 }
    675 }
    676 boundBox bb(localPoints, true);
    677
    • PackedBoolList isZonePoint(mesh.nPoints());
      
      forAll(cFaces, cFacei)
      {
          const face& f = mesh.faces()[cFaces[cFacei]]);
          forAll(f, fp)
          {
              if (isZonePoint.set(f[fp]))
              {
                  bb.add(mesh.points()[f[fp]]))
              }
           }
      }
      
      returnRedue(iszonePoint.count(), sumOp<label>());
    • Please register or sign in to reply
  • added 1 commit

    • cea20ae1 - ENH: checkMesh: output information about zones

    Compare with previous version

  • mentioned in commit a74d9a0a

  • Please register or sign in to reply
    Loading