itensor is broken, partly
- Subject: itensor is broken, partly
- From: Robert Dodier
- Date: Wed, 24 May 2006 08:24:16 -0600
For the record, here are the problems in share/tensor/test/itest.mac
which yield unexpected results:
********************** Problem 26 ***************
Input:
contract(expand(lc2kdt(levi_civita([i, j, k], []) a([], [j])
levi_civita([], [k, l, n]) b([l], []) c([n], []))))
Result:
a([], [n]) b([1], []) c([n], [])
- a([], [l]) c([1], []) b([l], [])
This differed from the expected result:
a([], [n]) b([i], []) c([n], [])
- a([], [l]) c([i], []) b([l], [])
********************** Problem 27 ***************
Input:
contract(expand(lc2kdt(levi_civita([i, j, k], []) a([], [j])
b([], [k]) levi_civita([], [i, l, n]) c([l], []) u([n], []))))
Result:
3 a([], [l]) b([], [n]) c([l], []) u([n], [])
- 3 b([], [l]) a([], [n]) c([l], []) u([n], [])
- a([], [1]) b([], [n]) c([1], []) u([n], [])
+ b([], [1]) a([], [n]) c([1], []) u([n], [])
+ a([], [1]) b([], [l]) u([1], []) c([l], [])
- b([], [1]) a([], [l]) u([1], []) c([l], [])
This differed from the expected result:
a([], [l]) b([], [n]) c([l], []) u([n], [])
- b([], [l]) a([], [n]) c([l], []) u([n], [])
********************** Problem 28 ***************
Input:
canform(contract(expand(applyb1(levi_civita([], [i, j, k]) a([j])
a([k]), lc_l, lc_u))))
Result:
levi_civita([], [%1, %2, i]) a([%1], []) a([%2], [])
This differed from the expected result:
0
load (itensor); batch ("./itest.mac", test); executes the test.
FWIW
Robert