polyTg.h

#ifndef POLYTG_H
#define POLYTG_H
#include "Tg.h"
#include "../common/ISterm.h"

#include "../wrapperMMGD/polyWrapper.h"
#include "../wrapperMMGD/matrixWrapper.h"

namespace etvo
{
 class seriesTg;
class polyTg:public ISterm
{
public:
  static polyTg Epsilon();
  static polyTg Top();
  static polyTg E();

  polyTg();
  polyTg(bool TopNotE); //true Top false E
  polyTg(const Tg & m);
  polyTg(const std::vector<Tg> & v);

  void canon();
  bool isCanon() const;
  bool isE() const;
  polyTg operator+(const polyTg &p) const;
  polyTg oplus(const polyTg &p) const;
  polyTg oplusCD(const polyTg &p) const;
  polyTg operator+(const Tg & m)const;
  void add(const Tg & m);
  polyTg operator+=(const Tg & m);

  polyTg operator*(const polyTg & p) const;
  polyTg otimes(const polyTg &p) const;
  polyTg operator*(const Tg & m) const;
  static polyTg otimes(const Tg & m, const polyTg & p);
  polyTg otimesCD(const polyTg &p) const;
  
  polyTg inf(const polyTg & p) const;
  polyTg infCD(const polyTg & p) const;
  
  seriesTg star() const;
    
  polyTg lfrac(const polyTg & p) const;
  polyTg lfracCD(const polyTg & p) const;
  polyTg lfrac(const Tg & m) const;
  polyTg rfrac(const polyTg & p) const;
  polyTg rfracCD(const polyTg & p) const;
  polyTg rfrac(const Tg & m) const;
  bool operator==(const polyTg &)const;
  bool operator!=(const polyTg &)const;
  bool operator<=(const polyTg &) const;
  bool operator>=(const polyTg &) const;

  static polyTg toPolyTg(const poly & p);
  poly toPoly() const;
  static polyTg toCausal(const polyTg & p);
  Tg getFirstDif(const polyTg & p) const;  
  polyTg transientStar(int Tmax) const;
  void getMaxGain(unsigned int & vee, unsigned int & wedge) const;
  void getLcmGain(unsigned int & vee, unsigned int & wedge) const;
  std::pair<unsigned int, unsigned int> getPeriodicity() const;
  std::vector<Tg> getTerms() const;
  void removeTerm(unsigned idx);
  Tg operator[](unsigned idx) const;
  unsigned int size() const;
  std::string toString() const;
  std::string toStringAsDeltaVar() const;
   
  matrix<poly> getCore(unsigned ratio = 1) const;
  static polyTg coreToPolyTg(const matrix<poly>& core);

  static etvo::matrix<poly> getMatN(unsigned size);
  matrix<poly> getCoreMax(unsigned ratio = 1) const;
 
  void toPov(graphicPR::PovRay & pov, graphicPR::PovRay::Color c);

private:
  void sort();
  void simplify();
    std::vector<Tg> _poly;
};
std::ostream & operator<<(std::ostream & st,const polyTg & p);
}

#endif // POLYTG_H