src/dynamicMesh/polyMeshModifiers/slidingInterface/enrichedPatch/enrichedPatchPointMap.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     Enriched patch point map
  27
  28 Author
  29     Hrvoje Jasak, Wikki Ltd.  All rights reserved.  Copyright Hrvoje Jasak.
  30
  31 \*---------------------------------------------------------------------------*/
  32
  33 #include "enrichedPatch.H"
  34
  35 // * * * * * * * * * * * * * Private Member Functions  * * * * * * * * * * * //
  36
  37 void Foam::enrichedPatch::completePointMap() const
  38 {
  39     if (pointMapComplete_)
  40     {
  41         FatalErrorIn("void enrichedPatch::completePointMap() const")
  42             << "Point map already completed"
  43             << abort(FatalError);
  44     }
  45
  46     pointMapComplete_ = true;
  47
  48     const Map<label>& pmm = pointMergeMap();
  49
  50     // Get the mesh points for both patches.  If the point has not been
  51     // merged away, add it to the map
  52
  53     // Do master patch
  54     const labelList& masterMeshPoints = masterPatch_.meshPoints();
  55     const pointField& masterLocalPoints = masterPatch_.localPoints();
  56
  57     forAll (masterMeshPoints, pointI)
  58     {
  59         if
  60         (
  61             !pmm.found(masterMeshPoints[pointI])
  62          && !pointMap_.found(masterMeshPoints[pointI])
  63         )
  64         {
  65             pointMap_.insert
  66             (
  67                 masterMeshPoints[pointI],
  68                 masterLocalPoints[pointI]
  69             );
  70         }
  71     }
  72
  73     // Do slave patch
  74     const labelList& slaveMeshPoints = slavePatch_.meshPoints();
  75     const pointField& slaveLocalPoints = slavePatch_.localPoints();
  76
  77     forAll (slaveMeshPoints, pointI)
  78     {
  79         if
  80         (
  81             !pmm.found(slaveMeshPoints[pointI])
  82          && !pointMap_.found(slaveMeshPoints[pointI])
  83         )
  84         {
  85             pointMap_.insert
  86             (
  87                 slaveMeshPoints[pointI],
  88                 slaveLocalPoints[pointI]
  89             );
  90         }
  91     }
  92 }
  93
  94
  95 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
  96
  97
  98 Foam::Map<Foam::point>& Foam::enrichedPatch::pointMap()
  99 {
 100     return pointMap_;
 101 }
 102
 103
 104 const Foam::Map<Foam::point>& Foam::enrichedPatch::pointMap() const
 105 {
 106     return pointMap_;
 107 }
 108
 109
 110 // ************************************************************************* //