src/OpenFOAM/meshes/primitiveMesh/primitiveMeshPointCells.C

   1 /*---------------------------------------------------------------------------*\
   2   =========                 |
   3   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   4    \\    /   O peration     |
   5     \\  /    A nd           | Copyright (C) 1991-2008 OpenCFD Ltd.
   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 "cell.H"
  29
  30 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  31
  32 namespace Foam
  33 {
  34
  35 // * * * * * * * * * * * * * Private Member Functions  * * * * * * * * * * * //
  36
  37 void primitiveMesh::calcPointCells() const
  38 {
  39     // Loop through cells and mark up points
  40
  41     if (debug)
  42     {
  43         Pout<< "primitiveMesh::calcPointCells() : "
  44             << "calculating pointCells"
  45             << endl;
  46     }
  47
  48     // It is an error to attempt to recalculate pointCells
  49     // if the pointer is already set
  50     if (pcPtr_)
  51     {
  52         FatalErrorIn("primitiveMesh::calcPointCells() const")
  53             << "pointCells already calculated"
  54             << abort(FatalError);
  55     }
  56     else
  57     {
  58         const cellList& cf = cells();
  59
  60         // Count number of cells per point
  61
  62         labelList npc(nPoints(), 0);
  63
  64         forAll (cf, cellI)
  65         {
  66             const labelList curPoints = cf[cellI].labels(faces());
  67
  68             forAll (curPoints, pointI)
  69             {
  70                 label ptI = curPoints[pointI];
  71
  72                 npc[ptI]++;
  73             }
  74         }
  75
  76
  77         // Size and fill cells per point
  78
  79         pcPtr_ = new labelListList(npc.size());
  80         labelListList& pointCellAddr = *pcPtr_;
  81
  82         forAll (pointCellAddr, pointI)
  83         {
  84             pointCellAddr[pointI].setSize(npc[pointI]);
  85         }
  86         npc = 0;
  87
  88
  89         forAll (cf, cellI)
  90         {
  91             const labelList curPoints = cf[cellI].labels(faces());
  92
  93             forAll (curPoints, pointI)
  94             {
  95                 label ptI = curPoints[pointI];
  96
  97                 pointCellAddr[ptI][npc[ptI]++] = cellI;
  98             }
  99         }
 100     }
 101 }
 102
 103
 104 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
 105
 106 const labelListList& primitiveMesh::pointCells() const
 107 {
 108     if (!pcPtr_)
 109     {
 110         calcPointCells();
 111     }
 112
 113     return *pcPtr_;
 114 }
 115
 116
 117 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 118
 119 } // End namespace Foam
 120
 121 // ************************************************************************* //