| 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 |