Loading...
Searching...
No Matches
solutionControlTemplates.C
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-2013 OpenFOAM Foundation
9 Copyright (C) 2019-2022 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#include "GeometricField.H"
30#include "volMesh.H"
31#include "fvPatchField.H"
32
33// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
34
35template<class Type>
36void Foam::solutionControl::storePrevIter() const
37{
38 typedef GeometricField<Type, fvPatchField, volMesh> GeoField;
39
40 for (GeoField& fld : mesh_.objectRegistry::sorted<GeoField>())
41 {
42 const word& fldName = fld.name();
43
44 if (!fldName.contains("PrevIter") && mesh_.relaxField(fldName))
45 {
46 DebugInfo
47 << algorithmName_ << ": storing previous iter for "
48 << fldName << endl;
49
50 fld.storePrevIter();
51 }
52 }
53}
54
55
56// ************************************************************************* //
fld
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))
Foam::GeometricField
Generic GeometricField class.
Definition GeometricField.H:75
Foam::solutionControl::storePrevIter
void storePrevIter() const
Store previous iteration field for vol<Type>Fields.
Definition solutionControlTemplates.C:29
Foam::solutionControl::algorithmName_
const word algorithmName_
The dictionary name, e.g. SIMPLE, PIMPLE.
Definition solutionControl.H:95
Foam::solutionControl::mesh_
fvMesh & mesh_
Reference to the mesh database.
Definition solutionControl.H:85
Foam::string::contains
bool contains(char c) const noexcept
True if string contains given character (cf. C++23).
Definition string.H:412
Foam::word
A class for handling words, derived from Foam::string.
Definition word.H:66
DebugInfo
#define DebugInfo
Report an information message using Foam::Info.
Definition messageStream.H:552
Foam::endl
Ostream & endl(Ostream &os)
Add newline and flush stream.
Definition Ostream.H:519