- fixed bug in multiply_lcm() (see paranoia_check9)

[ginac.git] / ginac / power.h

diff --git a/ginac/power.h b/ginac/power.h
index a946a82f5311be8d368fced4227dd4e19fc75d2e..13dfd5ffbd4da8b413956991a74a1819f8695c9e 100644 (file)
--- a/ginac/power.h
+++ b/ginac/power.h
@@ -1,8 +1,9 @@
 /** @file power.h
  *
- *  Interface to GiNaC's symbolic exponentiation (basis^exponent).
- *
- *  GiNaC Copyright (C) 1999 Johannes Gutenberg University Mainz, Germany
+ *  Interface to GiNaC's symbolic exponentiation (basis^exponent). */
+
+/*
+ *  GiNaC Copyright (C) 1999-2000 Johannes Gutenberg University Mainz, Germany
  *
  *  This program is free software; you can redistribute it and/or modify
  *  it under the terms of the GNU General Public License as published by
@@ -22,7 +23,13 @@
 #ifndef __GINAC_POWER_H__
 #define __GINAC_POWER_H__
 
-class power;
+#include "basic.h"
+#include "ex.h"
+
+#ifndef NO_NAMESPACE_GINAC
+namespace GiNaC {
+#endif // ndef NO_NAMESPACE_GINAC
+
 class numeric;
 class add;
 
@@ -30,6 +37,8 @@ class add;
  *  representing exponentiation. */
 class power : public basic
 {
+    GINAC_DECLARE_REGISTERED_CLASS(power, basic)
+
     friend class mul;
 
 // member functions
@@ -38,39 +47,39 @@ class power : public basic
 public:
     power();
     ~power();
-    power(power const & other);
-    power const & operator=(power const & other);
+    power(const power & other);
+    const power & operator=(const power & other);
 protected:
-    void copy(power const & other);
+    void copy(const power & other);
     void destroy(bool call_parent);
 
     // other constructors
 public:
-    power(ex const & lh, ex const & rh);
-    power(ex const & lh, numeric const & rh);
+    power(const ex & lh, const ex & rh);
+    power(const ex & lh, const numeric & rh);
 
     // functions overriding virtual functions from bases classes
 public:
     basic * duplicate() const;
+    void print(ostream & os, unsigned upper_precedence=0) const;
     void printraw(ostream & os) const;
     void printtree(ostream & os, unsigned indent) const;
-    void print(ostream & os, unsigned upper_precedence=0) const;
     void printcsrc(ostream & os, unsigned type, unsigned upper_precedence=0) const;
     bool info(unsigned inf) const;
-    int nops() const;
-    ex & let_op(int const i);
-    int degree(symbol const & s) const;
-    int ldegree(symbol const & s) const;
-    ex coeff(symbol const & s, int const n=1) const;
+    unsigned nops() const;
+    ex & let_op(int i);
+    int degree(const symbol & s) const;
+    int ldegree(const symbol & s) const;
+    ex coeff(const symbol & s, int n=1) const;
     ex eval(int level=0) const;
     ex evalf(int level=0) const;
-    ex diff(symbol const & s) const;
-    ex series(symbol const & s, ex const & point, int order) const;
-    ex subs(lst const & ls, lst const & lr) const;
+    ex diff(const symbol & s) const;
+    ex series(const symbol & s, const ex & point, int order) const;
+    ex subs(const lst & ls, const lst & lr) const;
     ex normal(lst &sym_lst, lst &repl_lst, int level=0) const;
-    ex simplify_ncmul(exvector const & v) const;
+    ex simplify_ncmul(const exvector & v) const;
 protected:
-    int compare_same_type(basic const & other) const;
+    int compare_same_type(const basic & other) const;
     unsigned return_type(void) const;
     unsigned return_type_tinfo(void) const;
     ex expand(unsigned options=0) const;
@@ -80,12 +89,12 @@ protected:
     
     // non-virtual functions in this class
 protected:
-    ex expand_add(add const & a, int const n) const;
-    ex expand_add_2(add const & a) const;
-    ex expand_mul(mul const & m, numeric const & n) const;
-    //ex expand_commutative_3(ex const & basis, numeric const & exponent,
+    ex expand_add(const add & a, int n) const;
+    ex expand_add_2(const add & a) const;
+    ex expand_mul(const mul & m, const numeric & n) const;
+    //ex expand_commutative_3(const ex & basis, const numeric & exponent,
     //                         unsigned options) const;
-    // ex expand_noncommutative(ex const & basis, numeric const & exponent, unsigned options) const;
+    // ex expand_noncommutative(const ex & basis, const numeric & exponent, unsigned options) const;
 
 // member variables
 
@@ -98,9 +107,13 @@ protected:
 // global constants
 
 extern const power some_power;
-extern type_info const & typeid_power;
+extern const type_info & typeid_power;
 
-#define ex_to_power(X) static_cast<power const &>(*(X).bp)
+// utility functions
+inline const power &ex_to_power(const ex &e)
+{
+       return static_cast<const power &>(*e.bp);
+}
 
 // wrapper functions
 
@@ -108,13 +121,15 @@ extern type_info const & typeid_power;
  *
  *  @param b the basis expression
  *  @param e the exponent expression */
-inline ex pow(ex const & b, ex const & e)
+inline ex pow(const ex & b, const ex & e)
 { return power(b,e); }
 
 /** Square root expression.  Returns a power-object with exponent 1/2 as a new
  *  expression.  */
-inline ex sqrt(ex const & a)
-{ return power(a,exHALF()); }
+ex sqrt(const ex & a);
 
-#endif // ndef __GINAC_POWER_H__
+#ifndef NO_NAMESPACE_GINAC
+} // namespace GiNaC
+#endif // ndef NO_NAMESPACE_GINAC
 
+#endif // ndef __GINAC_POWER_H__