add operator- Index

src/include/ranges/crange.h

@@ -26,6 +26,7 @@ namespace CNORXZ
 	CIndex& operator--();
 	CIndex operator+(Int n) const;
 	CIndex operator-(Int n) const;
+	SizeT operator-(const CIndex& i) const;
 	CIndex& operator+=(Int n);
 	CIndex& operator-=(Int n);
 

src/include/ranges/dindex.h

@@ -32,6 +32,7 @@ namespace CNORXZ
 	DIndex& operator--();
 	DIndex operator+(Int n) const;
 	DIndex operator-(Int n) const;
+	SizeT operator-(const DIndex& i) const;
 	DIndex& operator+=(Int n);
 	DIndex& operator-=(Int n);
 

src/include/ranges/mrange.cc.h

@@ -299,6 +299,12 @@ namespace CNORXZ
 	return o -= n;
     }
 
+    template <class FormatT, class... Indices>
+    SizeT GMIndex<FormatT,Indices...>::operator-(const GMIndex& i) const
+    {
+	return lex() - i.lex();
+    }
+
     template <class FormatT, class... Indices>
     GMIndex<FormatT,Indices...>& GMIndex<FormatT,Indices...>::operator+=(Int n)
     {

src/include/ranges/mrange.h

@@ -47,6 +47,7 @@ namespace CNORXZ
 	GMIndex& operator--();
 	GMIndex operator+(Int n) const;
 	GMIndex operator-(Int n) const;
+	SizeT operator-(const GMIndex& i) const;
 	GMIndex& operator+=(Int n);
 	GMIndex& operator-=(Int n);
 

src/include/ranges/prange.cc.h

@@ -52,6 +52,12 @@ namespace CNORXZ
 	return PIndex(mRangePtr, IB::mPos + n);
     }
 
+    template <class IndexT>
+    SizeT PIndex<IndexT>::operator-(const PIndex& i) const
+    {
+	return lex() - i.lex();
+    }
+    
     template <class IndexT>
     PIndex<IndexT> PIndex<IndexT>::operator-(Int n) const
     {

src/include/ranges/prange.h

@@ -26,6 +26,7 @@ namespace CNORXZ
 	PIndex& operator--();
 	PIndex operator+(Int n) const;
 	PIndex operator-(Int n) const;
+	SizeT operator-(const PIndex& i) const;
 	PIndex& operator+=(Int n);
 	PIndex& operator-=(Int n);
 

src/include/ranges/srange.cc.h

@@ -54,6 +54,12 @@ namespace CNORXZ
 	return SIndex(mRangePtr, IB::mPos - n);
     }
 
+    template <typename MetaT, SizeT S>
+    SizeT SIndex<MetaT,S>::operator-(const SIndex& i) const
+    {
+	return lex() - i.lex();
+    }
+
     template <typename MetaT, SizeT S>
     SIndex<MetaT,S>& SIndex<MetaT,S>::operator+=(Int n)
     {

src/include/ranges/srange.h

@@ -27,6 +27,7 @@ namespace CNORXZ
 	SIndex& operator--();
 	SIndex operator+(Int n) const;
 	SIndex operator-(Int n) const;
+	SizeT operator-(const SIndex& i) const;
 	SIndex& operator+=(Int n);
 	SIndex& operator-=(Int n);
 

src/include/ranges/urange.cc.h

@@ -58,6 +58,12 @@ namespace CNORXZ
 	return UIndex(mRangePtr, IB::mPos - n);
     }
 
+    template <typename MetaT>
+    SizeT UIndex<MetaT>::operator-(const UIndex& i) const
+    {
+	return lex() - i.lex();
+    }
+    
     template <typename MetaT>
     UIndex<MetaT>& UIndex<MetaT>::operator+=(Int n)
     {

src/include/ranges/urange.h

@@ -26,14 +26,12 @@ namespace CNORXZ
 	DEFAULT_MEMBERS(UIndex);
 	UIndex(const RangePtr& range, SizeT pos = 0);
 
-	void swap(UIndex& i) {}; // !!!
-	
 	UIndex& operator=(SizeT lexpos);
 	UIndex& operator++();
 	UIndex& operator--();
 	UIndex operator+(Int n) const;
 	UIndex operator-(Int n) const;
-	SizeT operator-(const UIndex& i) const { return lex() - i.lex(); } // !!!
+	SizeT operator-(const UIndex& i) const;
 	UIndex& operator+=(Int n);
 	UIndex& operator-=(Int n);
 

src/include/ranges/xindex.cc.h

@@ -65,6 +65,14 @@ namespace CNORXZ
 	return std::make_shared<XIndex<Index,Meta>>(*mI - n);
     }
     
+    template <class Index, typename Meta>
+    SizeT XIndex<Index,Meta>::operator-(const XIndexBase& i) const
+    {
+	const XIndex<Index,Meta>* ip = dynamic_cast<const XIndex<Index,Meta>*>(&i);
+	CXZ_ASSERT(ip != nullptr, "bad index type");
+	return (*mI) - (*ip->mI);
+    }
+    
     template <class Index, typename Meta>
     XIndex<Index,Meta>& XIndex<Index,Meta>::operator+=(Int n)
     {

src/include/ranges/xindex.h

@@ -25,6 +25,7 @@ namespace CNORXZ
 	virtual XIndexBase& operator--() = 0;
 	virtual XIndexPtr operator+(Int n) const = 0;
 	virtual XIndexPtr operator-(Int n) const = 0;
+	virtual SizeT operator-(const XIndexBase& i) const = 0;
 	virtual XIndexBase& operator+=(Int n) = 0;
 	virtual XIndexBase& operator-=(Int n) = 0;
 
@@ -80,6 +81,7 @@ namespace CNORXZ
 	virtual XIndex& operator--() override final;
 	virtual XIndexPtr operator+(Int n) const override final;
 	virtual XIndexPtr operator-(Int n) const override final;
+	virtual SizeT operator-(const XIndexBase& i) const override final;
 	virtual XIndex& operator+=(Int n) override final;
 	virtual XIndex& operator-=(Int n) override final;
 

src/include/ranges/yrange.h

@@ -38,6 +38,7 @@ namespace CNORXZ
 	YIndex& operator--();
 	YIndex operator+(Int n) const; // equivalent to applying n times ++
 	YIndex operator-(Int n) const;
+	SizeT operator-(const YIndex& i) const;
 	YIndex& operator+=(Int n);
 	YIndex& operator-=(Int n);
 

src/lib/ranges/crange.cc

@@ -41,6 +41,11 @@ namespace CNORXZ
 	return CIndex(mRangePtr, IB::mPos - n);
     }
 
+    SizeT CIndex::operator-(const CIndex& i) const
+    {
+	return lex() - i.lex();
+    }
+    
     CIndex& CIndex::operator+=(Int n)
     {
 	IB::mPos += n;

src/lib/ranges/dindex.cc

@@ -75,6 +75,11 @@ namespace CNORXZ
 	return (*mI) - n;
     }
 
+    SizeT DIndex::operator-(const DIndex& i) const
+    {
+	return (*mI) - (*i.mI);
+    }
+    
     DIndex& DIndex::operator+=(Int n)
     {
 	(*mI) += n;

src/lib/ranges/yrange.cc

@@ -247,6 +247,11 @@ namespace CNORXZ
 	return o -= n;
     }
 
+    SizeT YIndex::operator-(const YIndex& i) const
+    {
+	return lex() - i.lex();
+    }
+    
     YIndex& YIndex::operator+=(Int n)
     {
 	if(static_cast<Int>(lex()) + n < 0){