Go to the documentation of this file.00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028
00029
00030
00031
00032
00033 #ifndef MapFvSurfaceField_H
00034 #define MapFvSurfaceField_H
00035
00036 #include <OpenFOAM/Field.H>
00037 #include <finiteVolume/surfaceMesh.H>
00038
00039
00040
00041 namespace Foam
00042 {
00043
00044 template<class Type, class MeshMapper>
00045 class MapInternalField<Type, MeshMapper, surfaceMesh>
00046 {
00047 public:
00048
00049 MapInternalField()
00050 {}
00051
00052 void operator()
00053 (
00054 Field<Type>& field,
00055 const MeshMapper& mapper
00056 ) const;
00057 };
00058
00059
00060 template<class Type, class MeshMapper>
00061 void MapInternalField<Type, MeshMapper, surfaceMesh>::operator()
00062 (
00063 Field<Type>& field,
00064 const MeshMapper& mapper
00065 ) const
00066 {
00067 if (field.size() != mapper.surfaceMap().sizeBeforeMapping())
00068 {
00069 FatalErrorIn
00070 (
00071 "void MapInternalField<Type, MeshMapper, surfaceMesh>::operator()\n"
00072 "(\n"
00073 " Field<Type>& field,\n"
00074 " const MeshMapper& mapper\n"
00075 ") const"
00076 ) << "Incompatible size before mapping. Field size: " << field.size()
00077 << " map size: " << mapper.surfaceMap().sizeBeforeMapping()
00078 << abort(FatalError);
00079 }
00080
00081 field.autoMap(mapper.surfaceMap());
00082
00083
00084 const labelList flipFaces = mapper.surfaceMap().flipFaceFlux().toc();
00085
00086 forAll (flipFaces, i)
00087 {
00088 if (flipFaces[i] < field.size())
00089 {
00090 field[flipFaces[i]] *= -1.0;
00091 }
00092 }
00093 }
00094
00095
00096
00097
00098 }
00099
00100
00101
00102 #endif
00103
00104