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

Experimenting with least square vector correction: off

Browse source
This commit is contained in:
Hrvoje Jasak 2015-04-24 14:11:11 +01:00
parent 7d149ebb2d
commit b3c099be07
1 changed files with 9 additions and 7 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

16
src/finiteVolume/finiteVolume/gradSchemes/leastSquaresGrad/leastSquaresVectors.C
View file

@ -30,7 +30,11 @@ License
// * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * // // * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
// * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * // namespace Foam
{
defineTypeNameAndDebug(leastSquaresVectors, 0);
}
const Foam::scalar Foam::leastSquaresVectors::smallDotProdTol_ const Foam::scalar Foam::leastSquaresVectors::smallDotProdTol_
( (
@ -38,12 +42,6 @@ const Foam::scalar Foam::leastSquaresVectors::smallDotProdTol_
); );
namespace Foam
{
defineTypeNameAndDebug(leastSquaresVectors, 0);
}
// * * * * * * * * * * * * * * * * Constructors * * * * * * * * * * * * * * // // * * * * * * * * * * * * * * * * Constructors * * * * * * * * * * * * * * //
Foam::leastSquaresVectors::leastSquaresVectors(const fvMesh& mesh) Foam::leastSquaresVectors::leastSquaresVectors(const fvMesh& mesh)
@ -221,6 +219,8 @@ void Foam::leastSquaresVectors::makeLeastSquaresVectors() const
} }
} }
# if 0
// Coefficient sign correction on least squares vectors // Coefficient sign correction on least squares vectors
// The sign of the coefficient must be positive to ensure correct // The sign of the coefficient must be positive to ensure correct
// behaviour in coupled systems. This is checked by evaluating // behaviour in coupled systems. This is checked by evaluating
@ -406,6 +406,8 @@ void Foam::leastSquaresVectors::makeLeastSquaresVectors() const
} }
} }
#endif
if (debug) if (debug)
{ {
Info<< "leastSquaresVectors::makeLeastSquaresVectors() :" Info<< "leastSquaresVectors::makeLeastSquaresVectors() :"