bernstein/src/polynomial.cpp

   1 /*
   2  *      Bernstein Expansion
   3  *
   4  *      - polynomial class
   5  */
   6
   7 #include "polynomial.h"
   8
   9 using namespace GiNaC;
  10
  11 polynomial::polynomial(void)
  12 {
  13 }
  14
  15 polynomial::polynomial(ex &e)
  16 {
  17         *this += e;
  18 }
  19
  20 polynomial::polynomial(const ex &e)
  21 {
  22         *this += e;
  23 }
  24
  25
  26 polynomial& polynomial::operator=(ex &e)
  27 {
  28         *this = polynomial::polynomial(e);
  29
  30         return *this;
  31 }
  32
  33
  34 // term i of the polynomial
  35 ex polynomial::term(unsigned int i)
  36 {
  37         ex e = *this;
  38         if(!is_a<add>(e)) {
  39                 return e;
  40         }
  41         if(this->nops() == 0) {
  42                 return e;
  43         }
  44         return this->op(i);
  45 }
  46
  47
  48 // number of terms in the polynomial
  49 unsigned int polynomial::nbTerms(void)
  50 {
  51         ex e = *this;
  52         if(!is_a<add>(e)) {
  53                 return 1;
  54         }
  55
  56         if(this->nops() == 0) {
  57                 return 1;
  58         }
  59
  60         return this->nops();
  61 }