polyEd.h

Go to the documentation of this file.

#ifndef POLYED_H
#define POLYED_H

#include "Ed.h"
#include "../wrapperMMGD/matrixWrapper.h"
#include "../wrapperMMGD/seriesWrapper.h"
#include "../common/ISterm.h"
#include "../wrapperMMGD/polyWrapper.h"

namespace etvo
{ 
 
 class seriesEd;  // anticipate declaration

class polyEd:public ISterm
{
public:
  static polyEd Epsilon();  
  static polyEd E();
  static polyEd Top();

  polyEd();
  polyEd(bool TopNotE); //true Top false E
  polyEd(const Ed & m);
  polyEd(const std::vector<Ed> & v);
    
  bool isCanon() const;
  void canon();
  
   bool isE() const;
  
  polyEd oplus(const polyEd &p) const;  
  polyEd operator+(const polyEd &p) const;
  polyEd oplusCD(const polyEd &p) const;
  polyEd operator+(const Ed & m)const;

  void add(const Ed & m);
  polyEd operator+=(const Ed & m);
  
  polyEd operator*(const polyEd & p) const;
  polyEd otimes(const polyEd &p) const;
  polyEd operator*(const Ed & m) const;

  static polyEd otimes(const Ed & m,const polyEd & p);
  polyEd otimesCD(const polyEd &p) const;
  
  polyEd inf(const polyEd & p) const;
  polyEd infCD(const polyEd & p) const;

  seriesEd star() const;

  polyEd lfrac(const polyEd &) const;
  polyEd lfracCD(const polyEd &) const;
  polyEd lfrac(const Ed & m) const;
  polyEd rfrac(const polyEd &) const;
  polyEd rfracCD(const polyEd &) const;
  polyEd rfrac(const Ed & m) const;
  
  bool operator==(const polyEd &) const;
  bool operator!=(const polyEd &) const;
  bool operator<=(const polyEd &) const;
  bool operator>=(const polyEd & ) const; 
  
  static polyEd toPolyEd(const poly & p);
  poly toPoly() const;
  static polyEd toCausal(const polyEd & p);
  Ed getFirstDif(const polyEd & p) const;

  polyEd transientStar(int Tmax) const;

  void getMaxGain(unsigned int & mu, unsigned int & beta) const;
  void getLcmGain(unsigned int & mu, unsigned int & beta) const;
  std::pair<unsigned int, unsigned int> getPeriodicity() const;

  std::vector<Ed> getTerms() const;
  void removeTerm(unsigned idx);

  Ed operator[](unsigned idx) const;
  unsigned int size() const;

  std::string toString() const;
  std::string toStringAsMuVar() const;
    
   

   matrix<poly> getCore(unsigned ratio=1) const; 

   static polyEd coreToPolyEd(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);   

  polyEd odot(const polyEd &p) const;
  polyEd osum(const polyEd &p) const;
  
private:
  void sort();
  void simplify();
  //void SimplifyUpTo(unsigned idx); //TODO to improve the sum with Ed
    std::vector<Ed> _poly;
};
std::ostream & operator<<(std::ostream & st,const polyEd & p);
}
#endif // POLYED_H