sampledPatchTemplates.C
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-2019 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 #include "sampledPatch.H"
27 
28 // * * * * * * * * * * * * * Private Member Functions * * * * * * * * * * * //
29 
30 template<class Type>
32 Foam::sampledSurfaces::patch::sampleField
33 (
35 ) const
36 {
37  // One value per face
38  tmp<Field<Type>> tvalues(new Field<Type>(patchFaceLabels_.size()));
39  Field<Type>& values = tvalues.ref();
40  forAll(patchFaceLabels_, i)
41  {
42  label patchi = patchIDs_[patchIndex_[i]];
43  const Field<Type>& bField = vField.boundaryField()[patchi];
44  values[i] = bField[patchFaceLabels_[i]];
45  }
46 
47  return tvalues;
48 }
49 
50 
51 template<class Type>
53 Foam::sampledSurfaces::patch::sampleField
54 (
56 ) const
57 {
58  // One value per face
59  tmp<Field<Type>> tvalues(new Field<Type>(patchFaceLabels_.size()));
60  Field<Type>& values = tvalues.ref();
61 
62  forAll(patchFaceLabels_, i)
63  {
64  label patchi = patchIDs_[patchIndex_[i]];
65  values[i] = sField.boundaryField()[patchi][patchFaceLabels_[i]];
66  }
67 
68  return tvalues;
69 }
70 
71 
72 template<class Type>
74 Foam::sampledSurfaces::patch::interpolateField
75 (
76  const interpolation<Type>& interpolator
77 ) const
78 {
79  // One value per vertex
80  tmp<Field<Type>> tvalues(new Field<Type>(points().size()));
81  Field<Type>& values = tvalues.ref();
82 
83  const labelList& own = mesh().faceOwner();
84 
85  boolList pointDone(points().size(), false);
86 
87  forAll(faces(), cutFacei)
88  {
89  label patchi = patchIDs_[patchIndex_[cutFacei]];
90  const polyPatch& pp = mesh().boundaryMesh()[patchi];
91  label patchFacei = patchFaceLabels()[cutFacei];
92  const face& f = faces()[cutFacei];
93 
94  forAll(f, faceVertI)
95  {
96  label pointi = f[faceVertI];
97 
98  if (!pointDone[pointi])
99  {
100  label facei = patchFacei + pp.start();
101  label celli = own[facei];
102 
103  values[pointi] = interpolator.interpolate
104  (
105  points()[pointi],
106  celli,
107  facei
108  );
109  pointDone[pointi] = true;
110  }
111  }
112  }
113 
114  return tvalues;
115 }
116 
117 
118 // ************************************************************************* //
#define forAll(list, i)
Loop across all elements in list.
Definition: UList.H:434
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
A face is a list of labels corresponding to mesh vertices.
Definition: face.H:75
const Boundary & boundaryField() const
Return const-reference to the boundary field.
T & ref() const
Return non-const reference or generate a fatal error.
Definition: tmpI.H:174
Generic GeometricField class.
virtual Type interpolate(const vector &position, const label celli, const label facei=-1) const =0
Interpolate field to the given point in the given cell.
dynamicFvMesh & mesh
const pointField & points
Pre-declare SubField and related Field type.
Definition: Field.H:56
labelList f(nPoints)
label patchi
label start() const
Return start label of this patch in the polyMesh face list.
Definition: polyPatch.H:303
Abstract base class for interpolation.
A class for managing temporary objects.
Definition: PtrList.H:53
A patch is a list of labels that address the faces in the global face list.
Definition: polyPatch.H:66