src/OpenFOAM/meshes/primitiveMesh/PrimitivePatch/PrimitivePatchLocalPointOrder.C

   1 /*---------------------------------------------------------------------------*\
   2   =========                 |
   3   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   4    \\    /   O peration     |
   5     \\  /    A nd           | Copyright (C) 2011 OpenFOAM Foundation
   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
  13     the Free Software Foundation, either version 3 of the License, or
  14     (at your 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, see <http://www.gnu.org/licenses/>.
  23
  24 Description
  25      Orders the local points on the patch for most efficient search
  26
  27 \*---------------------------------------------------------------------------*/
  28
  29 #include "SLList.H"
  30 #include "boolList.H"
  31
  32 // * * * * * * * * * * * * * Private Member Functions  * * * * * * * * * * * //
  33
  34 template
  35 <
  36     class Face,
  37     template<class> class FaceList,
  38     class PointField,
  39     class PointType
  40 >
  41 void
  42 Foam::PrimitivePatch<Face, FaceList, PointField, PointType>::
  43 calcLocalPointOrder() const
  44 {
  45     // Note: Cannot use bandCompressing as point-point addressing does
  46     // not exist and is not considered generally useful.
  47     //
  48
  49     if (debug)
  50     {
  51         Pout<< "PrimitivePatch<Face, FaceList, PointField, PointType>::"
  52             << "calcLocalPointOrder() : "
  53             << "calculating local point order"
  54             << endl;
  55     }
  56
  57     if (localPointOrderPtr_)
  58     {
  59         // it is considered an error to attempt to recalculate
  60         // if already allocated
  61         FatalErrorIn
  62         (
  63             "PrimitivePatch<Face, FaceList, PointField, PointType>::"
  64             "calcLocalPointOrder()"
  65         )   << "local point order already calculated"
  66             << abort(FatalError);
  67     }
  68
  69     const List<Face>& lf = localFaces();
  70
  71     const labelListList& ff = faceFaces();
  72
  73     boolList visitedFace(lf.size(), false);
  74
  75     localPointOrderPtr_ = new labelList(meshPoints().size(), -1);
  76
  77     labelList& pointOrder = *localPointOrderPtr_;
  78
  79     boolList visitedPoint(pointOrder.size(), false);
  80
  81     label nPoints = 0;
  82
  83     forAll(lf, faceI)
  84     {
  85         if (!visitedFace[faceI])
  86         {
  87             SLList<label> faceOrder(faceI);
  88
  89             do
  90             {
  91                 const label curFace = faceOrder.first();
  92
  93                 faceOrder.removeHead();
  94
  95                 if (!visitedFace[curFace])
  96                 {
  97                     visitedFace[curFace] = true;
  98
  99                     const labelList& curPoints = lf[curFace];
 100
 101                     // mark points
 102                     forAll(curPoints, pointI)
 103                     {
 104                         if (!visitedPoint[curPoints[pointI]])
 105                         {
 106                             visitedPoint[curPoints[pointI]] = true;
 107
 108                             pointOrder[nPoints] = curPoints[pointI];
 109
 110                             nPoints++;
 111                         }
 112                     }
 113
 114                     // add face neighbours to the list
 115                     const labelList& nbrs = ff[curFace];
 116
 117                     forAll(nbrs, nbrI)
 118                     {
 119                         if (!visitedFace[nbrs[nbrI]])
 120                         {
 121                             faceOrder.append(nbrs[nbrI]);
 122                         }
 123                     }
 124                 }
 125             } while (faceOrder.size());
 126         }
 127     }
 128
 129     if (debug)
 130     {
 131         Pout<< "PrimitivePatch<Face, FaceList, PointField, PointType>::"
 132             << "calcLocalPointOrder() "
 133             << "finished calculating local point order"
 134             << endl;
 135     }
 136 }
 137
 138
 139 // ************************************************************************* //