src/tetDecompositionFiniteElement/tetPolyMeshCellDecomp/calcTetPolyMeshCellDecompAddressing.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 Description
  26     The insertion of edges in the edge list is dictated by the upper
  27     triangular ordering. The current decomposition of a polyhedral cell into
  28     tetrahedra requires insertion of cell centres. The points
  29     are ordered in the following way:
  30         1) points of the polyMesh (supporting polyhedral cells)
  31         2) cell centres
  32
  33     The algorithm for owner-neighbour insertion first adds all the
  34     points the owner point shares the edge with (only the ones with
  35     the higher label than the owner point), followed by the
  36     pointCells. This is because the the cell centres are
  37     guaranteed to have a higher index than the internal vertices.
  38
  39     Note:
  40     It is assumed that the edges are constructed such that the start
  41     label is lower than the end label and that the pointCells list are
  42     ordered.
  43
  44 \*---------------------------------------------------------------------------*/
  45
  46 #include "tetPolyMeshCellDecomp.H"
  47 #include "tetPolyMeshLduAddressingCellDecomp.H"
  48
  49 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  50
  51 namespace Foam
  52 {
  53
  54 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
  55
  56 const lduAddressing& tetPolyMeshCellDecomp::lduAddr() const
  57 {
  58     if (!lduPtr_)
  59     {
  60         lduPtr_ = new tetPolyMeshLduAddressingCellDecomp(*this);;
  61     }
  62
  63     return *lduPtr_;
  64 }
  65
  66
  67 label tetPolyMeshCellDecomp::maxNPointsForCell() const
  68 {
  69     if (maxNPointsForCell_ < 0)
  70     {
  71         const faceList& faces = mesh_.faces();
  72         const cellList& polyCells = mesh_.cells();
  73
  74         forAll (polyCells, cellI)
  75         {
  76             maxNPointsForCell_ =
  77                 max
  78                 (
  79                     maxNPointsForCell_,
  80                     polyCells[cellI].labels(faces).size() + 1
  81                 );
  82         }
  83     }
  84
  85     return maxNPointsForCell_;
  86 }
  87
  88
  89 // Fill buffer with addressing for the cell
  90 label tetPolyMeshCellDecomp::addressing
  91 (
  92     const label cellID,
  93     labelList& localToGlobalBuffer,
  94     labelList& globalToLocalBuffer
  95 ) const
  96 {
  97     const unallocFaceList& meshFaces = mesh_.faces();
  98
  99     const labelList& cellFaces = mesh_.cells()[cellID];
 100
 101     label nextLocal = 0;
 102
 103     forAll (cellFaces, faceI)
 104     {
 105         const face& curFace = meshFaces[cellFaces[faceI]];
 106
 107         forAll (curFace, pointI)
 108         {
 109             // If the point has not been already inserted into the local
 110             // buffer, add it
 111             if (globalToLocalBuffer[curFace[pointI]] == -1)
 112             {
 113                 localToGlobalBuffer[nextLocal] = curFace[pointI];
 114                 globalToLocalBuffer[curFace[pointI]] = nextLocal;
 115                 nextLocal++;
 116             }
 117         }
 118     }
 119
 120     // Mark up the cell centre
 121     localToGlobalBuffer[nextLocal] = cellOffset() + cellID;
 122     globalToLocalBuffer[cellOffset() + cellID] = nextLocal;
 123     nextLocal++;
 124
 125     // Return size of addressing
 126     return nextLocal;
 127 }
 128
 129
 130 // Clear global to local addressing
 131 void tetPolyMeshCellDecomp::clearAddressing
 132 (
 133     const label cellID,
 134     const label nCellPoints,
 135     labelList& localToGlobalBuffer,
 136     labelList& globalToLocalBuffer
 137 ) const
 138 {
 139     // Only clear the places that have been used.  The rest of the buffer
 140     // is already initiated to -1
 141     for (label localI = 0; localI < nCellPoints; localI++)
 142     {
 143         globalToLocalBuffer[localToGlobalBuffer[localI]] = -1;
 144     }
 145 }
 146
 147
 148 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 149
 150 } // End namespace Foam
 151
 152 // ************************************************************************* //