diff --git a/src/OpenFOAM/primitives/Tensor/tensor/tensor.C b/src/OpenFOAM/primitives/Tensor/tensor/tensor.C index cadbef6a9b..2d8da5e71c 100644 --- a/src/OpenFOAM/primitives/Tensor/tensor/tensor.C +++ b/src/OpenFOAM/primitives/Tensor/tensor/tensor.C @@ -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);