Mult_Table.hpp

Go to the documentation of this file.

#pragma once
#include "Utility/OpenMP_Macros.hpp"
#include "Utility/General.hpp"
#include "Spaces/Point.hpp"
 
 
namespace mackey
{
    template<typename group_t>
    struct MultTableData {
        int level; 
        std::map<std::vector<int>, IDGenerators<typename group_t::rank_t>> NonZeroHomology;
        std::vector<std::vector<Green<group_t>>> Greens;            
        std::map<std::array<int, 3>, Green<group_t> > tripleGreens; 
        std::vector<std::vector<int>> degree;       
        std::map<std::vector<int>, int> antidegree; 
        Eigen::Matrix<int, -1, -1> index_product;   
        std::vector<int> minsphere;
        std::vector<int> maxsphere;
        std::vector<std::vector<int>> basicIrreducibles;
    };
 
    template<typename group_t>
    class MultTable : public MultTableData<group_t>{
        typedef typename group_t::rank_t rank_t;
        typedef typename group_t::diff_t diff_t;
 
    public:
 
        using MultTableData<group_t>::level;
        using MultTableData<group_t>::NonZeroHomology;
        using MultTableData<group_t>::Greens;
        using MultTableData<group_t>::tripleGreens;
        using MultTableData<group_t>::degree;
        using MultTableData<group_t>::antidegree;
        using MultTableData<group_t>::index_product;
        using MultTableData<group_t>::minsphere;
        using MultTableData<group_t>::maxsphere;
        using MultTableData<group_t>::basicIrreducibles;
 
        MultTable(const MultTableData<group_t>& MTD);
 
        MultTable(MultTableData<group_t>&& MTD);
 
        MultTable(int, const std::vector<int>&, const std::vector<int>&, const std::vector<std::vector<int>>&);
 
        int getdegreeindex(const std::vector<int>&) const;
 
        bool degreewithinrange(const std::vector<int>&) const;
 
        Green<group_t> triple_product(int i, int j1, int j2) const;
 
    protected:
 
        std::vector<int> getsphere(const std::vector<int>& deg) const;
 
 
    private:
        Hash<std::vector<int>, std::vector<int>> perfect_hasher;
        std::vector<Chains<rank_t, diff_t>> basicChains;
        std::map<int, Chains<rank_t, diff_t>> sphereChains;
        void getIrreducibleChains();
        void multiply_all_indices();
        void multiply(int, int);
        void make();
        void generators();
        void compute(const std::vector<std::pair<int, int>>&);
    };
}
#include "impl/Mult_Table.ipp"