[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