WIP: slicing + partial ranges

src/include/array/array_base.cc.h

@@ -64,12 +64,12 @@ namespace CNORXZ
 
     template <typename T>
     template <typename I, typename M>
-    Sptr<CArrayBase<T>> CArrayBase<T>::sl(const IndexInterface<I,M>& i) const
+    Sptr<CArrayBase<T>> CArrayBase<T>::sl(const IndexInterface<I,M>& begin,
+					  const IndexInterface<I,M>& end) const
     {
-	auto beg = std::make_shared<const_iterator>(this->begin());
+	auto ai = itLexSave(begin);
-	auto si = i.slice(beg);
+	auto aj = itLexSave(end);
-	auto it = *beg + i.lex();
+	return std::make_shared<CSlice>(subcube(ai,aj), this, ai.format(), ai.pos());
-	return std::make_shared<CSlice>(this, *si, it.pos());
     }
 
     template <typename T>

src/include/array/array_base.h

@@ -40,7 +40,8 @@ namespace CNORXZ
 	const T& at(const DPack& pack) const;
 
 	template <typename I, typename M>
-	Sptr<CArrayBase<T>> sl(const IndexInterface<I,M>& i) const;
+	Sptr<CArrayBase<T>> sl(const IndexInterface<I,M>& begin,
+			       const IndexInterface<I,M>& end) const;
 
 	template <class Index>
 	COpRoot<T,Index> operator()(const Sptr<Index>& i) const;

src/include/ranges/prange.cc.h

@@ -326,6 +326,23 @@ namespace CNORXZ
     {
 	return Vector<Uuid> { mRange->id() };
     }
+
+
+    /****************************
+     *   non-member functions   *
+     ****************************/
+    
+    template <class I, typename M>
+    RangePtr prange(const IndexInterface<I,M>& begin, const IndexInterface<I,M>& end)
+    {
+	Vector<SizeT> parts(end-begin);
+	const SizeT off = begin.pos();
+	for(auto i = begin.THIS(); i != end.THIS(); ++i){
+	    parts[i.pos()-off] = i.pos();
+	}
+	return begin.range()->partial(parts); // implement!!!!
+	//return PRangeFactory<typename I::RangeType>(begin.range(), parts).create();
+    }
 }
 
 #endif

src/include/ranges/prange.h

@@ -118,6 +118,9 @@ namespace CNORXZ
 	Vector<SizeT> mParts;
     };
 
+    template <class I, typename M>
+    RangePtr prange(const IndexInterface<I,M>& begin, const IndexInterface<I,M>& end);
+    
 } // namespace CNORXZ
 
 #endif

src/include/ranges/range_base.h

@@ -124,6 +124,7 @@ namespace CNORXZ
 
     RangePtr getSub(const RangePtr& r, SizeT num);
     MArray<RangePtr> getSub(const RangePtr& r);
+    
 }
 
 #endif

src/include/ranges/yrange.h

@@ -9,6 +9,7 @@
 #include "xpr/xpr.h"
 #include "index_format.h"
 #include "index_pack.h"
+#include "prange.h"
 
 namespace CNORXZ
 {
@@ -159,6 +160,10 @@ namespace CNORXZ
 	static Sptr<YRange> func(const RangePtr& r);
     };
 
+    template <>
+    RangePtr prange(const IndexInterface<YIndex,Vector<DType>>& begin,
+		    const IndexInterface<YIndex,Vector<DType>>& end);
+    
 }
 
 #endif

src/lib/ranges/yrange.cc

@@ -560,4 +560,18 @@ namespace CNORXZ
     {
 	return std::dynamic_pointer_cast<YRange>( YRangeFactory({r}).create() );
     }
+
+    template <>
+    RangePtr prange(const IndexInterface<YIndex,Vector<DType>>& begin,
+		    const IndexInterface<YIndex,Vector<DType>>& end)
+    {
+	CXZ_WARNING("YRange specialization"); // test
+	const SizeT dim = begin.range()->dim();
+	Vector<RangePtr> v(dim);
+	for(SizeT i = 0; i != dim; ++i){
+	    v[i] = prange( DIndex(begin.THIS().pack()[i]), DIndex(end.THIS().pack()[i]) );
+	}
+	return YRangeFactory(v).create();
+    }
+
 }