zhyang-dev/openfoam
1
0
Fork 0
mirror of https://develop.openfoam.com/Development/openfoam.git synced 2025-11-28 03:28:01 +00:00
Code Issues Actions Packages Projects Releases Wiki Activity

ENH: coupledPolyPatch: truncate tolerance calculation for zero area faces.

Browse Source 
Check duplicate anchor points
This commit is contained in:
mattijs
2012-08-20 11:44:04 +01:00
parent 0b97b9aeaf
commit b4da6c5b7e
1 changed files with 30 additions and 2 deletions
Download Patch File Download Diff File Expand all files Collapse all files

32
src/OpenFOAM/meshes/polyMesh/polyPatches/basic/coupled/coupledPolyPatch.C
View File

@ -172,7 +172,7 @@ Foam::scalarField Foam::coupledPolyPatch::calcFaceTol
maxLenSqr = max(maxLenSqr, magSqr(pt - cc)); maxLenSqr = max(maxLenSqr, magSqr(pt - cc));
maxCmpt = max(maxCmpt, cmptMax(cmptMag(pt))); maxCmpt = max(maxCmpt, cmptMax(cmptMag(pt)));
} }
tols[faceI] = max(SMALL*maxCmpt, Foam::sqrt(maxLenSqr)); tols[faceI] = max(SMALL, max(SMALL*maxCmpt, Foam::sqrt(maxLenSqr)));
} }
return tols; return tols;
} }
@ -200,12 +200,40 @@ Foam::label Foam::coupledPolyPatch::getRotation
} }
} }
if (anchorFp == -1 || mag(minDistSqr) > tol) if (anchorFp == -1 || Foam::sqrt(minDistSqr) > tol)
{ {
return -1; return -1;
} }
else else
{ {
// Check that anchor is unique.
forAll(f, fp)
{
scalar distSqr = magSqr(anchor - points[f[fp]]);
if (distSqr == minDistSqr && fp != anchorFp)
{
WarningIn
(
"label coupledPolyPatch::getRotation\n"
"(\n"
" const pointField&,\n"
" const face&,\n"
" const point&,\n"
" const scalar\n"
")"
) << "Cannot determine unique anchor point on face "
<< UIndirectList<point>(points, f)
<< endl
<< "Both at index " << anchorFp << " and " << fp
<< " the vertices have the same distance "
<< Foam::sqrt(minDistSqr)
<< " to the anchor " << anchor
<< ". Continuing but results might be wrong."
<< endl;
}
}
// Positive rotation // Positive rotation
return (f.size() - anchorFp) % f.size(); return (f.size() - anchorFp) % f.size();
} }