From b257ab4db399337a53c3c9e05509229eda76c358 Mon Sep 17 00:00:00 2001
From: Christian Zimmermann <christian.zimmermann@ur.de>
Date: Thu, 13 Sep 2018 17:42:16 +0200
Subject: [PATCH] some utilities

---
 src/include/helper_tools.h    | 18 ++++++++++++++
 src/include/type_operations.h | 45 +++++++++++++++++++++++++++++++++++
 2 files changed, 63 insertions(+)

diff --git a/src/include/helper_tools.h b/src/include/helper_tools.h
index 7a40a42..4e8c8c1 100644
--- a/src/include/helper_tools.h
+++ b/src/include/helper_tools.h
@@ -52,6 +52,12 @@ namespace MultiArrayTools
 
     template <class IndexType>
     inline void For(const std::shared_ptr<IndexType>& ind, const std::function<void(void)>& ll);
+
+    template <class Index>
+    inline auto mkOp(const std::shared_ptr<Index>& i)
+	-> decltype(std::declval<FunctionalMultiArray<typename Index::MetaType,
+		    identity<typename Index::MetaType>,typename Index::RangeType> >
+		    ().exec(i));
 }
 
 /* ========================= *
@@ -151,6 +157,18 @@ namespace MultiArrayTools
 	}
     }
 
+    template <class Index>
+    inline auto mkOp(const std::shared_ptr<Index>& i)
+	-> decltype(std::declval<FunctionalMultiArray<typename Index::MetaType,
+		    identity<typename Index::MetaType>,typename Index::RangeType> >
+		    ().exec(i))
+    {
+	FunctionalMultiArray<typename Index::MetaType,
+			     identity<typename Index::MetaType>,
+			     typename Index::RangeType> fma(i->range());
+	return fma.exec(i);
+    }
+
 }
 
 #endif
diff --git a/src/include/type_operations.h b/src/include/type_operations.h
index 27a65dd..24e634b 100644
--- a/src/include/type_operations.h
+++ b/src/include/type_operations.h
@@ -3,6 +3,8 @@
 #define __type_operations_h__
 
 #include <cstdlib>
+#include <vector>
+#include <algorithm>
 
 #include "base_def.h"
 #include "mbase_def.h"
@@ -16,6 +18,8 @@ namespace MultiArrayTools
 	using namespace MultiArrayHelper;
     }
 
+    // MultiArray
+    
     template <typename T, class... Ranges>
     class operate
     {
@@ -54,7 +58,48 @@ namespace MultiArrayTools
 
     };
     
+    // vector
 
+    template <typename T>
+    class getter
+    {
+    private:
+	size_t mPos;
+    public:
+	static constexpr bool FISSTATIC = false;
+
+	getter(size_t i) : mPos(i) {}
+
+	inline T operator()(const std::vector<T>& in)
+	{
+	    return in[mPos];
+	}
+    };
+
+    template <typename T>
+    std::vector<T>& operator+=(std::vector<T>& a, const std::vector<T>& b)
+    {
+	std::transform(a.begin(), a.end(), b.begin(), a.begin(), std::plus<T>());
+	return a;
+    }
+    
+    template <class OperationClass, typename T>
+    class OperationTemplate<std::vector<T>,OperationClass> : public OperationBase<std::vector<T>,OperationClass>
+    {
+    public:
+	typedef OperationBase<std::vector<T>,OperationClass> OB;
+
+	auto operator[](size_t i)
+	    -> Operation<T,getter<T>,OperationClass>
+	{
+	    std::shared_ptr<getter<T> > ff = std::make_shared<getter<T> >(i);
+	    return Operation<T,getter<T>,OperationClass>(ff,OB::THIS());
+	}
+	
+    private:
+	OperationTemplate() = default;
+	friend OperationClass;
+    };
 
 } // namespace MultiArrayTools