growthrate.h

Go to the documentation of this file.

/*
Copyright (C) 2013, 2014, 2015, 2016
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 <vector>
#include <pari/pari.h>
#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 iPerron;
    int iPisot;
    int iSalem;
    string strGrowthRate;
    bool bComputed;
};

class GrowthRate
{
    private:
        vector< GEN > t_POLfactors; 
        
        GEN gGrowthRate; 
        GEN gMaximalRoots; 
        long int iIndexMaximalRoot; 
        
        GEN gEpsilon; 
        long int iPariPrecision; 
        
    public:
        GrowthRate();
        ~GrowthRate();
        
        GrowthRate_Result grrComputations( vector< mpz_class > iPolynomial, const bool& bOnlyGrowthRate = false );
        
    private:
        void irreducibleFactors( const vector< mpz_class >& iPolynomial );
        
        void minimalRoot();

        long int iNumberRootsUnitCircle( GEN gPol );
};

#endif // GROWTHRATE_H