Changeset d403a1 for molecuilder/src/tesselation.cpp

Timestamp:

Jan 11, 2010, 3:06:30 AM (16 years ago)

Author:

Frederik Heber <heber@…>

Children:

91b1e79

Parents:

23f6ec

Message:

InsideOutside unit test of tesselation is working correctly.

• FIX: BoundaryTriangleSet::GetIntersectionInsideTriangle() - don't need helper, just check whether CrossPoint is returned (and true) for all of the three sides.
• FIX: Tesselation::IsInnerPoint() - projection onto plane and stuff was nonsense, just take the Point ans Intersection which is on the plane anyway.
• FIX: Vector::GetIntersectionOfTwoLinesOnPlane() - coefficient MUST be zero (then vectors are coplanar), but parallelity check was missing. Also, we have to check whether s is in [0,1] in order to see whether we are inside the triangle side or outside.

Signed-off-by: Frederik Heber <heber@tabletINS.(none)>

File:

• molecuilder/src/tesselation.cpp (modified) (4 diffs)

molecuilder/src/tesselation.cpp

@@ -420 +420 @@
     if (CrossPoint.GetIntersectionOfTwoLinesOnPlane(endpoints[i%3]->node->node, endpoints[(i+1)%3]->node->node, endpoints[(i+2)%3]->node->node, Intersection, &NormalVector)) {
       CrossPoint.SubtractVector(endpoints[i%3]->node->node);  // cross point was returned as absolute vector
-      helper.CopyVector(endpoints[(i+1)%3]->node->node);
-      helper.SubtractVector(endpoints[i%3]->node->node);
+      i++;
+    } else 
       break;
-    } else 
-      i++;
   } while (i<3);
   if (i==3) {
-    eLog() << Verbose(0) << "Could not find any cross points, something's utterly wrong here!" << endl;
-  }
-  Log() << Verbose(1) << "INFO: Crosspoint is " << CrossPoint << "." << endl;
-
-  // check whether intersection is inside or not by comparing length of intersection and length of cross point
-  if ((CrossPoint.NormSquared() - helper.NormSquared()) < MYEPSILON) { // inside
+    Log() << Verbose(1) << "INFO: Crosspoint " << CrossPoint << " inside of triangle." << endl;
     return true;
-  } else { // outside!
-    Intersection->Zero();
+  } else {
+    Log() << Verbose(1) << "INFO: Crosspoint " << CrossPoint << " outside of triangle." << endl;
     return false;
   }
@@ -3241 +3234 @@
 bool Tesselation::IsInnerPoint(const Vector &Point, const LinkedCell* const LC, const double epsilon) const
 {
-        Info FunctionInfo(__func__);
+  Info FunctionInfo(__func__);
   class BoundaryTriangleSet *result = FindClosestTriangleToVector(&Point, LC);
   Vector Center;
@@ -3252 +3245 @@
     return false;
   } else {
-    Log() << Verbose(1) << "INFO: Closest triangle found is " << *result << "." << endl;
+    Log() << Verbose(1) << "INFO: Closest triangle found is " << *result << " with normal vector " << result->NormalVector << "." << endl;
   }
 
@@ -3259 +3252 @@
   DistanceToCenter.CopyVector(&Center);
   DistanceToCenter.SubtractVector(&Point);
-  Log() << Verbose(2) << "INFO: Vector from point to test to center is " << Center << "." << endl;
+  Log() << Verbose(2) << "INFO: Vector from point to test to center is " << DistanceToCenter << "." << endl;
 
   // check whether we are on boundary
   if (fabs(DistanceToCenter.ScalarProduct(&result->NormalVector)) < MYEPSILON) {
     // calculate whether inside of triangle
-    DistanceToCenter.ProjectOntoPlane(&result->NormalVector);
-    DistanceToCenter.AddVector(&Center);
-    Center.CopyVector(&DistanceToCenter);
+    DistanceToCenter.CopyVector(&Point);
+    Center.CopyVector(&Point);
     Center.SubtractVector(&result->NormalVector); // points towards MolCenter
     DistanceToCenter.AddVector(&result->NormalVector); // points outside