Changes in src/vector.cpp [7b36fe:fee69b]

File:

• src/vector.cpp (modified) (5 diffs)

src/vector.cpp

@@ -222 +222 @@
  * \param *Origin first vector of line
  * \param *LineVector second vector of line
- * \return true -  \a this contains intersection point on return, false - line is parallel to plane (even if in-plane)
+ * \return true -  \a this contains intersection point on return, false - line is parallel to plane
  */
 bool Vector::GetIntersectionWithPlane(const Vector * const PlaneNormal, const Vector * const PlaneOffset, const Vector * const Origin, const Vector * const LineVector)
@@ -235 +235 @@
   Direction.Normalize();
   Log() << Verbose(1) << "INFO: Direction is " << Direction << "." << endl;
-  //Log() << Verbose(1) << "INFO: PlaneNormal is " << *PlaneNormal << " and PlaneOffset is " << *PlaneOffset << "." << endl;
   factor = Direction.ScalarProduct(PlaneNormal);
-  if (fabs(factor) < MYEPSILON) { // Uniqueness: line parallel to plane?
-    Log() << Verbose(1) << "BAD: Line is parallel to plane, no intersection." << endl;
+  if (factor < MYEPSILON) { // Uniqueness: line parallel to plane?
+    eLog() << Verbose(2) << "Line is parallel to plane, no intersection." << endl;
     return false;
   }
@@ -244 +243 @@
   helper.SubtractVector(Origin);
   factor = helper.ScalarProduct(PlaneNormal)/factor;
-  if (fabs(factor) < MYEPSILON) { // Origin is in-plane
-    Log() << Verbose(1) << "GOOD: Origin of line is in-plane." << endl;
+  if (factor < MYEPSILON) { // Origin is in-plane
+    Log() << Verbose(1) << "Origin of line is in-plane, simple." << endl;
     CopyVector(Origin);
     return true;
@@ -259 +258 @@
   helper.SubtractVector(PlaneOffset);
   if (helper.ScalarProduct(PlaneNormal) < MYEPSILON) {
-    Log() << Verbose(1) << "GOOD: Intersection is " << *this << "." << endl;
+    Log() << Verbose(1) << "INFO: Intersection at " << *this << " is good." << endl;
     return true;
   } else {
@@ -354 +353 @@
   Vector parallel;
   double factor = 0.;
+  double pfactor = 0.;
   if (fabs(a.ScalarProduct(&b)*a.ScalarProduct(&b)/a.NormSquared()/b.NormSquared() - 1.) < MYEPSILON) {
     parallel.CopyVector(Line1a);