Additive.hpp

Go to the documentation of this file.

#pragma once
#include "Utility/OpenMP_Macros.hpp"
#include "Spaces/Point.hpp"
#include <iomanip>
 
namespace mackey {
 
    template<typename group_t>
    class AdditiveStructure {
    public:
        typedef typename group_t::rank_t rank_t;
        typedef typename group_t::diff_t diff_t;
 
        std::vector<int> minsphere;
        std::vector<int> maxsphere;
 
        AdditiveStructure(const std::vector<int>& minsphere, const std::vector<int>& maxsphere);
 
        template<typename sphere_t>
        std::vector<MackeyFunctor<rank_t>> getMackey(const sphere_t& sphere) const;
 
        template<typename sphere_t>
        MackeyFunctor<rank_t> getMackey(int degree, const sphere_t& sphere) const;
 
        template<typename T>
        void print_unique(T& stream);
 
        template<typename T>
        void print_unknown(T& stream);
 
        const std::vector<MackeyFunctor<rank_t>>& identified();
 
        const std::vector<MackeyFunctor<rank_t>>& unknown();
 
    private:
 
        void identify();
 
        std::vector<int> true_unknowns;
        std::vector<MackeyFunctor<rank_t>> uniqueMackeys;
        std::vector<MackeyFunctor<rank_t>> unknownMackeys;
        std::vector<std::vector<MackeyFunctor<rank_t>>> allMackeys;
        std::vector<char> found_unknown;
        void compute();
        void identify_unique(MackeyFunctor<rank_t>&, bool);
        void identify_unknown(int);
 
        void check_unique_and_assign(const MackeyFunctor<rank_t>&);
        void check_unknown_and_assign(const std::vector<MackeyFunctor<rank_t>>&);
 
        template<typename _group>
        friend std::ostream& operator<<(std::ostream& os, const AdditiveStructure<_group>&);
    };
 
    template<typename group_t>
    std::ostream& operator<<(std::ostream& os, const AdditiveStructure<group_t>& A);
}
 
#include "impl/Additive.ipp"