MapConsistentVolFields.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-2017 OpenFOAM Foundation
9 Copyright (C) 2019-2021 OpenCFD Ltd.
10-------------------------------------------------------------------------------
11License
12 This file is part of OpenFOAM.
13
14 OpenFOAM is free software: you can redistribute it and/or modify it
15 under the terms of the GNU General Public License as published by
16 the Free Software Foundation, either version 3 of the License, or
17 (at your option) any later version.
18
19 OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
20 ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
21 FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
22 for more details.
23
24 You should have received a copy of the GNU General Public License
25 along with OpenFOAM. If not, see <http://www.gnu.org/licenses/>.
26
27\*---------------------------------------------------------------------------*/
28
29#ifndef MapConsistentVolFields_H
30#define MapConsistentVolFields_H
31
32#include "GeometricField.H"
33#include "meshToMesh0.H"
34#include "IOobjectList.H"
35
36// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
37
38namespace Foam
39{
40
41template<class Type, class CombineOp>
43(
44 const IOobjectList& objects,
45 const meshToMesh0& meshToMesh0Interp,
46 const meshToMesh0::order& mapOrder,
47 const CombineOp& cop
48)
49{
51
52 const fvMesh& meshSource = meshToMesh0Interp.fromMesh();
53 const fvMesh& meshTarget = meshToMesh0Interp.toMesh();
54
55 IOobjectList fields = objects.lookupClass(fieldType::typeName);
56
57 forAllConstIters(fields, fieldIter)
58 {
59 Info<< " interpolating " << (*fieldIter)->name()
60 << endl;
61
62 // Read field. Do not auto-load old-time field
63 fieldType fieldSource(*fieldIter(), meshSource, false);
64
65 IOobject fieldTargetIOobject
66 (
67 (*fieldIter)->name(),
68 meshTarget.time().timeName(),
69 meshTarget,
72 );
73
74 if (fieldTargetIOobject.typeHeaderOk<fieldType>(true))
75 {
76 // Read fieldTarget
77 fieldType fieldTarget(fieldTargetIOobject, meshTarget, false);
78
79 // Interpolate field
80 meshToMesh0Interp.interpolate
81 (
82 fieldTarget,
83 fieldSource,
84 mapOrder,
85 cop
86 );
87
88 // Write field
89 fieldTarget.write();
90 }
91 else
92 {
93 fieldTargetIOobject.readOpt(IOobject::NO_READ);
94
95 // Interpolate field
96 fieldType fieldTarget
97 (
98 fieldTargetIOobject,
99 meshToMesh0Interp.interpolate
100 (
101 fieldSource,
102 mapOrder,
103 cop
104 )
105 );
106
107 // Write field
108 fieldTarget.write();
109 }
110 }
111}
112
113
114// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
115
116} // End namespace Foam
117
118// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
119
120#endif
121
122// ************************************************************************* //
Generic GeometricField class.
List of IOobjects with searching and retrieving facilities.
Definition: IOobjectList.H:59
IOobjectList lookupClass(const char *clsName) const
The list of IOobjects with the given headerClassName.
Definition: IOobjectList.C:313
Defines the attributes of an object for which implicit objectRegistry management is supported,...
Definition: IOobject.H:170
readOption readOpt() const noexcept
The read option.
Definition: IOobjectI.H:164
bool typeHeaderOk(const bool checkType=true, const bool search=true, const bool verbose=true)
Read header (uses typeFilePath to find file) and check its info.
static word timeName(const scalar t, const int precision=precision_)
Definition: Time.C:780
Mesh data needed to do the Finite Volume discretisation.
Definition: fvMesh.H:91
const Time & time() const
Return the top-level database.
Definition: fvMesh.H:290
Serial mesh to mesh interpolation class.
Definition: meshToMesh0.H:66
const fvMesh & toMesh() const
Definition: meshToMesh0.H:238
void interpolate(VolumeField< Type > &, const VolumeField< Type > &, order=INTERPOLATE, const CombineOp &cop=eqOp< Type >()) const
Interpolate volume field.
order
Enumeration specifying required accuracy.
Definition: meshToMesh0.H:148
const fvMesh & fromMesh() const
Definition: meshToMesh0.H:233
Namespace for OpenFOAM.
messageStream Info
Information stream (stdout output on master, null elsewhere)
Ostream & endl(Ostream &os)
Add newline and flush stream.
Definition: Ostream.H:372
void MapConsistentVolFields(const IOobjectList &objects, const meshToMesh0 &meshToMesh0Interp, const meshToMesh0::order &mapOrder, const CombineOp &cop)
multivariateSurfaceInterpolationScheme< scalar >::fieldTable fields
Definition: createFields.H:97
#define forAllConstIters(container, iter)
Iterate across all elements of the container object with const access.
Definition: stdFoam.H:278