diff --git a/src/include/multi_array_operation.cc.h b/src/include/multi_array_operation.cc.h
index 723936d..221ee33 100644
--- a/src/include/multi_array_operation.cc.h
+++ b/src/include/multi_array_operation.cc.h
@@ -98,13 +98,13 @@ namespace MultiArrayTools
     }
 
     template <typename T, class OperationClass>
-    template <template <class> class AOp, class... Indices>
-    auto OperationBase<T,OperationClass>::oto(const std::shared_ptr<Indices>&... inds) const
+    template <class... Indices>
+    auto OperationBase<T,OperationClass>::addto(const std::shared_ptr<Indices>&... inds) const
 	-> MultiArray<T,typename Indices::RangeType...>
     {
 	MultiArray<T,typename Indices::RangeType...> out(inds->range()...,
                                                          static_cast<T>(0));
-        AOp<decltype(out(inds...)),decltype(THIS())>::selfApply(out(inds...), THIS());
+        out(inds...) += THIS();
 	return out;
     }
 
diff --git a/src/include/multi_array_operation.h b/src/include/multi_array_operation.h
index 7be4d8a..c8f1734 100644
--- a/src/include/multi_array_operation.h
+++ b/src/include/multi_array_operation.h
@@ -69,8 +69,8 @@ namespace MultiArrayTools
 	auto to(const std::shared_ptr<Indices>&... inds) const
 	    -> MultiArray<T,typename Indices::RangeType...>;
 
-        template <template <class> class AOp, class... Indices>
-	auto oto(const std::shared_ptr<Indices>&... inds) const
+        template <class... Indices>
+	auto addto(const std::shared_ptr<Indices>&... inds) const
 	    -> MultiArray<T,typename Indices::RangeType...>;
 
         template <typename R, class... Args> // Args = Operation Classes