zhyang-dev/OpenFOAM-12
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

fvMesh: Construct changers after geometry pointers have been set

Browse Source 
This allows changers to access geometry from the fvMesh on construction.
It also prevents a crash that occurs when a changer loads a field in
which there is a cyclic boundary condition.
This commit is contained in:
Will Bainbridge
2022-02-16 15:51:27 +00:00
parent fd50abcc85
commit 1950a11974
1 changed files with 20 additions and 18 deletions
Download Patch File Download Diff File Expand all files Collapse all files

38
src/finiteVolume/fvMesh/fvMesh.C
View File

@ -260,24 +260,9 @@ Foam::fvMesh::fvMesh(const IOobject& io, const bool changers)
fvSolution(static_cast<const objectRegistry&>(*this)),
data(static_cast<const objectRegistry&>(*this)),
boundary_(*this, boundaryMesh()),
topoChanger_
(
changers
? fvMeshTopoChanger::New(*this)
: autoPtr<fvMeshTopoChanger>(nullptr)
),
distributor_
(
changers
? fvMeshDistributor::New(*this)
: autoPtr<fvMeshDistributor>(nullptr)
),
mover_
(
changers
? fvMeshMover::New(*this)
: autoPtr<fvMeshMover>(nullptr)
),
topoChanger_(nullptr),
distributor_(nullptr),
mover_(nullptr),
lduPtr_(nullptr),
curTimeIndex_(time().timeIndex()),
VPtr_(nullptr),
@ -294,6 +279,14 @@ Foam::fvMesh::fvMesh(const IOobject& io, const bool changers)
Pout<< FUNCTION_NAME << "Constructing fvMesh from IOobject" << endl;
}
// Construct changers
if (changers)
{
topoChanger_.set(fvMeshTopoChanger::New(*this).ptr());
distributor_.set(fvMeshDistributor::New(*this).ptr());
mover_.set(fvMeshMover::New(*this).ptr());
}
// Check the existence of the cell volumes and read if present
// and set the storage of V00
if (fileHandler().isFile(time().timePath()/"V0"))
@ -365,6 +358,9 @@ Foam::fvMesh::fvMesh
fvSolution(static_cast<const objectRegistry&>(*this)),
data(static_cast<const objectRegistry&>(*this)),
boundary_(*this, boundaryMesh()),
topoChanger_(nullptr),
distributor_(nullptr),
mover_(nullptr),
lduPtr_(nullptr),
curTimeIndex_(time().timeIndex()),
VPtr_(nullptr),
@ -407,6 +403,9 @@ Foam::fvMesh::fvMesh
fvSolution(static_cast<const objectRegistry&>(*this)),
data(static_cast<const objectRegistry&>(*this)),
boundary_(*this, boundaryMesh()),
topoChanger_(nullptr),
distributor_(nullptr),
mover_(nullptr),
lduPtr_(nullptr),
curTimeIndex_(time().timeIndex()),
VPtr_(nullptr),
@ -440,6 +439,9 @@ Foam::fvMesh::fvMesh
fvSolution(static_cast<const objectRegistry&>(*this)),
data(static_cast<const objectRegistry&>(*this)),
boundary_(*this),
topoChanger_(nullptr),
distributor_(nullptr),
mover_(nullptr),
lduPtr_(nullptr),
curTimeIndex_(time().timeIndex()),
VPtr_(nullptr),