diff --git a/src/include/multi_array_operation.cc.h b/src/include/multi_array_operation.cc.h
index 26ef0cf..4b1abed 100644
--- a/src/include/multi_array_operation.cc.h
+++ b/src/include/multi_array_operation.cc.h
@@ -89,10 +89,10 @@ namespace MultiArrayTools
 
     template <typename T, class OperationClass>
     template <class... Indices>
-    auto OperationBase<T,OperationClass>::operator()(const std::shared_ptr<Indices>&... inds) const
+    auto OperationBase<T,OperationClass>::to(const std::shared_ptr<Indices>&... inds) const
 	-> MultiArray<T,typename Indices::RangeType...>
     {
-	MultiArray<T,typename Indices::RangeType...> out(inds->range()... , static_cast<T>(0));
+	MultiArray<T,typename Indices::RangeType...> out(inds->range()...);
 	out(inds...) = THIS();
 	return out;
     }
diff --git a/src/include/multi_array_operation.h b/src/include/multi_array_operation.h
index ec6fb5d..a5a3dfe 100644
--- a/src/include/multi_array_operation.h
+++ b/src/include/multi_array_operation.h
@@ -66,7 +66,7 @@ namespace MultiArrayTools
 	    -> ConstOperationRoot<T,typename Indices::RangeType...>;
 
 	template <class... Indices>
-	auto operator()(const std::shared_ptr<Indices>&... inds) const
+	auto to(const std::shared_ptr<Indices>&... inds) const
 	    -> MultiArray<T,typename Indices::RangeType...>;
 
         template <typename R, class... Args> // Args = Operation Classes