zhyang-dev/OpenFOAM-6
1
0
Fork 0
mirror of https://github.com/OpenFOAM/OpenFOAM-6.git synced 2025-12-08 06:57:46 +00:00
Code Issues Actions Packages Projects Releases Wiki Activity

blockMesh: Stabilised arc-cosine argument for 180 degree arcs

Browse Source
This commit is contained in:
Will Bainbridge
2018-05-01 11:26:29 +01:00
parent 87e32ab499
commit 06ade28f93
1 changed files with 2 additions and 1 deletions
Download Patch File Download Diff File Expand all files Collapse all files

3
src/mesh/blockMesh/blockEdges/arcEdge/arcEdge.C
View File

@ -72,7 +72,8 @@ Foam::cylindricalCS Foam::blockEdges::arcEdge::calcAngle()
vector r3(p3_ - centre);
// find angles
angle_ = radToDeg(acos((r3 & r1)/(mag(r3) * mag(r1))));
const scalar cosAngle = (r3 & r1)/(mag(r3) * mag(r1));
angle_ = radToDeg(acos(max(-1, min(cosAngle, 1))));
// check if the vectors define an exterior or an interior arcEdge
if (((r1 ^ r2) & (r1 ^ r3)) < 0.0)