src/OpenFOAM/meshes/primitiveMesh/primitiveMeshEdgeFaces.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 "primitiveMesh.H"
  28 #include "ListOps.H"
  29
  30 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
  31
  32 const Foam::labelListList& Foam::primitiveMesh::edgeFaces() const
  33 {
  34     if (!efPtr_)
  35     {
  36         if (debug)
  37         {
  38             Pout<< "primitiveMesh::edgeFaces() : calculating edgeFaces"
  39                 << endl;
  40
  41             if (debug == -1)
  42             {
  43                 // For checking calls:abort so we can quickly hunt down
  44                 // origin of call
  45                 FatalErrorIn("primitiveMesh::edgeFaces()")
  46                     << abort(FatalError);
  47             }
  48         }
  49
  50         // Invert faceEdges
  51         efPtr_ = new labelListList(nEdges());
  52         invertManyToMany(nEdges(), faceEdges(), *efPtr_);
  53     }
  54
  55     return *efPtr_;
  56 }
  57
  58
  59 const Foam::labelList& Foam::primitiveMesh::edgeFaces
  60 (
  61     const label edgeI,
  62     DynamicList<label>& storage
  63 ) const
  64 {
  65     if (hasEdgeFaces())
  66     {
  67         return edgeFaces()[edgeI];
  68     }
  69     else
  70     {
  71         // Use the fact that pointEdges are sorted in incrementing edge order
  72         const edge& e = edges()[edgeI];
  73         const labelList& pFaces0 = pointFaces()[e[0]];
  74         const labelList& pFaces1 = pointFaces()[e[1]];
  75
  76         label i0 = 0;
  77         label i1 = 0;
  78
  79         storage.clear();
  80
  81         while (i0 < pFaces0.size() && i1 < pFaces1.size())
  82         {
  83             if (pFaces0[i0] < pFaces1[i1])
  84             {
  85                 ++i0;
  86             }
  87             else if (pFaces0[i0] > pFaces1[i1])
  88             {
  89                 ++i1;
  90             }
  91             else
  92             {
  93                 // Equal. Append.
  94                 storage.append(pFaces0[i0]);
  95                 ++i0;
  96                 ++i1;
  97             }
  98         }
  99
 100         return storage;
 101     }
 102 }
 103
 104
 105 const Foam::labelList& Foam::primitiveMesh::edgeFaces(const label edgeI) const
 106 {
 107     return edgeFaces(edgeI, labels_);
 108 }
 109
 110
 111 // ************************************************************************* //