zhyang-dev/openfoam
1
0
Fork 0
mirror of https://develop.openfoam.com/Development/openfoam.git synced 2025-11-28 03:28:01 +00:00
Code Issues Actions Packages Projects Releases Wiki Activity

ENH: tensor: stabilise eigenvectors

Browse Source
This commit is contained in:
mattijs
2014-07-07 11:43:08 +01:00
committed by Andrew Heather
parent 2ba526869a
commit 866ca15d2f
1 changed files with 2 additions and 2 deletions
Download Patch File Download Diff File Expand all files Collapse all files

4
src/OpenFOAM/primitives/Tensor/tensor/tensor.C
View File

@ -150,7 +150,7 @@ Foam::vector Foam::eigenValues(const tensor& t)
if (R2 < Q3) if (R2 < Q3)
{ {
scalar sqrtQ = sqrt(Q); scalar sqrtQ = sqrt(Q);
scalar theta = acos(R/(Q*sqrtQ)); scalar theta = acos(min(1.0, max(-1.0, R/(Q*sqrtQ))));
scalar m2SqrtQ = -2*sqrtQ; scalar m2SqrtQ = -2*sqrtQ;
scalar aBy3 = a/3; scalar aBy3 = a/3;
@ -345,7 +345,7 @@ Foam::vector Foam::eigenValues(const symmTensor& t)
if (R2 < Q3) if (R2 < Q3)
{ {
scalar sqrtQ = sqrt(Q); scalar sqrtQ = sqrt(Q);
scalar theta = acos(R/(Q*sqrtQ)); scalar theta = acos(min(1.0, max(-1.0, R/(Q*sqrtQ))));
scalar m2SqrtQ = -2*sqrtQ; scalar m2SqrtQ = -2*sqrtQ;
scalar aBy3 = a/3; scalar aBy3 = a/3;