[GiNaC-devel] release 1.3.6
Chris Dams
Chris.Dams at mi.infn.it
Wed Dec 13 13:13:54 CET 2006
Dear Alexei,
On Wed, 13 Dec 2006, Sheplyakov Alexei wrote:
> ex H1 = color_f(a, c, e)*color_f(b, d, e)
> + color_f(a, d, e)*color_f(b, c, e);
>
> ex test = (H1*H1).expand(expand_options::expand_indexed);
> cout << test << endl;
> // 2*f.e.b.d*f.e.b.c*f.e.a.c*f.e.a.d+f.symbol8.b.c*f.symbol8.a.d*f.e.b.c*f.e.a.d+f.symbol7.b.d*f.symbol7.a.c*f.e.b.d*f.e.a.c
> // inconsistent indices!
>
> cout << test.simplify_indexed() << endl;
> // prints 288, should be 216
At my place this prints 144-2*f.c.a.symbol8*f.a.d.e*f.d.symbol8.b*f.c.e.b.
Strange???
@dummy index renaming: yes, this has been much improved in HEAD. In
principle I think it would be possible to move the class
make_flat_inserter into the *.cpp files. But that would be code
duplication. Not so nice. Besides that, new functions were added to
indexed.cpp and these need to be called from indexed.h. Fixing all this in
a binary compatible way would get way too messy IMO. Maybe we should
release 1.4 and claim proudly that we now can rename dummy indices. Good
thing that we don't keep an accurate NEWS file. That way we can can
announce it like it is new. :-o ;-)
Best wishes,
Chris
More information about the GiNaC-devel
mailing list