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
00034
00035
00036
00037
00038
00039
00040 #ifndef faceEdge_H
00041 #define faceEdge_H
00042
00043 #include "label.H"
00044 #include "primitiveMesh.H"
00045
00046
00047
00048 namespace Foam
00049 {
00050
00051
00052
00053
00054
00055
00056
00057 class faceEdge
00058 {
00059
00060
00061
00062 label faceLabel_;
00063
00064
00065 label faceIndex_;
00066
00067 public:
00068
00069
00070
00071
00072 class faceEdgeHash
00073 {
00074
00075 public:
00076
00077 faceEdgeHash()
00078 {}
00079
00080
00081 label operator()(const faceEdge& fe, const label tableSize)
00082 const
00083 {
00084
00085
00086 return
00087 mag
00088 (
00089 fe.faceLabel()
00090 + (fe.faceLabel() * fe.faceLabel())
00091 + fe.faceIndex()
00092 + (fe.faceIndex() * fe.faceIndex())
00093 )
00094 % tableSize;
00095 }
00096 };
00097
00098
00099
00100
00101
00102 inline faceEdge()
00103 :
00104 faceLabel_(-1),
00105 faceIndex_(-1)
00106 {}
00107
00108
00109 inline faceEdge(const label faceLabel, const label faceIndex)
00110 :
00111 faceLabel_(faceLabel),
00112 faceIndex_(faceIndex)
00113 {}
00114
00115
00116
00117 label faceLabel() const
00118 {
00119 return faceLabel_;
00120 }
00121
00122 label faceIndex() const
00123 {
00124 return faceIndex_;
00125 }
00126
00127 template <class T>
00128 T interpolate
00129 (
00130 const primitiveMesh& mesh,
00131 const Field<T>& faceField,
00132 const Field<T>& vertField,
00133 const scalar weight
00134 ) const
00135 {
00136 const face& f = mesh.faces()[faceLabel_];
00137
00138 return
00139 (1-weight)*vertField[f[faceIndex_]]
00140 + weight*faceField[faceLabel_];
00141 }
00142
00143
00144 point coord(const primitiveMesh& mesh, const scalar weight) const
00145 {
00146 return interpolate(mesh, mesh.faceCentres(), mesh.points(), weight);
00147 }
00148
00149
00150
00151
00152 bool operator==(const faceEdge& fe) const
00153 {
00154 return
00155 (faceLabel() == fe.faceLabel())
00156 && (faceIndex() == fe.faceIndex());
00157 }
00158
00159
00160
00161
00162 inline friend Ostream& operator<<(Ostream& os, const faceEdge& fe)
00163 {
00164 os << token::BEGIN_LIST
00165 << fe.faceLabel_ << token::SPACE
00166 << fe.faceIndex_
00167 << token::END_LIST;
00168
00169
00170 os.check("Ostream& operator<<(Ostream&, const faceEdge&)");
00171
00172 return os;
00173 }
00174 };
00175
00176
00177
00178
00179 }
00180
00181
00182
00183 #endif
00184
00185