OpenFOAM logo
Open Source CFD Toolkit

MapTetFemPointField.H

Go to the documentation of this file.
00001 /*---------------------------------------------------------------------------*\
00002   =========                 |
00003   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
00004    \\    /   O peration     |
00005     \\  /    A nd           | Copyright (C) 1991-2005 OpenCFD Ltd.
00006      \\/     M anipulation  |
00007 -------------------------------------------------------------------------------
00008 License
00009     This file is part of OpenFOAM.
00010 
00011     OpenFOAM is free software; you can redistribute it and/or modify it
00012     under the terms of the GNU General Public License as published by the
00013     Free Software Foundation; either version 2 of the License, or (at your
00014     option) any later version.
00015 
00016     OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
00017     ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
00018     FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License
00019     for more details.
00020 
00021     You should have received a copy of the GNU General Public License
00022     along with OpenFOAM; if not, write to the Free Software Foundation,
00023     Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA
00024 
00025 Class
00026     MapTetFemPointField
00027 
00028 Description
00029     Map tetPoint internal field on topology change.  This is a partial
00030     template specialisation, see MapGeometricFields.
00031 
00032 \*---------------------------------------------------------------------------*/
00033 
00034 #ifndef MapTetFemPointField_H
00035 #define MapTetFemPointField_H
00036 
00037 #include "Field.H"
00038 #include "tetPointMesh.H"
00039 
00040 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
00041 
00042 namespace Foam
00043 {
00044 
00045 template<class Type, class MeshMapper>
00046 class MapInternalField<Type, MeshMapper, tetPointMesh>
00047 {
00048 public:
00049 
00050     MapInternalField()
00051     {}
00052 
00053     void operator()
00054     (
00055         Field<Type>& field,
00056         const MeshMapper& mapper
00057     ) const;
00058 };
00059 
00060 
00061 template<class Type, class MeshMapper>
00062 void MapInternalField<Type, MeshMapper, tetPointMesh>::operator()
00063 (
00064     Field<Type>& field,
00065     const MeshMapper& mapper
00066 ) const
00067 {
00068     if (field.size() != mapper.pointMap().sizeBeforeMapping())
00069     {
00070         FatalErrorIn
00071         (
00072             "void MapInternalField<Type, MeshMapper, tetPointMesh>"
00073             "::operator()\n"
00074             "(\n"
00075             "    Field<Type>& field,\n"
00076             "    const MeshMapper& mapper\n"
00077             ") const"
00078         )  << "Incompatible size before mapping.  Field size: " << field.size()
00079            << " map size: " << mapper.pointMap().sizeBeforeMapping()
00080            << abort(FatalError);
00081     }
00082 
00083     field.autoMap(mapper.pointMap());
00084 }
00085 
00086 
00087 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
00088 
00089 } // End namespace Foam
00090 
00091 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
00092 
00093 #endif
00094 
00095 // ************************************************************************* //
For further information go to www.openfoam.org