OpenFOAM logo
Open Source CFD Toolkit

MapTetFemElementField.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     MapTetFemElementField
00027 
00028 Description
00029     Map element internal field on topology change.  This is a partial
00030     template specialisation, see MapGeometricFields.
00031 
00032 \*---------------------------------------------------------------------------*/
00033 
00034 #ifndef MapTetFemElementField_H
00035 #define MapTetFemElementField_H
00036 
00037 #include "Field.H"
00038 #include "elementMesh.H"
00039 
00040 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
00041 
00042 namespace Foam
00043 {
00044 
00045 template<class Type, class MeshMapper>
00046 class MapInternalField<Type, MeshMapper, elementMesh>
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, elementMesh>::operator()
00063 (
00064     Field<Type>& field,
00065     const MeshMapper& mapper
00066 ) const
00067 {
00068     if (field.size() != mapper.elementMap().sizeBeforeMapping())
00069     {
00070         FatalErrorIn
00071         (
00072             "void MapInternalField<Type, MeshMapper, elementMesh>::operator()\n"
00073             "(\n"
00074             "    Field<Type>& field,\n"
00075             "    const MeshMapper& mapper\n"
00076             ") const"
00077         )  << "Incompatible size before mapping.  Field size: " << field.size()
00078            << " map size: " << mapper.elementMap().sizeBeforeMapping()
00079            << abort(FatalError);
00080     }
00081 
00082     field.autoMap(mapper.elementMap());
00083 }
00084 
00085 
00086 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
00087 
00088 } // End namespace Foam
00089 
00090 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
00091 
00092 #endif
00093 
00094 // ************************************************************************* //
For further information go to www.openfoam.org