src/foam/meshes/pointMesh/pointMesh.C

   1 /*---------------------------------------------------------------------------*\
   2   =========                 |
   3   \\      /  F ield         | foam-extend: Open Source CFD
   4    \\    /   O peration     | Version:     3.2
   5     \\  /    A nd           | Web:         http://www.foam-extend.org
   6      \\/     M anipulation  | For copyright notice see file Copyright
   7 -------------------------------------------------------------------------------
   8 License
   9     This file is part of foam-extend.
  10
  11     foam-extend 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 3 of the License, or (at your
  14     option) any later version.
  15
  16     foam-extend is distributed in the hope that it will be useful, but
  17     WITHOUT ANY WARRANTY; without even the implied warranty of
  18     MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  19     General Public License for more details.
  20
  21     You should have received a copy of the GNU General Public License
  22     along with foam-extend.  If not, see <http://www.gnu.org/licenses/>.
  23
  24 \*---------------------------------------------------------------------------*/
  25
  26 #include "pointMesh.H"
  27 #include "globalMeshData.H"
  28 #include "globalPointPatch.H"
  29 #include "pointMeshMapper.H"
  30 #include "pointFields.H"
  31 #include "MapGeometricFields.H"
  32 #include "MapPointField.H"
  33
  34 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  35
  36 defineTypeNameAndDebug(Foam::pointMesh, 0);
  37
  38
  39 // * * * * * * * * * * * * * Private Member Functions  * * * * * * * * * * * //
  40
  41 void Foam::pointMesh::mapFields(const mapPolyMesh& mpm) const
  42 {
  43     // Create a mapper
  44     const pointMeshMapper m(*this, mpm);
  45
  46     MapGeometricFields<scalar, pointPatchField, pointMeshMapper, pointMesh>(m);
  47     MapGeometricFields<vector, pointPatchField, pointMeshMapper, pointMesh>(m);
  48     MapGeometricFields
  49     <
  50         sphericalTensor,
  51         pointPatchField,
  52         pointMeshMapper,
  53         pointMesh
  54     >(m);
  55
  56     MapGeometricFields<symmTensor, pointPatchField, pointMeshMapper, pointMesh>
  57     (m);
  58
  59     MapGeometricFields
  60     <
  61         symmTensor4thOrder,
  62         pointPatchField,
  63         pointMeshMapper,
  64         pointMesh
  65     >(m);
  66
  67     MapGeometricFields<diagTensor, pointPatchField, pointMeshMapper, pointMesh>
  68     (m);
  69
  70     MapGeometricFields<tensor, pointPatchField, pointMeshMapper, pointMesh>(m);
  71 }
  72
  73
  74 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
  75
  76 Foam::pointMesh::pointMesh
  77 (
  78     const polyMesh& pMesh,
  79     bool alwaysConstructGlobalPatch
  80 )
  81 :
  82     MeshObject<polyMesh, pointMesh>(pMesh),
  83     GeoMesh<polyMesh>(pMesh),
  84     boundary_(*this, pMesh.boundaryMesh())
  85 {
  86     // Add the globalPointPatch if there are global points
  87     if
  88     (
  89         alwaysConstructGlobalPatch
  90      || GeoMesh<polyMesh>::mesh_.globalData().nGlobalPoints()
  91     )
  92     {
  93         boundary_.setSize(boundary_.size() + 1);
  94
  95         boundary_.set
  96         (
  97             boundary_.size() - 1,
  98             new globalPointPatch
  99             (
 100                 boundary_,
 101                 boundary_.size() - 1
 102             )
 103         );
 104     }
 105
 106     // Calculate the geometry for the patches (transformation tensors etc.)
 107     boundary_.calcGeometry();
 108 }
 109
 110
 111 bool Foam::pointMesh::movePoints() const
 112 {
 113     // Casting const-ness to answer the interface of meshObject
 114     // HJ, 30/Aug/2010
 115     const_cast<pointBoundaryMesh&>(boundary_).movePoints();
 116
 117     return true;
 118 }
 119
 120
 121 bool Foam::pointMesh::updateMesh(const mapPolyMesh& mpm) const
 122 {
 123     // Casting const-ness to answer the interface of meshObject
 124     // HJ, 30/Aug/2010
 125     const_cast<pointBoundaryMesh&>(boundary_).updateMesh();
 126
 127     // Map all registered point fields
 128     mapFields(mpm);
 129
 130     return true;
 131 }
 132
 133
 134 // ************************************************************************* //