1 /*---------------------------------------------------------------------------*\
3 \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
5 \\ / A nd | Copyright (C) 1991-2009 OpenCFD Ltd.
7 -------------------------------------------------------------------------------
9 This file is part of OpenFOAM.
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.
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
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
26 Update the polyMesh corresponding to the given map.
28 \*---------------------------------------------------------------------------*/
31 #include "mapPolyMesh.H"
33 #include "globalMeshData.H"
34 #include "pointMesh.H"
36 // * * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * //
38 void Foam::polyMesh::updateMesh(const mapPolyMesh& mpm)
40 // Update boundaryMesh (note that patches themselves already ok)
41 boundary_.updateMesh();
44 pointZones_.clearAddressing();
45 faceZones_.clearAddressing();
46 cellZones_.clearAddressing();
48 // Update parallel data
49 if (globalMeshDataPtr_)
51 globalMeshDataPtr_->updateMesh();
54 setInstance(time().timeName());
56 // Map the old motion points if present
59 // Make a copy of the original points
60 pointField oldMotionPoints = *oldPointsPtr_;
62 pointField& newMotionPoints = *oldPointsPtr_;
64 // Resize the list to new size
65 newMotionPoints.setSize(points_.size());
68 newMotionPoints.map(oldMotionPoints, mpm.pointMap());
71 // Reset valid directions (could change by faces put into empty patches)
72 geometricD_ = Vector<label>::zero;
73 solutionD_ = Vector<label>::zero;
76 // Hack until proper callbacks. Below are all the polyMesh-MeshObjects.
79 if (thisDb().foundObject<pointMesh>(pointMesh::typeName))
81 const_cast<pointMesh&>
83 thisDb().lookupObject<pointMesh>
92 // ************************************************************************* //