Fix the bounding sphere algorithm

diff --git a/source/geometry/boundingsphere.h b/source/geometry/boundingsphere.h
index dd20d96b6a4b701206f9e27b56a16cfc1402c158..b27e3a1bea7813bc3454f919471f8a1b79c09ca3 100644 (file)
--- a/source/geometry/boundingsphere.h
+++ b/source/geometry/boundingsphere.h
@@ -83,8 +83,8 @@ BoundingSphere<T, N> BoundingSphere<T, N>::from_point_cloud(const Iter &begin, c
                if(d>sqdist)
                {
                        d = sqrt(d);
-                       bsphere.center += v*(1-bsphere.radius/d);
-                       bsphere.radius += d/2;
+                       bsphere.center += v*((T(1)-bsphere.radius/d)/T(2));
+                       bsphere.radius = (bsphere.radius+d)/T(2);
                        sqdist = bsphere.radius*bsphere.radius;
                }
        }