Sympy избегает неправильного вывода при упрощении суммы дельты Кронекера с помощью Piecewise ⇐ Python

[Anonymous]

Код: Выделить всё

from sympy.abc import x, i, j, n
from sympy import *

x = IndexedBase('x')
U = -Sum(Piecewise((KroneckerDelta(i, j, (1, n)), x[i] > 0), (0, True)), (i, 1, n))
print(simplify(U))
# Output: Piecewise((-1, x[i] > 0), (0, True))
# Which is wrong because i should not be a free variable here
# because it wasn't a free variable in the input.
# answer should be Piecewise((-1, x[j] > 0), (0, True))

Это также работает, только если U начинается с «-Sum». Если я опущу отрицательный знак, это не упростится автоматически, но я могу упростить его с помощью Piecewise_fold и doit:

Код: Выделить всё

U1 = Sum(Piecewise((KroneckerDelta(i, j, (1, n)), x[i] > 0), (0, True)), (i, 1, n))
print(simplify(piecewise_fold(U1)).doit())
# Output: Piecewise((1, x[i] > 0), (0, True))
# (wrong for the same reason as before)

Так это просто ошибка в SymPy? Есть ли обходной путь, позволяющий избежать неправильного отображения свободной переменной i в выходных данных?
Связано:

• Sympy упростить сумма дельты Кронекера с помощью Piecewise
• Sympy упрощает сумму дельты Кронекера

Подробнее здесь: https://stackoverflow.com/questions/792 ... ith-piecew