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

SlicedGeometricField: Include empty patches in the splice function

Browse Source 
for consistency with slicing.
This commit is contained in:
Henry Weller
2020-03-03 09:07:15 +00:00
parent 779ed4de1f
commit d98efdbcb3
1 changed files with 25 additions and 9 deletions
Download Patch File Download Diff File Expand all files Collapse all files

34
src/OpenFOAM/fields/GeometricFields/SlicedGeometricField/SlicedGeometricField.C
View File

@ -416,9 +416,9 @@ splice() const
label completeSize = GeoMesh::size(mesh);
forAll(mesh.boundary(), patchi)
forAll(mesh.boundaryMesh(), patchi)
{
completeSize += mesh.boundary()[patchi].size();
completeSize += mesh.boundaryMesh()[patchi].size();
}
tmp<Field<Type>> tCompleteField(new Field<Type>(completeSize));
@ -429,16 +429,32 @@ splice() const
label start = GeoMesh::size(mesh);
forAll(mesh.boundary(), patchi)
forAll(mesh.boundaryMesh(), patchi)
{
typename Field<Type>::subField
if
(
completeField,
mesh.boundary()[patchi].size(),
start
) = this->boundaryField()[patchi];
mesh.boundary()[patchi].size()
== mesh.boundaryMesh()[patchi].size()
)
{
typename Field<Type>::subField
(
completeField,
mesh.boundary()[patchi].size(),
start
) = this->boundaryField()[patchi];
}
else
{
typename Field<Type>::subField
(
completeField,
mesh.boundaryMesh()[patchi].size(),
start
) = Zero;
}
start += mesh.boundary()[patchi].size();
start += mesh.boundaryMesh()[patchi].size();
}
return tCompleteField;