diff --git a/cwVectOps.h b/cwVectOps.h
index 8aa52e4..8aae087 100644
--- a/cwVectOps.h
+++ b/cwVectOps.h
@@ -389,6 +389,16 @@ namespace cw
       return y;
     }
 
+    template< typename T >
+      T abs_sum( const T* v, unsigned n )
+    {
+      T y = 0;
+      for(unsigned i=0; i<n; ++i)
+        y += std::abs(v[i]);
+
+      return y;
+    }
+    
     template< typename T >
       T prod( const T* v, unsigned n )
     {