infinitensymetries.h

Go to the documentation of this file.

/*
Copyright (C) 2014, 2015, 2016
Rafael Guglielmetti, rafael.guglielmetti@unifr.ch
*/

/*
This file is part of AlVin.

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 AlVin. If not, see <http://www.gnu.org/licenses/>.
*/

#ifndef INFINITENSYMETRIES_H
#define INFINITENSYMETRIES_H

#include "alvin.h"
#include "CoxIter/lib/numbers/rational.h"

#include <eigen3/Eigen/Dense>
#include <igraph/igraph.h>

#ifdef _OPENMP
#include <omp.h>
#endif

using namespace Eigen;

struct GraphInvolution
{
    vector< unsigned int > iVertices; 
    vector< unsigned int > iPermutation; 
    vector< unsigned int > iLinearlyIndependantVertices; 
    unsigned int iOrbitCount; 
};

class InfiniteNSymetries
{
    protected:
        AlVin* alvin;
        vector< AlgebraicInteger* > aiQF;
        
        vector< vector< unsigned int > > iGraphMatrix; 
        vector< vector< unsigned int > > iCoxeterMatrix;
        
        unsigned int iVectorsCount; 
        unsigned int iDimension; 
        const unsigned int iVectorSize; 
        
        unsigned int iFixedPointsDimension; 
        bool bFinished; 
        
        vector< GraphInvolution > usefulInvolutions;
        
    public:
        InfiniteNSymetries( AlVin* alvin );
        bool Run( const unsigned int& iNRMin, const unsigned int& iNRMax );
        
        unsigned int get_iFixedPointsDimension() const;
        virtual void print_basisFixedPoints( const string& strSpacer = "" ) const = 0;
        
        virtual ~InfiniteNSymetries();
        
        vector< GraphInvolution > get_usefulInvolutions() const;
        
    protected:
        virtual bool bDottedSameWeight( const unsigned int& v1, const unsigned int& v2, const unsigned int& w1, const unsigned int& w2 ) const = 0;
        virtual bool FindIntegralSymmetryFromSubgraph( const vector< unsigned int >& iVertices ) = 0;
};

#endif // INFINITENSYMETRIES_H