src/dynamicMesh/attachDetach/attachDetachPointMatchMap.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 "attachDetach.H"
  28 #include "polyMesh.H"
  29 #include "primitiveMesh.H"
  30 #include "primitiveFacePatch.H"
  31 #include "polyTopoChanger.H"
  32
  33 // * * * * * * * * * * * * * Private Member Functions  * * * * * * * * * * * //
  34
  35 const Foam::Map<Foam::label>&
  36 Foam::attachDetach::pointMatchMap() const
  37 {
  38     if (!pointMatchMapPtr_)
  39     {
  40         calcPointMatchMap();
  41     }
  42
  43     return *pointMatchMapPtr_;
  44 }
  45
  46
  47 void Foam::attachDetach::calcPointMatchMap() const
  48 {
  49     if (debug)
  50     {
  51         Pout<< "void attachDetach::calcPointMatchMap() const "
  52             << " for object " << name() << " : "
  53             << "Calculating point matching" << endl;
  54     }
  55
  56     if (pointMatchMapPtr_)
  57     {
  58         FatalErrorIn
  59         (
  60             "void attachDetach::calcPointMatchMap() const"
  61         )   << "Point match map already calculated for object " << name()
  62             << abort(FatalError);
  63     }
  64
  65     const polyMesh& mesh = topoChanger().mesh();
  66     const faceList& faces = mesh.faces();
  67
  68     const polyPatch& masterPatch = mesh.boundaryMesh()[masterPatchID_.index()];
  69     const polyPatch& slavePatch = mesh.boundaryMesh()[slavePatchID_.index()];
  70
  71     // Create the reverse patch out of the slave patch
  72     primitiveFacePatch reverseSlavePatch
  73     (
  74         faceList(slavePatch.size()),
  75         mesh.points()
  76     );
  77
  78     const label slavePatchStart = slavePatch.start();
  79
  80     forAll (reverseSlavePatch, faceI)
  81     {
  82         reverseSlavePatch[faceI] =
  83             faces[slavePatchStart + faceI].reverseFace();
  84     }
  85
  86     // Create point merge list and remove merged points
  87     const labelList& masterMeshPoints = masterPatch.meshPoints();
  88     const labelList& slaveMeshPoints = reverseSlavePatch.meshPoints();
  89
  90     const faceList& masterLocalFaces = masterPatch.localFaces();
  91     const faceList& slaveLocalFaces = reverseSlavePatch.localFaces();
  92
  93     pointMatchMapPtr_ = new Map<label>(2*slaveMeshPoints.size());
  94     Map<label>& removedPointMap = *pointMatchMapPtr_;
  95
  96     forAll (masterLocalFaces, faceI)
  97     {
  98         const face& curMasterPoints = masterLocalFaces[faceI];
  99         const face& curSlavePoints = slaveLocalFaces[faceI];
 100
 101         forAll (curMasterPoints, pointI)
 102         {
 103             // If the master and slave point labels are the same, the
 104             // point remains.  Otherwise, the slave point is removed and
 105             // replaced by the master
 106             if
 107             (
 108                 masterMeshPoints[curMasterPoints[pointI]]
 109              != slaveMeshPoints[curSlavePoints[pointI]]
 110             )
 111             {
 112 // Pout << "Matching slave point " << slaveMeshPoints[curSlavePoints[pointI]] << " with " << masterMeshPoints[curMasterPoints[pointI]] << endl;
 113
 114                 // Grab the addressing
 115                 removedPointMap.insert
 116                 (
 117                     slaveMeshPoints[curSlavePoints[pointI]],
 118                     masterMeshPoints[curMasterPoints[pointI]]
 119                 );
 120             }
 121         }
 122     }
 123
 124     if (debug)
 125     {
 126         Pout<< "void attachDetach::calcPointMatchMap() const "
 127             << " for object " << name() << " : "
 128             << "Finished calculating point matching" << endl;
 129     }
 130 }
 131
 132
 133 // ************************************************************************* //