Homology.hpp

Go to the documentation of this file.

#pragma once
#include "Chain_Complexes/Chains.hpp"
#include "Utility/General.hpp"
#include "Smith.hpp"
#include "Abelian.hpp"
 
 
namespace mackey {
 
    //
    template<typename rank_t, typename diff_t>
    class Homology {
 
        typedef std::conditional_t<SFINAE::is_finite_cyclic<scalar_t<diff_t>>::value, scalar_t<diff_t>, int64_t> HScalar; //scalar with extended accuracy for Z coefficients (no change for Z/N coefficients)
        typedef std::conditional_t<SFINAE::is_Dense<diff_t>::value, Eigen::Matrix<HScalar, -1, -1>, Eigen::SparseMatrix<HScalar, 0, storage_t<diff_t>>> diff_t_C; //diff_t column major with the scalar above
        typedef std::conditional_t<SFINAE::is_Dense<diff_t>::value, Eigen::Matrix<HScalar, -1, -1, 1>, Eigen::SparseMatrix<HScalar, 1, storage_t<diff_t>>> diff_t_R; //diff_t row major with the scalar above
    public:
 
        typedef diff_t_C Gens_t;
        typedef col_vector_t<diff_t_C> gen_t;
 
        AbelianGroup<rank_t> group;
        Gens_t generators;
        bool isZero;
        
        Homology()=default;
 
        Homology(const Junction<rank_t, diff_t>& J, bool getQ=0);
 
 
        rank_t basis(const gen_t&) const;
 
        gen_t boundary(const gen_t&) const;
 
    private:
        diff_t_R Out_Qi, In_P_full, In_P_reduced;
        std::vector<typename diff_t::StorageIndex> dontModOut;
        diff_t_C In_Q;
        row_vector_t<diff_t_C> diagonal;
        typename diff_t::StorageIndex M;
        diff_t_C getKernel(diff_t_C&);
        void KernelModImage(diff_t_C&, diff_t_C&, bool);
 
    };
}
#include "impl/Homology.ipp"