method 'subs()'

8 May 2002

      Hi all,
I don't understand the behavior of method 'subs()' in the following
case:
...
subs(a^x*b^x, $1^$2*$3^$2==($1*$3)^$2);
(a*b)^x
subs(a^5*b^5, $1^$2*$3^$2==($1*$3)^$2);
b^5*a^5
I also check it with a little program, that I attach, and the output is

Expression 1: b^x*a^x
List of match: {$1==b,$2==x,$3==a}
Expression 1 after subs: (b*a)^x

Expression 2: b^5*a^5
List of match: {$1==b,$2==5,$3==a}
Expression 2 after subs: b^5*a^5

Why is different the behavior of subs() when it is applied on
b^x*a^x or on b^5*a^5?

Best regards

     Tatiana Zolo

#include <iostream>
#include <ginac/ginac.h>

using namespace std;
using namespace GiNaC;

typedef GiNaC::ex GExpr;
typedef GiNaC::symbol GSymbol;
typedef GiNaC::numeric GNumber;
typedef GiNaC::lst GList;

int main() {
  GSymbol x("x");
  GSymbol a("a");
  GSymbol b("b");
  GList l_1, l_2;

  GExpr p = pow(a, x) * pow(b, x);
  cout << endl << "Expression 1: " << p << endl;
  match(p, pow(wild(1), wild(2)) * pow(wild(3), wild(2)), l_1);
  cout << "List of match: " << l_1 << endl;
  p = p.subs(pow(wild(1), wild(2)) * pow(wild(3), wild(2))
	     == pow(wild(1)*wild(3), wild(2)));
  cout << "Expression 1 after subs: " << p << endl << endl;

  GExpr q = pow(a, 5) * pow(b, 5);
  cout << "Expression 2: " << q << endl;
  match(q, pow(wild(1), wild(2)) * pow(wild(3), wild(2)), l_2);
  cout << "List of match: " << l_2 << endl;
  q = q.subs(pow(wild(1), wild(2)) * pow(wild(3), wild(2))
	     == pow(wild(1)*wild(3), wild(2)));
  cout << "Expression 2 after subs: " << q << endl << endl;

  return 0;
}