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

tensor: corrected eigenvector ordering function to handle cases when sub-determinants are equal

Browse Source
This commit is contained in:
Henry
2014-02-07 09:53:59 +00:00
committed by Andrew Heather
parent ef52265131
commit 945e4fe938
1 changed files with 5 additions and 5 deletions
Download Patch File Download Diff File Expand all files Collapse all files

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

@ -2,7 +2,7 @@
========= | ========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | \\ / O peration |
\\ / A nd | Copyright (C) 2011-2013 OpenFOAM Foundation \\ / A nd | Copyright (C) 2011-2014 OpenFOAM Foundation
\\/ M anipulation | \\/ M anipulation |
------------------------------------------------------------------------------- -------------------------------------------------------------------------------
License License
@ -223,7 +223,7 @@ Foam::vector Foam::eigenVector(const tensor& t, const scalar lambda)
scalar magSd2 = mag(sd2); scalar magSd2 = mag(sd2);
// Evaluate the eigenvector using the largest sub-determinant // Evaluate the eigenvector using the largest sub-determinant
if (magSd0 > magSd1 && magSd0 > magSd2 && magSd0 > SMALL) if (magSd0 >= magSd1 && magSd0 >= magSd2 && magSd0 > SMALL)
{ {
vector ev vector ev
( (
@ -235,7 +235,7 @@ Foam::vector Foam::eigenVector(const tensor& t, const scalar lambda)
return ev; return ev;
} }
else if (magSd1 > magSd2 && magSd1 > SMALL) else if (magSd1 >= magSd2 && magSd1 > SMALL)
{ {
vector ev vector ev
( (
@ -418,7 +418,7 @@ Foam::vector Foam::eigenVector(const symmTensor& t, const scalar lambda)
scalar magSd2 = mag(sd2); scalar magSd2 = mag(sd2);
// Evaluate the eigenvector using the largest sub-determinant // Evaluate the eigenvector using the largest sub-determinant
if (magSd0 > magSd1 && magSd0 > magSd2 && magSd0 > SMALL) if (magSd0 >= magSd1 && magSd0 >= magSd2 && magSd0 > SMALL)
{ {
vector ev vector ev
( (
@ -430,7 +430,7 @@ Foam::vector Foam::eigenVector(const symmTensor& t, const scalar lambda)
return ev; return ev;
} }
else if (magSd1 > magSd2 && magSd1 > SMALL) else if (magSd1 >= magSd2 && magSd1 > SMALL)
{ {
vector ev vector ev
( (