[GiNaC-list] subs does not work
Zhao Li
zhaoli at ihep.ac.cn
Thu Feb 21 04:42:14 CET 2013
Here I am trying to make a substitution. For testing, the function "has" works as shown below. But the "subs" just does not implement the substitution.
cout << has(temp, JF1(wild(1), GA(wild(5)), wild(2))*JF1(wild(3), GA(wild(5)), wild(4)) , has_options::algebraic ) << endl;
cout << temp.subs( JF1(wild(1),GA(wild(5)),wild(2))*JF1(wild(3),GA(wild(5)),wild(4)) == 2, subs_options::algebraic ) << endl;
The functions "JF1" and "GA" are both declared using DECLARE_FUNCTION and commutative.
DECLARE_FUNCTION_1P( GA );
REGISTER_FUNCTION( GA, dummy() );
Best regards,
Zhao
Theory Division
Institute of High Energy Physics
Chinese Academy of Sciences
More information about the GiNaC-list
mailing list