dynamic index multiplication

src/include/ranges/index_mul.cc.h

@@ -84,7 +84,29 @@ namespace CNORXZ
     constexpr decltype(auto) operator*(const IndexInterface<I1,Meta1>& a,
 				       const IndexInterface<I2,Meta2>& b)
     {
-	// special operations for DIndex / YIndex
+	if constexpr(std::is_same<I1,DIndex>::value){
+	    if constexpr(std::is_same<I2,DIndex>::value){
+		return YIndex({ a.THIS().xptr(), b.THIS().xptr() });
+	    }
+	    else if constexpr(std::is_same<I2,YIndex>::value){
+		auto p = b.THIS().pack();
+		p.insert(0, a.THIS().xptr());
+		return YIndex(p);
+	    }
+	}
+	else if constexpr(std::is_same<I1,YIndex>::value){
+	    if constexpr(std::is_same<I2,DIndex>::value){
+		auto p = a.THIS().pack();
+		p.push_back(b.THIS().xptr());
+		return YIndex(p);
+	    }
+	    else if constexpr(std::is_same<I2,YIndex>::value){
+		auto p = a.THIS().pack();
+		p.insert(p.end(), b.THIS().pack().begin(), b.THIS().pack().end());
+		return YIndex(p);
+	    }
+	}
+	else {
 	    constexpr SizeT I1D = index_dim<I1>::value;
 	    constexpr SizeT I2D = index_dim<I2>::value;
 	    if constexpr(I1D == 1){
@@ -106,6 +128,7 @@ namespace CNORXZ
 		}
 	    }
 	}
+    }
 
     /***************
      *   iptrMul   *
@@ -114,21 +137,47 @@ namespace CNORXZ
     template <class I1, class I2>
     decltype(auto) iptrMul(const Sptr<I1>& a, const Sptr<I2>& b)
     {
-	// special operations for DIndex / YIndex
+	if constexpr(std::is_same<I1,DIndex>::value){
-	if constexpr(index_dim<I1>::value == 1){
+	    if constexpr(std::is_same<I2,DIndex>::value){
+		return std::make_shared<YIndex>({ a->xptr(), b->xptr() });
+	    }
+	    else if constexpr(std::is_same<I2,YIndex>::value){
+		auto p = b->pack();
+		p.insert(0, a->xptr());
+		return std::make_shared<YIndex>(p);
+	    }
+	}
+	else if constexpr(std::is_same<I1,YIndex>::value){
+	    if constexpr(std::is_same<I2,DIndex>::value){
+		auto p = a->pack();
+		p.push_back(b->xptr());
+		return std::make_shared<YIndex>(p);
+	    }
+	    else if constexpr(std::is_same<I2,YIndex>::value){
+		auto p = a->pack();
+		p.insert(p.end(), b->pack().begin(), b->pack().end());
+		return std::make_shared<YIndex>(p);
+	    }
+	}
+	else {
+	    constexpr SizeT I1D = index_dim<I1>::value;
+	    constexpr SizeT I2D = index_dim<I2>::value;
+	    if constexpr(I1D == 1){
 		if constexpr(index_dim<I2>::value == 1){
 		    return std::make_shared<MIndex<I1,I2>>(a, b);
 		}
 		else {
-		return MIndexSptrMul::evalXM(a, b);
+		    return MIndexSptrMul::evalXM(a, b, std::make_index_sequence<I1D>{});
 		}
 	    }
 	    else {
 		if constexpr(index_dim<I2>::value == 1){
-		return MIndexSptrMul::evalMX(a, b);
+		    return MIndexSptrMul::evalMX(a, b, std::make_index_sequence<I2D>{});
 		}
 		else {
-		return MIndexSptrMul::evalMM(a, b);
+		    return MIndexSptrMul::evalMM(a, b, std::make_index_sequence<I1D>{},
+						 std::make_index_sequence<I2D>{});
+		}
 	    }
 	}
     }

src/tests/range_unit_test.cc

@@ -430,6 +430,10 @@ namespace
 	    for(auto ui = ur->begin(); ui != ur->end(); ++ui){
 		const SizeT p = ci.lex()*s1 + ui.lex();
 		EXPECT_EQ((ci*ui).lex(), p);
+		for(auto ci2 = cr->begin(); ci2 != cr->end(); ++ci2){
+		    const SizeT p2 = ci.lex()*s1*s2 + ui.lex()*s2 + ci2.lex();
+		    EXPECT_EQ((ci*ui*ci2).lex(), p2);
+		}
 	    }
 	}
     }