21 Feb
2013
21 Feb
'13
4:42 a.m.
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