nondim tolerance

This commit is contained in:
mattijs
2008-08-14 09:26:47 +01:00
parent 64bc136106
commit 3e7a5b230f
2 changed files with 3 additions and 13 deletions

View File

@ -34,12 +34,6 @@ Description
#include "scalar.H" #include "scalar.H"
// * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
//- Not VSMALL - but smaller than SMALL
const Foam::scalar Foam::triangleFuncs::NVSMALL = Foam::sqrt(Foam::VSMALL);
// * * * * * * * * * * * * * Private Member Functions * * * * * * * * * * * // // * * * * * * * * * * * * * Private Member Functions * * * * * * * * * * * //
void Foam::triangleFuncs::setIntersection void Foam::triangleFuncs::setIntersection
@ -120,7 +114,8 @@ bool Foam::triangleFuncs::intersectAxesBundle
scalar det = v2*u1 - u2*v1; scalar det = v2*u1 - u2*v1;
if (Foam::mag(det) < SMALL) // Fix for V0:(-31.71428 0 -15.10714) V10:(-1.285715 8.99165e-16 -1.142858) V20:(0 0 -1.678573) i0:0
if (Foam::mag(det)/max(max(mag(V10),mag(V20)),1) < SMALL)
{ {
// Triangle parallel to dir // Triangle parallel to dir
return false; return false;
@ -522,7 +517,7 @@ bool Foam::triangleFuncs::classify
bool hit = false; bool hit = false;
if (Foam::mag(u1) < NVSMALL) if (Foam::mag(u1) < ROOTVSMALL)
{ {
beta = u0/u2; beta = u0/u2;

View File

@ -71,11 +71,6 @@ public:
private: private:
// Private data
//- Not VSMALL - but smaller than SMALL
static const scalar NVSMALL;
// Private Member Functions // Private Member Functions
//- Helper function for intersect. Sets pt to be anywhere on the edge //- Helper function for intersect. Sets pt to be anywhere on the edge