interpolationCellPointWallModifiedI.H
Go to the documentation of this file.
1 /*---------------------------------------------------------------------------*\
2  ========= |
3  \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
4  \\ / O peration |
5  \\ / A nd | Copyright (C) 2011-2017 OpenFOAM Foundation
6  \\/ M anipulation |
7 -------------------------------------------------------------------------------
8 License
9  This file is part of OpenFOAM.
10 
11  OpenFOAM is free software: you can redistribute it and/or modify it
12  under the terms of the GNU General Public License as published by
13  the Free Software Foundation, either version 3 of the License, or
14  (at your option) any later version.
15 
16  OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
17  ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
18  FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
19  for more details.
20 
21  You should have received a copy of the GNU General Public License
22  along with OpenFOAM. If not, see <http://www.gnu.org/licenses/>.
23 
24 \*---------------------------------------------------------------------------*/
25 
26 // * * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * //
27 
28 template<class Type>
30 (
32 ) const
33 {
34  const barycentric& weights = cpw.weights();
35  const triFace& faceVertices = cpw.faceVertices();
36 
37  Type t = this->psi_[cpw.cell()]*weights[0];
38  t += this->psip_[faceVertices[0]]*weights[1];
39  t += this->psip_[faceVertices[1]]*weights[2];
40  t += this->psip_[faceVertices[2]]*weights[3];
41 
42  return t;
43 }
44 
45 
46 template<class Type>
48 (
49  const vector& position,
50  const label celli,
51  const label facei
52 ) const
53 {
54  return interpolate
55  (
57  (
58  this->pMesh_,
59  position,
60  celli,
61  facei
62  )
63  );
64 }
65 
66 
67 template<class Type>
69 (
70  const barycentric& coordinates,
71  const tetIndices& tetIs,
72  const label facei
73 ) const
74 {
75  if (facei >= 0)
76  {
77  if (facei != tetIs.face())
78  {
80  << "specified face " << facei << " inconsistent with the face "
81  << "stored by tetIndices: " << tetIs.face()
82  << exit(FatalError);
83  }
84 
85  const polyBoundaryMesh& bm = this->pMesh_.boundaryMesh();
86  label patchi = bm.whichPatch(facei);
87 
88  if (patchi != -1)
89  {
90  if (isA<wallPolyPatch>(bm[patchi]))
91  {
92  Type t = this->psi_[tetIs.cell()];
93 
94  return t;
95  }
96  }
97  }
98 
99  // If the wall face selection did not return, then use the normal
100  // interpolate method
101 
103  (
104  coordinates,
105  tetIs,
106  facei
107  );
108 }
109 
110 
111 // ************************************************************************* //
intWM_LABEL_SIZE_t label
A label is an int32_t or int64_t as specified by the pre-processor macro WM_LABEL_SIZE.
Definition: label.H:59
errorManipArg< error, int > exit(error &err, const int errNo=1)
Definition: errorManip.H:124
error FatalError
#define FatalErrorInFunction
Report an error message using Foam::FatalError.
Definition: error.H:319
const triFace & faceVertices() const
Interpolation addressing for points on face.
const barycentric & weights() const
Interpolation weights.
A triangular face using a FixedList of labels corresponding to mesh vertices.
Definition: triFace.H:68
label cell() const
Return the cell.
Definition: tetIndicesI.H:28
Storage and named access for the indices of a tet which is part of the decomposition of a cell...
Definition: tetIndices.H:81
label face() const
Return the face.
Definition: tetIndicesI.H:40
Foam::polyBoundaryMesh.
static tmp< GeometricField< Type, fvsPatchField, surfaceMesh > > interpolate(const GeometricField< Type, fvPatchField, volMesh > &tvf, const surfaceScalarField &faceFlux, Istream &schemeData)
Interpolate field onto faces using scheme given by Istream.
label cell() const
Cell index.
Type interpolate(const cellPointWeightWallModified &cpw) const
Interpolate field for the given cellPointWeight.
label patchi
label whichPatch(const label faceIndex) const
Return patch index for a given face label.
Foam::cellPointWeightWallModified.