Take care not to count initialization values when finding min/max

[libs/math.git] / source / geometry / transformedshape.h

diff --git a/source/geometry/transformedshape.h b/source/geometry/transformedshape.h
index 74912aa0b904cf1786c2efb7d78f6ec37c56d6d7..f214e4bb1aa9c578308088f2475879debc726378 100644 (file)
--- a/source/geometry/transformedshape.h
+++ b/source/geometry/transformedshape.h
@@ -88,10 +88,18 @@ inline BoundingBox<T, D> TransformedShape<T, D>::get_axis_aligned_bounding_box()
 
                point = transformation.transform(point);
 
-               for(unsigned j=0; j<D; ++j)
+               if(i==0)
+               {
+                       min_pt = point;
+                       max_pt = point;
+               }
+               else
                {
-                       min_pt[j] = std::min(min_pt[j], point[j]);
-                       max_pt[j] = std::max(min_pt[j], point[j]);
+                       for(unsigned j=0; j<D; ++j)
+                       {
+                               min_pt[j] = std::min(min_pt[j], point[j]);
+                               max_pt[j] = std::max(max_pt[j], point[j]);
+                       }
                }
        }