CarpetLib: Add ipow function to vect class

darcs-hash:20041225183147-891bb-84276981b1deab37af4c3b541f993e62d63200e6.gz

1 files changed, 8 insertions, 0 deletions

diff --git a/Carpet/CarpetLib/src/vect.hh b/Carpet/CarpetLib/src/vect.hh
index 45599057b..d6e84ec46 100644
--- a/Carpet/CarpetLib/src/vect.hh
+++ b/Carpet/CarpetLib/src/vect.hh
@@ -564,6 +564,14 @@ inline vect<T,D> pow (const vect<T,D>& a, const vect<U,D>& b) {
   return r;
 }
 
+/** Return the element-wise integer power of two vectors.  */
+template<class T,int D>
+inline vect<T,D> ipow (const vect<T,D>& a, const vect<int,D>& b) {
+  vect<T,D> r;
+  for (int d=0; d<D; ++d) r[d]=ipow(a[d],b[d]);
+  return r;
+}
+
 
 
 // Reduction operators