applications/solvers/solidMechanics/elasticNonLinULSolidFoam/findGlobalFaceZones.H

   1 labelList globalFaceZones;
   2
   3 {
   4     SLList<label> globalFaceZonesSet;
   5
   6     const faceZoneMesh& faceZones = mesh.faceZones();
   7
   8     forAll(faceZones, zoneI)
   9     {
  10         const faceZone& curFaceZone = faceZones[zoneI];
  11
  12         forAll(curFaceZone, faceI)
  13         {
  14             // if unused face exist
  15             if (curFaceZone[faceI] >= mesh.nFaces())
  16             {
  17                 globalFaceZonesSet.insert(zoneI);
  18                 break;
  19             }
  20         }
  21     }
  22
  23     globalFaceZones = labelList(globalFaceZonesSet);
  24 }
  25
  26 labelListList globalToLocalFaceZonePointMap
  27 (
  28     globalFaceZones.size()
  29 );
  30
  31 forAll(globalFaceZones, zoneI)
  32 {
  33     label curZoneID = globalFaceZones[zoneI];
  34
  35     labelList curMap(mesh.faceZones()[curZoneID]().nPoints(), -1);
  36
  37     vectorField fzGlobalPoints =
  38         mesh.faceZones()[curZoneID]().localPoints();
  39
  40     //- set all slave points to zero because only the master order is used
  41     if(!Pstream::master())
  42     {
  43         fzGlobalPoints *= 0.0;
  44     }
  45
  46     //- pass points to all procs
  47     reduce(fzGlobalPoints, sumOp<vectorField>());
  48
  49     //- now every proc has the master's list of FZ points
  50     //- every proc must now find the mapping from their local FZ points to
  51     //- the global FZ points
  52
  53     const vectorField& fzLocalPoints =
  54         mesh.faceZones()[curZoneID]().localPoints();
  55
  56     const edgeList& fzLocalEdges =
  57         mesh.faceZones()[curZoneID]().edges();
  58
  59     const labelListList& fzPointEdges =
  60         mesh.faceZones()[curZoneID]().pointEdges();
  61
  62     scalarField minEdgeLength(fzLocalPoints.size(), GREAT);
  63
  64     forAll(minEdgeLength, pI)
  65     {
  66         const labelList& curPointEdges = fzPointEdges[pI];
  67
  68         forAll(curPointEdges, eI)
  69         {
  70             scalar Le = fzLocalEdges[curPointEdges[eI]].mag(fzLocalPoints);
  71             if (Le < minEdgeLength[pI])
  72             {
  73                 minEdgeLength[pI] = Le;
  74             }
  75         }
  76     }
  77
  78     forAll(fzGlobalPoints, globalPointI)
  79     {
  80 //         scalar minDist = GREAT;
  81
  82         forAll(fzLocalPoints, procPointI)
  83         {
  84             scalar curDist =
  85                 mag
  86                 (
  87                     fzLocalPoints[procPointI]
  88                   - fzGlobalPoints[globalPointI]
  89                 );
  90
  91 //             if (curDist < minDist)
  92 //             {
  93 //                 minDist = curDist;
  94 //             }
  95
  96             if (curDist < 1e-4*minEdgeLength[procPointI])
  97             {
  98                 curMap[globalPointI] = procPointI;
  99                 break;
 100             }
 101         }
 102
 103 //         if (curMap[globalPointI] == -1)
 104 //         {
 105 //             Pout << "minDist: " << minDist << endl;
 106 //         }
 107     }
 108
 109     forAll(curMap, globalPointI)
 110     {
 111         if (curMap[globalPointI] == -1)
 112         {
 113             FatalErrorIn(args.executable())
 114                 << "local to global face zone point map is not correct"
 115                     << abort(FatalError);
 116         }
 117     }
 118
 119     globalToLocalFaceZonePointMap[zoneI] = curMap;
 120 }