Loading...
Searching...
No Matches
extendedCentredFaceToCellStencil.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 | www.openfoam.com
6 \\/ M anipulation |
7-------------------------------------------------------------------------------
8 Copyright (C) 2011-2016 OpenFOAM Foundation
9-------------------------------------------------------------------------------
10License
11 This file is part of OpenFOAM.
12
13 OpenFOAM is free software: you can redistribute it and/or modify it
14 under the terms of the GNU General Public License as published by
15 the Free Software Foundation, either version 3 of the License, or
16 (at your option) any later version.
17
18 OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
19 ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
20 FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
21 for more details.
22
23 You should have received a copy of the GNU General Public License
24 along with OpenFOAM. If not, see <http://www.gnu.org/licenses/>.
25
26Class
27 Foam::extendedCentredFaceToCellStencil
28
29Description
30
31SourceFiles
32 extendedCentredFaceToCellStencil.C
33
34\*---------------------------------------------------------------------------*/
35
36#ifndef extendedCentredFaceToCellStencil_H
37#define extendedCentredFaceToCellStencil_H
38
40
41// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
42
43namespace Foam
44{
45
48/*---------------------------------------------------------------------------*\
49 Class extendedCentredFaceToCellStencil Declaration
50\*---------------------------------------------------------------------------*/
51
52class extendedCentredFaceToCellStencil
53:
55{
56 // Private data
57
58 //- Swap map for getting neighbouring data
60
61 //- Per face the stencil.
62 labelListList stencil_;
63
64
65 // Private Member Functions
66
67 //- No copy construct
68 extendedCentredFaceToCellStencil
69 (
70 const extendedCentredFaceToCellStencil&
71 ) = delete;
72
73 //- No copy assignment
74 void operator=(const extendedCentredFaceToCellStencil&) = delete;
75
76
77public:
78
79 // Constructors
80
81 //- Construct from uncompacted face stencil
82 explicit extendedCentredFaceToCellStencil(const faceToCellStencil&);
83
84
85 // Member Functions
86
87 //- Return reference to the parallel distribution map
88 const mapDistribute& map() const
89 {
90 return *mapPtr_;
91 }
92
93 //- Return reference to the stencil
94 const labelListList& stencil() const
95 {
96 return stencil_;
97 }
98
99 //- After removing elements from the stencil adapt the schedule (map).
100 void compact();
101
102 //- Use map to get the data into stencil order
103 template<class T>
104 void collectData
105 (
107 List<List<T>>& stencilFld
108 ) const
109 {
111 (
112 map(),
113 stencil(),
114 fld,
115 stencilFld
116 );
118
119 //- Sum surface field contributions to create cell values
120 template<class Type>
122 (
124 const List<List<scalar>>& stencilWeights
125 ) const
126 {
128 (
129 map(),
130 stencil(),
131 fld,
132 stencilWeights
133 );
134 }
135
137
138
139// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
140
141} // End namespace Foam
142
143// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
144
145#endif
146
147// ************************************************************************* //
Info<< nl;Info<< "Write faMesh in vtk format:"<< nl;{ vtk::uindirectPatchWriter writer(aMesh.patch(), fileName(aMesh.time().globalPath()/vtkBaseFileName));writer.writeGeometry();globalIndex procAddr(aMesh.nFaces());labelList cellIDs;if(UPstream::master()) { cellIDs.resize(procAddr.totalSize());for(const labelRange &range :procAddr.ranges()) { auto slice=cellIDs.slice(range);slice=identity(range);} } writer.beginCellData(4);writer.writeProcIDs();writer.write("cellID", cellIDs);writer.write("area", aMesh.S().field());writer.write("normal", aMesh.faceAreaNormals());writer.beginPointData(1);writer.write("normal", aMesh.pointAreaNormals());Info<< " "<< writer.output().name()<< nl;}{ vtk::lineWriter writer(aMesh.points(), aMesh.edges(), fileName(aMesh.time().globalPath()/(vtkBaseFileName+"-edges")));writer.writeGeometry();writer.beginCellData(4);writer.writeProcIDs();{ Field< scalar > fld(faMeshTools::flattenEdgeField(aMesh.magLe(), true))
Generic GeometricField class.
A 1D array of objects of type <T>, where the size of the vector is known and used for subscript bound...
Definition List.H:72
Pointer management similar to std::unique_ptr, with some additional methods and type checking.
Definition autoPtr.H:65
void compact()
After removing elements from the stencil adapt the schedule (map).
void collectData(const GeometricField< T, fvsPatchField, surfaceMesh > &fld, List< List< T > > &stencilFld) const
Use map to get the data into stencil order.
const mapDistribute & map() const
Return reference to the parallel distribution map.
const labelListList & stencil() const
Return reference to the stencil.
tmp< GeometricField< Type, fvPatchField, volMesh > > weightedSum(const GeometricField< Type, fvsPatchField, surfaceMesh > &fld, const List< List< scalar > > &stencilWeights) const
Sum surface field contributions to create cell values.
static void collectData(const mapDistribute &map, const labelListList &stencil, const GeometricField< T, fvsPatchField, surfaceMesh > &fld, List< List< T > > &stencilFld)
Use map to get the data into stencil order.
static tmp< GeometricField< Type, fvPatchField, volMesh > > weightedSum(const mapDistribute &map, const labelListList &stencil, const GeometricField< Type, fvsPatchField, surfaceMesh > &fld, const List< List< scalar > > &stencilWeights)
Sum surface field contributions to create cell values.
extendedFaceToCellStencil(const polyMesh &)
Construct from mesh.
baseclass for extended cell centred addressing. Contains per cell a list of neighbouring faces in glo...
Class containing processor-to-processor mapping information.
A class for managing temporary objects.
Definition tmp.H:75
Namespace for OpenFOAM.
List< labelList > labelListList
List of labelList.
Definition labelList.H:38