Всегда думал что тензорные расчёты для maple не составят проблем. Но вдруг оказалось, что мэмпл даёт бредовый результат. Очень хотелось бы узнать что это всего лишь моя ошибка, а не косячность пакета physics.
Давайте рассмотрим пример:
Код:
Simplify(LeviCivita[4, sigma, lambda, rho]*LeviCivita[4, xi, eta, mu]*g_[rho, mu]*qp[sigma]*q[lambda]*qp[xi]*q[eta])
Ответ мэмл даёт ответ таким:
Но это неверно, так как индексы должны бежать 1,2,3 а здесь мэпл считает что и 4 индекс пробегается.
Для подтверждения этих слов можете сделать SumOverRepeatedIndices(%).
Моя приамбула такова:
Код:
with(Physics)
Setup(mathematicalnotation = true);
Coordinates(X);
Setup(spaceindices = lowercaselatin)
Setup(tensors = q[mu](X))
PDEtools:-declare(q(X))
Setup(tensors = qp[mu](X))
PDEtools:-declare(qp(X))