Changes in src/Line.cpp [112b09:f932b7]

File:

• src/Line.cpp (modified) (1 diff)

src/Line.cpp

@@ -215 +215 @@
 }
 
+std::vector<Vector> Line::getSphereIntersections() const{
+  std::vector<Vector> res;
+
+  // line is kept in normalized form, so we can skip a lot of calculations
+  double discriminant = 1-origin->NormSquared();
+  // we might have 2, 1 or 0 solutions, depending on discriminant
+  if(discriminant>=0){
+    if(discriminant==0){
+      res.push_back(*origin);
+    }
+    else{
+      Vector helper = sqrt(discriminant)*(*direction);
+      res.push_back(*origin+helper);
+      res.push_back(*origin-helper);
+    }
+  }
+  return res;
+}
+
 Line makeLineThrough(const Vector &x1, const Vector &x2){
   if(x1==x2){