Massey.hpp

Go to the documentation of this file.

#pragma once
#include "Green.hpp"
 
 
namespace mackey {
 
    template<typename group_t>
    class Massey {
        typedef typename group_t::rank_t rank_t;
        typedef typename group_t::diff_t diff_t;
 
    public:
        rank_t basis; 
        rank_t normalBasis; 
        AbelianGroup<rank_t> group; 
        std::array<rank_t, 2> indeterminacy; 
        IDGenerators<rank_t> boxID; 
        bool exists; 
        bool noIndeterminacy; 
        bool isZero; 
 
        Massey();
 
        Massey(const chains_t<group_t>&, const chains_t<group_t>&, const chains_t<group_t>&, int, int, int, int, int, int, int, bool);
    };
 
 
    namespace internal
    {
        template<typename group_t>
        class MasseyCompute {
            typedef typename group_t::rank_t rank_t;
            typedef typename group_t::diff_t diff_t;
 
            Massey<group_t> Mass;
 
            const chains_t<group_t>& C, D, E;
            gen_t<rank_t, diff_t> resgenC, resgenD, resgenE, boundaryCD, boundaryDE;
 
            chains_t<group_t> BoxCD, BoxDE;
            Junction<rank_t, diff_t> Box;
 
            IDGeneratorCompute<group_t> ID;
 
            std::vector<std::vector<int64_t>> CD_detailedrank, DE_detailedrank;
            std::vector<int64_t> CD_E_detailedrank, C_DE_detailedrank;
 
            const int level, degreeC, degreeD, degreeE, selectC, selectD, selectE;
            Eigen::PermutationMatrix<-1, -1, int64_t> BoxDE_to_BoxCD;
 
            typedef typename mackey::Tensor<chains_t<group_t>, std::vector<std::vector<int64_t>>> tensortype;
            gen_t<rank_t, diff_t> box_boundary(const chains_t<group_t>&, const chains_t<group_t>&, const tensortype&, const gen_t<rank_t, diff_t>&, const gen_t<rank_t, diff_t>&, int, int);
            gen_t<rank_t, diff_t> product_bottom(const chains_t<group_t>&, const chains_t<group_t>&, const std::vector<int64_t>&, const gen_t<rank_t, diff_t>&, const gen_t<rank_t, diff_t>&, int, int);
 
            MasseyCompute(const chains_t<group_t>&, const chains_t<group_t>&, const chains_t<group_t>&, int, int, int, int, int, int, int, bool);
 
            void getGens();
            void box();
            void compute();
            void indeterminacy();
 
            template<typename>
            friend class mackey::Massey;
        };
    }
}
#include "impl/Massey.ipp"