griefith/test2
1
0
Fork 0
Code Issues Pull Requests Actions 1 Packages Projects Releases Wiki Activity

Reverting new epsilons since they had bad global impact - instead just remove the problematic normalize() call

Browse Source
This commit is contained in:
Daniel Genrich
2008-04-20 09:57:28 +00:00
parent f265eeed8a
commit b43f24a944
2 changed files with 3 additions and 3 deletions
Download Patch File Download Diff File Expand all files Collapse all files

4
intern/elbeem/intern/ntl_vector3dim.h
View 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 */

2
intern/elbeem/intern/solver_util.cpp
View 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;