diff --git a/src/include/pack_num.h b/src/include/pack_num.h
index f197c6e..a9d4adb 100644
--- a/src/include/pack_num.h
+++ b/src/include/pack_num.h
@@ -37,6 +37,14 @@ namespace MultiArrayHelper
     template <size_t N>
     struct PackNum
     {
+	template <class MA, class ITuple, class... IPtrs>
+	static auto mkElemOperation(const MA& ma, const ITuple& ituple, IPtrs... iptrs)
+	    -> decltype(PackNum<N-1>::mkElemOperation(ma, ituple, std::get<N>(ituple), iptrs...))
+	{
+	    return PackNum<N-1>::mkElemOperation(ma, ituple, std::get<N>(ituple), iptrs...);
+	}
+
+	
 	template <typename... T>
 	static void printTuple(std::ostream& out, const std::tuple<T...>& tp)
 	{
@@ -94,6 +102,13 @@ namespace MultiArrayHelper
     template<>
     struct PackNum<0>
     {
+	template <class MA, class ITuple, class... IPtrs>
+	static auto mkElemOperation(const MA& ma, const ITuple& ituple, IPtrs... iptrs)
+	    -> decltype(ma(iptrs...))
+	{
+	    return ma(iptrs...);
+	}
+	
 	template <typename... T>
 	static void printTuple(std::ostream& out, const std::tuple<T...>& tp)
 	{	    
diff --git a/src/include/type_operations.h b/src/include/type_operations.h
index 3491378..27a65dd 100644
--- a/src/include/type_operations.h
+++ b/src/include/type_operations.h
@@ -27,7 +27,7 @@ namespace MultiArrayTools
 	inline auto apply(const MultiArray<T,Ranges...>& ma)
 	    -> OperationRoot<T,Ranges...>
 	{
-	    return mkElemOperation(ma, ituple); // -> pack_num
+	    return PackNum<sizeof...(Ranges)-1>::mkElemOperation(ma, ituple);
 	}
 
     private: