make first operation test run (index copying not supported at the moment)

src/container_range.cc

@@ -14,7 +14,7 @@ namespace MultiArrayTools
     /**********************
      *   ContainerIndex   *
      **********************/
-
+    /*
     template <class... Indices>
     ContainerIndex<Indices...>::ContainerIndex(const ContainerIndex& in) :
 	IndexInterface<std::tuple<typename Indices::MetaType...> >(in)
@@ -31,7 +31,7 @@ namespace MultiArrayTools
 	IB::mPos = PackNum<sizeof...(Indices)-1>::makePos(mIPack);
 	return *this;
     }
-
+    */
     template <class... Indices>
     template <class MRange>
     ContainerIndex<Indices...>::ContainerIndex(const std::shared_ptr<MRange>& range) :
@@ -145,7 +145,7 @@ namespace MultiArrayTools
     }
 
     template <class... Indices>
-    const std::shared_ptr<typename ContainerIndex<Indices...>::RangeType>& ContainerIndex<Indices...>::range() const;
+    std::shared_ptr<typename ContainerIndex<Indices...>::RangeType> ContainerIndex<Indices...>::range() const
     {
 	return std::dynamic_pointer_cast<RangeType>( IB::mRangePtr );
     }

src/container_range.h

@@ -32,8 +32,8 @@ namespace MultiArrayTools
     public:
 	ContainerIndex() = delete;
 
-	ContainerIndex(const ContainerIndex& in);
+	//ContainerIndex(const ContainerIndex& in);
-	ContainerIndex& operator=(const ContainerIndex& in);
+	//ContainerIndex& operator=(const ContainerIndex& in);
 
 	template <class MRange>
 	ContainerIndex(const std::shared_ptr<MRange>& range);
@@ -62,7 +62,7 @@ namespace MultiArrayTools
 
 	ContainerIndex& operator()(); // -> sync; just to shorten the code
 
-	const std::shared_ptr<RangeType>& range() const;
+	std::shared_ptr<RangeType> range() const;
     };
 
     

src/multi_array.cc

@@ -544,7 +544,7 @@ namespace MultiArrayTools
     template <typename T, class CRange, class Function>
     FunctionalMultiArray<T,CRange,Function>::FunctionalMultiArray(const std::shared_ptr<CRange>& range,
 								  const Function& func) :
-	MultiArrayBase<T,CRange>(range), mFunc(func) { /*CHECK;*/ }
+	MultiArrayBase<T,CRange>(range), mFunc(func) {}
 
     template <typename T, class CRange, class Function>
     const T& FunctionalMultiArray<T,CRange,Function>::operator[](const typename CRange::IndexType& i) const

src/multi_array_operation.cc

@@ -39,14 +39,13 @@ namespace MultiArrayTools
 	mSecond(second), mArrayRef(ma), mIndex()
     {
 	MultiRangeFactory<Ranges...> mrf( index->range() );
-	mIndex = mrf.create();	
+	std::shared_ptr<MultiRange<Ranges...> > mr =
-	CHECK;
+	    std::dynamic_pointer_cast<MultiRange<Ranges...> >( mrf.create() );
+	mIndex = std::make_shared<IndexType>( mr->begin() );
 	(*mIndex) = *index;
-	CHECK;
 	for(*mIndex = 0; mIndex->pos() != mIndex->max(); ++(*mIndex)){
 	    get() = mSecond.get();
 	}
-	CHECK;
     }
     
     template <typename T, class... Ranges>
@@ -72,15 +71,13 @@ namespace MultiArrayTools
 	OperationBase<T>(), OperationTemplate<ConstOperationRoot<T,Ranges...> >(this),
 	mArrayRef(ma), mIndex( std::make_shared<IndexType>( mArrayRef.range() ) )
     {
-	CHECK;
 	mIndex(indices...);
-	CHECK;
     }
 
     template <typename T, class... Ranges>
     const T& ConstOperationRoot<T,Ranges...>::get() const
     {
-	return mArrayRef[ *mIndex ];
+	return mArrayRef[ (*mIndex)() ];
     }
     
     /***********************
@@ -94,9 +91,7 @@ namespace MultiArrayTools
 	MutableOperationBase<T>(), OperationTemplate<OperationRoot<T,Ranges...> >(this),
 	mArrayRef(ma), mIndex( std::make_shared<IndexType>( mArrayRef.range() ) )
     {
-	CHECK;
 	(*mIndex)(indices...);
-	CHECK;
     }
 
     template <typename T, class... Ranges>
@@ -108,13 +103,13 @@ namespace MultiArrayTools
     template <typename T, class... Ranges>
     const T& OperationRoot<T,Ranges...>::get() const
     {
-	return mArrayRef[ *mIndex ];
+	return mArrayRef[ (*mIndex)() ];
     }
 
     template <typename T, class... Ranges>
     T& OperationRoot<T,Ranges...>::get()
     {
-	return mArrayRef[ *mIndex ];
+	return mArrayRef[ (*mIndex)() ];
     }
 
     /***********************

src/multi_range.cc

@@ -14,6 +14,7 @@ namespace MultiArrayTools
      *   MultiIndex   *
      ******************/
 
+    /*
     template <class... Indices>
     MultiIndex<Indices...>::MultiIndex(const MultiIndex<Indices...>& in) :
 	IndexInterface<std::tuple<typename Indices::MetaType...> >(in)
@@ -30,15 +31,12 @@ namespace MultiArrayTools
 	IB::mPos = PackNum<sizeof...(Indices)-1>::makePos(mIPack);
 	return *this;
     }
-
+    */
     template <class... Indices>
     MultiIndex<Indices...>& MultiIndex<Indices...>::operator=(ContainerIndex<Indices...>& ci)
     {	
-	CHECK;
 	PackNum<sizeof...(Indices)-1>::copyInst(mIPack, ci);
-	CHECK;
 	IB::mPos = PackNum<sizeof...(Indices)-1>::makePos(mIPack);
-	CHECK;
 	return *this;
     }
     
@@ -155,7 +153,7 @@ namespace MultiArrayTools
     }
 
     template <class... Indices>
-    const std::shared_ptr<typename MultiIndex<Indices...>::RangeType>& MultiIndex<Indices...>::range() const
+    std::shared_ptr<typename MultiIndex<Indices...>::RangeType> MultiIndex<Indices...>::range() const
     {
 	return std::dynamic_pointer_cast<RangeType>( IB::mRangePtr );
     }

src/multi_range.h

@@ -32,8 +32,8 @@ namespace MultiArrayTools
 	MultiIndex() = delete;
 	// NO DEFAULT HERE !!!
 	// ( have to assign sub-indices (ptr!) correctly )
-	MultiIndex(const MultiIndex& in);
+	//MultiIndex(const MultiIndex& in);
-	MultiIndex& operator=(const MultiIndex& in);
+	//MultiIndex& operator=(const MultiIndex& in);
 	MultiIndex& operator=(ContainerIndex<Indices...>& ci);
 
 	template <class MRange>
@@ -65,7 +65,7 @@ namespace MultiArrayTools
 	
 	virtual size_t dim() const override;
 
-	const std::shared_ptr<RangeType>& range() const;
+	std::shared_ptr<RangeType> range() const;
     };
 
     /*************************

src/op_unit_test.cc

@@ -66,8 +66,12 @@ namespace {
 	MultiArray<double,CRange> res(crptr);
 
 	auto i = std::dynamic_pointer_cast<SRange::IndexType>( srptr->index() );
-	CHECK;
 	res.operator()<SRange>(i) = ma1.operator()<SRange>(i) + ma2.operator()<SRange>(i);
+
+	EXPECT_EQ( fabs( res.at('a') - (2.917+8.870) ) < 0.0001, true);
+	EXPECT_EQ( fabs( res.at('l') - (9.436+4.790) ) < 0.0001, true );
+	EXPECT_EQ( fabs( res.at('f') - (0.373+8.215) ) < 0.0001, true );
+	EXPECT_EQ( fabs( res.at('g') - (7.192+5.063) ) < 0.0001, true );
     }
 
 } // anonymous namspace

src/single_range.h

@@ -23,7 +23,7 @@ namespace MultiArrayTools
 	typedef U MetaType;
 	typedef SingleRange<U,TYPE> RangeType;
 	
-	DEFAULT_MEMBERS_X(SingleIndex);
+	//DEFAULT_MEMBERS_X(SingleIndex);
 	
 	SingleIndex(const std::shared_ptr<SingleRange<U,TYPE> >& range);