[GiNaC-devel] infinite loop in simplify_indexed()

Jens Vollinga jensv at nikhef.nl
Sat May 21 01:43:40 CEST 2011


Hi,

Am 18.05.2011 00:33, schrieb PG CLARK:
> But what do you do with it next? add all the bits up and try to cancel
> things out treating the indexes as static labels? That would be a sure
> fire way to get a result. But is that all it's used for? If so there'd
> have been no need for it to be called for the test program where there
> was no symmetry to consider. Or is it also used to account for pseudo
> symmetry that may exist due to the equivalence of dummy indexes? So if
> dummy indexes in a term have the same dimension you treat them as a
> symmetric permutation and sum over all of them too?
>
> There has to be a better way!

I think I fixed the problem (and hopefully didn't mess up shortly before 
the next release ...).

simplify_indexed() now doesn't do a complete symmetrization anymore if 
no non-symmetric indices are involved, since a vanishing expression 
cannot be expected then.

Regards,
Jens


More information about the GiNaC-devel mailing list