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 64b1dbc341
commit 3ae0021235
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 there is a zero root
if (mag(c) < SMALL) 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; i = 0;
ii = -0.5*a + q; ii = -0.5*a + q;
iii = -0.5*a - q; iii = -0.5*a - q;
}
else
{
FatalErrorIn("eigenValues(const tensor&)")
<< "zero and complex eigenvalues in tensor: " << t
<< abort(FatalError);
}
} }
else else
{ {