[GiNaC-list] subs does not work

[Richard B. Kreckel]

On 02/21/2013 04:42 AM, Zhao Li wrote:
> 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() );

$ cat zhaoli.cc
#include <iostream>
#include <ginac/ginac.h>
using namespace std;
using namespace GiNaC;

DECLARE_FUNCTION_1P( GA );
REGISTER_FUNCTION( GA, dummy() );

DECLARE_FUNCTION_3P( JF1 );
REGISTER_FUNCTION( JF1, dummy() );

int main()
{
     symbol a("a");
     symbol b("b");
     symbol t("t");
     symbol x("x");
     symbol y("y");
     ex temp = JF1(a, GA(t), b)*JF1(x, GA(t), y);
     cout << temp << endl;
     cout << temp.has( 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;
}
$ g++ zhaoli.cc -lginac
$./a.out
JF1(a,GA(t),b)*JF1(x,GA(t),y)
1
2

Cheers
   -richy.
-- 
Richard B. Kreckel