interpolationCellPointI.H
Go to the documentation of this file.
1 /*---------------------------------------------------------------------------*\
2  ========= |
3  \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
4  \\ / O peration | Website: https://openfoam.org
5  \\ / A nd | Copyright (C) 2011-2020 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 (
31  const cellPointWeight& cpw
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 += psip_[faceVertices[0]]*weights[1];
39  t += psip_[faceVertices[1]]*weights[2];
40  t += 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(cellPointWeight(this->mesh_, position, celli, facei));
55 }
56 
57 
58 template<class Type>
60 (
61  const barycentric& coordinates,
62  const tetIndices& tetIs,
63  const label facei
64 ) const
65 {
66  // Assumes that the position is consistent with the supplied
67  // tetIndices. Does not pay attention to whether or not facei is
68  // supplied or not - the result will be essentially the same.
69  // Performs a consistency check, however.
70 
71  if (facei >= 0)
72  {
73  if (facei != tetIs.face())
74  {
76  << "specified face " << facei << " inconsistent with the face "
77  << "stored by tetIndices: " << tetIs.face()
78  << exit(FatalError);
79  }
80  }
81 
82  const triFace triIs = tetIs.faceTriIs(this->mesh_);
83 
84  return
85  this->psi_[tetIs.cell()]*coordinates[0]
86  + psip_[triIs[0]]*coordinates[1]
87  + psip_[triIs[1]]*coordinates[2]
88  + psip_[triIs[2]]*coordinates[3];
89 }
90 
91 
92 // ************************************************************************* //
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:306
const triFace & faceVertices() const
Interpolation addressing for points on face.
triFace faceTriIs(const polyMesh &mesh) const
Return the indices corresponding to the tri on the face for.
Definition: tetIndicesI.H:67
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:31
Type interpolate(const cellPointWeight &cpw) const
Interpolate field for the given cellPointWeight.
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:43
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.
Foam::cellPointWeight.