added xfmBox to match xfmVector and xfmPoint etc

Author: ingowald

cuBQL/math/affine.h

@@ -181,6 +181,25 @@ namespace cuBQL {
                                               const typename AffineSpaceT<L>::vector_t& n)
   { return xfmNormal(m.l,n); }
 
+  template<typename T>
+  box_t<T,3> xfmBox(const AffineSpaceT<LinearSpace3<vec_t<T,3>>> &xfm,
+                    const box_t<T,3> &box)
+  {
+    box_t<T,3> out;
+    if (in.empty()) return out;
+    auto l = in.lower;
+    auto u = in.upper;
+    out.extend(xfmPoint(xfm,vec_t<T,3>(l.x,l.y,l.z)));
+    out.extend(xfmPoint(xfm,vec_t<T,3>(l.x,l.y,u.z)));
+    out.extend(xfmPoint(xfm,vec_t<T,3>(l.x,u.y,l.z)));
+    out.extend(xfmPoint(xfm,vec_t<T,3>(l.x,u.y,u.z)));
+    out.extend(xfmPoint(xfm,vec_t<T,3>(u.x,l.y,l.z)));
+    out.extend(xfmPoint(xfm,vec_t<T,3>(u.x,l.y,u.z)));
+    out.extend(xfmPoint(xfm,vec_t<T,3>(u.x,u.y,l.z)));
+    out.extend(xfmPoint(xfm,vec_t<T,3>(u.x,u.y,u.z)));
+    return out;
+  }
+
 
   ////////////////////////////////////////////////////////////////////////////////
   /// Comparison Operators