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

Limitting cosTheta to be inside [-1,1] interval

Browse source
This commit is contained in:
Vanja Skuric 2018-02-07 13:31:47 +01:00
parent f500917045
commit 9c11730d7c
1 changed files with 6 additions and 2 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

8
src/foam/primitives/Tensor/tensor/tensor.C
View file

@ -151,7 +151,9 @@ vector eigenValues(const tensor& t)
if (R2 < Q3)
{
scalar sqrtQ = sqrt(Q);
scalar theta = acos(R/(Q*sqrtQ));
scalar cosTheta = R/(Q*sqrtQ);
cosTheta += (neg(cosTheta) - pos(cosTheta))*SMALL;
scalar theta = acos(cosTheta);
scalar m2SqrtQ = -2*sqrtQ;
scalar aBy3 = a/3;
@ -406,7 +408,9 @@ vector eigenValues(const symmTensor& t)
if (R2 < Q3)
{
scalar sqrtQ = sqrt(Q);
scalar theta = acos(R/(Q*sqrtQ));
scalar cosTheta = R/(Q*sqrtQ);
cosTheta += (neg(cosTheta) - pos(cosTheta))*SMALL;
scalar theta = acos(cosTheta);
scalar m2SqrtQ = -2*sqrtQ;
scalar aBy3 = a/3;