research/OpenFOAM-4.x
1
0
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity Actions

tensor::eigenValues: Test the magnitude of 'QQ' before dividing by it.

Browse source 
Resolves bug-report http://bugs.openfoam.org/view.php?id=2351
This commit is contained in:
Henry Weller 2016-11-26 20:53:19 +00:00
parent fcb789a267
commit 422f1770e0
1 changed files with 1 additions and 1 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

2
src/OpenFOAM/primitives/Tensor/tensor/tensor.C
View file

@ -123,7 +123,7 @@ Foam::vector Foam::eigenValues(const tensor& t)
}
// Two identical roots and one distinct root
else if (mag(PPP/QQ - 1) < SMALL)
else if (mag(QQ) > SMALL && mag(PPP/QQ - 1) < SMALL)
{
scalar sqrtP = sqrt(P);
scalar signQ = sign(Q);