Implement maxval and minval reduction operations

1 files changed, 24 insertions, 0 deletions

diff --git a/vec_pseudo.h b/vec_pseudo.h
index b8aa808..0ee5a60 100644
--- a/vec_pseudo.h
+++ b/vec_pseudo.h
@@ -656,6 +656,18 @@ namespace vecmathlib {
       return res /= x;
     }
     
+    real_t maxval() const
+    {
+      real_t res = v[0];
+      for (int d=1; d<size; ++d) res = std::fmax(res, v[d]);
+      return res;
+    }
+    real_t minval() const
+    {
+      real_t res = v[0];
+      for (int d=1; d<size; ++d) res = std::fmin(res, v[d]);
+      return res;
+    }
     real_t prod() const
     {
       real_t res = R(1.0);
@@ -1072,6 +1084,18 @@ namespace vecmathlib {
   }
   
   template<typename real_t, int size>
+  inline real_t maxval(realpseudovec<real_t, size> x)
+  {
+    return x.maxval();
+  }
+  
+  template<typename real_t, int size>
+  inline real_t minval(realpseudovec<real_t, size> x)
+  {
+    return x.minval();
+  }
+  
+  template<typename real_t, int size>
   inline real_t prod(realpseudovec<real_t, size> x)
   {
     return x.prod();