From: Christian Bauer <Christian.Bauer@uni-mainz.de>
Date: Wed, 28 Jan 2004 19:34:55 +0000 (+0000)
Subject: expand_indexed expands the base expression
X-Git-Tag: release_1-2-0~14
X-Git-Url: https://www.ginac.de/ginac.git//ginac.git?p=ginac.git;a=commitdiff_plain;h=2053dbc953da03b57dfd7da01fcdc196e8b8e912

expand_indexed expands the base expression
---

diff --git a/ginac/indexed.cpp b/ginac/indexed.cpp
index 0bb34d8a..40a6232c 100644
--- a/ginac/indexed.cpp
+++ b/ginac/indexed.cpp
@@ -325,20 +325,24 @@ ex indexed::expand(unsigned options) const
 {
 	GINAC_ASSERT(seq.size() > 0);
 
-	if ((options & expand_options::expand_indexed) && is_exactly_a<add>(seq[0])) {
-
-		// expand_indexed expands (a+b).i -> a.i + b.i
-		const ex & base = seq[0];
-		ex sum = _ex0;
-		for (size_t i=0; i<base.nops(); i++) {
+	if (options & expand_options::expand_indexed) {
+		ex newbase = seq[0].expand(options);
+		if (is_exactly_a<add>(newbase)) {
+			ex sum = _ex0;
+			for (size_t i=0; i<newbase.nops(); i++) {
+				exvector s = seq;
+				s[0] = newbase.op(i);
+				sum += thiscontainer(s).expand(options);
+			}
+			return sum;
+		}
+		if (!are_ex_trivially_equal(newbase, seq[0])) {
 			exvector s = seq;
-			s[0] = base.op(i);
-			sum += thiscontainer(s).expand();
+			s[0] = newbase;
+			return ex_to<indexed>(thiscontainer(s)).inherited::expand(options);
 		}
-		return sum;
-
-	} else
-		return inherited::expand(options);
+	}
+	return inherited::expand(options);
 }
 
 //////////