src/tetDecompositionFiniteElement/tetPolyMeshCellDecomp/tetPolyPatches/constraint/processor/processorTetPolyPatchCellDecomp.C

   1 /*---------------------------------------------------------------------------*\
   2   =========                 |
   3   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   4    \\    /   O peration     |
   5     \\  /    A nd           | Copyright held by original author
   6      \\/     M anipulation  |
   7 -------------------------------------------------------------------------------
   8 License
   9     This file is part of OpenFOAM.
  10
  11     OpenFOAM is free software; you can redistribute it and/or modify it
  12     under the terms of the GNU General Public License as published by the
  13     Free Software Foundation; either version 2 of the License, or (at your
  14     option) any later version.
  15
  16     OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
  17     ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
  18     FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License
  19     for more details.
  20
  21     You should have received a copy of the GNU General Public License
  22     along with OpenFOAM; if not, write to the Free Software Foundation,
  23     Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA
  24
  25 \*---------------------------------------------------------------------------*/
  26
  27 #include "processorTetPolyPatchCellDecomp.H"
  28 #include "tetPolyMeshCellDecomp.H"
  29 #include "globalTetPolyPatchCellDecomp.H"
  30
  31 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  32
  33 namespace Foam
  34 {
  35
  36 // * * * * * * * * * * * * * Private Member Functions  * * * * * * * * * * * //
  37
  38 template<template<class> class FaceList>
  39 labelList processorTetPolyPatchCellDecomp::calcProcLocalEdgesIndices
  40 (
  41     const PrimitivePatch<face, FaceList, const pointField&>& p
  42 ) const
  43 {
  44     if (debug)
  45     {
  46         Info<< "labelList processorTetPolyPatchCellDecomp::"
  47             << "calcProcLocalEdgesIndices(const primitivePatch& p) const : "
  48             << "calculating local edge indices"
  49             << endl;
  50     }
  51
  52     // Count number of edges in the patch
  53
  54     // Get reference to the mesh
  55     const tetPolyMeshCellDecomp& mesh = boundaryMesh().mesh();
  56
  57     // Get reference to edges of the primitive patch
  58     const edgeList& patchEdges = p.edges();
  59
  60     // Get reference to faces of the primitive patch
  61     const faceList& patchFaces = p;
  62
  63     // Get reference to faces of the primitive patch
  64     const faceList& patchLocalFaces = p.localFaces();
  65
  66     // Get reference to mesh points of the primitive patch
  67     const labelList& patchMeshPoints = p.meshPoints();
  68
  69     // Make a point rejection list to remove shared processor points
  70
  71     // Get reference to shared processor points
  72     const labelList& sharedPoints =
  73         boundaryMesh().globalPatch().meshPoints();
  74
  75     boolList acceptPoint(patchMeshPoints.size(), true);
  76
  77     forAll (patchMeshPoints, pointI)
  78     {
  79         label curP = patchMeshPoints[pointI];
  80
  81         forAll (sharedPoints, sharedI)
  82         {
  83             if (sharedPoints[sharedI] == curP)
  84             {
  85                 acceptPoint[pointI] = false;
  86                 break;
  87             }
  88         }
  89     }
  90
  91     // edges of the polyPatch
  92     label maxNEdgesInPatch = patchEdges.size();
  93
  94     // diagonal edges across faces
  95     forAll (patchFaces, faceI)
  96     {
  97         maxNEdgesInPatch += patchFaces[faceI].size() - 3;
  98     }
  99
 100     labelList localEdgeInd(maxNEdgesInPatch, -1);
 101
 102     label nEdges = 0;
 103
 104     const lduAddressing& lduAddr = mesh.lduAddr();
 105
 106     // First do the edges of the primitive patch
 107     forAll (patchEdges, edgeI)
 108     {
 109         if
 110         (
 111             acceptPoint[patchEdges[edgeI].start()]
 112          || acceptPoint[patchEdges[edgeI].end()]
 113         )
 114         {
 115             localEdgeInd[nEdges] =
 116                 lduAddr.triIndex
 117                 (
 118                     patchMeshPoints[patchEdges[edgeI].start()],
 119                     patchMeshPoints[patchEdges[edgeI].end()]
 120                 );
 121
 122             nEdges++;
 123         }
 124     }
 125
 126     // Now do the diagonal edges
 127     forAll (patchFaces, faceI)
 128     {
 129         const face& curFace = patchFaces[faceI];
 130         const face& curLocalFace = patchLocalFaces[faceI];
 131
 132         for (label pointI = 2; pointI < curFace.size() - 1; pointI++)
 133         {
 134             if
 135             (
 136                 acceptPoint[curLocalFace[0]]
 137              || acceptPoint[curLocalFace[pointI]]
 138             )
 139             {
 140                 localEdgeInd[nEdges] =
 141                     lduAddr.triIndex(curFace[0], curFace[pointI]);
 142
 143                 nEdges++;
 144             }
 145         }
 146     }
 147
 148     // Reset the size of the list
 149     localEdgeInd.setSize(nEdges);
 150
 151 #   ifdef DEBUGtetFemMatrix
 152     if (min(localEdgeInd) < 0)
 153     {
 154         FatalErrorIn
 155         (
 156             "labelList processorTetPolyPatchCellDecomp::"
 157             "calcProcLocalEdgesIndices(const primitivePatch& p) const"
 158         )   << "Problem in local edge addressing"
 159             << abort(FatalError);
 160     }
 161 #   endif
 162
 163     if (debug)
 164     {
 165         Info<< "labelList processorTetPolyPatchCellDecomp::"
 166             << "calcProcLocalEdgesIndices(const primitivePatch& p) const : "
 167             << "finished calculating local edge indices"
 168             << endl;
 169     }
 170
 171     return localEdgeInd;
 172 }
 173
 174
 175 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 176
 177 } // End namespace Foam
 178
 179 // ************************************************************************* //