Add ceil function for our Vector class

1 files changed, 11 insertions, 0 deletions

diff --git a/src/utilities/vectorHelpers.h b/src/utilities/vectorHelpers.h
index 17453dd..e4fc958 100644
--- a/src/utilities/vectorHelpers.h
+++ b/src/utilities/vectorHelpers.h
@@ -188,6 +188,17 @@ Vector<T,Size> ceil(const Vector<T,Size>& a)
   return out;
 }
 
+/// applies ceil to each component of a vector
+template <typename T, unsigned Size>
+Vector<T,Size> ceil(const Vector<T,Size>& a)
+{
+  Vector<T,Size> out;
+  for (unsigned int iDim=0; iDim < Size; ++iDim) {
+    out[iDim] = std::ceil(a[iDim]);
+  }
+  return out;
+}
+
 /*
 /// algorithm by Möller–Trumbore (TODO add ref), implemented by Lucas Cruz and Mathias J. Krause
 /// returns true if there is an intersection of a triangle given by (point0, point1, point1) and a ray given by its origin and direction and computes the distance