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 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

Best regards,

Mithun Bairagi