Changeset 0e9394 for molecuilder/src/Plane.cpp

Timestamp:

Apr 30, 2010, 9:03:57 AM (16 years ago)

Author:

Tillmann Crueger <crueger@…>

Children:

7a8319

Parents:

c53e0b

git-author:

Tillmann Crueger <crueger@…> (04/30/10 08:26:28)

git-committer:

Tillmann Crueger <crueger@…> (04/30/10 09:03:57)

Message:

Added a class Space to represent general subspaces of R^{3 (planes and lines)}

File:

• molecuilder/src/Plane.cpp (modified) (3 diffs)

molecuilder/src/Plane.cpp

@@ -8 +8 @@
 #include "Plane.hpp"
 #include "vector.hpp"
+#include "defs.hpp"
 #include "Exceptions/LinearDependenceException.hpp"
 #include "info.hpp"
@@ -13 +14 @@
 #include "verbose.hpp"
 #include "Helpers/Assert.hpp"
+#include <cmath>
 
 /**
@@ -163 +165 @@
   return res;
 };
+
+/************ Methods inherited from Space ****************/
+
+double Plane::distance(Vector &point){
+  double res = point.ScalarProduct(*normalVector)-offset;
+  return fabs(res);
+}
+
+Vector Plane::getClosestPoint(Vector &point){
+  Vector difference = distance(point) * (*normalVector);
+  if(difference.IsZero()){
+    // the point itself lies on the plane
+    return point;
+  }
+  // get the direction this vector is pointing
+  double sign = difference.ScalarProduct(*normalVector);
+  // sign cannot be zero, since normalVector and difference are both != zero
+  sign = sign/fabs(sign);
+  return (point - (sign * difference));
+}
+
+bool Plane::isContained(Vector &point){
+  return (point.ScalarProduct(*normalVector) - offset) < MYEPSILON;
+}