ianhinder / Kranc

A Mathematica package for generating code for solving time dependent partial differential equations
http://kranccode.org
GNU General Public License v2.0
28 stars 10 forks source link

TestTensorTools.nb, contractions in partial derivatives fail #87

Closed eschnett closed 11 years ago

eschnett commented 11 years ago

The test PD[u[ua],la] fails.

This is also a real-world problem; I have to work around this issue by contracting via KroneckerDelta outside of PD in McLachlan.

eschnett commented 11 years ago

This version of the test succeeds:

test[MakeExplicit[ KD[ua, lb] PD[u[ub], la]], {PD[u1, 1] + PD[u2, 2] + PD[u3, 3]}]

Note that I needed to modify both input and output. I believe the expected output for the test is wrong as well.

ianhinder commented 11 years ago

I believe that this was fixed in 4ac728ac483c7b1a488474f87ff097905466d7eb.. 321220cf5b4ef8de0fc1a57cb679f0d2f3ff599e. Please reopen if not.