[GiNaC-devel] Simplifying powers

[Vladimir V. Kisil]

	Dear All,

	The method power::eval() contains the following rules:

// ^(*(...,x;c1),c2) -> *(^(*(...,x;1),c2),c1^c2)  (c1, c2 numeric(), c1>0)
// ^(*(...,x;c1),c2) -> *(^(*(...,x;-1),c2),(-c1)^c2)  (c1, c2 numeric(), c1<0)

  Is there a good reason why c1 here shall be numeric? Cannot we simply check the
  positive/negative info flag for c1 for that simplification? 
  Shall I prepare a patch for this?

  Currently we have:

  possymbol  n("n"), m("m");
  cout << "e=" << pow(pow(n,2),half) << endl;  
  // -> n
  cout << "e=" << pow(pow(n*m,2),half) << endl;  
  // -> sqrt(m^2*n^2)

  Best wishes,
  Vladimir
-- 
Vladimir V. Kisil     email: kisilv at maths.leeds.ac.uk
                        www: http://www.maths.leeds.ac.uk/~kisilv/
Book: Geometry of Mobius Transformations
http://www.worldscientific.com/worldscibooks/10.1142/p835