Changeset 4087cb for src/base/helper.cpp

Timestamp:

Oct 15, 2013, 2:55:12 PM (12 years ago)

Author:

Frederik Heber <heber@…>

Children:

fb3b4c

Parents:

b4ac80

git-author:

Frederik Heber <heber@…> (08/30/13 20:45:35)

git-committer:

Frederik Heber <heber@…> (10/15/13 14:55:12)

Message:

Increased numerical tolerance of assert in LevelOperatorFAS::Restrict().

• modified VMG::Vector::IsComponentWise..(), e.g. ..LessOr..() to accept equality within interval of width tol.

File:

• src/base/helper.cpp (modified) (1 diff)

src/base/helper.cpp

@@ -93 +93 @@
 bool Helper::AssertVectorsEqual(const Vector& pos_1, const Vector& pos_2, const vmg_float& tol)
 {
-  bool equal = (pos_1 - pos_2).Abs().IsComponentwiseLessOrEqual(tol);
+  // create factor with similar magnitude as pos1 but always greater equal 1
+  const VMG::Vector factor( abs(pos_1.X())+1, abs(pos_1.Y())+1, abs(pos_1.Z())+1);
+  const VMG::Vector temp = (pos_1 - pos_2).Abs() / factor;
+  bool equal = temp.IsComponentwiseLessOrEqual(tol);
   assert(equal);
   return equal;