exaFOAM/foam-extend4.1-coherent-io
Archived
1
1
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity

Bugfix: Round-off stabilisation in sqrt

Browse source
This commit is contained in:
Hrvoje Jasak 2016-03-17 10:13:30 +00:00
parent 5673cec54b
commit 1a7654b0a6
1 changed files with 3 additions and 1 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

4
src/ODE/sixDOF/finiteRotation/finiteRotation.C
View file

@ -82,7 +82,9 @@ Foam::vector Foam::finiteRotation::eulerAngles(const tensor& rotT)
// Calculate roll angle // Calculate roll angle
rollAngle = atan2(rotT.yz(), rotT.zz()); rollAngle = atan2(rotT.yz(), rotT.zz());
const scalar c2 = sqrt(rotT.xx() + rotT.xy()); // Use mag to avoid negative value due to round-off
// HJ, 24/Feb/2016
const scalar c2 = sqrt(Foam::max(0, rotT.xx() + rotT.xy()));
// Calculate pitch angle // Calculate pitch angle
pitchAngle = atan2(-rotT.xz(), c2); pitchAngle = atan2(-rotT.xz(), c2);