src/OpenFOAM/meshes/polyMesh/polyMeshUpdate.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     Update the polyMesh corresponding to the given map.
  26
  27 \*---------------------------------------------------------------------------*/
  28
  29 #include "polyMesh.H"
  30 #include "mapPolyMesh.H"
  31 #include "Time.H"
  32 #include "globalMeshData.H"
  33 #include "pointMesh.H"
  34
  35 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
  36
  37 void Foam::polyMesh::updateMesh(const mapPolyMesh& mpm)
  38 {
  39     // Update boundaryMesh (note that patches themselves already ok)
  40     boundary_.updateMesh();
  41
  42     // Update zones
  43     pointZones_.clearAddressing();
  44     faceZones_.clearAddressing();
  45     cellZones_.clearAddressing();
  46
  47     // Update parallel data
  48     if (globalMeshDataPtr_)
  49     {
  50         globalMeshDataPtr_->updateMesh();
  51     }
  52
  53     setInstance(time().timeName());
  54
  55     // Map the old motion points if present
  56     if (oldPointsPtr_)
  57     {
  58         // Make a copy of the original points
  59         pointField oldMotionPoints = *oldPointsPtr_;
  60
  61         pointField& newMotionPoints = *oldPointsPtr_;
  62
  63         // Resize the list to new size
  64         newMotionPoints.setSize(points_.size());
  65
  66         // Map the list
  67         newMotionPoints.map(oldMotionPoints, mpm.pointMap());
  68
  69         // Any points created out-of-nothing get set to the current coordinate
  70         // for lack of anything better.
  71         forAll(mpm.pointMap(), newPointI)
  72         {
  73             if (mpm.pointMap()[newPointI] == -1)
  74             {
  75                 newMotionPoints[newPointI] = points_[newPointI];
  76             }
  77         }
  78     }
  79
  80     // Reset valid directions (could change by faces put into empty patches)
  81     geometricD_ = Vector<label>::zero;
  82     solutionD_ = Vector<label>::zero;
  83
  84
  85     // Hack until proper callbacks. Below are all the polyMesh-MeshObjects.
  86
  87     // pointMesh
  88     if (thisDb().foundObject<pointMesh>(pointMesh::typeName))
  89     {
  90         const_cast<pointMesh&>
  91         (
  92             thisDb().lookupObject<pointMesh>
  93             (
  94                 pointMesh::typeName
  95             )
  96         ).updateMesh(mpm);
  97     }
  98 }
  99
 100
 101 // ************************************************************************* //