| 53 | 
  | 
#include <list> | 
| 54 | 
  | 
#include <map> | 
| 55 | 
  | 
#include <utility> | 
| 56 | 
< | 
 | 
| 56 | 
> | 
#include "config.h" | 
| 57 | 
  | 
namespace oopse { | 
| 58 | 
  | 
 | 
| 59 | 
  | 
  template<typename ElemType> ElemType pow(ElemType x, int N) { | 
| 189 | 
  | 
      return polyPairMap_.size(); | 
| 190 | 
  | 
    } | 
| 191 | 
  | 
 | 
| 192 | 
+ | 
    PolynomialType& operator = (const PolynomialType& p) { | 
| 193 | 
+ | 
 | 
| 194 | 
+ | 
      if (this != &p)  // protect against invalid self-assignment | 
| 195 | 
+ | 
      { | 
| 196 | 
+ | 
        typename Polynomial<ElemType>::const_iterator i; | 
| 197 | 
+ | 
 | 
| 198 | 
+ | 
        polyPairMap_.clear();  // clear out the old map | 
| 199 | 
+ | 
       | 
| 200 | 
+ | 
        for (i =  p.begin(); i != p.end(); ++i) { | 
| 201 | 
+ | 
          this->setCoefficient(i->first, i->second); | 
| 202 | 
+ | 
        } | 
| 203 | 
+ | 
      } | 
| 204 | 
+ | 
      // by convention, always return *this | 
| 205 | 
+ | 
      return *this;  | 
| 206 | 
+ | 
    } | 
| 207 | 
+ | 
 | 
| 208 | 
  | 
    PolynomialType& operator += (const PolynomialType& p) { | 
| 209 | 
  | 
        typename Polynomial<ElemType>::const_iterator i; | 
| 210 | 
  | 
 | 
| 220 | 
  | 
        for (i =  p.begin(); i  != p.end(); ++i) { | 
| 221 | 
  | 
          this->addCoefficient(i->first, -i->second); | 
| 222 | 
  | 
        }         | 
| 223 | 
+ | 
        return *this; | 
| 224 | 
  | 
    } | 
| 225 | 
  | 
 | 
| 226 | 
  | 
    PolynomialType& operator *= (const PolynomialType& p) { | 
| 327 | 
  | 
 | 
| 328 | 
  | 
  /** | 
| 329 | 
  | 
   * Tests if two polynomial have the same exponents | 
| 330 | 
< | 
   * @return true if these all of the exponents in these Polynomial are identical | 
| 330 | 
> | 
   * @return true if all of the exponents in these Polynomial are identical | 
| 331 | 
  | 
   * @param p1 the first polynomial | 
| 332 | 
  | 
   * @param p2 the second polynomial | 
| 333 | 
  | 
   * @note this function does not compare the coefficient | 
| 351 | 
  | 
    return true; | 
| 352 | 
  | 
  } | 
| 353 | 
  | 
 | 
| 354 | 
< | 
  typedef Polynomial<double> DoublePolynomial; | 
| 354 | 
> | 
  typedef Polynomial<RealType> DoublePolynomial; | 
| 355 | 
  | 
 | 
| 356 | 
  | 
} //end namespace oopse | 
| 357 | 
  | 
#endif //MATH_POLYNOMIAL_HPP |