rationalinteger_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 RATIONALINTEGER_INFINITENSYMETRIES_H
#define RATIONALINTEGER_INFINITENSYMETRIES_H

#include "infinitensymetries.h"
#include "rationalinteger.h"

namespace Eigen 
{
    template<class> struct NumTraits;
    template<> struct NumTraits<Rational<RationalInteger>> 
    {
        typedef Rational<RationalInteger> Real;
        typedef Rational<RationalInteger> NonInteger;
        typedef Rational<RationalInteger> Nested;
        typedef Rational<RationalInteger> Literal;
        
        static inline Real epsilon() { return 0; }
        static inline Real dummy_precision() { return 0; }
        enum {
            IsInteger = 0,
            IsSigned = 1,
            IsComplex = 0,
            RequireInitialization = 1,
            ReadCost = Eigen::HugeCost,
            AddCost = Eigen::HugeCost,
            MulCost = Eigen::HugeCost
        };
    };
    
    namespace internal {
        template<>
        struct significant_decimals_impl<Rational<RationalInteger>> 
        {
            // Infinite precision when printing
            static inline int run() { return 0; }
        };
    }
};

class RationalInteger_InfiniteNSymetries : public InfiniteNSymetries
{
    private:
        vector< Matrix< Rational<RationalInteger>, Dynamic, 1 > > rriVectorsC; 
        vector< vector< Rational<RationalInteger> > > riDottedWeights; 
        
        Matrix< Rational<RationalInteger>, Dynamic, Dynamic > rriBasisFixedPoints; 
        
        vector< Rational<RationalInteger> > rriQF; 
        vector< vector< RationalInteger > > riVectorsProducts; 

    public:
        RationalInteger_InfiniteNSymetries( AlVin* alvin );
        ~RationalInteger_InfiniteNSymetries();
        
        Rational<RationalInteger> rriVectorNorm( const Matrix< Rational<RationalInteger>, Dynamic, 1 >& rriV ) const;
        Rational<RationalInteger> rriVectorsProduct( const Matrix< Rational<RationalInteger>, Dynamic, 1 >& rriV1, const Matrix< Rational<RationalInteger>, Dynamic, 1 >& rriV2 ) const;
        bool FindIntegralSymmetryFromSubgraph( const vector< unsigned int >& iVertices );
        
        virtual void print_basisFixedPoints( const string& strSpacer = "" ) const;
        
    private:
        bool bDottedSameWeight( const unsigned int& v1, const unsigned int& v2, const unsigned int& w1, const unsigned int& w2 ) const;
        void WorkWithIsomorphisms( const vector< unsigned int >& iVertices, const vector< GraphInvolution >& iIsomorphisms );
        vector< GraphInvolution > FindIsomorphismsInSubgraph( const vector< unsigned int >& iVertices );
};


#endif // RATIONALINTEGER_INFINITENSYMETRIES_H