growthrate.h

Go to the documentation of this file.

/*
Copyright (C) 2013-2017
Rafael Guglielmetti, rafael.guglielmetti@unifr.ch
*/
 
/*
This file is part of CoxIter.
 
CoxIter is free software: you can redistribute it and/or modify
it under the terms of the GNU General Public License as
published by the Free Software Foundation, either version 3 of the
License, or (at your option) any later version.
 
CoxIter is distributed in the hope that it will be useful,
but WITHOUT ANY WARRANTY; without even the implied warranty of
MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
GNU General Public License for more details.
 
You should have received a copy of the GNU General Public License
along with CoxIter. If not, see <http://www.gnu.org/licenses/>.
*/
 
#ifndef GROWTHRATE_H
#define GROWTHRATE_H
 
#include <pari/pari.h>
#include <vector>
#ifdef _USE_LOCAL_GMP_
#include "gmpxx.h"
#else
#include <gmpxx.h>
#endif
#include <iostream>
 
#include "lib/paripolynomials.h"
 
using namespace std;
using namespace PariPolynomials;
 
struct GrowthRate_Result {
  int perron;
  int pisot;
  int salem;
  string growthRate;
  bool isComputed;
};
 
class GrowthRate {
private:
  vector<GEN>
      t_POLfactors; 
 
  GEN gGrowthRate;   
  GEN gMaximalRoots; 
  long int indexMaximalRoot; 
 
  GEN gEpsilon;           
  long int pariPrecision; 
 
public:
  GrowthRate();
  ~GrowthRate();
 
  GrowthRate_Result grrComputations(vector<mpz_class> polynomial,
                                    const bool &onlyGrowthRate = false);
 
private:
  void irreducibleFactors(const vector<mpz_class> &polynomial);
 
  void minimalRoot();
 
  long int numberRootsUnitCircle(GEN gPol);
};
 
#endif // GROWTHRATE_H