Reverting new epsilons since they had bad global impact - instead just remove the...
authorDaniel Genrich <daniel.genrich@gmx.net>
Sun, 20 Apr 2008 09:57:28 +0000 (09:57 +0000)
committerDaniel Genrich <daniel.genrich@gmx.net>
Sun, 20 Apr 2008 09:57:28 +0000 (09:57 +0000)
intern/elbeem/intern/ntl_vector3dim.h
intern/elbeem/intern/solver_util.cpp

index 35a3ad59ec916cc5b171db800e6c79fec79fe73d..6f17f9f5e7b864c559cfa4a0c4f138630168e4bb 100644 (file)
@@ -833,7 +833,7 @@ typedef float gfxReal;
 //#define vecGfx2F(x) (x)
 //#define vecD2Gfx(x) vecD2F(x)
 //#define vecGfx2D(x) vecF2D(x)
-#define VECTOR_EPSILON (1.192092896e-07F)
+#define VECTOR_EPSILON (1e-5f)
 #else
 typedef double gfxReal;
 #define GFX_REAL_MAX __DBL_MAX__
@@ -841,7 +841,7 @@ typedef double gfxReal;
 //#define vecGfx2F(x) vecF2D(x)
 //#define vecD2Gfx(x) (x)
 //#define vecGfx2D(x) (x)
-#define VECTOR_EPSILON (2.2204460492503131e-016)
+#define VECTOR_EPSILON (1e-10)
 #endif
 
 /* fixed double prec. type, for epxlicitly double values */
index e74875a0dbfbd74d755f1e33cf830ba7261ceeca..25b8b98b3c6f7a8723d8b563e81e9ab64512da25 100644 (file)
@@ -40,7 +40,7 @@ void LbmFsgrSolver::prepareVisualization( void ) {
        int mainGravDir=6; // if normalizing fails, we asume z-direction gravity
        LbmFloat mainGravLen = 0.;
        FORDF1{
-               LbmFloat thisGravLen = dot(LbmVec(dfVecX[l],dfVecY[l],dfVecZ[l]), getNormalized(mLevel[mMaxRefine].gravity) );  
+               LbmFloat thisGravLen = dot(LbmVec(dfVecX[l],dfVecY[l],dfVecZ[l]), mLevel[mMaxRefine].gravity ); 
                if(thisGravLen>mainGravLen) {
                        mainGravLen = thisGravLen;
                        mainGravDir = l;