[GiNaC-list] bugs in indexed object
Dr. Mithun Bairagi
bairagirasulpur at gmail.com
Sun Apr 14 10:01:15 CEST 2024
Dear all ginac lovers, I have found simplify_indexed() is not working as expected.
I have compiled and run the following c++ program using ginac 1.87 (mingw64 compiler, windows 10 os),
#include <ginac/ginac.h>
using namespace GiNaC;
using namespace std;
int main()
{
symbol A("A"), B("B"), C("C"), i_sym("i"),j_sym("j");
varidx i(i_sym, 3);
varidx ic(i_sym, 3,true);
varidx j(j_sym, 3);
varidx jc(j_sym, 3,true);
ex e = indexed(A, {i,i}) * indexed(B, {ic,ic})-indexed(A, {j,j}) * indexed(B, {jc,jc});
cout << (e).simplify_indexed() << endl;
return 0;
}
the above program shows the following result:
A.i.i*B~symbol16.symbol16-A.i.i*B~symbol17.symbol17
But I think the result should be zero. Please tell me what is the problem here?
Best regards,
Mithun Bairagi
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://www.ginac.de/pipermail/ginac-list/attachments/20240414/7f73994c/attachment.htm>
More information about the GiNaC-list
mailing list