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

Robustness improvement

Browse source
This commit is contained in:
Hrvoje Jasak 2011-02-27 02:19:31 +00:00
parent 5296a9fe1f
commit 3eb0abe907
1 changed files with 5 additions and 15 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

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

@ -317,24 +317,14 @@ vector eigenValues(const symmTensor& t)
// If there is a zero root
if (mag(c) < SMALL)
// if (mag(c) < 1.0e-100)
{
scalar disc = sqr(a) - 4*b;
const scalar disc = Foam::max(sqr(a) - 4*b, 0.0);
if (disc >= -SMALL)
{
scalar q = -0.5*sqrt(max(scalar(0), disc));
scalar q = -0.5*sqrt(max(scalar(0), disc));
i = 0;
ii = -0.5*a + q;
iii = -0.5*a - q;
}
else
{
FatalErrorIn("eigenValues(const tensor&)")
<< "zero and complex eigenvalues in tensor: " << t
<< abort(FatalError);
}
i = 0;
ii = -0.5*a + q;
iii = -0.5*a - q;
}
else
{