Scene::CalcBounds(), including Box::Extend(Box&), Box::Clear()

[MicroTrace.git] / Sphere.cxx

diff --git a/Sphere.cxx b/Sphere.cxx
index 5706cc3..137ce8f 100644 (file)
--- a/Sphere.cxx
+++ b/Sphere.cxx
@@ -26,29 +26,41 @@ Sphere::Intersect(Ray& ray)
   if(t > ray.t())
     return false;
 
-  if( t < 1e-6 )
+  if( t < Epsilon )
     {
       t = (-B+root)/(2.0f*A);
-      if( t < 1e-6 || t > ray.t())
+      if( t < Epsilon || t > ray.t())
        return false;
     }
+
   ray.setT(t);
+  ray.setHit(this);
+
   return true;
 }
 
 Vec3f
 Sphere::GetNormal(Ray& ray)
 {
-  // We don't want to modify the ray (probably this is not needed, but I am
-  // too lazy to think about it now ...)
-  Ray tempRay = ray;
-
-  // Surface normal is the difference between intersection and center point
-  Intersect(tempRay);
-  // intersection
-  Vec3f i = tempRay.origin() + tempRay.direction() * tempRay.t();
-  // normal
-  Vec3f n = (i - m_center);
-  n.normalize();
-  return n;
+  Vec3f p = ray.origin()+ray.direction()*ray.t();
+  Vec3f N = (p-m_center);
+  N.normalize();
+
+  return N;
+}
+
+Box
+Sphere::CalcBounds()
+{
+  Vec3f min(m_center[0] - m_radius/2, m_center[1] - m_radius/2,
+    m_center[2] - m_radius/2);
+  Vec3f max(m_center[0] + m_radius/2, m_center[1] + m_radius/2,
+    m_center[2] + m_radius/2);
+  return Box(min, max);
+}
+
+bool
+Sphere::InVoxel(const Box& box)
+{
+  return CalcBounds().Overlaps(box);
 }