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

ENH: pointConstraints: allow non-facePointPatch. Fixes #3085

Browse Source
This commit is contained in:
mattijs
2024-01-15 15:23:09 +00:00
parent e3f0521b19
commit 26904bf6df
1 changed files with 24 additions and 8 deletions
Download Patch File Download Diff File Expand all files Collapse all files

32
src/finiteVolume/interpolation/volPointInterpolation/pointConstraints.C
View File

@ -56,8 +56,6 @@ void Foam::pointConstraints::makePatchPatchAddressing()
const polyMesh& mesh = pMesh();
const pointBoundaryMesh& pbm = pMesh.boundary();
const polyBoundaryMesh& bm = mesh.boundaryMesh();
// first count the total number of patch-patch points
@ -65,9 +63,15 @@ void Foam::pointConstraints::makePatchPatchAddressing()
forAll(pbm, patchi)
{
if (!isA<emptyPointPatch>(pbm[patchi]) && !pbm[patchi].coupled())
const auto* fpp = isA<facePointPatch>(pbm[patchi]);
if
(
fpp
&& !isA<emptyPointPatch>(pbm[patchi])
&& !pbm[patchi].coupled()
)
{
const labelList& bp = bm[patchi].boundaryPoints();
const labelList& bp = fpp->patch().boundaryPoints();
nPatchPatchPoints += bp.size();
@ -103,9 +107,15 @@ void Foam::pointConstraints::makePatchPatchAddressing()
forAll(pbm, patchi)
{
if (!isA<emptyPointPatch>(pbm[patchi]) && !pbm[patchi].coupled())
const auto* fpp = isA<facePointPatch>(pbm[patchi]);
if
(
fpp
&& !isA<emptyPointPatch>(pbm[patchi])
&& !pbm[patchi].coupled()
)
{
const labelList& bp = bm[patchi].boundaryPoints();
const labelList& bp = fpp->patch().boundaryPoints();
const labelList& meshPoints = pbm[patchi].meshPoints();
forAll(bp, pointi)
@ -163,9 +173,15 @@ void Foam::pointConstraints::makePatchPatchAddressing()
// Copy from patchPatch constraints into coupledConstraints.
forAll(pbm, patchi)
{
if (!isA<emptyPointPatch>(pbm[patchi]) && !pbm[patchi].coupled())
const auto* fpp = isA<facePointPatch>(pbm[patchi]);
if
(
fpp
&& !isA<emptyPointPatch>(pbm[patchi])
&& !pbm[patchi].coupled()
)
{
const labelList& bp = bm[patchi].boundaryPoints();
const labelList& bp = fpp->patch().boundaryPoints();
const labelList& meshPoints = pbm[patchi].meshPoints();
forAll(bp, pointi)